CN112329220A - Visual dynamics simulation method for fan blade structure - Google Patents

Visual dynamics simulation method for fan blade structure Download PDF

Info

Publication number
CN112329220A
CN112329220A CN202011166629.5A CN202011166629A CN112329220A CN 112329220 A CN112329220 A CN 112329220A CN 202011166629 A CN202011166629 A CN 202011166629A CN 112329220 A CN112329220 A CN 112329220A
Authority
CN
China
Prior art keywords
blade structure
excitation
simulation
response
fan blade
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.)
Pending
Application number
CN202011166629.5A
Other languages
Chinese (zh)
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 Jiaotong University
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN202011166629.5A priority Critical patent/CN112329220A/en
Publication of CN112329220A publication Critical patent/CN112329220A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D17/00Monitoring or testing of wind motors, e.g. diagnostics
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/12Geometric CAD characterised by design entry means specially adapted for CAD, e.g. graphical user interfaces [GUI] specially adapted for CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/20Configuration CAD, e.g. designing by assembling or positioning modules selected from libraries of predesigned modules

Abstract

The invention discloses a visual dynamics simulation method of a fan blade structure, which relates to the technical field of dynamics simulation and comprises the following steps: establishing a visual vibration monitoring experiment system of a fan blade structure; applying an exciting force to a fan blade structure in a video in a man-machine interaction mode; comprehensively generating dynamic response of a new blade structure under excitation; generating simulation animation of dynamic response prediction of the blade structure in real time; response prediction accuracy is calculated. By implementing the method, the dynamic response of the blade structure under different excitation forces can be simulated, other complex geometric modeling and finite element modeling are not needed, a complex model correction process is also not needed, the simulation calculation speed is high, and the efficiency is high.

Description

