CN109657196B - MOM-KA mixing method based on airspace strong coupling path acceleration - Google Patents

MOM-KA mixing method based on airspace strong coupling path acceleration Download PDF

Info

Publication number
CN109657196B
CN109657196B CN201811429153.2A CN201811429153A CN109657196B CN 109657196 B CN109657196 B CN 109657196B CN 201811429153 A CN201811429153 A CN 201811429153A CN 109657196 B CN109657196 B CN 109657196B
Authority
CN
China
Prior art keywords
target
environment
mom
strong coupling
coupling
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
CN201811429153.2A
Other languages
Chinese (zh)
Other versions
CN109657196A (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.)
Shanghai Radio Equipment Research Institute
Original Assignee
Shanghai Radio Equipment Research 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 Shanghai Radio Equipment Research Institute filed Critical Shanghai Radio Equipment Research Institute
Priority to CN201811429153.2A priority Critical patent/CN109657196B/en
Publication of CN109657196A publication Critical patent/CN109657196A/en
Application granted granted Critical
Publication of CN109657196B publication Critical patent/CN109657196B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Computing Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention relates to an MOM-KA mixing method based on airspace strong coupling path acceleration, which comprises the following steps of: s1, calculating target initial current; s2, acquiring a strong coupling path between the target and the environment, and establishing a strong coupling area set of the target and the environment; s3, calculating mutual coupling between the target and the environment according to the strong coupling area set established in the step S2; and S4, calculating the target current after the target and the environment are coupled, judging whether iteration converges according to the current change rate, if so, ending the iteration process, and if not, repeating the steps S3-S4. According to the invention, a coupling path between the target and the rough surface is provided, a strong coupling area set between the target and the rough surface is established, only the coupling between the target with strong mutual coupling and the rough surface area is calculated, and the calculation efficiency can be greatly improved on the premise of ensuring the solving precision.

Description

