CN107169217B - A kind of equivalent method of turbulent boundary layer load model - Google Patents
A kind of equivalent method of turbulent boundary layer load model Download PDFInfo
- Publication number
- CN107169217B CN107169217B CN201710378734.7A CN201710378734A CN107169217B CN 107169217 B CN107169217 B CN 107169217B CN 201710378734 A CN201710378734 A CN 201710378734A CN 107169217 B CN107169217 B CN 107169217B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- equivalent
- msup
- load model
- 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
Classifications
-
- 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
- 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)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of equivalent method of turbulent boundary layer load model, include the following steps:(1) turbulent boundary layer load model forms equivalent completely random face pressure load model after equivalent;(2) magnitude of the equivalent correlation function of the equivalent completely random face pressure load model is determined;(3) structural bending wavelength and turbulent boundary layer load characteristic wavelength are determined, and then calculates consistent resistant frequency, determines the applicable frequency range of equivalent random face pressure load model on this basis.A kind of equivalent method of turbulent boundary layer load model provided by the invention, it is a kind of technology that turbulent boundary layer load model is equivalent to completely random face pressure load model, the technology can effectively reduce the calculation amount of turbulent boundary layer load effect lower structure dynamic response analysis, shorten the design cycle, save design cost.
Description
Technical field
The present invention relates to random face pressure load model equivalent method field, and in particular to a kind of turbulent boundary layer load model
Equivalent method.
Background technology
As spacecraft develops to high flying speed, it faces the environment such as the random noise of sternness in duty cycle, this
It may cause structural failure or precision instrument, instrument malfunction.Therefore, in the design process of spacecraft, mechanical oscillation need to be considered
With the influence of noise.The sound of something astir of test method, theoretical method and numerical method forecasting system under noise excitation can be used
Should.Wherein, test method can obtain it is reliable as a result, but carry out analysis of experiments cost it is higher, design cycle length;Theoretical method
It is only applicable to single system, it is difficult to solve the problems, such as the dynamic response indication of complication system;Numerical method can save design cost, shorten
Design cycle, is effective supplementary means of analysis of experiments.
At present in a kind of generally acknowledged turbulent boundary layer load model, the coherence length of low-frequency range is longer, and high band is concerned with
Length is shorter.The random response that mode superposition method analysis structure in using FInite Element is encouraged in turbulent boundary layer load
When, with the rise of analysis frequency, the coherence length of turbulent boundary layer load shortens, it is desirable to and the size of finite element grid diminishes,
This causes calculation amount to increase by geometric progression.Therefore, in higher frequency band, need to adopt an effective measure to solve above-mentioned turbulent boundary layer
The problem of load model analysis efficiency is low, and then shorten the design cycle, save design cost.
The content of the invention
Goal of the invention:In order to overcome the deficiencies in the prior art, for a kind of existing turbulent boundary layer load mould
Type in the application there are the problem of, the present invention provides a kind of equivalent method of the turbulent boundary layer load model, which can
Effectively improve the efficiency of turbulent boundary layer load excitation lower structure dynamic response simulation analysis.
Technical solution:To achieve the above object, the technical solution adopted by the present invention is:
A kind of equivalent method of turbulent boundary layer load model, comprises the following steps:
Step 1:Turbulent boundary layer load model forms equivalent completely random face pressure load model after equivalent;
Step 2:Determine the magnitude of the equivalent correlation function of the equivalent completely random face pressure load model;
Step 3:The equivalent random face pressure load model is determined according to structural model and turbulent boundary layer load model
It is applicable in frequency range.
Further, the turbulent boundary layer load model in the step 1 is:
Wherein ξxFor 2 points of distances in the direction of the x axis, ξyFor 2 points of distances in the y-axis direction, ω is angular frequency, S0
For the magnitude of load power spectrum, Dx=αx/kc、Dy=αy/kcThe respectively coherence length of downbeam and cross-wind direction, dimensionless
Constant αx=8, αy=1.2, kc=ω/UcFor convection current wave number, Uc=0.7U is convection velocity, and U is speed of incoming flow.
Further, the equivalent completely random face pressure load model in the step 1 is:
Spp(ξx,ξy, ω) and=S0Ceq(ω)δ(ξx)δ(ξy) (2)
Wherein Ceq(ω) is the magnitude of equivalent correlation function, and function δ (ξ) is Kronecker function:
Further, in the step 2 the equivalent correlation function of equivalent completely random face pressure load model magnitude Ceq
(ω) meets following formula:
Further, the magnitude C of the equivalent correlation function of the equivalent completely random face pressure load modeleq(ω) is:
Further, the applicable frequency range of equivalent completely random face pressure load model is f >=f in the step 3crit,
fcritFor critical frequency.
Further, the critical frequency is:
fcrit=4fc (6)
Wherein fcFor consistent resistant frequency.
Further, the consistent resistant frequency fcTo make structural bending wavelength XB(ω) and reverberation field load characteristic wavelength XT
(ω) equal i.e. λB(ω)=λTConsistent resistant frequency when (ω):
Wherein, E is elasticity modulus of materials, and ρ is density of material, and ν is material Poisson's ratio, and h is thick for body structure surface plate class member
Degree.
Beneficial effect:A kind of equivalent method of turbulent boundary layer load model provided by the invention is a kind of by turbulence edge
Interlayer load model is equivalent to the technology of completely random face pressure load model, which can effectively reduce turbulent boundary layer load and swash
The calculation amount of lower structure dynamic response analysis is encouraged, shortens the design cycle, saves design cost.
Brief description of the drawings
Fig. 1 is the logical procedure diagram of the present invention;
Fig. 2 is the schematic diagram of a rectangle simply supported slab;
Fig. 3 is the dynamic respond power spectral density schematic diagram at point A on rectangle simply supported slab.
Embodiment
The present invention is further described below in conjunction with the accompanying drawings.
Be a kind of logical procedure diagram of the equivalent method of turbulent boundary layer load model as shown in Figure 1, it is main include with
Lower step:
Step (1) turbulent boundary layer load model forms equivalent completely random face pressure load model after equivalent;
(1.1) turbulent boundary layer load model, its cross-spectrum spatially at any two points between face pressure load are:
Wherein ξxFor 2 points of distances in the direction of the x axis, ξyFor 2 points of distances in the y-axis direction, ω is angular frequency, S0
For the magnitude of load power spectrum, Dx=αx/kc、Dy=αy/kcThe respectively coherence length of downbeam and cross-wind direction, dimensionless
Constant αx=8, αy=1.2, kc=ω/UcFor convection current wave number, Uc=0.7U is convection velocity, and U is speed of incoming flow.
(1.2) equivalent completely random face pressure load model, spatially the cross-spectrum of face pressure is at any two points for it:
Spp(ξx,ξy, ω) and=S0Ceq(ω)δ(ξx)δ(ξy) (2)
Wherein Ceq(ω) is the magnitude of equivalent correlation function, and function δ (ξ) is Kronecker function:
Step (2)) determine the equivalent completely random face pressure load model equivalent correlation function magnitude Ceq(ω), into
And determine the equivalent completely random face pressure load model;
The magnitude C of the equivalent correlation function of equivalent completely random face pressure load modeleq(ω) meets following formula:
Solution formula (4) obtains the magnitude C of the equivalent correlation function of equivalent completely random face pressure load modeleq(ω) is:
Step (3) determines the equivalent random face pressure load model according to structural model and turbulent boundary layer load model
It is applicable in frequency range;Specifically include:
(3.1) bending wavelength of structure is determined:
Wherein, E is elasticity modulus of materials, and ρ is density of material, and ν is material Poisson's ratio, and h is thick for body structure surface plate class member
Degree.
(3.2) characteristic wavelength of turbulent boundary layer load is determined:
λTThe π U of (ω)=2c/ω (7)
(3.3) when calculating makes structural bending wavelength and equal turbulent boundary layer load characteristic wavelength, i.e. λB(ω)=λD(ω)
When consistent resistant frequency:
(3.4) critical frequency that equivalent completely random face pressure load model is applicable in is calculated:
fcrit=4fc (9)
(3.5) the applicable frequency range for determining the equivalent completely random face pressure load model in step (2) is f >=fcrit。
Embodiment
As shown in Fig. 2, by taking a rectangle simply supported slab as an example, consistent resistant frequency is calculated.The size of rectangle simply supported slab is:X-axis
To length Lx=1m, y-axis is to length Ly=1m, thickness h=0.005m.The parameter of rectangle simply supported slab material therefor is:Elasticity modulus
E=120GPa, density of material ρ=7800kg/m3, Poisson's ratio υ=0.3.As speed of incoming flow U=100m/s, by each parameter
Value substitutes into formula (8) and obtains fc=102Hz.
The critical frequency being applicable in by the equivalent completely random face pressure load model of step (3.4) calculating is fcrit=408Hz.
Determine that the applicable frequency range of equivalent completely random face pressure load model is f >=f by step (3.5)crit, that is, work as
When analyzing frequency f >=408Hz, in this example, formula can be replaced as the equivalent completely random face pressure load model shown in formula (2)
(1) the turbulent boundary layer load model shown in.
The equivalent completely random face pressure load obtained as above-mentioned steps is put in the Simply-Supported Rectangular Plates shown in Fig. 2, is counted
Calculate obtain point A (0.3m, 0.2m) place dynamic respond power spectral density (in units of dB, reference value 1m2Hz-1), such as Fig. 3 institutes
Show.In Fig. 3 the result shows that, in this example, as f >=fcrit, i.e. during f >=408Hz, equivalent completely random face that above-mentioned steps obtain
Ballast model can effectively represent turbulent boundary layer load model.
The effect explanation that the present embodiment finally obtains, method proposed by the invention can be effectively by turbulent boundary layer load
Model conversion improves the efficiency that subsequent response is analyzed into equivalent completely random face pressure load model.
The above is only the preferred embodiment of the present invention, it should be pointed out that:For the ordinary skill people of the art
For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications also should
It is considered as protection scope of the present invention.
Claims (4)
- A kind of 1. equivalent method of turbulent boundary layer load model, it is characterised in that:Comprise the following steps:Step 1:Turbulent boundary layer load model forms equivalent completely random face pressure load model after equivalent;Step 2:Determine the magnitude of the equivalent correlation function of the equivalent completely random face pressure load model;Step 3:The equivalent completely random face pressure load model is determined according to structural model and turbulent boundary layer load model It is applicable in frequency range;Wherein:Turbulent boundary layer load model in the step 1 is:<mrow> <msub> <mi>S</mi> <mrow> <mi>p</mi> <mi>p</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>&xi;</mi> <mi>x</mi> </msub> <mo>,</mo> <msub> <mi>&xi;</mi> <mi>y</mi> </msub> <mo>,</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>S</mi> <mn>0</mn> </msub> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>|</mo> <msub> <mi>&xi;</mi> <mi>x</mi> </msub> <mo>|</mo> </mrow> <mo>/</mo> <msub> <mi>D</mi> <mi>x</mi> </msub> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>|</mo> <msub> <mi>&xi;</mi> <mi>y</mi> </msub> <mo>|</mo> </mrow> <mo>/</mo> <msub> <mi>D</mi> <mi>y</mi> </msub> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>jk</mi> <mi>c</mi> </msub> <msub> <mi>&xi;</mi> <mi>x</mi> </msub> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>Wherein ξxFor 2 points of distances in the direction of the x axis, ξyFor 2 points of distances in the y-axis direction, ω is angular frequency, S0To carry The magnitude of lotus power spectrum, Dx=αx/kc、Dy=αy/kcThe respectively coherence length of downbeam and cross-wind direction, dimensionless constant αx=8, αy=1.2, kc=ω/UcFor convection current wave number, Uc=0.7U is convection velocity, and U is speed of incoming flow;Equivalent completely random face pressure load model is described in the step 1:Spp(ξx,ξy, ω) and=S0Ceq(ω)δ(ξx)δ(ξy) (2)Wherein Ceq(ω) is the magnitude of equivalent correlation function, and function δ (ξ) is Kronecker function:The magnitude C of the equivalent correlation function of equivalent completely random face pressure load model in the step 2eq(ω) meets following formula:<mrow> <msub> <mo>&Integral;</mo> <mi>&infin;</mi> </msub> <msub> <mo>&Integral;</mo> <mi>&infin;</mi> </msub> <msup> <mi>C</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msup> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> <mi>&delta;</mi> <mrow> <mo>(</mo> <msub> <mi>&xi;</mi> <mi>x</mi> </msub> <mo>)</mo> </mrow> <mi>&delta;</mi> <mrow> <mo>(</mo> <msub> <mi>&xi;</mi> <mi>y</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>d&xi;</mi> <mi>x</mi> </msub> <msub> <mi>d&xi;</mi> <mi>y</mi> </msub> <mo>=</mo> <msub> <mo>&Integral;</mo> <mi>&infin;</mi> </msub> <msub> <mo>&Integral;</mo> <mi>&infin;</mi> </msub> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>|</mo> <msub> <mi>&xi;</mi> <mi>x</mi> </msub> <mo>|</mo> </mrow> <mo>/</mo> <msub> <mi>D</mi> <mi>x</mi> </msub> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>|</mo> <msub> <mi>&xi;</mi> <mi>y</mi> </msub> <mo>|</mo> </mrow> <mo>/</mo> <msub> <mi>D</mi> <mi>y</mi> </msub> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>jk</mi> <mi>c</mi> </msub> <msub> <mi>&xi;</mi> <mi>x</mi> </msub> </mrow> </msup> <msub> <mi>d&xi;</mi> <mi>x</mi> </msub> <msub> <mi>d&xi;</mi> <mi>y</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>The applicable frequency range of equivalent completely random face pressure load model in the step 3 is f >=fcrit, fcritFor critical frequency Rate.
- 2. the equivalent method of turbulent boundary layer load model according to claim 1, it is characterised in that:It is described equivalent complete The magnitude C of the equivalent correlation function of random face pressure load modeleq(ω) is:<mrow> <msup> <mi>C</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msup> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mn>4</mn> <msub> <mi>D</mi> <mi>x</mi> </msub> <msub> <mi>D</mi> <mi>y</mi> </msub> </mrow> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>D</mi> <mi>x</mi> </msub> <msub> <mi>k</mi> <mi>c</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
- 3. the equivalent method of turbulent boundary layer load model according to claim 1, it is characterised in that:The critical frequency For:fcrit=4fc (6)Wherein fcFor consistent resistant frequency.
- 4. the equivalent method of the random face pressure load model of reverberation field according to claim 3, it is characterised in that:It is described consistent Resistant frequency fcTo make structural bending wavelength XB(ω) and reverberation field load characteristic wavelength XT(ω) equal i.e. λB(ω)=λTWhen (ω) Consistent resistant frequency:<mrow> <msub> <mi>f</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <msubsup> <mi>U</mi> <mi>c</mi> <mn>2</mn> </msubsup> <mrow> <mn>2</mn> <mi>&pi;</mi> </mrow> </mfrac> <msqrt> <mfrac> <mrow> <mn>12</mn> <mi>&rho;</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>&upsi;</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow> <mrow> <msup> <mi>Eh</mi> <mn>2</mn> </msup> </mrow> </mfrac> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>Wherein, E is elasticity modulus of materials, and ρ is density of material, and ν is material Poisson's ratio, and h is body structure surface plate class member thickness.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710378734.7A CN107169217B (en) | 2017-05-25 | 2017-05-25 | A kind of equivalent method of turbulent boundary layer load model |
PCT/CN2018/082580 WO2018214657A1 (en) | 2017-05-25 | 2018-04-10 | Equalisation technique for turbulent boundary layer load model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710378734.7A CN107169217B (en) | 2017-05-25 | 2017-05-25 | A kind of equivalent method of turbulent boundary layer load model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107169217A CN107169217A (en) | 2017-09-15 |
CN107169217B true CN107169217B (en) | 2018-04-24 |
Family
ID=59821574
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710378734.7A Active CN107169217B (en) | 2017-05-25 | 2017-05-25 | A kind of equivalent method of turbulent boundary layer load model |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN107169217B (en) |
WO (1) | WO2018214657A1 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107169217B (en) * | 2017-05-25 | 2018-04-24 | 东南大学 | A kind of equivalent method of turbulent boundary layer load model |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106055798A (en) * | 2016-06-02 | 2016-10-26 | 东南大学 | Acoustic vibration response analysis method under non-stable random dynamic loads |
CN106227947A (en) * | 2016-07-26 | 2016-12-14 | 南京航空航天大学 | A kind of cooling tower inner surface Equivalent Wind Load obtaining value method |
CN106484952A (en) * | 2016-09-14 | 2017-03-08 | 东南大学 | A kind of equivalence techniques of the random face pressure load model of reverberation field |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6363789B1 (en) * | 2000-05-31 | 2002-04-02 | The Boeing Company | Acoustic pressure load conversion method to vibration spectra |
GB201104413D0 (en) * | 2011-03-16 | 2011-04-27 | Airbus Operations Ltd | Stresses induced by random loading |
CN107169217B (en) * | 2017-05-25 | 2018-04-24 | 东南大学 | A kind of equivalent method of turbulent boundary layer load model |
-
2017
- 2017-05-25 CN CN201710378734.7A patent/CN107169217B/en active Active
-
2018
- 2018-04-10 WO PCT/CN2018/082580 patent/WO2018214657A1/en active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106055798A (en) * | 2016-06-02 | 2016-10-26 | 东南大学 | Acoustic vibration response analysis method under non-stable random dynamic loads |
CN106227947A (en) * | 2016-07-26 | 2016-12-14 | 南京航空航天大学 | A kind of cooling tower inner surface Equivalent Wind Load obtaining value method |
CN106484952A (en) * | 2016-09-14 | 2017-03-08 | 东南大学 | A kind of equivalence techniques of the random face pressure load model of reverberation field |
Non-Patent Citations (1)
Title |
---|
基于模态叠加法的面压载荷与基础激励等效方法;张鹏等;《东南大学学报》;20170331;第47卷(第2期);第331-336页 * |
Also Published As
Publication number | Publication date |
---|---|
WO2018214657A1 (en) | 2018-11-29 |
CN107169217A (en) | 2017-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wolf et al. | Trailing edge noise reduction of wind turbine blades by active flow control | |
De La Torre et al. | Experimental investigation of added mass effects on a hydrofoil under cavitation conditions | |
CN107133422B (en) | A kind of certainty acoustic power flow response predicting method based on antithesis modal equations | |
CN104897330B (en) | Film structure pre-tension measuring instrument and method based on static pressure deformation | |
Amin et al. | Characteristics of wind forces and responses of rectangular tall buildings | |
CN107368630A (en) | A kind of numerical method for obtaining coupling loss factor | |
Rainbird et al. | Blockage-tolerant wind tunnel measurements for a NACA 0012 at high angles of attack | |
CN107169217B (en) | A kind of equivalent method of turbulent boundary layer load model | |
Chen et al. | Unsteady aerodynamics on a tapered prism under forced excitation | |
Mayer et al. | Design of a Kevlar-walled test section with dynamic turntable and aeroacoustic investigation of an oscillating airfoil | |
Mayer et al. | Aeroacoustic investigation of an oscillating airfoil in the pre-and post-stall regime | |
CN109508490B (en) | Acoustic model equivalent method for hollow aluminum profile | |
Du et al. | Enabling a wind energy harvester based on ZnO thin film as the building skin | |
Bowen et al. | A thorough experimental investigation on airfoil turbulence interaction noise | |
Barone et al. | Aerodynamic and aeroacoustic properties of a flatback airfoil: an update | |
Poplingher et al. | Transonic Shock Buffet on the Benchmark Supercritical Wing | |
CN106484952B (en) | A kind of equivalence techniques of the random face pressure load model of reverberation field | |
Luschi et al. | A simple analytical model for the resonance frequency of perforated beams | |
CN204085595U (en) | A kind of gas flow sensor | |
CN103226944B (en) | The system of selection of medium and low frequency sound insulation vacuum clip laminate and panel material thereof | |
CN112541276B (en) | Dynamic response prediction equivalent method based on scaling model | |
Berntsen | The remote microphone probe for unsteady surface pressure measurement | |
Cheng et al. | Rotor broadband noise due to surface roughness during ice accretion | |
Kawashima et al. | Measurements of Unsteady Force Response on Airfoils with Arbitrarily Shaped Thickness due to Incident Large-Scale Turbulence | |
D’Alessio | Flow past a slippery cylinder: part 1—circular cylinder |
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 |