Visual dynamics simulation method for fan blade structure
Technical Field
The invention relates to the technical field of dynamics simulation, in particular to a visual dynamics simulation method for a fan blade structure.
Background
The blade structure of the wind driven generator is an important component of a fan, and the research on the dynamic characteristics of the blade structure is important for the working performance of the fan. The dynamic simulation of the fan blade structure is an important means for identifying structural modal parameters and analyzing the dynamic response of the structure, is beneficial to analyzing the dynamic response characteristics of the blade structure under the action of external load, and has important significance for the safe service, stability and reliability of the wind driven generator.
The traditional dynamic simulation method is generally based on finite element software, and a complex geometric model and a finite element model of a blade structure are established to analyze the dynamic characteristics of the blade. The method needs professional personnel, is complex to operate and difficult to model, and the established finite element model is very different from the result of an actual experiment, so that a complex correction process is needed to obtain an available finite element model. In addition, the finite element analysis method requires long modeling and simulation calculation time, and the kinetic analysis period is long and the efficiency is low.
The vibration measurement based on vision is an important method for monitoring the dynamic response of the fan blade structure and evaluating the health condition of the structure. A camera is adopted to shoot a structural vibration video, and structural vibration information and structural modal parameters can be extracted from video data. However, most of the existing visual vibration analysis methods for the blade structure can only monitor and identify modal parameters from existing video data, and cannot simulate and predict the dynamic response of the fan blade under a new excitation force. In fact, video monitoring data of blade structure vibration, which contains abundant static and dynamic structural information, can be used to analyze and predict the dynamic response of the blade under new excitation force, and thus, a vision-based wind turbine blade structure dynamic simulation method needs to be developed.
Therefore, technical personnel in the field are dedicated to developing a visual dynamics simulation method of a fan blade structure, the problems that when finite element analysis is adopted in the prior art, professionals are required to perform complex and tedious geometric modeling and physical modeling on a target structure, time consumption is long, efficiency is low, and a complex model correction process is required to meet experimental measurement results are solved, and the visual dynamics simulation is adopted, so that the method is convenient, fast, efficient and high in calculation speed, and the period of dynamics analysis is greatly shortened.
Disclosure of Invention
In view of the above-mentioned drawbacks of the prior art, the technical problems to be solved by the present invention are: the traditional dynamics simulation method aims at solving the problems of complex geometric modeling and finite element modeling of the fan blade structure, low simulation efficiency and complex correction caused by inconsistency of a finite element model and an experimental result.
In order to achieve the purpose, the invention provides a visual dynamics simulation method of a fan blade structure, which comprises the following steps:
step 1, establishing a visual vibration monitoring experiment system of a fan blade structure;
step 2, applying an exciting force to a fan blade structure in a video in a man-machine interaction mode;
step 3, comprehensively generating a dynamic response of a new blade structure under excitation;
step 4, generating a simulation animation of dynamic response prediction of the blade structure in real time;
and 5, calculating response prediction accuracy.
Further, the step 1 further comprises: obtaining video data of a blade structure from an experiment, identifying full-field modal parameters, and obtaining a modal shape phi and a natural frequency omegadAnd modal damping ξ.
Further, the applying of the excitation force in the human-computer interaction manner in step 2 specifically includes: manually moving the mouse to any point (x, y) on the leaf structure, pressing the left button of the mouse and simultaneously moving the mouse a certain distance in a certain direction, and then at a new position (x)e,ye) The left mouse button is released.
Further, the action point of the applied exciting force is (x, y), the applied direction of the exciting force is the moving direction of the mouse, and the magnitude of the exciting force is proportional to the moving distance of the mouse.
Further, the excitation force is:
Fu(x,y,t)=λ(xe-x)F(t)
Fv(x,y,t)=λ(ye-y)F(t)
wherein, Fu(x, y, t) and Fv(x, y, t) is the component of the applied excitation force F (x, y, t) in the horizontal and vertical directions; λ is a scale factor; f (t) isA function of the time variation, representing the type of excitation.
Further, the f (t) may be a pulsed excitation, a random excitation, or a swept excitation.
Further, the dynamic response expression in step 3 is:
Figure BDA0002746005590000021
wherein U (x, y, t) is the dynamic response of the blade structure, and comprises a horizontal displacement response U (x, y, t) and a vertical displacement response v (x, y, t); q is the truncated modal order; phi is aiIs the ith order mode shape; q. q.si(t) is the ith order modal coordinate; q. q.si(t) is:
Figure BDA0002746005590000022
wherein, (.)TAnd
Figure BDA0002746005590000023
respectively representing transpose and convolution operations; h isi(t) is a unit impulse response function, and the expression is as follows:
Figure BDA0002746005590000024
wherein m isiAnd ωniRespectively, the ith order modal mass and the undamped natural frequency.
Further, in step 4, each frame of the simulation animation is subjected to motion amplification, and the animation frame at time t is:
I(x+αu,y+αv,t)=I(x,y,0)
where α is a motion amplification factor.
Further, the step 5 further comprises: and calculating the cross-correlation coefficient of the simulation result and the experimental result by adopting a cross-correlation method.
Further, the step 5 further comprises: the response actual measurement result obtained by the experiment is G (x, y, t), the visual simulation response prediction result is U (x, y, t), and the specific expression of the cross-correlation coefficient is as follows:
Figure BDA0002746005590000031
wherein r (G, U) is the cross-correlation coefficient of the two, Cov (G, U) is the covariance of G and U, and Var [ G ] and Var [ U ] are the variance of G and U.
Compared with the prior art, the invention at least has the following beneficial technical effects:
1. the invention establishes a visual dynamics simulation model of the fan blade structure, can simulate the dynamics response of the blade structure under different excitation forces, and intuitively, conveniently and efficiently realizes the dynamics analysis of the fan blade structure;
2. establishing a dynamic visual simulation model of the fan blade structure based on the video data without other complicated geometric modeling and finite element modeling and a complicated model correction process;
3. the visual dynamics simulation calculation speed is high, the efficiency is high, and the effect of real-time simulation calculation can be achieved;
4. and the correlation degree of the simulation result and the experimental actual measurement result is accurately evaluated by adopting the cross-correlation coefficient, so that the response prediction accuracy of the visual simulation can be evaluated.
The conception, the specific structure and the technical effects of the present invention will be further described with reference to the accompanying drawings to fully understand the objects, the features and the effects of the present invention.
Drawings
FIG. 1 is a schematic view of a visual vibration monitoring experiment system of a fan blade structure of the present invention;
FIG. 2 is a first eight-order modal shape of a fan blade structure extracted from video data according to the present invention;
FIG. 3 is a schematic diagram of the excitation point A and vibration pick-up point B of the present invention;
FIG. 4 is a waveform of the response prediction at vibration pick-up point B of the present invention;
FIG. 5 is a diagram illustrating the accuracy of response prediction evaluation according to the present invention.
Detailed Description
The technical contents of the preferred embodiments of the present invention will be more clearly and easily understood by referring to the drawings attached to the specification. The present invention may be embodied in many different forms of embodiments and the scope of the invention is not limited to the embodiments set forth herein.
Step 1: as shown in fig. 1, firstly, a visual vibration monitoring experiment system of a fan blade structure is set up, video data of the blade structure is obtained from the experiment, full-field modal parameters are identified, and modal shape phi and natural frequency omega are obtaineddAnd modal damping ξ. This is the data preparation phase of the present invention. The extraction of modality parameters from video data is shown in fig. 2.
Step 2: and (4) applying exciting force by the visual dynamics simulation of the fan blade structure.
As shown in fig. 3, excitation force is applied to a fan blade structure in a video in a man-machine interaction manner. And applying an excitation force in a man-machine interaction mode, namely manually moving the mouse to any point (x, y) on the blade structure, pressing a left mouse button and simultaneously moving the mouse for a certain distance in a certain direction, then releasing the left mouse button at a new position (xe, ye), moving the mouse to a point A on the fan blade structure in the process, pressing the left mouse button and moving the left mouse button for a certain distance in the horizontal direction, and releasing the mouse to finish the application of the excitation force. And selecting the point B as a vibration pickup point for comparing the experimental result with the visual simulation result. The applied direction of the exciting force is the moving direction of the mouse, the size is in direct proportion to the moving distance, the exciting force is divided into the horizontal direction and the vertical direction, and then the specific expression of the exciting force is as follows:
Fu(x,y,t)=λ(xe-x)F(t)
Fv(x,y,t)=λ(ye-y)F(t)
wherein, Fu(x, y, t) and Fv(x, y, t) is the component of the applied excitation force F (x, y, t) in the horizontal and vertical directions; λ is a scale factor to control the magnitude of the force; f (t) is timeAnd the function of the chemical equation represents the excitation type, such as pulse excitation, random excitation, sweep excitation and the like.
And step 3: and (5) a fan blade structure visual dynamics simulation result. The response prediction waveform at point B is shown in fig. 4. And comprehensively generating the dynamic response of the blade structure under the new excitation. The dynamic response expression is:
Figure BDA0002746005590000041
wherein U (x, y, t) is the dynamic response of the blade structure, and comprises a horizontal displacement response U (x, y, t) and a vertical displacement response v (x, y, t); q is the truncated modal order; phi is aiIs the ith order mode shape; q. q.siAnd (t) is the ith order modal coordinate. q. q.si(t) can be expressed as:
Figure BDA0002746005590000042
wherein, (.)TAnd
Figure BDA0002746005590000043
respectively representing transpose and convolution operations; h isi(t) is a unit impulse response function, and the expression is as follows:
Figure BDA0002746005590000044
wherein m isiAnd ωniRespectively, the ith order modal mass and the undamped natural frequency.
After excitation is applied, the dynamic response of the fan blade structure under the excitation force can be automatically simulated and predicted, and the response result can be displayed in a video in real time. From B point vibration pickup, the consistency of the experimental result and the simulation result of blade response prediction is verified, and the fact that the response prediction waveforms of the two are very close to each other and the cross correlation coefficient is as high as 0.85 can be found.
And 4, step 4: and generating simulation animation of dynamic response prediction of the blade structure in real time. In the steps, displacement responses u (x, y, t) and v (x, y, t) of each pixel point in the horizontal and vertical directions after excitation is applied are obtained, so that simulation animation of dynamic response of the blade structure can be generated in real time, and each frame of the animation is subjected to motion amplification processing to enable tiny motion to be visualized. The animation frame at time t may be expressed as:
I(x+αu,y+αv,t)=I(x,y,0)
where α is a motion amplification factor. For amplifying micro-movements (u, v) in horizontal and vertical directions. Note that the step is mainly assignment operation, the calculation amount is small and far lower than that of the traditional complex gold tower algorithm, and therefore the real-time animation simulation can be realized.
And 5: and evaluating the accuracy of the simulation response prediction of the visual dynamics of the fan blade structure. And calculating the cross-correlation coefficient of the simulation result and the experimental result by adopting a cross-correlation method, wherein the response actual measurement result obtained by the experiment is G (x, y, t), the visual simulation response prediction result is U (x, y, t), and the specific expression of the cross-correlation coefficient is as follows:
Figure BDA0002746005590000051
wherein r (G, U) is the cross-correlation coefficient of the two, Cov (G, U) is the covariance of G and U, and Var [ G ] and Var [ U ] are the variance of G and U.
And (3) calculating the cross-correlation coefficient between the experimental measurement result of all the pixel points on the fan blade structure and the visual simulation result, as shown in fig. 5, finding that the correlation of the response prediction result is only weak at the root position of the blade because the root displacement is almost 0. In addition, the experimental actual measurement and the visual simulation response prediction results have good correlation at most positions, which shows that the response prediction accuracy of the visual dynamics simulation of the fan blade structure is high.
According to the visual dynamics simulation method for the fan blade structure, the camera is used for shooting the structural vibration video, modal parameters can be extracted from video data, the dynamic response of the blade structure under different excitation forces is simulated, other complex geometric modeling and finite element modeling are not needed, a complex model correction process is not needed, the simulation calculation speed is high, and the efficiency is high.
The foregoing detailed description of the preferred embodiments of the invention has been presented. It should be understood that numerous modifications and variations could be devised by those skilled in the art in light of the present teachings without departing from the inventive concepts. Therefore, the technical solutions available to those skilled in the art through logic analysis, reasoning and limited experiments based on the prior art according to the concept of the present invention should be within the scope of protection defined by the claims.

Claims (10)

1. A visual dynamics simulation method of a fan blade structure is characterized by comprising the following steps:
step 1, establishing a visual vibration monitoring experiment system of a fan blade structure;
step 2, applying an exciting force to a fan blade structure in a video in a man-machine interaction mode;
step 3, comprehensively generating a dynamic response of a new blade structure under excitation;
step 4, generating a simulation animation of dynamic response prediction of the blade structure in real time;
and 5, calculating response prediction accuracy.
2. The method of claim 1, wherein step 1 further comprises: obtaining video data of a blade structure from an experiment, identifying full-field modal parameters, and obtaining a modal shape phi and a natural frequency omegadAnd modal damping ξ.
3. The method of claim 1, wherein the applying the excitation force in a human-computer interaction manner in step 2 specifically comprises: manually moving the mouse to any point (x, y) on the leaf structure, pressing the left button of the mouse and simultaneously moving the mouse a certain distance in a certain direction, and then at a new position (x)e,ye) The left mouse button is released.
4. The method of claim 3, wherein the point of application of the excitation force is (x, y), the direction of application of the excitation force is the direction of mouse movement, and the magnitude of the excitation force is proportional to the distance of mouse movement.
5. The method of claim 4, wherein the energizing force is:
Fu(x,y,t)=λ(xe-x)F(t)
Fv(x,y,t)=λ(ye-y)F(t)
wherein, Fu(x, y, t) and Fv(x, y, t) is the component of the applied excitation force F (x, y, t) in the horizontal and vertical directions; λ is a scale factor; f (t) is a function of time, representing the excitation type.
6. The method of claim 5, wherein F (t) can be pulsed excitation, random excitation, or swept excitation.
7. The method of claim 1, wherein the dynamic response expression in step 3 is:
Figure FDA0002746005580000011
wherein U (x, y, t) is the dynamic response of the blade structure, and comprises a horizontal displacement response U (x, y, t) and a vertical displacement response v (x, y, t); q is the truncated modal order; phi is aiIs the ith order mode shape; q. q.si(t) is the ith order modal coordinate; q. q.si(t) is:
Figure FDA0002746005580000021
wherein, (.)TAnd
Figure FDA0002746005580000022
respectively representing transpose and convolution operations; h isi(t) is a unit impulse response function, and the expression is as follows:
Figure FDA0002746005580000023
wherein m isiAnd ωniRespectively, the ith order modal mass and the undamped natural frequency.
8. The method of claim 1, wherein in step 4, each frame of the simulation animation is subjected to motion amplification, and the animation frames at time t are:
I(x+αu,y+αv,t)=I(x,y,0)
where α is a motion amplification factor.
9. The method of claim 1, wherein step 5 further comprises: and calculating the cross-correlation coefficient of the simulation result and the experimental result by adopting a cross-correlation method.
10. The method of claim 9, wherein step 5 further comprises: the response actual measurement result obtained by the experiment is G (x, y, t), the visual simulation response prediction result is U (x, y, t), and the specific expression of the cross-correlation coefficient is as follows:
Figure FDA0002746005580000024
wherein r (G, U) is the cross-correlation coefficient of the two, Cov (G, U) is the covariance of G and U, and Var [ G ] and Var [ U ] are the variance of G and U.
CN202011166629.5A 2020-10-27 2020-10-27 Visual dynamics simulation method for fan blade structure Pending CN112329220A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011166629.5A CN112329220A (en) 2020-10-27 2020-10-27 Visual dynamics simulation method for fan blade structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011166629.5A CN112329220A (en) 2020-10-27 2020-10-27 Visual dynamics simulation method for fan blade structure

