CN113125262A - Method for quickly forecasting deformation of breakable calcareous sand in loading process - Google Patents
Method for quickly forecasting deformation of breakable calcareous sand in loading process Download PDFInfo
- Publication number
- CN113125262A CN113125262A CN202110413272.4A CN202110413272A CN113125262A CN 113125262 A CN113125262 A CN 113125262A CN 202110413272 A CN202110413272 A CN 202110413272A CN 113125262 A CN113125262 A CN 113125262A
- Authority
- CN
- China
- Prior art keywords
- strain
- epsilon
- stress
- deformation
- shear
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/24—Investigating strength properties of solid materials by application of mechanical stress by applying steady shearing forces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Pathology (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Computational Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention relates to a method for quickly forecasting deformation of breakable calcareous sand in a loading process, which is used for decoupling the calculation of stress ratio and body strain, and can calculate the change of the stress ratio along with shear strain independently and then calculate the change of the body strain along with the shear strain without interactively calculating the stress ratio and the body strain increment in each increment step, thereby simplifying the calculation process and improving the calculation efficiency.
Description
Technical Field
The invention belongs to the field of geotechnical engineering research, and particularly relates to a rapid prediction method for deformation of breakable calcareous sand in a loading process.
Background
The calcareous sand used for the marine artificial hydraulic reclamation island has obvious particle crushing property after being loaded. For the breakable calcareous sand, a large number of triaxial test researches are carried out to analyze the stress-strain characteristics of the breakable calcareous sand, but a method is not available at present, namely, the deformation of the breakable calcareous sand in the loaded breaking process can be quickly forecasted, for example, when the vertical direction has continuous compressive strain, the changes of the forecasting stress and the volume deformation are decoupled. In the existing method, the stress and the volume deformation need to be interactively calculated at each time step, the calculation of the stress and the volume deformation is not decoupled, and the calculation program is more complex.
Disclosure of Invention
The invention provides a rapid forecasting method for deformation of breakable calcareous sand in a loading process, which aims to decouple and forecast the stress and the volume deformation of the breakable calcareous sand in the loading process and carry out rapid forecasting on the deformation of the breakable calcareous sand in the loading process without interactively calculating the stress and the volume deformation of a soil sample at each time step.
The invention relates to some abbreviations and symbols, the following are notes:
σ1and σ'1: vertical stress, σ ', to which the aggregate of particles is subjected'1For effective stress, where σ1And σ'1Taking the same value;
σ2and σ3,σ′2And σ'3: horizontal stress, σ, to which the assembly of particles is subjected2And σ3Is perpendicular to the direction of'2And σ'3For effective stress, where σ2And σ'2Taking the same value, σ3And σ'3Taking the same value;
ε1、ε2and ε3: strain and respectively stress sigma1、σ2And σ3The directions are the same;
εv: bulk strain, epsilonv=ε1+ε2+ε3;
t0,t1,t2,…,ti,…,tn: the recorded starting time in the loading process is t0The time recorded later is t from small to large1,t2,…,ti,…,tnWhere i is more than or equal to 1 and less than or equal to n, and n +1 is the number of recorded time points;
(εv)i: t thiBody strain epsilon corresponding to timev;
(εs)i: t thiShear strain epsilon corresponding to times;
(△εv)i: increase in bulk strain, (. DELTA.. epsilon.)v)i=(εv)i-(εv)i-1;
(△εs)i: increase in shear strain, (. DELTA.. epsilon.)s)i=(εs)i-(εs)i-1;
M: material parameters and equal to critical stress ratio;
k2: a parameter;
pa: atmospheric pressure;
p′0: an initial average effective stress;
ξ1and xi2: calculating a beta parameter;
α: material parameter,. alpha. - (C)2+1)exp(-C1εs)+(C2-1)exp(-εs)+1;
C1And C2: calculating an alpha parameter;
the technical scheme of the invention is as follows: a method for quickly forecasting the deformation of breakable calcareous sand in a loading process comprises the following steps:
step 1: let the vertical stress of the aggregate of sand particles be sigma1Stress on the horizontal plane is respectively sigma2And σ3Where σ is2And σ3Is perpendicular to the direction of the particle assembly and the strain of the particle assembly is epsilon1、ε2And ε3In which strain epsilon1、ε2And ε3Respectively in the direction of the stress sigma1、σ2And σ3The directions are the same; defining average effective stress p', shear stress q, stress ratio eta, and bulk strain epsilonvAnd shear strain epsilons:
εv=ε1+ε2+ε3 (4)
Step 2: calculating the stress ratio:
the stress ratio eta is calculated by the following formula:
where M is a material parameter and is equal to the critical stress ratio, k2Is a parameter, β is a parameter;
in the above formula paIs largeGas pressure, p'0Is the initial mean effective stress, ξ1And xi2Is a parameter;
thus resulting from shear strain epsilonsThe stress ratio eta can be directly calculated;
and step 3: calculating the volume deformation:
let the start time recorded during loading be t0The time recorded later is t from small to large1,t2,…,ti,…,tnWhere 1 ≦ i ≦ n, n +1 is the number of recorded time points, let tiBody strain epsilon corresponding to timevIs (epsilon)v)iLet a tiShear strain epsilon corresponding to timesIs (epsilon)s)iLet the shear strain increment (Delta epsilon) generated by two adjacent time differencess)iEqual, define the bulk strain increment ([ Delta ] [ epsilon ]v)iAnd increase in shear strain (DELTA ε)s)i:
(△εv)i=(εv)i-(εv)i-1 (8)
(△εs)i=(εs)i-(εs)i-1 (9)
Where the strain (epsilon) is cut at every moments)iAnd increase in shear strain (DELTA ε)s)iFor a known set value, the volume strain increment (Δ ε) is calculated as followsv)i:
In the above formula, α is a material parameter, and a relationship of α with a change in shear strain is defined as:
α=-(C2+1)exp(-C1εs)+(C2-1)exp(-εs)+1 (11)
in the above formula C1And C2For the material parameter, each t can be calculated from the equation (10)i(i is more than or equal to 1 and less than or equal to n) corresponding to the momentv)iBy increasing these bulk strains by an amount (. DELTA.. di-elect cons.)v)iWhen added together, the body strain (epsilon) at each moment is obtainedv)iI.e., (. epsilon.)v)i=(εv)i-1+(△εv)iAnd the initial t is known0Time body strain (Delta epsilon)v)iIs 0.
In step 2, calculating to obtain a stress ratio eta, and calculating to obtain an effective stress sigma 'in the horizontal direction'2And σ'3When the stress is constant and equal, the vertical effective stress sigma 'can be obtained'1And average effective and shear stresses:
in step 3, when the strain is horizontal2And ε3When the two strains are equal, the vertical strain epsilon can be obtained1And strain in horizontal direction epsilon3:
The invention has the beneficial effects that: stress ratio and body strain are calculated and decoupled, stress ratio changing along with shear strain can be calculated independently, then body strain changing along with shear strain is calculated, and stress ratio and body strain increment do not need to be calculated in each increment step in an interactive mode, so that the calculation process is simplified, and the calculation efficiency is improved. In addition, the calculation formula of the ratio of the volume strain increment and the shear strain increment given by the formula (10) along with the stress ratio is more in line with the test curve of the breakable calcareous sand compared with the existing formula.
Drawings
FIG. 1 is a schematic view of an assembly of sand particles stressed vertically and horizontally;
FIG. 2 is a graph showing the variation of stress ratio and bulk strain with shear strain;
FIG. 3 is a graph of the ratio of bulk strain increase to shear strain increase versus stress ratio.
Fig. 1. assembly of sand particles.
Detailed Description
In order to make the technical means, innovative features, objectives and effects of the present invention apparent, the present invention will be further described with reference to the following detailed drawings.
With the aggregate 1 of sand particles shown in FIG. 1, stress σ is applied in the vertical direction1Is stressed horizontally by a stress2And σ3,σ2And σ3Is vertical. Carrying out a triaxial test with a confining pressure sigma2=σ3And a change curve of the stress ratio with the shear strain and a change curve of the body strain with the shear strain can be drawn. The confining pressures in the three-axis test of the breakable calcareous sand are 100kPa, 200kPa, 300kPa, 500kPa, 700kPa and 900kPa, and the corresponding stress ratio curve and the corresponding body strain curve are shown in FIG. 2. FIG. 3 shows σ2=σ3The ratio of the increase in bulk strain to the increase in shear strain at 900kPa is plotted against the stress ratio.
The invention relates to some abbreviations and symbols, the following are notes:
σ1and σ'1: vertical stress, σ ', to which the aggregate of particles is subjected'1For effective stress, where σ1And σ'1Take the same value
σ2And σ3,σ′2And σ'3: horizontal stress, σ, to which the assembly of particles is subjected2And σ3Is perpendicular to the direction of'2And σ'3For effective stress, where σ2And σ'2Taking the same value, σ3And σ'3Take the same value
ε1、ε2And ε3: strain and respectively stress sigma1、σ2And σ3Same direction
εv: bulk strain, epsilonv=ε1+ε2+ε3
t0,t1,t2,…,ti,…,tn: the recorded starting time in the loading process is t0The time recorded later is t from small to large1,t2,…,ti,…,tnWhere 1 ≦ i ≦ n, n +1 is the number of recorded time points
(εv)i: t thiBody strain epsilon corresponding to timev
(εs)i: t thiShear strain epsilon corresponding to times
(△εv)i: increase in bulk strain, (. DELTA.. epsilon.)v)i=(εv)i-(εv)i-1
(△εs)i: increase in shear strain, (. DELTA.. epsilon.)s)i=(εs)i-(εs)i-1
M: material parameter equal to critical stress ratio
k2: parameter(s)
pa: atmospheric pressure
p′0: initial mean effective stress
ξ1And xi2: calculating beta parameter
α: material parameter,. alpha. - (C)2+1)exp(-C1εs)+(C2-1)exp(-εs)+1
C1And C2: calculating alpha parameter
The technical scheme of the invention is as follows: a method for quickly forecasting the deformation of breakable calcareous sand in a loading process comprises the following steps:
step 1: let the vertical stress of the aggregate of sand particles be sigma1Stress on the horizontal plane is respectively sigma2And σ3Where σ is2And σ3Is perpendicular to the direction of the particle assembly and the strain of the particle assembly is epsilon1、ε2And ε3In which strain epsilon1、ε2And ε3Respectively in the direction of the stress sigma1、σ2And σ3The directions are the same; defining average effective stress p', shear stress q, stress ratio eta, and bulk strain epsilonvAnd shear strain epsilons:
εv=ε1+ε2+ε3 (4)
Step 2: calculating the stress ratio:
the stress ratio eta is calculated by the following formula:
where M is a material parameter and is equal to the critical stress ratio, k2As parameters, β is:
in the above formula paIs atmospheric pressure, p'0Is the initial mean effective stress, ξ1And xi2Is a parameter such that the shear strain εsThe stress ratio eta can be directly calculated;
in this example, the horizontal effective stress σ'2And σ'3Constant and equal, so that the vertical effective stress sigma 'can be obtained'1And average effective and shear stresses:
and step 3: calculating the volume deformation:
let the start time recorded during loading be t0The time recorded later is t from small to large1,t2,…,ti,…,tnWhere 1 ≦ i ≦ n, n +1 is the number of recorded time points, let tiBody strain epsilon corresponding to timevIs (epsilon)v)iLet a tiShear strain epsilon corresponding to timesIs (epsilon)s)iLet the shear strain increment (Delta epsilon) generated by two adjacent time differencess)iEquality, defining a volume strain increaseAmount (. DELTA.. di-elect cons.)v)iAnd increase in shear strain (DELTA ε)s)i:
(△εv)i=(εv)i-(εv)i-1 (8)
(△εs)i=(εs)i-(εs)i-1 (9)
Where the strain (epsilon) is cut at every moments)iAnd increase in shear strain (DELTA ε)s)iFor a known set value, the volume strain increment (Δ ε) is calculated as followsv)i:
In the above formula, α is a material parameter, and a relationship of α with a change in shear strain is defined as:
α=-(C2+1)exp(-C1εs)+(C2-1)exp(-εs)+1 (11)
in the above formula C1And C2For the material parameter, each t can be calculated from the equation (10)i(i is more than or equal to 1 and less than or equal to n) corresponding to the momentv)iBy increasing these bulk strains by an amount (. DELTA.. di-elect cons.)v)iWhen added together, the body strain (epsilon) at each moment is obtainedv)iI.e., (. epsilon.)v)i=(εv)i-1+(△εv)iAnd the initial t is known0Time body strain (Delta epsilon)v)iIs 0;
k2is a parameter, C1And C2It is calculated in the following manner,C2=-b1(p′0/pa)+b2,a1、a2、b1、b2are all parameters.
In this example, when the strain is horizontal ε2And ε3When the two strains are equal, the vertical strain epsilon can be obtained1And strain in horizontal direction epsilon3:
The values of the parameters are as follows:
TABLE 1 parameter values
The predicted shear stress curve and body strain curve are shown in fig. 2, and the predicted change curve of the ratio of the increment of body strain to the increment of shear strain along with the stress ratio is shown in fig. 3. In fig. 3, the prediction curves given by the conventional Cam-clay model are also given, and it can be seen that the conventional Cam-clay model does not describe the initial bending segment and the right-side hooking segment, and the prediction method given herein is more consistent with the actual experiment.
Claims (6)
1. A method for quickly forecasting the deformation of breakable calcareous sand in a loading process is characterized in that:
which comprises the following steps:
step 1: let the vertical stress of the aggregate of sand particles be sigma1Stress on the horizontal plane is respectively sigma2And σ3Where σ is2And σ3Is perpendicular to the direction of the particle assembly and the strain of the particle assembly is epsilon1、ε2And ε3In which strain epsilon1、ε2And ε3Respectively in the direction of the stress sigma1、σ2And σ3The directions are the same; defining average effective stress p', shear stress q, stress ratio eta, and bulk strain epsilonvAnd shear strain epsilons:
εv=ε1+ε2+ε3 (4)
Step 2: calculating the stress ratio:
the stress ratio eta is calculated by the following formula:
where M is a material parameter and is equal to the critical stress ratio, k2Is a parameter, β is a parameter;
thus resulting from shear strain epsilonsThe stress ratio eta can be directly calculated;
and step 3: calculating the volume deformation:
let the start time recorded during loading be t0The time recorded later is t from small to large1,t2,…,ti,…,tnWhere 1 ≦ i ≦ n, n +1 is the number of recorded time points, let tiBody strain epsilon corresponding to timevIs (epsilon)v)iLet a tiShear strain epsilon corresponding to timesIs (epsilon)s)iLet the shear strain increment (Delta epsilon) generated by two adjacent time differencess)iEqual, define the bulk strain increment ([ Delta ] [ epsilon ]v)iAnd increase in shear strain (DELTA ε)s)i:
(△εv)i=(εv)i-(εv)i-1 (8)
(△εs)i=(εs)i-(εs)i-1 (9)
Where the strain (epsilon) is cut at every moments)iAnd increase in shear strain (DELTA ε)s)iFor a known set value, the volume strain increment (Δ ε) is calculated as followsv)i:
In the above formula, α is a material parameter, and each t can be calculated from the formula (10)iTime-dependent increase in body strain (Delta epsilon)v)iBy increasing these bulk strains by an amount (. DELTA.. di-elect cons.)v)iWhen added together, the body strain (epsilon) at each moment is obtainedv)iI.e., (. epsilon.)v)i=(εv)i-1+(△εv)iAnd the initial t is known0Time body strain (Delta epsilon)v)iIs 0.
2. The method for rapidly forecasting the deformation of the breakable calcareous sand in the loading process according to claim 1, characterized in that: in the step 2, the step of the method is carried out,in the above formula paIs atmospheric pressure, p'0Is the initial mean effective stress, ξ1And xi2Is a parameter.
3. The method for rapidly forecasting the deformation of the breakable calcareous sand in the loading process according to claim 1, characterized in that: in step 3, the relationship of α with the change in shear strain is defined as α ═ C2+1)exp(-C1εs)+(C2-1)exp(-εs)+1,C1And C2Is a material parameter.
4. The method for rapidly forecasting the deformation of the breakable calcareous sand in the loading process according to claim 1, characterized in that: in step 2, calculating to obtain a stress ratio eta, and calculating to obtain an effective stress sigma 'in the horizontal direction'2And σ3'constant and equal,' the vertical effective stress sigma 'can be obtained'1And average effective and shear stresses:
5. the method for rapidly forecasting the deformation of the breakable calcareous sand in the loading process according to claim 1, characterized in that: in step 3, when the strain is horizontal2And ε3When the two strains are equal, the vertical strain epsilon can be obtained1And strain in horizontal direction epsilon3:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110413272.4A CN113125262B (en) | 2021-04-16 | 2021-04-16 | Method for quickly forecasting deformation of breakable calcareous sand in loading process |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110413272.4A CN113125262B (en) | 2021-04-16 | 2021-04-16 | Method for quickly forecasting deformation of breakable calcareous sand in loading process |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113125262A true CN113125262A (en) | 2021-07-16 |
CN113125262B CN113125262B (en) | 2022-09-09 |
Family
ID=76777082
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110413272.4A Active CN113125262B (en) | 2021-04-16 | 2021-04-16 | Method for quickly forecasting deformation of breakable calcareous sand in loading process |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113125262B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4854175A (en) * | 1988-02-29 | 1989-08-08 | The Research Foundation Of State University Of New York | Simple shear device for testing earthen materials and powders |
CN110411804A (en) * | 2019-09-02 | 2019-11-05 | 上海交通大学 | A kind of contact face mechanical property test sample, preparation method and the test method of the soil body and structure |
CN111024486A (en) * | 2019-12-19 | 2020-04-17 | 南京航空航天大学 | Creep behavior prediction method for unidirectional ceramic matrix composite |
CN112014242A (en) * | 2020-09-04 | 2020-12-01 | 长沙理工大学 | Three-dimensional strain failure criterion-based asphalt pavement load checking method |
-
2021
- 2021-04-16 CN CN202110413272.4A patent/CN113125262B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4854175A (en) * | 1988-02-29 | 1989-08-08 | The Research Foundation Of State University Of New York | Simple shear device for testing earthen materials and powders |
CN110411804A (en) * | 2019-09-02 | 2019-11-05 | 上海交通大学 | A kind of contact face mechanical property test sample, preparation method and the test method of the soil body and structure |
CN111024486A (en) * | 2019-12-19 | 2020-04-17 | 南京航空航天大学 | Creep behavior prediction method for unidirectional ceramic matrix composite |
CN112014242A (en) * | 2020-09-04 | 2020-12-01 | 长沙理工大学 | Three-dimensional strain failure criterion-based asphalt pavement load checking method |
Non-Patent Citations (4)
Title |
---|
孙吉主 等: "三轴压缩条件下钙质砂的颗粒破裂过程研究", 《岩土力学》 * |
王年香 等: "软粘土的归一化应力应变模式", 《水利水运工程学报》 * |
邓陈艳 等: "纳米硅溶胶加固标准砂三轴试验研究", 《科技通报》 * |
金炜枫 等: "引入流体方程的离散颗粒–连续土体耦合方法研究", 《岩石力学与工程学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113125262B (en) | 2022-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sawaguchi et al. | Designing Fe–Mn–Si alloys with improved low-cycle fatigue lives | |
Fan et al. | Internal-state-variable based self-consistent constitutive modeling for hot working of two-phase titanium alloys coupling microstructure evolution | |
Lapin et al. | Room temperature mechanical behaviour of cast in-situ TiAl matrix composite reinforced with carbide particles | |
Koyama et al. | Influence of dislocation separation on dynamic strain aging in a Fe–Mn–C austenitic steel | |
Bao et al. | On plane stress–plane strain interactions in fatigue crack growth | |
Voyiadjis et al. | Damage model for concrete using bounding surface concept | |
CN113125262B (en) | Method for quickly forecasting deformation of breakable calcareous sand in loading process | |
Wu et al. | Enhanced strength-ductility synergy in fiber-like structural titanium matrix composites by controlling TiB content | |
Xiao et al. | Microstructural development under interrupted hot deformation and the mechanical properties of a cast Mg–Gd–Y–Zr alloy | |
Maeda et al. | Experimental analysis of aluminum yield surface for binary Al Mg alloy sheet samples | |
Krauß et al. | Revealing the local fatigue behavior of bimodal copper laminates by micropillar fatigue tests | |
Wang et al. | Hot deformation behavior and microstructure evolution of flaky Ti3AlC2 particles reinforced pure Al composite | |
CN116432284A (en) | Method for predicting carbonization depth of sprayed concrete based on initial damage degree | |
Pan et al. | Rate-change instrumented indentation for measuring strain rate sensitivity | |
Huang et al. | Modeling texture evolution during rolling process of AZ31 magnesium alloy with elasto-plastic self consistent model | |
Doucet et al. | Plane-strain work hardening and transient behavior of interstitial-free steel | |
Ju et al. | Anisotropic microstructure and mechanical properties of as-forged (Ti, Nb) B/Ti2AlNb composites | |
Bolesta et al. | Molecular dynamics simulation of sphere indentation in a thin copper film | |
Drozd et al. | Degradation of the mechanical properties of a Mg–Li–Al composite at elevated temperatures studied by the stress relaxation technique | |
Greeff et al. | Modeling dynamic phase transitions in Ti and Zr | |
CN113155612B (en) | Deformation prediction method for microfiber mixed silica sol solidified calcareous sand | |
Chuvil'Deev et al. | Dispersion Limit upon Equal-Channel Angular Pressing. Temperature Effect. | |
CN113125263B (en) | Forecasting method for stress deformation of silica sol cured non-breakable sandy soil | |
Graves et al. | Evaluation of A Ti-22ai-23nb “Orthorhombic” Alloy for use as the Matrix in A High Temperature Ti-Based Composite | |
Champagne et al. | Some Material Characteristics of Cold‐Sprayed 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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |