CN117034645B - YIG microstrip nonlinear amplitude limiting device joint design method based on numerical calculation and three-dimensional electromagnetic field simulation - Google Patents

YIG microstrip nonlinear amplitude limiting device joint design method based on numerical calculation and three-dimensional electromagnetic field simulation Download PDF

Info

Publication number
CN117034645B
CN117034645B CN202311074134.3A CN202311074134A CN117034645B CN 117034645 B CN117034645 B CN 117034645B CN 202311074134 A CN202311074134 A CN 202311074134A CN 117034645 B CN117034645 B CN 117034645B
Authority
CN
China
Prior art keywords
amplitude limiting
microstrip
simulation
yig
resonance
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
CN202311074134.3A
Other languages
Chinese (zh)
Other versions
CN117034645A (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.)
Chengdu Normal University
Original Assignee
Chengdu Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chengdu Normal University filed Critical Chengdu Normal University
Priority to CN202311074134.3A priority Critical patent/CN117034645B/en
Publication of CN117034645A publication Critical patent/CN117034645A/en
Application granted granted Critical
Publication of CN117034645B publication Critical patent/CN117034645B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Motors That Do Not Use Commutators (AREA)

Abstract

The invention belongs to the technical field of microwaves and provides a combined design method of a YIG microstrip nonlinear amplitude limiting device based on numerical calculation and three-dimensional electromagnetic field simulation, which comprises the steps of firstly selecting a non-resonance region and a resonance region, then dividing the non-resonance and resonance region in turn according to a parallel coupling microstrip structure, then obtaining an initial radio frequency field response of the YIG microstrip nonlinear amplitude limiting device under excitation of a certain microwave power by using simulation software, and finally calculating attenuation of each divided region of the YIG microstrip nonlinear amplitude limiting device on the input microwave power according to equivalent simulation, thereby obtaining the integral amplitude limiting isolation response of the device; according to the invention, the COMSOL iterative simulation flow is controlled by the MATLAB script so as to realize automatic simulation of the amplitude limiting isolation performance of the amplitude limiting device based on the ferrite material, and the Y I G microstrip amplitude limiting device amplitude limiting isolation response performance designed based on the gyromagnetic nonlinear loss mechanism can be efficiently and accurately obtained by using the design method.

Description