Publications (1)

Publication Number Publication Date
CN112329220A true CN112329220A (en) 2021-02-05

Family

ID=74296675

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011166629.5A Pending CN112329220A (en) 2020-10-27 2020-10-27 Visual dynamics simulation method for fan blade structure

Country Status (1)

Country Link
CN (1) CN112329220A (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160267664A1 (en) * 2015-03-11 2016-09-15 Massachusetts Institute Of Technology Methods and apparatus for modeling deformations of an object
US20160300341A1 (en) * 2014-12-11 2016-10-13 Jeffrey R. Hay Apparatus and Method for Visualizing Periodic Motions in Mechanical Components
CN107111755A (en) * 2015-01-13 2017-08-29 摩福公司 The video personation detection method and system evaluated based on liveness
US20180084195A1 (en) * 2016-09-19 2018-03-22 Los Alamos National Security, Llc System and method for automated extraction of high resolution structural dynamics from video
US20200065957A1 (en) * 2018-05-31 2020-02-27 Rdi Technologies, Inc. Monitoring of objects based on frequency spectrum of motion and frequency filtering
CN111784647A (en) * 2020-06-19 2020-10-16 西北工业大学 High-precision structural modal testing method based on video vibration amplification

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160300341A1 (en) * 2014-12-11 2016-10-13 Jeffrey R. Hay Apparatus and Method for Visualizing Periodic Motions in Mechanical Components
CN107111755A (en) * 2015-01-13 2017-08-29 摩福公司 The video personation detection method and system evaluated based on liveness
US20160267664A1 (en) * 2015-03-11 2016-09-15 Massachusetts Institute Of Technology Methods and apparatus for modeling deformations of an object
US20180084195A1 (en) * 2016-09-19 2018-03-22 Los Alamos National Security, Llc System and method for automated extraction of high resolution structural dynamics from video
US20200065957A1 (en) * 2018-05-31 2020-02-27 Rdi Technologies, Inc. Monitoring of objects based on frequency spectrum of motion and frequency filtering
CN111784647A (en) * 2020-06-19 2020-10-16 西北工业大学 High-precision structural modal testing method based on video vibration amplification

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ABE DAVIS ET AL.: "Image-Space Modal Bases for Plausible Manipulation of Objects in Video", 《ACM TRANSACTIONS ON GRAPHICS》 *
YONGCHAO YANG ET AL.: "Spatiotemporal video-domain high-fidelity simulation and realistic visualization of full-field dynamic responses of structures by a combination of high-spatial-resolution modal model and video motion manipulations", 《STRUCT CONTROL HEALTH MONIT》 *
刘万锋等: "汽车驱动桥壳动态应变场的模态分析", 《西安公路交通大学学报》 *
张茜平: "基于区域视频图像序列的微小细节增强技术研究", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 *

Similar Documents

Publication Publication Date Title
Abdeljaber et al. 1-D CNNs for structural damage detection: Verification on a structural health monitoring benchmark data
Ozbek et al. Challenges in testing and monitoring the in-operation vibration characteristics of wind turbines
JP5561787B2 (en) Method and apparatus for correlating a simulation model with a physical device based on a correlation metric
Tamura et al. Proper orthogonal decomposition of random wind pressure field
Bhowmick et al. Spatiotemporal compressive sensing of full-field Lagrangian continuous displacement response from optical flow of edge: Identification of full-field dynamic modes
CN104200236A (en) Quick target detection method based on DPM (deformable part model)
Avci et al. Efficiency validation of one dimensional convolutional neural networks for structural damage detection using a SHM benchmark data
Zhang et al. A flutter prediction method with low cost and low risk from test data
CN107729592A (en) Traced back the Time variable structure Modal Parameters Identification of track based on broad sense subspace
CN105892304B (en) Curve data automatic interpretation method
Karmakar et al. Conditional simulation of non-Gaussian wind velocity profiles: Application to buffeting response of Vincent Thomas suspension bridge
CN111488670B (en) Nonlinear mass point spring soft tissue deformation simulation method
Awadallah et al. Automated multiclass structural damage detection and quantification using augmented reality
CN116861544B (en) Building abnormal vibration source positioning method based on edge cloud cooperation and related equipment
Butler et al. Physical and numerical modeling of downburst generated gust fronts
CN112329220A (en) Visual dynamics simulation method for fan blade structure
CN105960613A (en) System identification device
CN114066108A (en) Method and device for calculating turbulence space representativeness of anemometer tower
JP2021108031A (en) Model generation device, vehicle simulation system, model generation method, vehicle simulation method and computer program
JP2021056707A (en) Vehicle simulation system, vehicle simulation method and computer program
KR101624615B1 (en) Method and apparatus for modeling dynamic muscle stiffness and method tereof
de Meneses et al. Trajectory analysis for sport and video surveillance
JP7301008B2 (en) Hammering sound inspection system and method by AI
CN103020374B (en) Water pump Flow Field Calculation real-time emulation system and its development approach based on virtual reality technology
Popescu et al. Blind source separation: A preprocessing tool for monitoring of structures

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210205

RJ01 Rejection of invention patent application after publication