CN105243185A - Monte Carlo method based tensegrity structure form-finding method - Google Patents

Monte Carlo method based tensegrity structure form-finding method Download PDF

Info

Publication number
CN105243185A
CN105243185A CN201510573837.XA CN201510573837A CN105243185A CN 105243185 A CN105243185 A CN 105243185A CN 201510573837 A CN201510573837 A CN 201510573837A CN 105243185 A CN105243185 A CN 105243185A
Authority
CN
China
Prior art keywords
configuration
unit
node
monte carlo
new configuration
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
Application number
CN201510573837.XA
Other languages
Chinese (zh)
Other versions
CN105243185B (en
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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201510573837.XA priority Critical patent/CN105243185B/en
Publication of CN105243185A publication Critical patent/CN105243185A/en
Application granted granted Critical
Publication of CN105243185B publication Critical patent/CN105243185B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Complex Calculations (AREA)

Abstract

The present invention provides a Monte Carlo method based tensegrity structure form-finding method, and belongs to the field of tensegrity structure statics analysis. Aiming at the problems of lack of corresponding programs, a relatively low calculation speed, and lack of definitions for basic concepts and a processing method for particular cases, the steps adopted by the method of the present invention are as follows: step 1. setting an initial condition; step 2. randomly generating an initial configuration; step 3. randomly generating a new configuration and determining whether to accept the new configuration; step 4. checking whether a system is balanced; and step 5. obtaining a form-finding result. According to the Monte Carlo method based tensegrity structure form-finding method provided by the present invention, the particular case and technical problem in form-finding calculation of a tensegrity structure are solved, so that a speed of form-finding calculation is quickened; and a form-finding calculation process of the tensegrity structure by a Monte Carlo method is implemented by using MATLAB programming, so that a self-balanced configuration can be obtained rapidly and accurately, and the effects are significant.

Description

A kind of tension integral structure based on Monte Carlo method looks for shape method
Technical field
The invention belongs to tension integral structure statics Analysis field, what relate to a kind of tension integral structure looks for shape method.
Background technology
The seif-citing rate cable net structure that tension integral structure is made up of continuous stay cable and individual stamp, its concept is proposed by Fuller the earliest, owing to having graceful profile, and take full advantage of properties of materials, can weight reduction save material, to get the attention in recent years and developing rapidly, and achieving in many fields such as civil engineering work, aircraft deployable mechanism, biological cell engineering and materialogies and apply widely.
After the topological structure of a tension integral structure and element length are designated, need to determine its self-balance structure, be namely called and look for shape.Can shape be looked for be related to structure normally work, the research for tension integral structure be very important.Researchers, to the large quantifier elimination of having looked for shape method to carry out of tension integral structure, can use analytical method to carry out looking for shape for the simple tension integral structure with high symmetry.But for structure more generally, need to use numerical method.Force density method, Dynamic Relaxation, evolution algorithm, genetic algorithm etc. are all seen in report.However, still to lack more flexibly, more powerful, can be used for large-size or geometry is irregular and asymmetric tension integral structure form looks for the method for shape.For this reason, researchers propose Monte Carlo method introduced tension integral structure look for shape.
But these study the method not utilizing specific procedure to realize it, and being therefore difficult to is applied in engineering actual computation and guides the solution of technical matters.Meanwhile, method still has lower computing velocity, lacks the definition of basic concepts and the disposal route to special circumstances in computation process simultaneously.
Summary of the invention
For lacking corresponding program, lower and the problem lacked for the definition of key concept and the disposal route of special circumstances of computing velocity, the present invention proposes a kind of tension integral structure based on Monte Carlo method and look for shape method, look for shape to calculate to tension integral structure, thus obtain self-balance structure in engineering.
The object of the invention is to be achieved through the following technical solutions:
Tension integral structure based on Monte Carlo method looks for a shape method, comprises the steps:
Step one, setting starting condition:
In a program the geometric configuration of tension integral structure and material properties are set, the rigidity of unit is calculated, proposes the concept of connection matrix simultaneously and the unit order of connection is specified;
Step 2, stochastic generation initial configuration:
Utilize program stochastic generation initial configuration, obtain the coordinate vector matrix of each unit, then calculate the energy of configuration, and set node random file radius;
The new configuration of step 3, stochastic generation also judges whether to accept:
First node or last node are fixed, thus raises the efficiency and do not affect and look for shape result.Select a node at random again, with selected node for the centre of sphere, r is the reposition of interior random this node of setting of the ball of radius, thus generates new configuration.Constantly judge whether that drag-line occurs lax and deals with to it in the process, if without cable slack, carried out the calculating of new configuration energy, if there is cable slack, retightens node and carry out next round calculating.
Calculate the energy of new configuration, if be less than the energy of previous step configuration, accept new configuration, and replace original configuration with new configuration, if be greater than the energy of previous step, do not replace, carry out next round calculating.
When iterations reaches the integral multiple of 100, then calculate the ratio accepting new configuration in this 100 step, if ratio is lower, then reduced radius.
Whether step 4, check system balance:
Whether the balance of system is balanced by node and judges, if the internal force meeting each unit is 0, and system balancing.In fact want to reach complete equipilibrium due to error and just need more iterations, and reduce gradually along with the increase nodal displacement radius of iterations, more iterations does not just have too large practical significance.Therefore can limit according to error specification error, setting approximate equilibrium condition.
Step 5, obtain looking for shape result:
Shape result is looked in output, i.e. each node coordinate, and draws tension integral structure self-equilibrating configuration according to coordinate.
The present invention has the following advantages:
1, the look for special circumstances in shape calculating and technical matters for tension integral structure solve, thus accelerate the speed of looking for shape to calculate;
2, utilize MATLAB programming realization Monte Carlo method to look for shape computation process to tension integral structure, quick and precisely can obtain its self-equilibrating configuration, Be very effective.
Accompanying drawing explanation
Fig. 1 is the process flow diagram that Monte Carlo method looks for shape to calculate;
Fig. 2 is that shape computation process energy curves figure looked for by T prism;
Fig. 3 is that the initial configuration generated in shape computation process looked for by T prism;
Fig. 4 is that the configuration obtained after 10000 steps calculate in shape computation process looked for by T prism;
Fig. 5 is that the configuration obtained after 30000 steps calculate in shape computation process looked for by T prism;
Fig. 6 is the final configuration meeting the setting of approximate equilibrium condition that T prism looks for shape to calculate.
Embodiment
Below in conjunction with accompanying drawing, technical scheme of the present invention is further described; but be not limited thereto; everyly technical solution of the present invention modified or equivalent to replace, and not departing from the spirit and scope of technical solution of the present invention, all should be encompassed in protection scope of the present invention.
Embodiment one: present embodiments provide for parameter definition and establishing method that a kind of tension integral structure based on Monte Carlo method looks for shape method, as shown in Figure 1, specifically comprise the steps:
Step one, setting starting condition, can be divided into the following steps to carry out:
1.1, remember that unit adds up to m, nodes is n, and depression bar number is e b, drag-line number is e s, wherein e b+ e s=m.
1.2, material properties: (depression bar and drag-line are designated as A respectively to need the material properties of setting to comprise the sectional area A of unit b, A s) and the elastic modulus E (E of unit b, E s, the material of general homogeneous unit is all identical) and the former long L of unit o(L ob, L osif the length of unit is different, respectively the former long matrix L of row input cell o).
1.3, the rigidity S of computing unit: the rigidity of a kth unit is according to formula S k=E ka k/ L okcalculate.
1.4, the connection matrix C of unit: can according to from down to up for starting point and terminal, counterclockwise mode specifies the order of connection of unit, and starting point is 1, and terminal is-1.
Step 2, stochastic generation initial configuration, can be divided into the following steps to carry out:
2.1, stochastic generation initial configuration P:p is the coordinates matrix of a node, and dimension is n × 3, and the coordinate vector matrix g of so each unit is g=Cp.
2.2, the ENERGY E of configuration is calculated: gross energy is according to formula calculate, wherein, S k, l k, l okrepresent a kth element stiffness respectively, length and former length.
2.3, node random file radius r is set: the half being chosen to be the longest unit in structure.
The new configuration of step 3, stochastic generation also judges whether to accept, and the following steps can be divided into carry out:
3.1, fix a node (as first node or last node) and look for shape result to improve counting yield and not affect.
3.2, random selected a certain node, with selected node for the centre of sphere, r is the reposition of interior random this node of setting of the ball of radius, thus generates new configuration P '.
3.3, judged whether that drag-line occurs lax, if having, repeated step 3.1, if without cable slack, carry out 3.4.
3.4, the ENERGY E of new configuration P ' is calculated ', if E ' < E, accept new configuration, and replace original configuration and P=P ', E=E with new configuration '; If E ' >=E, does not replace, and repeat step 3.1.
If 3.5 iterationses reach the integral multiple of 100, then calculate the ratio accepting new configuration in this 100 step, if ratio is lower than 50%, then reduced radius makes r=0.95r, if ratio is not less than 50%, then radius is constant.
Whether step 4, check system balance.
Whether the balance of system is balanced by node and judges, the internal force of each unit is f k=S k(l k-l ok), then force density q kfor if meet then system balancing.Balance equation can be rewritten as C according to connection matrix C tdiag (q) CP=0.
Within the scope of 10 ~ 100N, set that one is less according to scantling and character and numerical value is suitable limits of error ε obtains approximate equilibrium condition, can by limiting the maximal value F of out-of-balance force matrix in result of calculation on the left of balance equation equation maxjudge whether system balances, if F max< ε then carries out step 5, if F max> ε then carries out step 3.
Step 5, obtain looking for shape result.
Export each node coordinate, draw tension integral structure self-equilibrating configuration according to coordinate, obtain final self-equilibrating configuration, and obtain the situation of change of other process variable in computation process.
Embodiment two: present embodiments provide for the MATLAB programming realization that a kind of T prism tension integral structure based on Monte Carlo method looks for shape method, comprise the steps:
Step one, setting starting condition, can be divided into the following steps to carry out:
1.1, T prism geometric parameter setting, only provides unit number setting parameter statement.
m=12。
1.2, material properties: set respectively the material properties of depression bar and drag-line, only provides the material properties setting program section of depression bar.
Ab=200;
Eb=2.06e+5;
Lob=2.7806e+3。
1.3, the rigidity S of computing unit, only provides depression bar element stiffness component (matrix the first row).
Sb=Eb*Ab/Lob; The element stiffness of % depression bar
Sr1=[Sb; Sb; Sb]; % depression bar stiffness matrix
1.4, the connection matrix C (the capable n row of m) of unit, front eb=3 behavior depression bar, rear es=9 behavior drag-line (only providing the first row depression bar L1 connection matrix):
Cr1=[1,0,0,0,0 ,-1]; % depression bar L1, node 1 is starting point 6 is terminal
Step 2, stochastic generation initial configuration, can be divided into the following steps to carry out:
2.1, stochastic generation initial configuration P.
P=randi ([-2000,2000], n, 3); % [imin, imax] uniform random number group
g=C*P:
2.2, the ENERGY E of configuration is calculated.
E=1/2*sum (S.* (L '-Lo) .^2); The gross energy of % system initial configuration
2.3, node random file radius r is set.
R=max (L)/2; % is set as the half of the longest unit
The new configuration of step 3, stochastic generation also judges whether to accept, and the following steps can be divided into carry out:
3.1, a node is fixed.
J=randi ([1, n-1]); % is selected node 1 ~ (n-1) at random, and node n fixes
3.2, random selected a certain node, with selected node for the centre of sphere, r is the reposition of interior random this node of setting of the ball of radius, thus generates new configuration T, only provides the statement of setting X-coordinate.
T (j, 1)=2*r* (rand (1,1)-0.5)+P (j, 1); % sets X-coordinate at random
3.3, judged whether that drag-line occurs lax.
ifmin(L(eb+1:m)′-Lo(eb+1:m))<0
3.4, the ENERGY E of new configuration P ' is calculated ', and determine whether accept new configuration, provide the program segment accepting new configuration.
W=1/2*sum(S.*(L′-Lo).^2);
ifW<E
E=W;
P=T;
3.5, calculate the ratio accepting new configuration in 100 steps, according to the ratio of acceptance judge whether change radius, only provide the ratio of acceptance lower than 50% program segment.
ifk/100<0.5
r=0.95*r;
end
Whether step 4, check system balance.
Step 5, obtain looking for shape result, and draw configuration with line order.
Embodiment three: present embodiment for example with T prism, gives process and result that the tension integral structure based on Monte Carlo method looks for shape to calculate, comprises the steps:
Step one, setting starting condition.
T prism geometric parameter sets: unit number is m=12, e b=3 depression bars, e s=9 drag-lines, nodes n=6.The material properties setting of pull bar and depression bar.
Step 2, stochastic generation initial configuration.
Stochastic generation initial configuration P, as shown in Figure 3, and calculates the ENERGY E of configuration.Often to walk in iterative computation all calculate obtain the energy of configuration to understand its variation tendency, as Fig. 2.
The new configuration of step 3, stochastic generation also judges whether to accept.
Calculate the ENERGY E of new configuration P ' ', and determine whether accept new configuration.T prism looks for the configuration that obtains after 10000 steps and 30000 steps calculate in shape computation process as illustrated in figures 4-5.Calculate the ratio accepting new configuration in every 100 steps, judge whether to change radius.
Whether step 4, check system balance.
Step 5, obtain looking for shape result.What T prism looked for shape to calculate meets the final configuration of approximate equilibrium condition setting as shown in Figure 6.

Claims (5)

1. the tension integral structure based on Monte Carlo method looks for a shape method, it is characterized in that described method step is as follows:
Step one, setting starting condition:
(1) remember that unit adds up to m, nodes is n, and depression bar number is e b, drag-line number is e s, wherein e b+ e s=m;
(2) material properties: the former long L of the sectional area A of unit, the elastic modulus E of unit and unit o;
(3) the rigidity S of computing unit;
(4) the connection matrix C of unit: for starting point and terminal according to from down to up, counterclockwise mode specifies the order of connection of unit, and starting point is 1, and terminal is-1;
Step 2, stochastic generation initial configuration:
(1) stochastic generation initial configuration P:p is the coordinates matrix of each node, and dimension is n × 3, then the coordinate vector matrix g of each unit is: g=Cp;
(2) ENERGY E of configuration is calculated;
(3) node random file radius r is set;
The new configuration of step 3, stochastic generation also judges whether to accept:
(1) node is fixed;
(2) random selected a certain node, with selected node for the centre of sphere, the node random file r set in step 2 is the reposition of interior random this node of setting of the ball of radius, thus generates new configuration P ';
(3) judged whether that drag-line occurs lax, if having, repeated step (1), if without cable slack, carry out (4);
(4) ENERGY E of new configuration P ' is calculated ', if E ' <E, accept new configuration, and replace original configuration and P=P ', E=E with new configuration '; If E ' >=E, does not replace, and repeat step (1);
(5) if iterations reaches the integral multiple of 100, then calculate the ratio accepting new configuration in this 100 step, if ratio is lower than 50%, then reduced radius makes r=0.95r, if ratio is not less than 50%, then radius is constant;
Whether step 4, check system balance:
Whether the balance of system is balanced by node and judges, the internal force of each unit is f k=S k(l k-l ok), wherein, S k, l k, l okrepresent a kth element stiffness, length and former length, then force density q respectively kfor if meet then system balancing, is rewritten as C according to connection matrix C balance equation tdiag (q) CP=0;
Set a limits of error ε and obtain approximate equilibrium condition, by limiting the maximal value F of out-of-balance force matrix in result of calculation on the left of balance equation equation maxjudge whether system balances, if F max< ε then carries out step 5, if F max> ε then carries out step 3;
Step 5, obtain looking for shape result:
Export each node coordinate, draw tension integral structure self-equilibrating configuration according to coordinate, obtain final self-equilibrating configuration, and obtain the situation of change of other process variable in computation process.
2. the tension integral structure based on Monte Carlo method according to claim 1 looks for shape method, and when it is characterized in that the rigidity S of described computing unit, the rigidity of a kth unit is according to formula S k=E ka k/ L okcalculate.
3. the tension integral structure based on Monte Carlo method according to claim 1 looks for shape method, it is characterized in that described ENERGY E is according to formula calculate.
4. the tension integral structure based on Monte Carlo method according to claim 1 looks for shape method, it is characterized in that described displacement radius r is the half of the longest unit in structure.
5. the tension integral structure based on Monte Carlo method according to claim 1 looks for shape method, it is characterized in that described limits of error ε is within the scope of 10 ~ 100N.
CN201510573837.XA 2015-09-10 2015-09-10 A kind of tension integral structure based on Monte Carlo method looks for shape method Expired - Fee Related CN105243185B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510573837.XA CN105243185B (en) 2015-09-10 2015-09-10 A kind of tension integral structure based on Monte Carlo method looks for shape method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510573837.XA CN105243185B (en) 2015-09-10 2015-09-10 A kind of tension integral structure based on Monte Carlo method looks for shape method

Publications (2)

Publication Number Publication Date
CN105243185A true CN105243185A (en) 2016-01-13
CN105243185B CN105243185B (en) 2018-06-26

Family

ID=55040833

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510573837.XA Expired - Fee Related CN105243185B (en) 2015-09-10 2015-09-10 A kind of tension integral structure based on Monte Carlo method looks for shape method

Country Status (1)

Country Link
CN (1) CN105243185B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106354906A (en) * 2016-08-19 2017-01-25 西安电子科技大学 Mesh-shaped extendible antenna shape-finding method based on triangular patch errors
CN106934155A (en) * 2017-03-13 2017-07-07 北京市建筑设计研究院有限公司 A kind of cable-truss structure looks for shape method
CN107194080A (en) * 2017-05-25 2017-09-22 同济大学建筑设计研究院(集团)有限公司 A kind of structure determination methodology of radial tensioning system
CN108229054A (en) * 2018-01-26 2018-06-29 东南大学 A kind of symmetrical tension integral structure based on the group theory looks for shape method
CN108345729A (en) * 2018-01-26 2018-07-31 东南大学 A kind of symmetrical cable-rod structure based on the group theory looks for shape method
CN110792173A (en) * 2019-10-28 2020-02-14 中冶建筑研究总院有限公司 Tensioning integral triangular prism structure and shape finding method
US11413553B1 (en) 2021-02-08 2022-08-16 Marcio Sequeira De Oliveira Structural parts kit for formation of architectural and structural didactic models

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101582095A (en) * 2009-06-19 2009-11-18 东南大学 Non-linear dynamic finite element method for determining cable-strut system static balancing state
CN102622485A (en) * 2012-03-16 2012-08-01 东南大学 Morphological analysis inverse iteration method for determining space structural lofting states of cable-strut-beam system
CN104765925A (en) * 2015-04-14 2015-07-08 温州大学 Cable-strut pretension structure stretch forming optimal selection method based on error sensitivity analysis

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101582095A (en) * 2009-06-19 2009-11-18 东南大学 Non-linear dynamic finite element method for determining cable-strut system static balancing state
CN102622485A (en) * 2012-03-16 2012-08-01 东南大学 Morphological analysis inverse iteration method for determining space structural lofting states of cable-strut-beam system
CN104765925A (en) * 2015-04-14 2015-07-08 温州大学 Cable-strut pretension structure stretch forming optimal selection method based on error sensitivity analysis

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李悦 等: "张拉整体结构的蒙特卡罗找形方法", 《损伤、断裂与微纳米力学学术研讨会》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106354906A (en) * 2016-08-19 2017-01-25 西安电子科技大学 Mesh-shaped extendible antenna shape-finding method based on triangular patch errors
CN106354906B (en) * 2016-08-19 2019-08-02 西安电子科技大学 A kind of netted deployable antenna based on triangle surface error looks for shape method
CN106934155A (en) * 2017-03-13 2017-07-07 北京市建筑设计研究院有限公司 A kind of cable-truss structure looks for shape method
CN106934155B (en) * 2017-03-13 2020-04-17 北京市建筑设计研究院有限公司 Shape-finding method of cable truss structure
CN107194080A (en) * 2017-05-25 2017-09-22 同济大学建筑设计研究院(集团)有限公司 A kind of structure determination methodology of radial tensioning system
CN107194080B (en) * 2017-05-25 2024-04-05 同济大学建筑设计研究院(集团)有限公司 Structure determining method of spoke type tensioning system
CN108229054A (en) * 2018-01-26 2018-06-29 东南大学 A kind of symmetrical tension integral structure based on the group theory looks for shape method
CN108345729A (en) * 2018-01-26 2018-07-31 东南大学 A kind of symmetrical cable-rod structure based on the group theory looks for shape method
CN108345729B (en) * 2018-01-26 2021-07-27 东南大学 Symmetrical cable-strut structure shape finding method based on group theory
CN110792173A (en) * 2019-10-28 2020-02-14 中冶建筑研究总院有限公司 Tensioning integral triangular prism structure and shape finding method
CN110792173B (en) * 2019-10-28 2020-12-01 中冶建筑研究总院有限公司 Tensioning integral triangular prism structure and shape finding method
US11413553B1 (en) 2021-02-08 2022-08-16 Marcio Sequeira De Oliveira Structural parts kit for formation of architectural and structural didactic models

Also Published As

Publication number Publication date
CN105243185B (en) 2018-06-26

Similar Documents

Publication Publication Date Title
CN105243185A (en) Monte Carlo method based tensegrity structure form-finding method
Tran et al. An edge‐based smoothed finite element method for primal–dual shakedown analysis of structures
Zhang et al. Analytical form-finding of tensegrities using determinant of force-density matrix
CN106503368A (en) A kind of multi-state System Reliability modeling method related based on failure mechanism
CN104182209B (en) A kind of GCRO DR algorithm method for parallel processing based on PETSc
Castonguay et al. Application of high-order energy stable flux reconstruction schemes to the Euler equations
Chen et al. A two-layered parallel static security assessment for large-scale grids based on GPU
CN103246541A (en) Method for evaluating auto-parallelization and multistage parallelization cost
Moxey et al. Optimising the performance of the spectral/hp element method with collective linear algebra operations
Wu et al. Robust fuzzy structural safety assessment using mathematical programming approach
Pham et al. Optimal truss sizing by modified Rao algorithm combined with feasible boundary search method
Pourbakhshian et al. Shape optimization of arch dams using sensitivity analysis
Li et al. Fast time integration of Navier-Stokes equations with an exponential-integrator scheme
Zavala et al. Integrating a multi-objective optimization framework into a structural design software
CN105608267A (en) Multivariable global optimization algorithm
CN105701291B (en) Finite element fraction analysis apparatus and information acquisition method, sytem matrix parallel generation method
Huang et al. Parallel performance studies for a three-species application problem on the cluster maya
CN105550479B (en) A kind of load-bearing part stability loading prediction method considering kinematic hardening behavior
Huang et al. Parallel performance studies for a three-species application problem on maya 2013
Shoid et al. A new conjugate gradient method with exact line search
Shen et al. A principal component analysis-based dimension reduction method for parametric power flow
Djinevski et al. Network simulator tools and GPU parallel systems
Ahrari et al. Efficient simulation for optimization of topology, shape and size of modular truss structures
Li et al. Prestress Iterative Design and Precision Control Method for Cable-strut Structures using Enhanced Differential Evolution
Shabash et al. Improving the portability and performance of jViz. RNA–a dynamic RNA visualization software

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Zheng Zhiwei

Inventor after: He Jingfeng

Inventor after: Zhao Jianying

Inventor after: Wang Yang

Inventor before: Zheng Zhiwei

Inventor before: Zhao Jianying

Inventor before: Wang Yang

GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180626

Termination date: 20200910