YIG microstrip nonlinear amplitude limiting device joint design method based on numerical calculation and three-dimensional electromagnetic field simulation
Technical Field
The invention belongs to the technical field of microwaves, and particularly relates to a YIG microstrip nonlinear amplitude limiting device joint design method based on numerical calculation and three-dimensional electromagnetic field simulation.
Background
Ferrite is a metal oxide with ferromagnetism. The method can be generally divided into permanent magnetic ferrite, soft magnetic ferrite and gyromagnetic ferrite, and gyromagnetic nonlinear loss effect generated by the ferrite under the excitation of high-power microwave signals can be applied to the technical field of high-power microwave amplitude limiting protection at the front end of a microwave receiver;
at present, the existing commercial three-dimensional electromagnetic simulation tools such as high-frequency structure simulator HFSS, three-dimensional electromagnetic field simulation software CST, multi-physical field simulation software COMSOL and the like can only simulate and solve small signal linear response generated by the action of ferrite materials and electromagnetic waves, and gyromagnetic nonlinear loss response generated by ferrite under large signal excitation cannot be simulated; in addition, some open-source micromagnetic simulation software, such as OOMMF and MuMax, can simulate nonlinear response only on the nano-or micro-scale by simulating quantum mechanical mechanisms.
However, when the simulation object is a three-dimensional model with a millimeter magnitude, the simulation object is limited by the complex property of the magnetic material, the limited computing resources and other factors, so that the micromagnetic simulation software faces great challenges in simulating an electrically large-size ferrite amplitude limiting device;
therefore, a combined design method of YIG microstrip nonlinear amplitude limiting devices based on numerical calculation and three-dimensional electromagnetic field simulation is provided by a person skilled in the art to solve the problems provided by the background art.
Disclosure of Invention
In order to solve the technical problems, the invention provides a combined design method of a YIG microstrip nonlinear amplitude limiting device based on numerical calculation and three-dimensional electromagnetic field simulation, which aims to solve the problems that in the prior art, when a simulation object is a three-dimensional model with a millimeter magnitude, the simulation object is limited by factors such as complex properties of magnetic materials and limited computing resources, so that micromagnetism simulation software faces great challenges in simulating an electrically large-size ferrite amplitude limiting device.
A YIG microstrip nonlinear amplitude limiting device joint design method based on numerical calculation and three-dimensional electromagnetic field simulation comprises the following steps:
s1, firstly, selecting a non-resonance region and a resonance region;
s2, sequentially dividing non-resonance and resonance regions of the parallel coupling microstrip structure according to the parallel coupling microstrip structure;
s3, obtaining initial radio frequency field response of the YIG microstrip nonlinear amplitude limiting device under excitation of certain microwave power by using simulation software;
and S4, finally, calculating attenuation of each divided area of the YIG microstrip nonlinear amplitude limiting device on the input microwave power according to equivalent simulation, so as to obtain the integral amplitude limiting isolation response of the device.
Preferably, in S1, the non-resonant area is a non-resonant area formed by input and output ends with impedance matching of 50Ω; the resonance area is a resonance unit area formed by parallel coupling microstrip structures.
Preferably, in S2, the divided resonance areas are marked as resonance unit #1, resonance units #2, …, resonance unit #n, and the like.
Preferably, the S3 specifically includes: firstly, a three-dimensional electromagnetic field simulation model of a YIG microstrip nonlinear amplitude limiting device with a parallel coupling structure is established in COMSOL electromagnetic simulation software, basic simulation parameters such as an externally applied bias magnetic field, an input microwave frequency, a microwave power and the like are set, and linear response simulation is executed; then calculating and acquiring initial radio frequency magnetic field distribution of the YIG substrate of the device under the excitation of a microwave signal with a certain input power Pin through simulation softwareThen setting critical threshold field h for exciting gyromagnetic nonlinear effect c When the initial RF field is greater than the critical threshold field h for exciting gyromagnetic nonlinear effect c During the process, the radio frequency magnetic field distribution value obtained by the simulation of the input end area is extracted from simulation software>And a critical threshold magnetic field h calculated according to the pulse width parameter tau of the excitation microwave signal c (tau) simultaneous->Form substitution into equivalent magnetic loss tangent parameterAnd will->Equivalent tensor permeability of substrate material introduced into input region>And performing one-time equivalent loss numerical simulation calculation, thereby obtaining the radio frequency magnetic field distribution of microwave power attenuated by the YIG substrate in the input end area.
Preferably, the S4 specifically includes: firstly, the radio frequency magnetic field value of a first parallel coupling resonance region which is adjacent to an input end region is extracted, and then h is combined c (tau) calculation of equivalent magnetic loss tangent parameter of resonant cell #1Equivalent tensor permeability of substrate material brought there +.>In the imaginary part, continuing to execute numerical simulation calculation to obtain radio frequency magnetic field distribution attenuated by the resonance unit # 1; and the other areas are similar, and finally, the final amplitude limiting response result of the parallel coupling YIG microstrip amplitude limiting device is obtained after the attenuation of the output end substrate to the rest microwave power is calculated.
Preferably, the YIG microstrip nonlinear amplitude limiting device comprises parallel coupling resonance structures, the parallel coupling resonance structures take polycrystalline YIG materials as substrates, the YIG microstrip nonlinear amplitude limiting device comprises coupling microstrip units, and radio frequency magnetic fields between the coupling microstrip units are strong coupling resonance fields with extremely uneven distribution.
Preferably, the non-uniform field strength distribution is not only related to the input microwave frequency, power, substrate thickness, material properties; but also on the dimensions of the parallel coupled microstrip line elements, the gap parameters and the relative geometrical positions in the poly YI G substrate.
Compared with the prior art, the invention has the following beneficial effects:
according to the invention, the COMSOL iterative simulation flow is controlled by the MATLAB script so as to realize the automatic simulation of the amplitude limiting isolation performance of the amplitude limiting device based on the ferrite material, and the YIG microstrip amplitude limiting device amplitude limiting isolation response performance designed based on the gyromagnetic nonlinear loss mechanism can be efficiently and accurately obtained by using the design method.
Drawings
Fig. 1 is a schematic diagram of dividing YIG microstrip nonlinear limiter with parallel coupling structure in the embodiment of the invention;
FIG. 2 is a flow chart of a numerical calculation and three-dimensional electromagnetic field joint simulation design in an embodiment of the invention;
FIG. 3 is a schematic diagram of a COMSOL three-dimensional electromagnetic simulation model of a parallel coupling YIG microstrip nonlinear limiting device in an embodiment of the invention, wherein 3-1 is an input end 50Ω impedance matching microstrip line, 3-2 is an output end 50Ω impedance matching microstrip line, 3-3 is a coupling microstrip resonance unit, and 3-4 is a polycrystalline Yttrium Iron Garnet (YIG) substrate;
FIG. 4 is a diagram showing the distribution of the RF magnetic field before and after the microwave loss is input by the joint simulation calculation in the embodiment of the present invention;
fig. 5 is a graph of simulation results of amplitude limiting isolation performance of the YI G microstrip nonlinear amplitude limiting device in an embodiment of the present invention.
Detailed Description
Embodiments of the present invention are described in further detail below with reference to the accompanying drawings and examples. The following examples are illustrative of the invention but are not intended to limit the scope of the invention.
As shown in fig. 1 to 5:
the invention provides a YIG microstrip nonlinear amplitude limiting device joint design method based on numerical calculation and three-dimensional electromagnetic field simulation, which comprises the following steps:
s1, firstly, selecting a non-resonance region and a resonance region;
s2, sequentially dividing non-resonance and resonance areas of the parallel coupling microstrip structure according to the parallel coupling microstrip structure, wherein the divided resonance areas are marked as resonance units #1, resonance units #2, …, resonance units #n and the like, and the division can take into consideration the loss of microwave power caused by gyromagnetic nonlinear loss effect generated by YIG substrate materials when the input microwave power passes through each area.
S3, firstly, establishing three-dimensional electricity of the YIG microstrip nonlinear amplitude limiting device with a parallel coupling structure in COMSOL electromagnetic simulation softwareSetting basic simulation parameters such as an externally applied bias magnetic field, an input microwave frequency, a microwave power and the like, and executing linear response simulation; then calculating and acquiring initial radio frequency magnetic field distribution of the YIG substrate of the device under the excitation of a microwave signal with a certain input power Pin through simulation softwareThen setting critical threshold field h for exciting gyromagnetic nonlinear effect c When the initial RF field is greater than the critical threshold field h for exciting gyromagnetic nonlinear effect c During the process, the radio frequency magnetic field distribution value obtained by the simulation of the input end area is extracted from simulation software>And a critical threshold magnetic field h calculated according to the pulse width parameter tau of the excitation microwave signal c (tau) simultaneous->The form is substituted into the equivalent magnetic loss tangent parameter>And will->Equivalent tensor permeability of substrate material introduced into input region>Performing one-time equivalent loss numerical simulation calculation, thereby obtaining radio frequency magnetic field distribution of microwave power attenuated by the YIG substrate in the input end area;
s4, firstly extracting the radio frequency magnetic field value of a first parallel coupling resonance region which is close to the input end region, and then combining h c (tau) calculation of equivalent magnetic loss tangent parameter of resonant cell #1Equivalent tensor of substrate material brought therePermeability->In the imaginary part, continuing to execute numerical simulation calculation to obtain radio frequency magnetic field distribution attenuated by the resonance unit # 1; and the other areas are similar, and finally, the final amplitude limiting response result of the parallel coupling YIG microstrip amplitude limiting device is obtained after the attenuation of the output end substrate to the rest microwave power is calculated.
Specifically, in S1, the non-resonant area is a non-resonant area formed by input and output ends with impedance matching of 50Ω; the resonance area is a resonance unit area formed by parallel coupling microstrip structures.
Specifically, the YIG microstrip nonlinear amplitude limiting device comprises a parallel coupling resonance structure, the parallel coupling resonance structure takes polycrystalline YIG material as a substrate, the YIG microstrip nonlinear amplitude limiting device comprises coupling microstrip units, a radio frequency magnetic field between the coupling microstrip units is a strong coupling resonance field with extremely uneven distribution, and the uneven field intensity distribution is not only related to the input microwave frequency, power, substrate thickness and material properties; but also on the dimensions of the parallel coupled microstrip line elements, the gap parameters and the relative geometrical positions in the poly YIG substrate.
Embodiments of the present invention will be described in further detail below with reference to the accompanying drawings and examples, to illustrate the practical effects of the present invention:
as shown in fig. 3, fig. 3 is an example of a COMSOL simulation model for a parallel coupled YI G microstrip nonlinear limiter according to the present invention; the frequency domain response of the model structure is bandpass filtering with passband center frequency of 9.6GHz and bandwidth of 400 MHz; the substrate of the model adopts polycrystalline YIG material with relative dielectric constant epsilon r =14.5, saturation magnetization M s =1800 Gauss. The size parameters of the polycrystalline YIG substrate are respectively 22mm in length, 12mm in width and 0.3mm in thickness; further, in the model, the width of a metal conductor of the 50 omega impedance matching input and output microstrip line is 0.17mm; further, the resonance areas in the model are formed into resonance single by two parallel gold-plated microstrip linesA meta-element; the length of the coupled resonance units is l from left to right in sequence 1 =2.2mm,l 2 =2.55mm,l 3 =2.25mm,l 4 =2.64mm,l 5 =2.31 mm; the interval between the resonance units is g from left to right 1 =0.05mm,g 2 =0.25mm,g 3 =0.43mm,g 4 =0.4mm,g 5 =0.05mm。
As shown in fig. 4, fig. 4 is a diagram showing the radio frequency magnetic field distribution of the YIG microstrip amplitude limiting device before and after the amplitude limiting occurs, the radio frequency magnetic field distribution of the YIG microstrip amplitude limiting device before and after the amplitude limiting is started and the finally obtained equivalent amplitude limiting isolation degree value response under the conditions that the microwave frequency is 9.6GHz and the microwave power is 200W, 600W and 1000W respectively.
Compared with the linear response without consumption, the distribution of the radio frequency magnetic field is weaker and weaker, and the loss of microwave power is further increased after microwaves sequentially pass through each region of the YIG microstrip amplitude limiting device in the nonlinear loss response simulation; in addition, as the input power increases, the rf magnetic field coupling between the resonant cells increases, which also results in an increase in nonlinear losses, thus reflecting an increase in the clipping isolation of the YIG microstrip nonlinear clipping device.
As shown in fig. 5, fig. 5 is a simulation result of clipping isolation of an example model of a parallel coupled YIG microstrip nonlinear clipping device according to a simulation design flow in the present invention; under the condition that the externally applied bias magnetic field is 900Gauss, the excitation microwave frequency is 9.6GHz, and the input microwave pulse width is 60ns, 100ns and 200ns respectively. Simulation results show that: the limiting performance of the YIG microstrip device on the microwave narrow pulse signal with the maximum input power of 1000W can reach 20.61dB, and the limiting threshold level is 0.1W at the minimum.
From the above, the method controls the COMSOL iteration simulation flow through the MATLAB script to realize the automatic simulation of the amplitude limiting isolation performance of the amplitude limiting device based on the ferrite material, thereby solving the problem that the commercial electromagnetic simulation software cannot perform simulation due to the lack of the gyromagnetic nonlinear loss physical model.
While embodiments of the present invention have been shown and described above for purposes of illustration and description, it will be understood that the above embodiments are illustrative and not to be construed as limiting the invention, and that variations, modifications, alternatives, and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the invention.

