CN110346719B - Rotor eccentricity analysis method based on accurate angle-preserving mapping - Google Patents
Rotor eccentricity analysis method based on accurate angle-preserving mapping Download PDFInfo
- Publication number
- CN110346719B CN110346719B CN201910554843.9A CN201910554843A CN110346719B CN 110346719 B CN110346719 B CN 110346719B CN 201910554843 A CN201910554843 A CN 201910554843A CN 110346719 B CN110346719 B CN 110346719B
- Authority
- CN
- China
- Prior art keywords
- air gap
- mapping
- gap region
- rotor
- motor
- 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
Links
- 238000013507 mapping Methods 0.000 title claims abstract description 64
- 238000004458 analytical method Methods 0.000 title claims abstract description 15
- 238000009826 distribution Methods 0.000 claims abstract description 29
- 230000006698 induction Effects 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 230000009466 transformation Effects 0.000 claims abstract description 5
- 238000004804 winding Methods 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 4
- 230000005672 electromagnetic field Effects 0.000 abstract description 6
- 230000003068 static effect Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/34—Testing dynamo-electric machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Manufacture Of Motors, Generators (AREA)
Abstract
The invention discloses a rotor eccentricity analysis method based on accurate angle-preserving mapping, which comprises the following steps: the method comprises the steps of mapping an air gap region of a surface-mounted motor with eccentric rotor to an annular air gap region with a tooth slot concentric with a stator and a rotor through bilinear mapping, mapping the annular air gap region with the tooth slot to a polygonal air gap region through logarithmic mapping, mapping the polygonal air gap region to a rectangular air gap region through Schwarz-Christoffel mapping, and finally mapping the rectangular air gap region to two air gap regions without concentric slots through exponential mapping. Calculating the electromagnetic field distribution of the air gap region of the non-groove concentric circle according to a black lattice equation; the air gap magnetic field distribution of the motor under the original rotor eccentricity can be obtained through the inverse transformation of the 4 conformal mappings, and then the induction voltage and the electromagnetic torque of the rotor eccentric motor are obtained. The time for calculating the dynamic eccentricity and the static eccentricity is shorter than the finite element calculation time, and the prediction result is completely consistent with the finite element result.
Description
Technical Field
The invention relates to the field of motor electromagnetic field analysis, in particular to a rotor eccentricity analysis method based on accurate angle-preserving mapping.
Background
The analysis of the air gap field of the permanent magnet motor is the basis of motor design and performance calculation, and is generally solved by adopting a numerical method or an analytical method. The finite element has the characteristics of high calculation precision and wide application range, but the calculation time is long, so that the calculation period is prolonged. The analytic method has the characteristic of high calculation speed, but the precision is relatively low compared with a finite element. In the initial optimization design process of the motor, modeling optimization is generally performed by using an analytical method. After a proper basic structure of the motor stator and the motor rotor is obtained, a finite element method is adopted for further optimization design and result verification.
The method for solving the motor electromagnetic field by using the analytic method mainly comprises a subfield method (typically shown in patent CN105005692A), a magnetic circuit method (typically shown in patent CN103823926A) and the like, but the methods cannot consider the rotor eccentricity of the motor. Due to inevitable influence factors in actual production and manufacturing, such as manufacturing process, structural deformation, bearing errors and the like, the rotor of the motor is often eccentric to the stator, so that the performance of the motor is changed. In order to accurately predict the electromagnetic performance of the rotor eccentric motor, the invention provides a rotor eccentricity analysis method based on accurate angle-preserving mapping, which can convert an eccentric motor into a concentric motor through angle-preserving mapping, thereby considering the influence of the rotor eccentricity and having the characteristics of high precision and high calculation speed.
Disclosure of Invention
The invention mainly solves the technical problem that the analytic analysis of the electromagnetic field under the eccentricity of the motor rotor is difficult.
In order to solve the technical problems, the invention adopts a technical scheme that:
a rotor eccentricity analysis method based on accurate conformal mapping comprises the following steps:
1. mapping an air gap area of the rotor eccentric motor to an annular air gap area with tooth grooves, which is concentric with the stator and the rotor, through bilinear mapping;
2. mapping the concentric annular air gap areas with tooth grooves of the stator and the rotor to polygonal air gap areas through logarithmic mapping;
3. mapping the polygonal airgap region to the rectangular airgap region by Schwarz-Christoffel mapping;
4. mapping the rectangular air gap area to two non-groove concentric circle air gap areas through exponential mapping, and calculating the magnetic induction intensity distribution of the non-groove concentric circle air gap areas according to a black grid equation;
5. and (4) obtaining the air gap magnetic field distribution of the rotor eccentric motor by utilizing the magnetic induction intensity distribution of the air gap area of the non-groove concentric circle obtained in the step (4) through the inverse transformation of the angle keeping mapping in the step (1-4), and further obtaining the induced voltage and the electromagnetic torque of the rotor eccentric motor.
Further, the formula of the bilinear mapping in step 1 is as follows:
d=c0ejβ
wherein S represents the air gap region of the eccentric motor, T represents the annular air gap region, and the stator inner diameter of the motor is RsThe outer diameter of the rotor is Rr,ρdRepresents OsOrDistance, beta, represents OsOrAngle to the horizontal axis.
Further, the formula of the logarithmic mapping in step 2 is Z ═ log (t), where Z represents a polygonal air gap region.
Further, the formula of the Schwarz-christofel mapping in step 3 is as follows:
where W represents the rectangular air gap area,Z0,Z1,w0,w,wk,αjis a parameter of the Schwarz-christofel mapping, derived from the values of the points of the polygon on the polygonal airgap region Z.
Further, the formula of the exponential mapping in step 4 is as follows:
where Ψ represents the circular air gap region, Δ x and Δ y are the length and width of the rectangular air gap region W, and j is an imaginary unit;
the black lattice equation is as follows:
wherein B isrψkIs the internal polar coordinate (r) of the air gap region psi of the circular ringψ,αψ) Distribution of radial magnetic induction in position, BtψkIs a tangential magnetic induction intensity distribution on polar coordinates in a psi annular air gap regionrψAnd RsψIs the inner and outer diameters, mu, of the annular air gap region psi0Is a vacuum permeability, IkThe magnitude of the current in polar coordinates in the circular air gap region Ψ.
Further, the step 5 specifically comprises the following steps:
the inverse of the conformal mapping is to map the radial direction at the annular air gap region ΨMagnetic induction intensity distribution BrψkAnd tangential magnetic induction distribution BtψkAir-gap field distribution B mapped to air-gap region S of eccentric motorskI.e. by
Bsk(r,α)=Bsrk(r,α)+jBstk(r,α)=λ*(Brψk(r,α)+jBtψk(r,α))
according to the air-gap magnetic field distribution B of the rotor eccentric motorskObtaining the induced voltage E of the rotor eccentric motorpAnd an electromagnetic torque T, the calculation formula being as follows:
wherein N iscFor the number of turns per slot of the motor, lefIs the effective length of the motor, alphaiIs the initial angle of the coil winding, τ is the winding span, and r is the radius of any circle within the air gap region S of the eccentric machine.
The invention has the beneficial effects that:
1. the electromagnetic performance of the eccentric motor is directly calculated by an analytical method, and the calculation time under the condition of calculating the dynamic eccentricity and the static eccentricity is shorter than that of a finite element;
2. the predicted result is completely consistent with the result calculated by the finite element.
Drawings
FIG. 1 is a diagram of a surface-mounted permanent magnet machine with rotor eccentricity in the S domain;
FIG. 2 is a diagram of the concentric air gap region of the stator and rotor in the T field;
FIG. 3 is a diagram of a polygonal air gap region in the Z domain;
FIG. 4 is a diagram of a rectangular air gap region under the W domain;
FIG. 5 is a diagram of a slotless concentric circular air gap region under the Ψ -domain;
fig. 6 is an induced voltage waveform of the rotor eccentric motor;
fig. 7 is an electromagnetic torque waveform of the rotor eccentric motor.
Detailed Description
The following detailed description of the preferred embodiments of the present invention is provided to enable those skilled in the art to more readily understand the advantages and features of the present invention, and to clearly and unequivocally define the scope of the present invention.
The embodiment of the invention is as follows:
step one, for the eccentric motor shown in figure 1, the center of a circle of a rotor is OrThe center of the stator is OsThe center of the rotating shaft is OaThe stator of the motor has an inner diameter of RsThe outer diameter of the rotor is Rr. When O is presentaAnd OrCoincident with but not bound to OsWhen the motor is overlapped, the motor is in static eccentricity; when O is presentaAnd OsCoincident with but not bound to OrWhen the motor is overlapped, the motor is in dynamic eccentricity; when O is presentaIs not reacted with OrCoincident with O and notsWhen the two phases are coincident, the motor is in mixed eccentricity. Mapping a surface-mounted motor air gap region S with eccentric rotor to an annular air gap region T with a tooth slot and concentric with a stator and a rotor through bilinear mapping, as shown in FIG. 2;
d=c0ejβ
wherein S represents the air gap region of the eccentric motor, T represents the annular air gap region, and the stator inner diameter of the motor is RsThe outer diameter of the rotor is Rr,ρdRepresents OsOrDistance, beta, represents OsOrAngle to the horizontal axis.
Step two, utilizing the annular air gap area T with the tooth grooves, which is obtained in the step one, of the stator and the rotor, and mapping the annular air gap area T with the tooth grooves, which is obtained in the step one, of the stator and the rotor to a polygonal air gap area Z through logarithmic mapping, as shown in FIG. 3;
Z=log(T)
step three, mapping the polygonal air gap region Z obtained in the step two to a rectangular air gap region W through Schwarz-Christoffel mapping by utilizing the polygonal air gap region Z, as shown in FIG. 4;
wherein: z0,Z1,w0,w,wk,αjThe parameters are parameters mapped by Schwarz-Christoffel, and are obtained according to the values of the points of the polygon on the complex domain;
step four, using the rectangular air gap region W obtained in step three, mapping the rectangular air gap region W to two concentric circular air gap regions psi without grooves by exponential mapping, as shown in fig. 5,
wherein: Δ x and Δ y are the length and width of the rectangular air gap region W, and j is an imaginary unit;
the electromagnetic field distribution of the air gap region psi of concentric circles without grooves can be calculated according to the equation of black grids:
wherein B isrψkIs the internal polar coordinate (r) of the air gap region psi of the circular ringψ,αψ) Distribution of radial magnetic induction in position, BtψkIs a tangential magnetic induction intensity distribution on polar coordinates in a psi annular air gap regionrψAnd RsψIs the inner and outer diameters, mu, of the annular air gap region psi0Is a vacuum permeability, IkThe magnitude of the current in polar coordinates in the circular air gap region Ψ.
Step five, utilizing the electromagnetic field distribution of the non-groove concentric circle region psi obtained in the step four, and obtaining the air gap magnetic field distribution of the motor under the original rotor eccentricity through inverse transformation of 4 conformal mappings from the step one to the step four, wherein the inverse transformation of the conformal mapping is to obtain the radial magnetic induction intensity distribution B under the annular air gap region psirψkAnd tangential magnetic induction distribution BtψkAir-gap field distribution B mapped to air-gap region S of eccentric motorskI.e. by
Bsk(r,α)=Bsrk(r,α)+jBstk(r,α)=λ*(Brψk(r,α)+jBtψk(r,α))
according to the air-gap magnetic field distribution B of the rotor eccentric motorskObtaining the induced voltage E of the rotor eccentric motorpAnd an electromagnetic torque T, the calculation formula being as follows:
wherein N iscFor the number of turns per slot of the motor, lefIs the effective length of the motor, alphaiIs the initial angle of the coil winding, τ is the winding span, and r is the radius of any circle within the air gap region S of the eccentric machine.
As shown in fig. 6-7, fig. 6 is the waveform of the induced voltage of the rotor eccentric motor, and fig. 7 is the waveform of the electromagnetic torque of the rotor eccentric motor, it can be seen that the induced voltage and the electromagnetic torque calculated by the method are very consistent with the calculation result of the finite element, and the calculation time by the method is shorter than that by the finite element method.
The above description is only an embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes performed by the present specification and drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.
Claims (5)
1. A rotor eccentricity analysis method based on accurate conformal mapping is characterized by comprising the following steps:
1) mapping an air gap area of the rotor eccentric motor to an annular air gap area with tooth grooves, which is concentric with the stator and the rotor, through bilinear mapping; step 1) the bilinear mapping formula is as follows:
d=c0ejβ
wherein S represents the air gap region of the eccentric motor, T represents the annular air gap region, and the stator inner diameter of the motor is RsThe outer diameter of the rotor is Rr,ρdRepresents OsOrDistance, beta, represents OsOrThe included angle between the horizontal shaft and the horizontal shaft;
2) mapping the concentric annular air gap areas with tooth grooves of the stator and the rotor to polygonal air gap areas through logarithmic mapping;
3) mapping the polygonal airgap region to the rectangular airgap region by Schwarz-Christoffel mapping;
4) mapping the rectangular air gap region to a circular air gap region through exponential mapping, and calculating the magnetic induction intensity distribution of the circular air gap region according to a black grid equation;
5) and (3) obtaining the air gap magnetic field distribution of the rotor eccentric motor by utilizing the magnetic induction intensity distribution of the annular air gap area obtained in the step 4) through inverse transformation of the angle-preserving mapping from the step 1) to the step 4), and further obtaining the induction voltage and the electromagnetic torque of the rotor eccentric motor.
2. The method for analyzing the eccentricity of the rotor based on the precise conformal mapping of claim 1, wherein the logarithmic mapping of step 2) has a formula of Z ═ log (T), wherein Z represents a polygonal air gap region.
3. The method for analyzing the rotor eccentricity based on the precise conformal mapping of claim 1, wherein the formula of the Schwarz-Christoffel mapping in step 3) is as follows:
wherein W represents a rectangular air gap region, Z0,Z1,w0,w,wk,αjIs a parameter of the Schwarz-christofel mapping, derived from the values of the points of the polygon on the polygonal airgap region Z.
4. The method for analyzing the eccentricity of a rotor based on a precise conformal mapping according to claim 1, wherein the formula of the exponential mapping of step 4) is as follows:
where Ψ represents the circular air gap region, Δ x and Δ y are the length and width of the rectangular air gap region W, and j is an imaginary unit;
the black lattice equation is as follows:
wherein B isrψkIs the internal polar coordinate (r) of the air gap region psi of the circular ringψ,αψ) Distribution of radial magnetic induction in position, BtψkIs a tangential magnetic induction intensity distribution on polar coordinates in a psi annular air gap regionrψAnd RsψIs the inner and outer diameters, mu, of the annular air gap region psi0Is a vacuum permeability, IkThe magnitude of the current in polar coordinates in the circular air gap region Ψ.
5. The method for analyzing the rotor eccentricity based on the precise conformal mapping according to claim 1, wherein the step 5) is as follows:
the inverse of said conformal mapping isThe radial magnetic induction intensity distribution B under the annular air gap region psi is distributedrψkAnd tangential magnetic induction distribution BtψkAir-gap field distribution B mapped to air-gap region S of eccentric motorskI.e. by
Bsk(r,α)=Bsrk(r,α)+jBstk(r,α)
=λ*(Brψk(r,α)+jBtψk(r,α))
according to the air-gap magnetic field distribution B of the rotor eccentric motorskObtaining the induced voltage E of the rotor eccentric motorpAnd an electromagnetic torque T, the calculation formula being as follows:
wherein N iscFor the number of turns per slot of the motor, lefIs the effective length of the motor, alphaiIs the initial angle of the coil winding, τ is the winding span, R is the radius of any circle in the air gap region S of the eccentric machine, RsIs the stator inner diameter of the motor, mu0Is a vacuum magnetic permeability.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910554843.9A CN110346719B (en) | 2019-06-25 | 2019-06-25 | Rotor eccentricity analysis method based on accurate angle-preserving mapping |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910554843.9A CN110346719B (en) | 2019-06-25 | 2019-06-25 | Rotor eccentricity analysis method based on accurate angle-preserving mapping |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110346719A CN110346719A (en) | 2019-10-18 |
CN110346719B true CN110346719B (en) | 2020-12-01 |
Family
ID=68182985
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910554843.9A Active CN110346719B (en) | 2019-06-25 | 2019-06-25 | Rotor eccentricity analysis method based on accurate angle-preserving mapping |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110346719B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112906259B (en) * | 2021-01-25 | 2022-05-03 | 浙江大学 | Method for analyzing electromagnetic performance of permanent magnet vernier motor based on accurate conformal mapping |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005692A (en) * | 2015-07-08 | 2015-10-28 | 三峡大学 | Analytical method based permanent magnet motor field analysis and torque calculation method |
CN108111088A (en) * | 2017-12-23 | 2018-06-01 | 西安交通大学 | A kind of accurate Forecasting Methodology of permanent magnetic linear synchronous motor thrust for considering air gap fluctuation |
CN108319768A (en) * | 2018-01-23 | 2018-07-24 | 湖北西浦电机科技有限责任公司 | A kind of computer based permanent magnet motor armature magnetic reaction fields prediction technique |
CN108563912A (en) * | 2018-05-11 | 2018-09-21 | 北京理工大学 | A kind of analytic method of durface mounted permanent magnet synchronous motor air-gap field |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7421379B2 (en) * | 2004-10-29 | 2008-09-02 | Synopsys, Inc. | Segmentation and interpolation of current waveforms |
KR102012502B1 (en) * | 2017-09-26 | 2019-10-21 | 경북대학교 산학협력단 | Designing method of dielectric cavity |
-
2019
- 2019-06-25 CN CN201910554843.9A patent/CN110346719B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105005692A (en) * | 2015-07-08 | 2015-10-28 | 三峡大学 | Analytical method based permanent magnet motor field analysis and torque calculation method |
CN108111088A (en) * | 2017-12-23 | 2018-06-01 | 西安交通大学 | A kind of accurate Forecasting Methodology of permanent magnetic linear synchronous motor thrust for considering air gap fluctuation |
CN108319768A (en) * | 2018-01-23 | 2018-07-24 | 湖北西浦电机科技有限责任公司 | A kind of computer based permanent magnet motor armature magnetic reaction fields prediction technique |
CN108563912A (en) * | 2018-05-11 | 2018-09-21 | 北京理工大学 | A kind of analytic method of durface mounted permanent magnet synchronous motor air-gap field |
Non-Patent Citations (2)
Title |
---|
Magnetic Field Analysis of External Rotor Permanent-Magnet Synchronous Motors Using Conformal Mapping;Kamel Bough et al.;《IEEE TRANSACTIONS ON MAGNETICS》;20100930;第3685-3688页 * |
Study of Eccentric Round-Rotor Synchronous Machines through Conformal Mapping. Part I:Inductances and Air-Gap Field;A. Tessaro et al.;《2016 Eleventh International Conference on Ecological Vehicles and Renewable Energies》;20160523;第2页-3页第1栏,第6页第2栏-7页第1栏,图1 * |
Also Published As
Publication number | Publication date |
---|---|
CN110346719A (en) | 2019-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109600006B (en) | Solving method for electromagnetic design of surface-mounted permanent magnet motor | |
US8847460B2 (en) | Asymmetric stator teeth in an electric motor | |
CN108563912B (en) | Method for analyzing air gap field of surface-mounted permanent magnet synchronous motor | |
Alipour-Sarabi et al. | Performance analysis of concentrated wound-rotor resolver for its applications in high pole number permanent magnet motors | |
Zhu et al. | Multi-objective optimization design of a magnetic planetary geared permanent magnet brushless machine by combined design of experiments and response surface methods | |
Nasiri-Gheidari et al. | Performance evaluation of disk type variable reluctance resolvers | |
CN204179778U (en) | The motor of cogging torque can be reduced | |
Li et al. | Analytical synthesis of air-gap field distribution in permanent magnet machines with rotor eccentricity by superposition method | |
CN204258458U (en) | A kind of permanent magnet synchronous motor stator punching | |
CN110346719B (en) | Rotor eccentricity analysis method based on accurate angle-preserving mapping | |
Simón-Sempere et al. | Cogging torque cancellation by magnet shaping in surface-mounted permanent-magnet motors | |
CN112906259B (en) | Method for analyzing electromagnetic performance of permanent magnet vernier motor based on accurate conformal mapping | |
Guodong et al. | Estimation of maximum angular operation range for permanent-magnet slotted limited-angle torque motor | |
KR20050084515A (en) | Method for controlling electric motor cogging and permanent magnet electric motor | |
WO2013156363A1 (en) | Device and method for monitoring an eccentricity or a smallest air gap of a rotating electric machine | |
Škofic et al. | Parametric Study of a Permanent-Magnet Stepper Motor’s Stepping Accuracy Potential | |
JP2013021829A (en) | Manufacturing method of electric motor and electric motor | |
Ma et al. | The effect of a single‐sided pole shoe and slot on reducing torque ripple in a switched reluctance motor | |
CA2737357A1 (en) | Ac permanent magnet synchronous electrical machine | |
CN110309519B (en) | Method for improving magneto-resistive type rotary transformation precision | |
Ma et al. | Characteristic analysis and direct measurement for air gap magnetic field of external rotor permanent magnet synchronous motors in electric vehicles | |
CN208140810U (en) | A kind of detection device of armature coil phase sequence | |
He et al. | Analysis and Simulation on UMP and EMT Characters of Turbogenerator under Axial Air‐Gap Eccentricity | |
Zheng et al. | Calculation and optimization of a permanent magnetic hybrid driven 3-DOF motor | |
Li et al. | Analysis of magnetic field and levitation force characteristics for 3-DOF deflection type PM motors |
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 |