CN112765735B - Optimization method for suspension parameters of virtual rail train - Google Patents
Optimization method for suspension parameters of virtual rail train Download PDFInfo
- Publication number
- CN112765735B CN112765735B CN202110115978.2A CN202110115978A CN112765735B CN 112765735 B CN112765735 B CN 112765735B CN 202110115978 A CN202110115978 A CN 202110115978A CN 112765735 B CN112765735 B CN 112765735B
- Authority
- CN
- China
- Prior art keywords
- rail train
- virtual rail
- suspension
- parameters
- target virtual
- 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.)
- Expired - Fee Related
Links
- 239000000725 suspension Substances 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000005457 optimization Methods 0.000 title claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 230000007246 mechanism Effects 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 230000002068 genetic effect Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 230000006835 compression Effects 0.000 claims description 2
- 238000007906 compression Methods 0.000 claims description 2
- 238000012937 correction Methods 0.000 claims description 2
- 238000013016 damping Methods 0.000 claims description 2
- 230000009347 mechanical transmission Effects 0.000 claims description 2
- 230000003595 spectral effect Effects 0.000 claims description 2
- 238000001228 spectrum Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000003915 air pollution Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008846 dynamic interplay Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000007474 system interaction Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
The invention discloses an optimization method of virtual rail train suspension parameters, which comprises the steps of obtaining structural parameters of a target virtual rail train and the irregularity grade of a road, and carrying out dynamic modeling on the target virtual rail train to obtain a dynamic model of the target virtual rail train; according to the rough degree grade of the road, establishing a rough degree condition of the road corresponding to the target virtual rail train; carrying out parameterized modeling on the suspension of the target virtual rail train according to the dynamic model and the irregularity condition to obtain a parameterized model of the suspension; carrying out optimization calculation on a parameterized model of the suspension by taking three objective functions of a train stability index, a comfort index and a tire wear index as objects to obtain optimized suspension parameters; the method can solve the problem that a method for evaluating and optimizing the suspension parameters is lacked in the prior art, and has the advantages of comprehensive coverage influence factors, reliable calculation and wide application range.
Description
Technical Field
The invention relates to the technical field of rail trains, in particular to a method for optimizing suspension parameters of a virtual rail train.
Background
In recent years, with the rapid development of urban traffic, many different types of urban traffic systems have been developed to solve the problems of severe traffic congestion, energy shortage, and air pollution. Compared with the traditional railway, the articulated train has the advantages of energy conservation, environmental protection, small turning radius, low running noise, large passenger capacity and the like. The method can be used as a supplement of a large-scale urban rail transit system and can also be used as a main traffic mode of medium and small cities. In recent years, china has proposed a new type of transportation system, named virtual rail train system, which has the ability to track automatically along expected lanes. The device consists of three vehicle groups, and rubber tires and hub motors are adopted, so that the device can run on the existing roads in cities. The virtual rail train has the characteristics of complex structure, strong maneuverability and small curve passing performance.
In order to optimize the structural design of the virtual rail train public transportation system and ensure the safety and reliability of the operation of the virtual rail train public transportation system, the evaluation of the dynamic performance of the virtual rail train public transportation system is necessary. At present, the research on the dynamic performance evaluation of the virtual rail train is very rare, and the research on the virtual rail train is on the design of a controller based on a simplified dynamic model, neglects certain vibration characteristics of a vehicle system and dynamic interaction among parts of a mechanical system, and does not consider the effect of road surface irregularity.
And the dynamic performance of the virtual rail train is greatly influenced by a suspension system, particularly the key parameters of the suspension, and the selection of the key parameters of the suspension is reasonable or not, so that the dynamic performance of the virtual rail train directly has great influence on the motion stability, safety and comfort of the train.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides the optimization method of the suspension parameters of the virtual rail train, which can solve the problem that the method for evaluating and optimizing the suspension parameters is lacked in the prior art.
In order to solve the technical problem, the invention adopts the following technical scheme:
the method for optimizing the suspension parameters of the virtual rail train comprises the following steps:
s1, acquiring structural parameters of the target virtual rail train and the irregularity grade of the road, and performing dynamic modeling on the target virtual rail train to obtain a dynamic model of the target virtual rail train;
s2, according to the rough level of the road, establishing rough conditions of the road corresponding to the target virtual rail train;
s3, carrying out parameterized modeling on the suspension of the target virtual rail train according to the dynamic model and the irregularity condition to obtain a parameterized model of the suspension;
and S4, carrying out optimization calculation on the parameterized model of the suspension by taking three objective functions of a train stability index, a comfort index and a tire wear index as objects to obtain optimized suspension parameters.
The method for optimizing the suspension parameters of the virtual rail train provided by the invention has the main beneficial effects that:
the invention realizes the parameterization establishment of the suspension of the virtual rail train to a space dynamics model by combining the interaction among the suspension system, the steering system and the road, fully ensures the full consideration of the interaction and the connection among all the structures, and ensures the accurate simulation of the suspension performance, thereby ensuring the reliability of the optimization result.
Drawings
Fig. 1 is a flow chart of a method for optimizing suspension parameters of a virtual rail train according to the present invention.
Fig. 2 is a model diagram of a spatial structure of a virtual rail train.
FIG. 3 is a diagram showing a comparison between road unevenness and a simulation structure.
Detailed Description
The invention will be further described with reference to the accompanying drawings in which:
fig. 1 is a flowchart of a method for optimizing suspension parameters of a virtual rail train according to the present invention.
The method for optimizing the suspension parameters of the virtual rail train comprises the following steps of:
s1, acquiring structural parameters of the target virtual rail train and the irregularity grade of the road, and performing dynamic modeling on the target virtual rail train to obtain a dynamic model of the target virtual rail train;
further, the virtual rail train steering systems in the scheme all adopt a disconnected steering trapezoidal mechanism driven by a gear rack, and the subsequent description is also carried out based on the structure.
Wherein, as shown in fig. 2, the structural parameter includes the length L of the knuckle arm in the steering mechanism3Length L of pull rod2And the distance L between adjacent end points of the two side pull rods1And also includes the relative angles between the components.
The grade of the unevenness of the road can be divided into five grades from good to bad.
Further, the method for obtaining the target virtual rail train dynamics model comprises the following steps:
s1-1, changing the structural parameters of the steering mechanism of the target virtual rail train into coordinated parameters;
specifically, the coordinate parameter is a planar two-dimensional coordinate converted based on the aforementioned steering mechanism size and relative angle, and includes:
and S1-2, calculating the steering angle of the target virtual rail train according to the coordinate parameters.
Further, it comprises the following steps:
s1-2-1, when the rack in the gear rack mechanism is transversely displaced by y in steering, B1The new coordinates of (a) are:
wherein the content of the first and second substances,andrespectively is changed B1Point abscissa and ordinate;
s1-2-2, according to B1Calculating A according to the new coordinate and mechanical transmission relation1New coordinates of (2):
wherein the content of the first and second substances,andrespectively is A after change1、C1Point abscissa and ordinate;
s1-2-3, the steering angle alpha of the left wheel is as follows:
s1-2-4, repeating the steps S1-2-1 to S1-2-3 to obtain the steering angle beta of the right wheel:
the steering angles of the left and right wheels are constrained and determined by the formulas of S1-2-3 and S1-2-4 by the lateral displacement of the rack. I.e., given a lateral rack displacement, the steering angle of the left and right wheels is dictated by the steering system.
And S1-3, calculating the steering constraint condition of the target virtual rail train according to the steering angle and the structure parameters.
Specifically, the steering constraint conditions are as follows:
in the formula, kt、ctVertical stiffness and damping, respectively, of the wheel, ztIs the vertical displacement of the tyre, zr0The unevenness of the road surface is the unevenness of the road surface,the radial displacement change rate of the tire is shown, and delta r is the radial compression of the tire; kappaxwAnd alphaywRespectively the amount and angle of slip of the wheel, BxAnd ByLongitudinal and transverse stiffness coefficients, CxAnd CyForm factors in the longitudinal and transverse directions, respectively, DxAnd DyPeak values in the longitudinal and transverse directions, respectively, ExAnd EyLongitudinal and transverse curvature factors, respectively.
In addition, the left suspension system, the right suspension system and the wheels are independent in movement and stress, so that the left suspension system, the right suspension system and the wheels are combined together in a model to effectively evaluate the dynamic characteristics of the train.
S2, establishing an irregularity condition of the road corresponding to the target virtual rail train according to the irregularity grade of the road;
further, the irregularity condition is:
s2-1, calculating the irregularity z of the unilateral road according to the irregularity1:
Where G (n) is the power spectral density function of spatial frequency n, which is:
in the formula, phiiTo represent the random number of the phase angle, let-nminAnd-n abovemaxThe frequency limit is 0.01m-1And 10m-1;
Since the virtual rail trains actually travel on two roads parallel to each other, the irregularity of the roads on both sides needs to be calculated separately.
S2-2, calculating the irregularity z of the road on the other side2:
Wherein beta isiIs a random phase angle, G, between 0 and 2 pixAnd (n) is a cross-power spectrum, A, B, C, D, E is a fitting coefficient of a polynomial, and all parameters are empirical parameters.
The scheme adopts different parameters of different levels of road irregularity, as shown in the following table 1:
TABLE 1 parameters of different road irregularities
Based on the formulas of S2-1 and S2-2, road irregularity on the left and right sides can be simulated. FIG. 3 is a class A road irregularity. It can be seen that the road irregularity on the left and right sides and the irregularity obtained by modeling keep good consistency, thereby verifying the effectiveness of the model.
And S3, carrying out parameterized modeling on the suspension of the target virtual rail train according to the dynamic model and the irregularity condition to obtain a parameterized model of the suspension.
Preferably, the suspension model is obtained by calling a subvar file in Simpack software through matlab software and substituting the adjusted structural parameters into the subvar file.
The subvar file is provided with a suspension model, and the requirements can be met only by adjusting corresponding parameters of the suspension model.
And S4, carrying out optimization calculation on the parameterized model of the suspension by taking three objective functions of the train stability index, the comfort index and the tire wear index as objects to obtain optimized suspension parameters.
The optimization calculation is obtained by NSGA-II multi-target genetic algorithm.
Further, the stationarity index W is:
wherein A represents a vibration acceleration, f represents a vibration frequency, and F (f) represents a frequency correction coefficient;
comfort index NMVComprises the following steps:
in the formula, ax、ay、azRespectively in the longitudinal and transverse directionsAnd effective value of vertical acceleration, WdRespectively representing the weight coefficients of the acceleration, XP representing the coverage area of the bottom surface of the vehicle, and 95 representing the distribution probability;
the tire wear index WI is:
where μ is the tire-ground friction coefficient, TxAnd TyRespectively representing longitudinal and transverse slip forces, vxAnd vyRespectively, the longitudinal and lateral slip rates, and S the average area of contact of each tire with the ground.
Optimizing suspension parameter variables of the virtual rail train by a multi-objective genetic algorithm NSGA-II, obtaining a series of pareto optimal solutions by multi-objective optimization, optimizing design variables and optimization targets tend to converge through iterative calculation of the optimization algorithm, optimizing suspension parameters, improving performance indexes of the virtual rail train, and finally selecting one optimized suspension parameter by a minimum distance method.
The above description is for the purpose of describing particular embodiments of the present invention in order to facilitate the understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the particular embodiments, and it will be apparent to those of ordinary skill in the art that various changes may be made therein without departing from the spirit and scope of the present invention as defined and defined in the appended claims, and all such changes are intended to be protected by the present invention.
Claims (8)
1. A method for optimizing suspension parameters of a virtual rail train is characterized by comprising the following steps:
s1, obtaining the structural parameters of the target virtual rail train and the rough grade of the road, wherein the structural parameters comprise the length L of a steering knuckle arm in a steering mechanism3Length L of pull rod2And the distance L between adjacent end points of the two side pull rods1Performing dynamic modeling on the target virtual rail train to obtain a target virtual railA dynamic model of a railroad train, comprising:
s1-1, changing the structural parameters of the steering mechanism of the target virtual rail train into coordinate parameters;
s1-2, calculating the steering angle of the target virtual rail train according to the coordinated parameters;
s1-3, calculating a steering constraint condition of the target virtual rail train according to the steering angle and the structural parameters;
s2, establishing an irregularity condition of the road corresponding to the target virtual rail train according to the irregularity grade of the road;
s3, carrying out parameterized modeling on the suspension of the target virtual rail train according to the dynamic model and the irregularity condition to obtain a parameterized model of the suspension;
and S4, carrying out optimization calculation on the parameterized model of the suspension by taking three objective functions of a train stability index, a comfort index and a tire wear index as objects to obtain optimized suspension parameters.
2. The method of optimizing virtual rail train suspension parameters of claim 1, wherein the coordinated parameters comprise:
3. the method for optimizing suspension parameters of a virtual rail train as claimed in claim 2, wherein the calculation method of the steering angle is:
s1-2-1, when the rack in the gear rack mechanism is transversely displaced by y in steering, B1The new coordinates of (c) are:
wherein the content of the first and second substances,andrespectively is changed B1Point abscissa and ordinate;
s1-2-2, according to B1New coordinate and mechanical transmission relation of (2), calculating A1New coordinates of (2):
wherein, the first and the second end of the pipe are connected with each other,andrespectively is A after change1、C1Point abscissa and ordinate;
s1-2-3, the steering angle alpha of the left wheel is as follows:
s1-2-4, repeating the steps S1-2-1 to S1-2-3 to obtain the steering angle beta of the right wheel:
4. the method of optimizing virtual rail train suspension parameters of claim 3, wherein the steering constraints are:
in the formula, kt、ctVertical stiffness and damping, respectively, of the wheel, ztIs the vertical displacement of the tyre, zr0In order to obtain the unevenness of the road surface,the tire radial displacement change rate is shown, and delta r is tire radial compression; kappaxwAnd alphaywRespectively the amount and angle of slip of the wheel, BxAnd ByLongitudinal and transverse stiffness coefficients, CxAnd CyForm factor, D, in the longitudinal and transverse directions, respectivelyxAnd DyPeak values in the longitudinal and transverse directions, respectively, ExAnd EyLongitudinal and transverse curvature factors, respectively.
5. The method of optimizing virtual rail train suspension parameters of claim 4, wherein the out-of-flatness condition is:
s2-1, calculating the irregularity z of the unilateral road according to the irregularity1:
Where G (n) is a power spectral density function of spatial frequency n, which is:
in the formula, phiiTo represent the random number of the phase angle, let-nminAnd-n abovemaxThe frequency limit is 0.01m-1And 10m-1;
S2-2, calculating the irregularity z of the road on the other side2:
Wherein beta isiIs a random phase angle, G, between 0 and 2 pix(n) is the cross-power spectrum and A, B, C, D, E is the fitting coefficient of the polynomial.
6. The method for optimizing the suspension parameters of the virtual rail train according to claim 5, wherein the suspension model is obtained by calling a subvar file in Simpack software through matlab software and substituting the adjusted structural parameters into the subvar file.
7. The method according to claim 6, wherein the optimization calculation is obtained by NSGA-I multi-objective genetic algorithm.
8. The method of optimizing virtual rail train suspension parameters of claim 7, wherein the stationarity index W is:
wherein A represents a vibration acceleration, f represents a vibration frequency, and F (f) represents a frequency correction coefficient;
comfort index NMVComprises the following steps:
in the formula, ax、ay、azRespectively representing effective values of longitudinal, transverse and vertical acceleration, WdRespectively representing the weight coefficients of the acceleration, XP representing the coverage area of the bottom surface of the vehicle, and 95 representing the distribution probability;
the tire wear index WI is:
where μ is the tire-ground friction coefficient, TxAnd TyRespectively representing longitudinal and transverse slip forces, vxAnd vyRespectively, the longitudinal and lateral slip rates, and S the contact area of each tire with the ground.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110115978.2A CN112765735B (en) | 2021-01-28 | 2021-01-28 | Optimization method for suspension parameters of virtual rail train |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110115978.2A CN112765735B (en) | 2021-01-28 | 2021-01-28 | Optimization method for suspension parameters of virtual rail train |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112765735A CN112765735A (en) | 2021-05-07 |
CN112765735B true CN112765735B (en) | 2022-06-17 |
Family
ID=75706341
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110115978.2A Expired - Fee Related CN112765735B (en) | 2021-01-28 | 2021-01-28 | Optimization method for suspension parameters of virtual rail train |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112765735B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113910856B (en) * | 2021-11-05 | 2023-01-03 | 福州大学 | Vehicle active suspension adjusting system considering three-dimensional linear shape and unevenness of road curved surface |
CN115577448A (en) * | 2022-09-30 | 2023-01-06 | 中车长春轨道客车股份有限公司 | Screening method, device, equipment and medium for suspension parameters of railway vehicle |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015053782A1 (en) * | 2013-10-11 | 2015-04-16 | Halliburton Energy Services Inc. | Control of drill path using smoothing |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8457841B2 (en) * | 2007-08-31 | 2013-06-04 | GM Global Technology Operations LLC | Suspension system with optimized damping response |
US9256791B2 (en) * | 2012-12-04 | 2016-02-09 | Mobileye Vision Technologies Ltd. | Road vertical contour detection |
CN106250592A (en) * | 2016-07-21 | 2016-12-21 | 江苏大学 | A kind of vehicle inertial matter optimization of suspension parameters method for designing |
CN111523173B (en) * | 2020-03-05 | 2022-12-20 | 东南大学 | Virtual rail train hub motor driving independent suspension mechanism and design method thereof |
-
2021
- 2021-01-28 CN CN202110115978.2A patent/CN112765735B/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015053782A1 (en) * | 2013-10-11 | 2015-04-16 | Halliburton Energy Services Inc. | Control of drill path using smoothing |
Non-Patent Citations (1)
Title |
---|
基于双管液压减振器控制的车辆悬架振动研究;闫云敬等;《中国工程机械学报》;20200815(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112765735A (en) | 2021-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112765735B (en) | Optimization method for suspension parameters of virtual rail train | |
KR102331677B1 (en) | Inertia adjustment active suspension system and control method based on vehicle attitude deviation | |
CN102729760B (en) | Real-time optimal damping control algorithm of automobile semi-active suspension system | |
CN107391911A (en) | The determination method and device of the dynamic equivalent taper of rail vehicle | |
Yu et al. | Optimal design and control of 4-IWD electric vehicles based on a 14-DOF vehicle model | |
CN114491821B (en) | Linear parameter optimization method for flat longitudinal section of suspension type monorail line | |
CN103942392A (en) | Automotive chassis technical parameter robust design method based on full life circle | |
CN101419724A (en) | Three-dimensional pavement generating method applying to virtual prototype system | |
CN103921786A (en) | Nonlinear model prediction control method of regenerative braking of electric vehicle | |
CN111523173B (en) | Virtual rail train hub motor driving independent suspension mechanism and design method thereof | |
CN112836296B (en) | Modeling method of virtual rail train space dynamics model | |
Gong et al. | Robust control strategy of heavy vehicle active suspension based on road level estimation | |
CN110598275A (en) | Wheel profile optimization method based on response surface modeling and improved particle swarm optimization | |
CN105117546B (en) | City expressway Alignment Design optimization method under icy conditions | |
CN112364426B (en) | Railway pier damage assessment method and system based on driving safety and dynamic response | |
Munoz et al. | Conceptual design of a hybrid electric off-road vehicle | |
Wang et al. | Study on rail profile optimization based on the nonlinear relationship between profile and wear rate | |
Zheng et al. | Evaluation on braking stability of autonomous vehicles running along curved sections based on asphalt pavement adhesion properties | |
CN112987574B (en) | Control method of cloud control intelligent chassis system based on multiple intelligent agents | |
CN113761658A (en) | Optimized matching method for stiffness damping characteristics of vehicle air suspension | |
Liu | Research on the Damage of Heavy Vehicles to the Pavement | |
CN116882311B (en) | Computational fluid dynamics determination method for normalized working attack angle of lift wing of high-speed train | |
Wu et al. | Dynamic Performance Optimization of Monorail Rapid Transit Formation Vehicles Based on Improved Artificial Potential Field Algorithm | |
Anyakwo et al. | Modelling rail vehicle dynamics using a novel 2D wheel-rail contact model | |
CN113246749B (en) | Multi-wheel/multi-shaft independent driving electric automobile multi-target dynamics control method |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220617 |