Claims (2)

1. A YIG microstrip nonlinear amplitude limiting device joint design method based on numerical calculation and three-dimensional electromagnetic field simulation is characterized in that: the method comprises the following steps:
s1, firstly, selecting a non-resonance region and a resonance region;
s2, sequentially dividing non-resonance and resonance areas according to the parallel coupling microstrip structure, wherein the divided resonance areas are marked as resonance units #1, resonance units #2, … and resonance units #n;
s3, obtaining initial radio frequency magnetic field response of the YIG microstrip nonlinear amplitude limiting device under excitation of certain microwave power by using simulation software, specifically: firstly, a three-dimensional electromagnetic field simulation model of a YIG microstrip nonlinear amplitude limiting device with a parallel coupling structure is established in COMSOL electromagnetic simulation software, basic simulation parameters such as an externally applied bias magnetic field, an input microwave frequency and a microwave power are set, and linear response simulation is executed; then calculating and acquiring initial radio frequency magnetic field distribution of the YIG substrate of the device under the excitation of a microwave signal with a certain input power Pin through simulation softwareThe method comprises the steps of carrying out a first treatment on the surface of the Then setting critical threshold field for exciting gyromagnetic nonlinear effecth c When the initial RF magnetic field is greater than the critical threshold field for exciting gyromagnetic nonlinear effecth c During the process, the radio frequency magnetic field distribution numerical value obtained by the simulation of the input end area is extracted from the simulation softwareAnd a critical threshold magnetic field calculated according to the pulse width parameter tau of the excitation microwave signalh c (τ) Simultaneous use of two-purpose devicesForm substitution into equivalent magnetic loss tangent parameterAnd willEquivalent tensor permeability of substrate material introduced into input regionPerforming one-time equivalent loss numerical simulation calculation, thereby obtaining radio frequency magnetic field distribution of microwave power attenuated by the YIG substrate in the input end area;
s4, finally, calculating attenuation of each divided area of the YIG microstrip nonlinear amplitude limiting device on input microwave power according to equivalent simulation, so as to obtain the integral amplitude limiting isolation response of the device, and specifically: first, the RF magnetic field value of the first parallel coupled resonance region adjacent to the input end region is extracted, and then combinedh c (τ) Calculating equivalent magnetic loss tangent parameter of resonant cell #1Equivalent tensor permeability of substrate material brought thereIn the imaginary part, continuing to execute numerical simulation calculation to obtain radio frequency magnetic field distribution attenuated by the resonance unit # 1; and the other areas are similar, and finally, the final amplitude limiting response result of the parallel coupling YIG microstrip amplitude limiting device is obtained after the attenuation of the output end substrate to the rest microwave power is calculated.
2. The combined design method of YIG microstrip nonlinear amplitude limiting device based on numerical calculation and three-dimensional electromagnetic field simulation as set forth in claim 1, wherein the method is characterized in that: in S1, the non-resonance area is formed by input and output ends with 50 omega impedance matching; the resonance area is a resonance unit area formed by parallel coupling microstrip structures.
CN202311074134.3A 2023-08-24 2023-08-24 YIG microstrip nonlinear amplitude limiting device joint design method based on numerical calculation and three-dimensional electromagnetic field simulation Active CN117034645B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311074134.3A CN117034645B (en) 2023-08-24 2023-08-24 YIG microstrip nonlinear amplitude limiting device joint design method based on numerical calculation and three-dimensional electromagnetic field simulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311074134.3A CN117034645B (en) 2023-08-24 2023-08-24 YIG microstrip nonlinear amplitude limiting device joint design method based on numerical calculation and three-dimensional electromagnetic field simulation