MOM-KA mixing method based on airspace strong coupling path acceleration
Technical Field
The invention relates to a high-frequency and low-frequency mixing method, in particular to an MOM-KA mixing method based on airspace strong coupling path acceleration, and belongs to the technical field of target and environmental characteristics.
Background
The electromagnetic scattering property of the target in the background of the sea surface has a basic supporting effect on technologies such as radar detection and identification of the target. Therefore, the modeling research of the electromagnetic scattering characteristics of the target under the rough surface environment has important significance. The moment method (MOM) is a strict numerical method, can simulate the electromagnetic scattering of targets with any shapes, and has the characteristics of accurate simulation and no limitation on the shapes of the targets. The analysis method of rough surface scattering calculation, such as Kirchhoff Approximation (KA), is combined with the MOM numerical method of the body target, and the coupling calculation of the surface-body target scattering interaction is completed, so that the solving speed of the composite scattering of the large-scale rough surface and the complex target can be essentially increased, and meanwhile, a very clear physical explanation is given by the analysis solution.
In the MOM-KA mixing method based on the iteration method adopted in the prior art, the influence of a KA region on the MOM region is used as the variable quantity of an excitation term on the right side of an integration equation of the MOM region in each iteration, and therefore the convergence speed is accelerated. However, the method does not accelerate the coupling calculation of the KA region and the MOM region, and single iteration calculation is time-consuming.
In the MOM-KA hybrid method based on fast multi-pole algorithm (MLFMA) acceleration adopted in the prior art, the MLFMA is used for accelerating the solution of a target MOM area, and meanwhile, the complexity of coupling calculation of the MOM area and the KA area is reduced by combining methods such as fast far-field approximation and the like. However, this method is still limited in its applicability due to the large area of the matte.
In the prior art, coupling areas of different patches under excitation of an incident electromagnetic field are completed by ray tracing, then impedance matrixes of nearby sub-areas and coupling areas are filled by a moment method, impedance matrixes of non-nearby sub-areas and non-coupling areas are arranged, the whole impedance matrix sparseness degree is high, and memory space and calculation time are greatly saved.
Disclosure of Invention
The invention aims to provide an MOM-KA mixing method based on airspace strong coupling path acceleration, according to a target and environment strong coupling path, only the mutual coupling of a target and an environment with a strong coupling area is calculated, the calculation efficiency can be greatly improved on the premise of ensuring the solving precision, and the defects and the limitations in the prior art are overcome.
In order to achieve the above object, the present invention provides a MOM-KA mixing method based on airspace strong coupling path acceleration, comprising the following steps:
s1, calculating target initial current;
s2, acquiring a strong coupling path between the target and the environment, and establishing a strong coupling area set of the target and the environment; s3, calculating mutual coupling between the target and the environment according to the strong coupling area set established in the step S2;
and S4, calculating a target current after the target is coupled with the environment, judging whether iteration is converged according to the current change rate, if so, ending the iteration process, and if not, repeating the steps S3-S4.
In the step S1, triangular mesh discretization is carried out on the target, RWG basis functions are defined, a matrix equation (1) obtained by discretization of a target Electric Field Integral Equation (EFIE) is solved, the right side of the equation is an initial excitation term, namely, plane wave irradiation is added, and initial current I is obtained 0，MOM 。
ZI 0，MOM ＝V 0 (1)
The step S2 specifically includes the following steps,
s21, partitioning the target and the environment, and decomposing the target into N 1 Blocks, respectively marked as
Figure BDA0001882308800000021
The environment is decomposed into N 2 Blocks, respectively marked as->
Figure BDA0001882308800000022
Ray tracing along the direction of electromagnetic wave incidence, recording ray paths and labeling as a sequence of target, environmental block numbers, e.g. M i ,P j ,…M k 8230, classifying rays according to paths;
s22, for each path, respectively extracting the grouping of the target with strong coupling to the environment according to the wave propagation direction
Figure BDA0001882308800000023
And a packet with a stronger coupling of the environment to the target>
Figure BDA0001882308800000024
All paths are traversed to obtain a total set U of the target with strong coupling to the environment MP And a total set U with stronger coupling to the target from the environment PM 。
The step S3 specifically includes the following steps,
s31, in the iterative process of the step i, according to the total set U which is obtained in the step S2 and has stronger coupling of the target to the environment MP Target faradaic coefficient I calculated from step I-1 i-1,MOM Compute set U MP Middle target block
Figure BDA0001882308800000025
On a rough surface block->
Figure BDA0001882308800000026
Induced current I generated by scattered field i,KA The method specifically comprises the following steps:
Figure BDA0001882308800000027
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0001882308800000028
is the target block>
Figure BDA0001882308800000029
RWG base function of->
Figure BDA00018823088000000210
Is a rough surface block>
Figure BDA00018823088000000211
Normal direction of the surface element, is greater than or equal to>
Figure BDA00018823088000000212
The expression is as follows
Figure BDA0001882308800000031
Where g (r, r') is a scalar gray function.
S32, in the iterative process of the step i, according to the total set U which is obtained in the step S2 and has stronger coupling to the target PM Rough surface current I obtained according to step S31 i,KA Updating the right excitation term of the target matrix equation (1), specifically:
Figure BDA0001882308800000032
Figure BDA0001882308800000033
wherein the content of the first and second substances,
Figure BDA0001882308800000034
is expressed as
Figure BDA0001882308800000035
Wherein j is an imaginary unit, ω is an angular frequency, μ is a magnetic permeability,
Figure BDA0001882308800000036
is a unit vector.
In the step S4, the excitation item V updated according to the step S33 i,MOM Solving the matrix equation (1) again to obtain a new target induced current I i,MOM (ii) a Defining an error decision function epsilon for the ith iteration i Comprises the following steps:
Figure BDA0001882308800000037
wherein, I i,MOM And I i-1,MOM Coefficients of the target induced current calculated for steps i and i-1, respectively; such as e i If the value is less than the set threshold value, the iteration converges; such as e i And if the threshold value is larger than the set threshold value, repeating the steps S3 to S4.
In summary, the MOM-KA hybrid method based on airspace strong coupling path acceleration provided by the invention extracts a strong coupling path between a target and a rough surface based on ray tracing and diversity; according to the space strong coupling path between the target and the rough surface, only the coupling between the target with strong mutual coupling and the rough surface area is calculated, and the calculation efficiency can be greatly improved on the premise of ensuring the solving precision.
Drawings
FIG. 1 is a flow chart of an MOM-KA hybrid method based on airspace strong coupling path acceleration in the present invention;
FIG. 2 is a schematic diagram of the strong coupling path between the target and the rough surface airspace in the present invention.
Detailed Description
The technical contents, construction features, achieved objects and effects of the present invention will be described in detail by preferred embodiments with reference to fig. 1 to 2.
As shown in fig. 1, the MOM-KA mixing method based on spatial domain strong coupling path acceleration provided by the present invention includes the following steps:
s1, calculating target initial current;
s2, acquiring a strong coupling path between the target and the environment, and establishing a strong coupling area set of the target and the environment;
s3, calculating mutual coupling between the target and the environment according to the strong coupling area set established in the step S2;
and S4, calculating the target current after the target and the environment are coupled, judging whether iteration converges according to the current change rate, if so, ending the iteration process, and if not, repeating the steps S3-S4.
In the step S1, triangular mesh discretization is carried out on the target, RWG basis functions are defined, a matrix equation (1) obtained by discretization of a target Electric Field Integral Equation (EFIE) is solved, the right side of the equation is an initial excitation term, namely, plane wave irradiation is added, and initial current I is obtained 0，MOM 。
ZI 0，MOM ＝V 0 (1)
Wherein Z is an impedance matrix, I 0，MOM The moment method region, i.e., the target initial current, is represented.
The step S2 specifically includes the following steps:
s21, partitioning the target and the environment, and decomposing the target into N 1 Blocks, respectively marked as
Figure BDA0001882308800000041
The environment is decomposed into N 2 Blocks, respectivelyIs marked as->
Figure BDA0001882308800000042
Ray tracing along the direction of electromagnetic wave incidence, recording ray paths and labeling as a sequence of target, environmental block numbers, e.g. M i ,P j ,…M k 8230, classifying rays according to paths;
s22, for each path, respectively extracting the grouping of the target with strong coupling to the environment according to the wave propagation direction
Figure BDA0001882308800000043
And packets for which the context has a stronger coupling to the target>
Figure BDA0001882308800000044
All paths are traversed to obtain a total set U of targets with strong coupling to the environment MP And a total set U with strong coupling to the target from the environment PM . As shown in FIG. 2, the packet in coupled path 1 with the stronger coupling of the target to the environment is (M) 1 ,P 1 ) (ii) a The packet in coupled path 2 with the stronger coupling of the target to the environment is
Figure BDA0001882308800000045
The packet whose environment has strong coupling to the target is (P) 3 ,M 2 )；
The step S3 specifically includes the following steps:
s31, in the iterative process of the step i, according to the total set U which is obtained in the step S2 and has stronger coupling of the target to the environment MP Target faradaic coefficient I calculated from step I-1 i-1,MOM Compute set U MP Middle target block
Figure BDA0001882308800000046
On a rough surface block->
Figure BDA0001882308800000047
Induced current I generated by scattered field i,KA The method specifically comprises the following steps:
Figure BDA0001882308800000048
wherein the content of the first and second substances,
Figure BDA0001882308800000051
is the target block>
Figure BDA0001882308800000052
RWG base function of->
Figure BDA0001882308800000053
For a block of rough surface>
Figure BDA0001882308800000054
Normal direction of the surface element, is greater than or equal to>
Figure BDA0001882308800000055
The expression is as follows
Figure BDA0001882308800000056
Where g (r, r') is a scalar gray function.
S32, in the iterative process of the step i, according to the total set U which is obtained in the step S2 and has stronger coupling to the target PM The rough surface current I obtained in step S31 i,KA Updating the right excitation term of the target matrix equation (1), specifically:
Figure BDA0001882308800000057
Figure BDA0001882308800000058
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0001882308800000059
is expressed as
Figure BDA00018823088000000510
Wherein j is an imaginary unit, ω is an angular frequency, μ is a magnetic permeability,
Figure BDA00018823088000000511
is a unit vector.
In the step S4, the excitation item V updated according to the step S32 i,MOM Solving the matrix equation (1) again to obtain a new target induced current I i,MOM (ii) a Defining an error decision function epsilon for the ith iteration i Comprises the following steps:
Figure BDA00018823088000000512
wherein, I i,MOM And I i-1,MOM Coefficients of the target induced current calculated for steps i and i-1, respectively; such as e i If the value is less than the set threshold value, the iteration converges; such as e i And if the threshold value is larger than the set threshold value, repeating the steps S3 to S4.
In summary, compared with the prior art, the MOM-KA mixing method based on airspace strong coupling path acceleration provided by the present invention has the following advantages and beneficial effects: based on ray tracing and diversity, extracting a strong coupling path between the target and the rough surface, and establishing a strong coupling region set between the target and the rough surface; according to the strong coupling area set between the target and the rough surface, only the coupling between the target with strong mutual coupling and the rough surface area is calculated, and the calculation efficiency can be greatly improved on the premise of ensuring the solving precision.
While the present invention has been described in detail with reference to the preferred embodiments, it should be understood that the above description should not be taken as limiting the invention. Various modifications and alterations to this invention will become apparent to those skilled in the art upon reading the foregoing description. Accordingly, the scope of the invention should be determined from the following claims.

Claims (4)

1. An MOM-KA mixing method based on airspace strong coupling path acceleration is characterized by comprising the following steps of:
s1, calculating a target initial current;
s2, extracting a strong coupling path between the target and the environment, and establishing a target and environment strong coupling area set, which further comprises the steps of;
s21, partitioning the target and the environment, and decomposing the target into N 1 Blocks, respectively marked as
Figure FDA0003955629430000011
The environment is decomposed into N 2 The number of the blocks is such that, are respectively marked as->
Figure FDA0003955629430000012
Ray tracing is carried out along the incident direction of the electromagnetic wave, and the ray path is recorded and marked as a series of sequences of target and environment block numbers: m i ,P j ,...M k The method includes classifying rays according to paths;
s22, for each path, respectively extracting the grouping of the target with strong coupling to the environment according to the wave propagation direction
Figure FDA0003955629430000013
And a packet with a stronger coupling of the environment to the target>
Figure FDA0003955629430000014
All paths are traversed to obtain a total set U of the target with strong coupling to the environment MP And a total set U with stronger coupling to the target from the environment PM ；
S3, calculating mutual coupling between the target and the environment according to the strong coupling area set established in the step S2;
and S4, calculating a target current after the target is coupled with the environment, judging whether iteration is converged according to the current change rate, if so, ending the iteration process, and if not, repeating the steps S3-S4.
2. The MOM-KA mixing method based on spatial domain strong coupling path acceleration as claimed in claim 1, wherein in step S1, the target is triangularly gridded and discretized, RWG basis functions are defined, matrix equation (1) obtained by discretization of target electric field integral equation is solved, the right side of the equation is initial excitation term, i.e. external plane wave irradiation, to obtain initial current I 0，MOM
ZI 0，MOM ＝V 0 (1)。
3. The method of claim 2, wherein the step S3 further comprises the steps of:
s31, in the iterative process of the step i, according to the total set U which is obtained in the step S2 and has stronger coupling of the target to the environment MP Target faradaic coefficient I calculated from step I-1 i-1,MOM Compute set U MP Middle target block
Figure FDA0003955629430000015
On a rough surface block->
Figure FDA0003955629430000016
Induced current I generated by scattered field i,KA The method specifically comprises the following steps:
Figure FDA0003955629430000017
wherein the content of the first and second substances,
Figure FDA0003955629430000021
is the target block>
Figure FDA0003955629430000022
RWG base function of->
Figure FDA0003955629430000023
For a block of rough surface>
Figure FDA0003955629430000024
Bin normal of>
Figure FDA0003955629430000025
The expression is as follows
Figure FDA0003955629430000026
Wherein g (r, r') is a scalar gray function;
s32, in the iterative process of the step i, according to the total set U which is obtained in the step S2 and has stronger coupling to the target PM The rough surface induced current I obtained in step S31 i,KA Updating the right excitation term of the target matrix equation (1), specifically:
Figure FDA0003955629430000027
Figure FDA0003955629430000028
wherein, the first and the second end of the pipe are connected with each other,
Figure FDA0003955629430000029
is expressed as
Figure FDA00039556294300000210
/>
Wherein j is an imaginary unit, ω is an angular frequency, μ is a magnetic permeability,
Figure FDA00039556294300000211
is a unit vector.
4. The MOM-KA mixing method based on spatial domain strong coupling path acceleration as claimed in claim 3, wherein in step S4, the updated excitation term V is obtained according to step S32 i,MOM Solving the matrix equation (1) again to obtain a new target induced current I i,MOM (ii) a Defining an error decision function epsilon for the ith iteration i Comprises the following steps:
Figure FDA00039556294300000212
wherein, I i,MOM And I i-1,MOM Coefficients of the target induced current calculated for steps i and i-1, respectively; such as e i If the value is smaller than the set threshold value, iterative convergence is carried out; such as e i And if the threshold value is larger than the set threshold value, repeating the steps S3 to S4.
CN201811429153.2A 2018-11-27 2018-11-27 MOM-KA mixing method based on airspace strong coupling path acceleration Active CN109657196B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811429153.2A CN109657196B (en) 2018-11-27 2018-11-27 MOM-KA mixing method based on airspace strong coupling path acceleration

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811429153.2A CN109657196B (en) 2018-11-27 2018-11-27 MOM-KA mixing method based on airspace strong coupling path acceleration

Publications (2)

Publication Number Publication Date
CN109657196A CN109657196A (en) 2019-04-19
CN109657196B true CN109657196B (en) 2023-03-28

Family

ID=66111610

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811429153.2A Active CN109657196B (en) 2018-11-27 2018-11-27 MOM-KA mixing method based on airspace strong coupling path acceleration

Country Status (1)

Country Link
CN (1) CN109657196B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113935157A (en) * 2021-09-27 2022-01-14 中国人民解放军32802部队 Cross-scale multi-mechanism fusion electromagnetic scattering numerical solving method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104573289A (en) * 2015-02-06 2015-04-29 西安电子科技大学 Electromagnetic scattering simulation method of one-dimensional conductor rough sea surface and two-dimensional conductor target
CN106294287A (en) * 2015-06-08 2017-01-04 西安电子科技大学 A kind of core exterior domain decomposition method based on parallel high order MoM
CN106650048A (en) * 2016-12-05 2017-05-10 中国舰船研究设计中心 Ship and sea mutual-coupling scattering prediction method based on slope distribution

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI474294B (en) * 2013-05-10 2015-02-21 Univ Yuan Ze Method to set a navigation environment

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104573289A (en) * 2015-02-06 2015-04-29 西安电子科技大学 Electromagnetic scattering simulation method of one-dimensional conductor rough sea surface and two-dimensional conductor target
CN106294287A (en) * 2015-06-08 2017-01-04 西安电子科技大学 A kind of core exterior domain decomposition method based on parallel high order MoM
CN106650048A (en) * 2016-12-05 2017-05-10 中国舰船研究设计中心 Ship and sea mutual-coupling scattering prediction method based on slope distribution

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
低空目标与环境复合电磁散射特性研究;陈博韬;《中国博士学位论文全文数据库》;20130315;第1-112页 *

Also Published As

Publication number Publication date
CN109657196A (en) 2019-04-19

Similar Documents

Publication Publication Date Title
CN113158527B (en) Method for calculating frequency domain electromagnetic field based on implicit FVFD
CN110208375A (en) A kind of detection method and terminal device of anchor rod anchored defect
CN104778151A (en) Electromagnetic scattering analysis method of target with cavity on the basis of moment method and parabolic equation
CN111046603A (en) Electromagnetic scattering characteristic analysis method based on GPU parallel acceleration characteristic basis function algorithm
Guo et al. Hierarchical matrices method and its application in electromagnetic integral equations
CN105786765A (en) Method for generating incentive irrelevant characteristic basis function rapidly in self-adaption mode
CN109657196B (en) MOM-KA mixing method based on airspace strong coupling path acceleration
He et al. A multihybrid FE-BI-KA technique for 3-D electromagnetic scattering from a coated object above a conductive rough surface
CN105630740B (en) Emi analysis method based on matrix Taylor series expansion
CN112799063B (en) Double-station scattering center modeling method based on bouncing ray method
CN112784459A (en) Electromagnetic simulation method based on compression type finite element tearing and butt joint method
Guo et al. A novel JMCFIE-DDM for analysis of EM scattering and radiation by composite objects
CN104778286B (en) Sea skimming device Electromagnetic Scattering Characteristics rapid simulation method
Gu et al. Integral equation‐based nonoverlapping domain decomposition method for electromagnetic scattering analysis of PEC targets in half‐space environment
CN105205299B (en) The quick Dimension Reduction Analysis method of TV university electromagnetic characteristic of scattering
CN114491992B (en) Efficient electromagnetic scattering method based on equivalent dipole moment and physical optical method
CN113567941B (en) Rapid simulation method for near-field Doppler echo of target
CN105760343A (en) Time-domain integration method for analyzing electromagnetic scattering characteristics of multiple complicated metal targets
Guo et al. Investigation of Adam for low-frequency electromagnetic problems
Lei et al. Study on MPI-based parallel modified conformal FDTD for 3-D electrically large coated targets by using effective parameters
Gu et al. Parallel Integral Equation Based Non-overlapping DDM for Fast Solving Electromagnetic Scattering Problems with Changeable Parts
Guan et al. A multi-solver framework for electromagnetic scattering from complex objects
Delgado et al. Application of a sparsity pattern and region clustering for near field sparse approximate inverse preconditioners in method of moments simulations
Guan et al. A CFIE-based electromagnetic solver for composite objects
Wei et al. Domain decomposition WLP-FDTD for efficient electromagnetic scattering analysis

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