CN109190328B - Multi-parameter structure dynamic response sensitivity analysis method based on mixed finite-complex difference - Google Patents
Multi-parameter structure dynamic response sensitivity analysis method based on mixed finite-complex difference Download PDFInfo
- Publication number
- CN109190328B CN109190328B CN201811424490.2A CN201811424490A CN109190328B CN 109190328 B CN109190328 B CN 109190328B CN 201811424490 A CN201811424490 A CN 201811424490A CN 109190328 B CN109190328 B CN 109190328B
- Authority
- CN
- China
- Prior art keywords
- response
- parameter
- parameters
- dynamic response
- perturbation
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000010206 sensitivity analysis Methods 0.000 title description 14
- 230000035945 sensitivity Effects 0.000 claims abstract description 28
- 238000004458 analytical method Methods 0.000 claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims description 29
- 230000001133 acceleration Effects 0.000 claims description 15
- 238000005457 optimization Methods 0.000 claims description 15
- 238000013016 damping Methods 0.000 claims description 12
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 230000010354 integration Effects 0.000 claims description 7
- 230000005284 excitation Effects 0.000 claims description 5
- 239000000284 extract Substances 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 2
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects 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/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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/06—Power analysis or power optimisation
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)
- Complex Calculations (AREA)
Abstract
The invention provides a method for analyzing the dynamic response sensitivity of a multi-parameter structure with mixed finite-complex difference, which aims at the problem of the optimized design of the structure dynamics with multiple parameters, expands design parameters from a real number domain to a complex number domain, performs real part perturbation and imaginary part perturbation on two different design parameters respectively, performs dynamic analysis on the structure, extracts the real part response and the imaginary part response of an analysis result, obtains the dynamic response sensitivity of the two design parameters simultaneously, realizes the analysis of the dynamic response sensitivity of the multi-parameter structure, and provides a quick and effective method for analyzing the dynamic response sensitivity of the multi-parameter structure for engineering application.
Description
Technical Field
The invention belongs to the field of structure dynamics optimization design, and particularly relates to a dynamic response sensitivity analysis method of a mixed finite-complex difference multi-parameter structure.
Background
The dynamic response sensitivity analysis of the structural parameters is used as an important link in the structural dynamics optimization design based on the gradient analysis, and the influence degree and the influence rule of the structural design parameters on the output response are represented. An optimization target can be constructed by utilizing the output response of the structure, and the structure dynamics optimization design is realized according to the sensitivity analysis result of the structure dynamic response.
The existing methods can be divided into a theoretical analytical method, a finite difference method, a complex variation difference method and a semi-analytical method. The finite difference method and the complex variation difference method are simple in theory and easy to apply. However, when analyzing the structural dynamic response sensitivity of multiple parameters, perturbation needs to be performed on each parameter to increase the calculation time of re-analysis, and the analysis efficiency is low. How to realize the sensitivity analysis of a plurality of parameters in one-time analysis, reduce the calculation times, improve the sensitivity analysis efficiency and the structural dynamics optimization efficiency becomes an actual engineering problem to be solved urgently.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the problems, the invention provides a method for analyzing the dynamic response sensitivity of a multi-parameter structure with mixed finite-complex difference, which can be used for simultaneously analyzing the sensitivity of the structure with a plurality of optimized parameters by simultaneously constructing the real part perturbation and the imaginary part perturbation of the optimized parameters, and can effectively improve the analysis efficiency of the dynamic response sensitivity of the multi-parameter structure.
The technical scheme is as follows: in order to realize the purpose of the invention, the technical scheme adopted by the invention is as follows: a method of mixed finite-complex differential structure dynamic response sensitivity analysis, the method comprising the steps of:
(1) finite element modeling is carried out on the engineering structure to obtain a structure dynamics control equation, and the equation is solved by utilizing a numerical analysis method of time domain integration to obtain dynamic response before parameter perturbation;
(2) constructing finite difference and complex variation difference shooting quantities for different optimization parameters at the same time, and substituting the parameters after shooting into a structural dynamics control equation to carry out dynamic response solution;
(3) and extracting a real part and an imaginary part in the structure dynamic response data after the parameters perturb, simultaneously calculating the dynamic response sensitivity of different parameters, and constructing a time-course curve of the structure dynamic response sensitivity.
Further, in the step (1), finite element modeling is carried out on the engineering structure to obtain a dynamic control equation of the structure, and the equation is solved by using a numerical analysis method of time domain integration to obtain the structural dynamic response before parameter perturbation, and the specific steps are as follows:
(1.1) carrying out finite element modeling on the engineering structure to obtain a dynamic control equation of the structure:
wherein M represents the mass matrix of the finite element model, C represents the damping matrix of the finite element model, K represents the stiffness matrix of the finite element model, F (t) represents the external excitation vector to which the finite element model is subjected, x (p, t) represents the displacement response vector of the structure,a velocity response vector representing the structure is generated,representing an acceleration response vector of the structure, and p represents a design parameter vector;
and (1.2) solving the dynamic control equation by using a numerical analysis method of time domain integration to obtain a structural dynamic response before parameter perturbation, namely a displacement response, a speed response or an acceleration response.
Further, in the step (2), a real part shooting amount and an imaginary part shooting amount are respectively constructed for different optimization parameters, and then a structure dynamics control equation after shooting is solved, and the specific steps are as follows:
(2.1) respectively constructing finite difference and complex variation difference shooting quantities for different optimization parameters:
wherein h and v respectively represent real part perturbation coefficient and imaginary part perturbation coefficient, p1And p2To design the components in the parameter vector p, representing two different parameters thereof,representing a parameter p1The real part amount of shooting of (a), i.e. the finite difference amount of shooting,representing a parameter p2The Real part and the imaginary part of the Real-time-based perturbation amount are respectively represented by the superscript Real and Imag, and the perturbed parameters are respectively as follows:
wherein,andrespectively represents the optimized parameters after finite difference perturbation and complex variation differential perturbation, i represents an imaginary number unit, i represents2Get perturbed parameter vector as-1
(2.2) the perturbed parameters can cause the mass matrix, the damping matrix and the rigidity matrix in the finite element structure to change, and a perturbed structure dynamics control equation is obtained:
wherein,a perturbed structural mass matrix is represented,a damping matrix of the structure after perturbation is represented,expressing the perturbed structural rigidity matrix, and solving a dynamic control equation (4) by using a time domain integral numerical analysis method to obtain a structural dynamic response after parameter perturbation:
further, in the step (3), real part and imaginary part results in the perturbed structural dynamic response data are extracted, and dynamic response sensitivities of corresponding parameters are respectively calculated, wherein the method comprises the following steps:
(3.1) based on the structural dynamic response obtained in (2.2), expanding the structural dynamic response according to a Taylor series:
wherein,andfirst and second partial derivatives of the response (.) to the parameter p, i denotes an imaginary unit, i2=-1,onRepresenting an n-th order term;
(3.2) neglecting a second order term and an n order term in the formulas (15) - (17), extracting real part and imaginary part results in the structural dynamic response data, and obtaining the dynamic response sensitivity of the corresponding parameters:
wherein Re (.) andrepresenting the real and imaginary parts of the response (#) respectively,andrespectively representing parameters p of structure displacement response1And parameter p2The sensitivity of (a) to (b) is,andrespectively representing parameters p of structure speed response1And parameter p2The sensitivity of (a) to (b) is,andrespectively representing structural acceleration response parameters p1And parameter p2The superscripts d, v, a represent displacement, velocity and acceleration, respectively.
Has the advantages that: compared with the prior art, the technical scheme of the invention has the following beneficial technical effects:
the invention provides a multi-parameter sensitivity analysis method of mixed finite-complex difference, which expands design parameters from a real number field to a complex number field and constructs real part shooting quantity and imaginary part shooting quantity of different parameters at the same time, thereby realizing the analysis of the structural dynamic response sensitivity of a plurality of parameters in one analysis.
Drawings
FIG. 1 is a model diagram of finite element analysis of a structure according to the present invention;
FIG. 2 is a graph of the calculated A-point acceleration response with respect to the parameter c1A graph comparing the results of the sensitivity analysis;
FIG. 3 is a graph of the calculated A-point acceleration response with respect to the parameter k1The sensitivity analysis of (2) is compared with the result graph.
Detailed Description
As shown in fig. 1, a method for analyzing the dynamic response sensitivity of a multi-parameter structure by mixing finite-complex differences comprises the following steps:
(1) carrying out finite element modeling on the structure shown in FIG. 1 to obtain a dynamic control equation of the structure, and solving the equation to obtain the dynamic response of the structure before parameter perturbation;
(2) respectively constructing real part shot amount and imaginary part shot amount for different optimization parameters according to the optimization requirements of the structural dynamics, and then substituting the perturbed parameters into a structural dynamics control equation to carry out dynamic response solution;
(3) and extracting real part and imaginary part results in the perturbed structure dynamic response data, and respectively calculating the structure dynamic response sensitivity time-course curves of the corresponding parameters.
A dynamic response sensitivity analysis method of a mixed finite-complex difference multi-parameter structure is verified by adopting a cantilever beam structure with an elastic support, as shown in figure 1, the parameters of the structure are respectively as follows: the cantilever beam structure has length of 0.2m, circular cross section and radius of 0.005m, and is made of aluminum material with elastic modulus of 70GPa and density of 2.7 × 103kg/m3The poisson ratio μ is 0.3. And the parameters can be set to different values according to different design requirements.
The design parameters are the rigidity coefficient and the damping coefficient of five elastic supports, and the total number of the design parameters is 10, and a formed design parameter vector p ═ c1,k1,c2,k2,c3,k3,c4,k4,c5,k5Specific values of the design parameters are listed in table 1:
TABLE 1 values of design parameters for elastic bearings
The specific operation is as follows:
(1.1) carrying out finite element modeling on the structure, dividing the finite element modeling into 10 beam units, as shown in figure 1, obtaining a mass matrix, a damping matrix and a rigidity matrix according to the materials and geometrical parameters of the structure and the support, and finally obtaining a structure dynamics control equation as follows:
wherein M represents the mass matrix of the finite element model, K represents the stiffness matrix of the finite element model, C represents the damping matrix of the finite element model, and the proportional damping is selected in the present example, wherein C is α M + β K, wherein α is β is 0.0001. F (t) represents the external excitation vector to which the finite element model is subjected, in the present example, 10sin (20t) external excitation is applied in the vertical direction of A point, as shown in FIG. 1, x (p, t) represents the displacement response vector of the structure,a velocity response vector representing the structure is generated,an acceleration response vector representing the structure, and p a design parameter vector.
(1.2) solving the equation to obtain the response of the structure before the parameter perturbation, wherein the calculation cut-off time of the structure dynamic response is 0.01 s.
(2) Respectively constructing real part shot amount and imaginary part shot amount for different optimization parameters according to the optimization requirements of the structural dynamics, and then substituting the perturbed parameters into a structural dynamics control equation to carry out dynamic response solution;
(2.1) respectively constructing real part shooting quantity and imaginary part shooting quantity for different optimization parameters:
wherein h and v respectively represent a real part perturbation coefficient and an imaginary part perturbation coefficient, and h-v-10 is selected in the calculation example-3,k1And c1The components in the parameter vector are designed for the structure, representing two different parameters thereof,andrespectively represent the parameter c1Real part of the amount of shooting, parameter k1The Real and imaginary parts are labeled Real and Imag, respectively. The parameters after perturbation are respectively as follows:
wherein,andrespectively representing the optimized parameters after real part perturbation and imaginary part perturbation, i represents an imaginary unit, i2Is-1. Perturbed design parameter vectorWhereinIn the case of a real number,is a plurality of numbers. To this end, the design parameters are expanded from the real-number domain to the complex-number domain.
(2.2) in the present embodiment, the perturbed parameters may cause the damping matrix and the stiffness matrix in the finite element structure to change, so as to obtain a perturbed structural dynamics control equation:
wherein,represents the perturbed structural mass matrix, which in this example is equal to M in equation (21), i.e. A damping matrix of the structure after perturbation is represented,and F (t) representing a structural rigidity matrix after perturbation, wherein the external excitation vector received by the finite element model is the same as the value before parameter perturbation. Solving a dynamic control equation (24) by using a numerical analysis method of time domain integration to obtain a perturbed structural dynamic response:
(3) and extracting real part and imaginary part results in the perturbed structure dynamic response data, and respectively calculating the dynamic response sensitivity time-course curves of the corresponding parameters.
(3.1) based on the structural dynamic response obtained in (2.2), expanding the structural dynamic response according to a Taylor series:
wherein,andrespectively representing responses (. alpha.). to parameter c1The first and second partial derivatives of (a),andrespectively representing responses (. alpha.). to the parameter k1I denotes the imaginary unit, i2=-1,onRepresenting an n-th order term.
(3.2) neglecting second-order terms and n-order terms of the equations (25) - (27), extracting real part and imaginary part results in the structural dynamic response data after A point perturbation, and respectively calculating the dynamic response sensitivity time-course curve of the corresponding parameter:
wherein Re (.) andrepresenting the real and imaginary parts of the response (#) respectively,andrespectively representing parameters p of structure displacement response1And parameter p2The sensitivity of (a) to (b) is,andrespectively representing parameters p of structure speed response1And parameter p2The sensitivity of (a) to (b) is,andrespectively representing structural acceleration response parameters p1And parameter p2The superscripts d, v, a represent displacement, velocity and acceleration, respectively.Respectively representing the displacement response, the velocity response and the acceleration response of the point A on the structure. FIG. 2 and FIG. 3 show the acceleration response at point A with respect to the parameter c1And k1The results of the method, the finite difference sensitivity analysis method and the complex variation difference sensitivity analysis method are respectively compared, and the results are consistent.
Claims (2)
1. A method for analyzing the sensitivity of a structure's dynamic response to a mixed finite-complex difference, the method comprising the steps of:
(1) finite element modeling is carried out on the engineering structure to obtain a structure dynamics control equation, and the equation is solved by utilizing a numerical analysis method of time domain integration to obtain dynamic response before parameter perturbation;
(2) constructing finite difference and complex variation difference shooting quantities for different optimization parameters at the same time, and substituting the parameters after shooting into a structural dynamics control equation to carry out dynamic response solution;
the method comprises the following specific steps:
(2.1) respectively constructing finite difference and complex variation difference shooting quantities for different optimization parameters:
wherein h and v respectively represent real part perturbation coefficient and imaginary part perturbation coefficient, p1And p2To design the components in the parameter vector p, representing two different parameters thereof,representing a parameter p1The real part amount of shooting of (a), i.e. the finite difference amount of shooting,representing a parameter p2The Real part and the imaginary part of the Real-time-based perturbation amount are respectively represented by the superscript Real and Imag, and the perturbed parameters are respectively as follows:
wherein,andrespectively represents the optimized parameters after finite difference perturbation and complex variation differential perturbation, i represents an imaginary number unit, i represents2Get perturbed parameter vector as-1
(2.2) the perturbed parameters can cause the mass matrix, the damping matrix and the rigidity matrix in the finite element structure to change, and a perturbed structure dynamics control equation is obtained:
wherein,a perturbed structural mass matrix is represented,a damping matrix of the structure after perturbation is represented,expressing the perturbed structural rigidity matrix, and solving a dynamic control equation (3) by using a time domain integral numerical analysis method to obtain a structural dynamic response after parameter perturbation:
(3) extracting a real part and an imaginary part in the structure dynamic response data after the parameter perturbation, simultaneously calculating the dynamic response sensitivity of different parameters, and constructing a time-course curve of the structure dynamic response sensitivity, wherein the method comprises the following steps:
(3.1) based on the structural dynamic response obtained in (2.2), expanding the structural dynamic response according to a Taylor series:
wherein,andfirst and second partial derivatives of the response (.) to the parameter p, i denotes an imaginary unit, i2=-1,onRepresenting an n-th order term;
(3.2) neglecting second-order terms and n-order terms in the formulas (4) - (6), extracting real part and imaginary part results in the structural dynamic response data, and obtaining the dynamic response sensitivity of the corresponding parameters:
wherein Re (.) andrepresenting the real and imaginary parts of the response (#) respectively,andrespectively representing parameters p of structure displacement response1And parameter p2The sensitivity of (a) to (b) is,andrespectively representing parameters p of structure speed response1And parameter p2The sensitivity of (a) to (b) is,andrespectively representing structural acceleration response parameters p1And parameter p2The superscripts d, v, a represent displacement, velocity and acceleration, respectively.
2. The method for analyzing the sensitivity of the structural dynamic response of the hybrid finite-complex difference as claimed in claim 1, wherein in the step (1), finite element modeling is performed on the engineering structure to obtain a structural dynamic control equation, and the equation is solved by using a numerical analysis method of time domain integration to obtain the structural dynamic response before parameter perturbation, and the specific steps are as follows:
(1.1) carrying out finite element modeling on the engineering structure to obtain a dynamic control equation of the structure:
wherein M represents the mass matrix of the finite element model, C represents the damping matrix of the finite element model, K represents the stiffness matrix of the finite element model, F (t) represents the external excitation vector to which the finite element model is subjected, x (p, t) represents the displacement response vector of the structure,a velocity response vector representing the structure is generated,representing an acceleration response vector of the structure, and p represents a design parameter vector;
and (1.2) solving the dynamic control equation by using a numerical analysis method of time domain integration to obtain a structural dynamic response before parameter perturbation, namely a displacement response, a speed response or an acceleration response.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811424490.2A CN109190328B (en) | 2018-11-27 | 2018-11-27 | Multi-parameter structure dynamic response sensitivity analysis method based on mixed finite-complex difference |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811424490.2A CN109190328B (en) | 2018-11-27 | 2018-11-27 | Multi-parameter structure dynamic response sensitivity analysis method based on mixed finite-complex difference |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109190328A CN109190328A (en) | 2019-01-11 |
CN109190328B true CN109190328B (en) | 2020-03-31 |
Family
ID=64938079
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811424490.2A Active CN109190328B (en) | 2018-11-27 | 2018-11-27 | Multi-parameter structure dynamic response sensitivity analysis method based on mixed finite-complex difference |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109190328B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109885896B (en) * | 2019-01-25 | 2020-04-24 | 东南大学 | Nonlinear structure finite element model correction method based on complex variation differential sensitivity |
CN109902357B (en) * | 2019-01-31 | 2020-06-02 | 东南大学 | Method for analyzing dynamic response sensitivity of complex-variant differential smooth nonlinear structure |
CN114722664B (en) * | 2022-03-30 | 2023-06-20 | 哈尔滨理工大学 | Dynamic stiffness topology optimization method for multi-material structure |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7844431B2 (en) * | 2004-02-20 | 2010-11-30 | The Mathworks, Inc. | Method and apparatus for integrated modeling, simulation and analysis of chemical and biochemical reactions |
CN105912508B (en) * | 2016-04-28 | 2019-05-10 | 北京航空航天大学 | A kind of Random Perturbation method of the improved repetition structural vibration characteristic value based on agent model |
CN107609221B (en) * | 2017-08-15 | 2018-07-17 | 东南大学 | It is a kind of that hinged structure nonlinear parameter recognition methods is contained based on genetic algorithm |
CN108334670A (en) * | 2018-01-12 | 2018-07-27 | 厦门大学 | A kind of complex model uncertainty correction method for finite element model based on perturbation method |
CN108319778A (en) * | 2018-01-31 | 2018-07-24 | 东南大学 | A kind of structural parameters Sensitivity Analysis Method based on Complex Function Method |
-
2018
- 2018-11-27 CN CN201811424490.2A patent/CN109190328B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN109190328A (en) | 2019-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109190328B (en) | Multi-parameter structure dynamic response sensitivity analysis method based on mixed finite-complex difference | |
Atif et al. | Essentially entropic lattice Boltzmann model | |
CN109902404B (en) | Unified recursion calculation method for structural time-course response integral of different damping forms | |
CN104268342A (en) | Method for analyzing vibration characteristics of bearing retainer based on finite elements | |
CN105631125A (en) | Aerodynamic-thermal-structural coupling analysis method based on reduced-order model | |
CN108121865B (en) | A kind of structure section response propagation analysis method containing more uncertain parameters based on adjoint variable method | |
Akhmedzyanov et al. | Computer-aided design and construction development of the main elements of aviation engines | |
Rasquin et al. | Scalable fully implicit finite element flow solver with application to high-fidelity flow control simulations on a realistic wing design | |
CN109902357B (en) | Method for analyzing dynamic response sensitivity of complex-variant differential smooth nonlinear structure | |
CN107563102A (en) | A kind of power transmission skeleton method for visualizing of bearing structure | |
Tang et al. | Unsteady flow simulation with mesh adaptation | |
Ahmad et al. | Design and modal analysis of a quadcopter propeller through finite element analysis | |
CN107291992A (en) | A kind of applicable desert area electronics combined environment test analogue system and method | |
CN111563297A (en) | Supersonic aircraft thermal environment calculation method based on BP network | |
Štoller et al. | Design and assessment of shape of protective structure by usage of CFD software environment Ansys Fluent | |
Lamberson et al. | Aeroelastic simulations with modal and finite-element structural solvers using CREATE-AV/Kestrel v5 | |
Ahmad et al. | FEA based frequency analysis of unmanned aerial vehicle (UAV) | |
Di Zhang et al. | An improved least square similitude method based on Lagrange energy for estimating scaling laws and eliminating coupling effects | |
Ghoreyshi et al. | CFD calculation of aerodynamic indicial functions for a generic fighter configuration | |
Silva et al. | Computational Results for the KTH-NASA Wind-Tunnel Model Used for Acquisition of Transonic Nonlinear Aeroelastic Data | |
Ghoreyshi et al. | Indicial aerodynamic modeling for arbitrary gust responses | |
Zuo et al. | Efficient aero-structural design optimization: Coupling based on reverse iteration of structural model | |
Holmes et al. | Simulating bird strikes using smoothed particle hydrodynamics for improved aircraft safety | |
Chaaban et al. | A multiscale study of the retention behavior and hydraulic anisotropy in deformable porous media | |
Abu-Zurayk et al. | Development of the adjoint approach for aeroelastic wing optimization |
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 |