Publications (2)

Publication Number Publication Date
CN117034645A CN117034645A (en) 2023-11-10
CN117034645B true CN117034645B (en) 2024-02-23

Family

ID=88637220

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311074134.3A Active CN117034645B (en) 2023-08-24 2023-08-24 YIG microstrip nonlinear amplitude limiting device joint design method based on numerical calculation and three-dimensional electromagnetic field simulation

Country Status (1)

Country Link
CN (1) CN117034645B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115663427A (en) * 2022-10-31 2023-01-31 电子科技大学 High-power microwave self-adaptive protection device based on gyromagnetic material

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115663427A (en) * 2022-10-31 2023-01-31 电子科技大学 High-power microwave self-adaptive protection device based on gyromagnetic material

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
X-ray detected magnetic resonance of YIG thin films in the nonlinear regime of spin waves;J. Goulon等;Journal of Magnetism and Magnetic Materials;第2308-2329页 *
基于旋磁非线性传输线的小型化强电磁脉冲源的仿真研究;朱丹妮等;电子与信息学报;第737-744页 *

Also Published As

Publication number Publication date
CN117034645A (en) 2023-11-10

Similar Documents

Publication Publication Date Title
Vlaminck et al. Spin-wave transduction at the submicrometer scale: Experiment and modeling
Apaydin et al. Nonreciprocal leaky-wave antenna based on coupled microstrip lines on a non-uniformly biased ferrite substrate
Goto et al. One-dimensional magnonic crystal with Cu stripes for forward volume spin waves
Marzall et al. Microstrip ferrite circulator design with control of magnetization distribution
CN115663427A (en) High-power microwave self-adaptive protection device based on gyromagnetic material
CN117034645B (en) YIG microstrip nonlinear amplitude limiting device joint design method based on numerical calculation and three-dimensional electromagnetic field simulation
Dou et al. 3D FDTD method for arbitrary anisotropic materials
Abaei et al. Two dimensional multi-port method for analysis of propagation characteristics of substrate integrated waveguide
Gil et al. Left‐handed and right‐handed transmission properties of microstrip lines loaded with complementary split rings resonators
Nelin et al. Effective electromagnetocrystalline inhomogeneities
CN113206362B (en) Non-reciprocal ferrite phase shifter based on artificial surface plasmon transmission line
Salahun et al. Application of ferromagnetic composite in different planar tunable microwave devices
Kayano et al. Embedded F-SIR type transmission line with open-stub for negative group delay characteristic
Khalaj-Amirhosseini Analysis of longitudinally inhomogeneous waveguides using Taylor's series expansion
Peng et al. Design of microstrip Y-junction circulator based on ferrite thin films
Andhita et al. Miniaturization of Circular Resonators Waveguide at SHF Frequency using Floral Foam Electromagnetic Engineering
Reddy et al. Design of Microstrip Based Dual Junction Four-Port Circulator for 8.2 GHz X-band
Popov et al. Design of the Ferrite Band-pass Filter with Dual Tunability
Chen et al. Relate the GSTCs method and ECM through surface impedance in studying periodic structures
CN101909419B (en) Transmission line electromagnetic noise suppressor and preparation method thereof
Aranganadin et al. Effects of Realistic Magnetic Field in Ferrite on the Waveguide Circulator for Industrial Applications
Prabhakar et al. Auto-oscillation thresholds at the main resonance in ferrimagnetic films
Kim et al. Transmission-line analysis of an epsilon near zero tunneling circuit using a double ridge rectangular waveguide
Chao et al. Gyro-magnetically induced transparency and opaqueness at microwave frequency
Hanin et al. Performance of patch antenna loaded of ferrites films

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