CN111651840A - Method for detecting arch state on track slab based on deep learning technology - Google Patents

Method for detecting arch state on track slab based on deep learning technology Download PDF

Info

Publication number
CN111651840A
CN111651840A CN202010275388.1A CN202010275388A CN111651840A CN 111651840 A CN111651840 A CN 111651840A CN 202010275388 A CN202010275388 A CN 202010275388A CN 111651840 A CN111651840 A CN 111651840A
Authority
CN
China
Prior art keywords
track slab
track
slab
data
mortar
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
CN202010275388.1A
Other languages
Chinese (zh)
Other versions
CN111651840B (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.)
Beijing Jiaotong University
China State Railway Group Co Ltd
Original Assignee
Beijing Jiaotong University
China State Railway Group Co Ltd
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 Beijing Jiaotong University, China State Railway Group Co Ltd filed Critical Beijing Jiaotong University
Priority to CN202010275388.1A priority Critical patent/CN111651840B/en
Publication of CN111651840A publication Critical patent/CN111651840A/en
Application granted granted Critical
Publication of CN111651840B publication Critical patent/CN111651840B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B61RAILWAYS
    • B61KAUXILIARY EQUIPMENT SPECIALLY ADAPTED FOR RAILWAYS, NOT OTHERWISE PROVIDED FOR
    • B61K9/00Railway vehicle profile gauges; Detecting or indicating overheating of components; Apparatus on locomotives or cars to indicate bad track sections; General design of track recording vehicles
    • B61K9/08Measuring installations for surveying permanent way
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/241Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2415Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Machines For Laying And Maintaining Railways (AREA)

Abstract

The invention belongs to the technical field of track slab arch-up detection, and particularly relates to a track slab arch-up state detection method based on a deep learning technology. The method of the invention converts the classification problem into a time sequence classification problem, namely inputting the time sequence classification problem into a series of track slab displacement data signal segments and outputting the time sequence classification problem into a state classification corresponding to the track slab; aiming at the characteristics of the identification task of the arch state on the track slab, the identification process is divided into two stages: a characteristic extraction stage and a classification stage; considering that the noise of the original signal is more, particularly the signal including the track slab irregularity runs through all the time, the electroencephalogram signal is represented by multiple types of features and divided into two stages to simplify the calculation amount of the classification process and improve the calculation speed.

Description

Method for detecting arch state on track slab based on deep learning technology
Technical Field
The invention belongs to the technical field of track slab arch-up detection, and particularly relates to a track slab arch-up state detection method based on a deep learning technology.
Background
The novel high-speed track structure CRTS II type slab ballastless track is developed by studying and innovating after the introduction of the German Bo lattice slab ballastless track. The Jingjin high-speed railway is the first high-speed railway with the speed per hour of 300-350 km/h in China, the total length of a high-speed line is 1318km, the high-speed railway is the highest-standard high-speed railway with the longest construction mileage, the largest investment and the highest standard after the establishment of new China, ballastless tracks are laid on the whole line, and the CRTS II type plate ballastless track technology is adopted. At present, CRTS II plate-type ballastless tracks are already used in a plurality of passenger special lines such as Jingjin, Jinghu, HuHangzhou and the like. The CRTS II type plate ballastless track on the roadbed mainly comprises components such as a steel rail, a fastener system, a track plate, a Cement Asphalt Mortar (CA Mortar for short) layer, a concrete supporting layer and the like. As the CRTS II type plate ballastless track is influenced by factors such as train impact, environment temperature and the like, various diseases inevitably occur to the structure below the ballastless track, wherein the occurrence of gap of the plate ballastless track and the upward arching of the track plate are typical diseases.
At present, the field detection and maintenance of the open joint and the upper arch of the slab ballastless track mainly adopts methods of visual inspection, steel ruler insertion measurement, field slab uncovering and track inspection. The visual inspection and steel ruler insertion measurement method has the defects that the accurate distribution conditions of the middle local open joints and the open joints cannot be detected; the plate uncovering method has the defects that the plate uncovering method is only suitable for building a railway, and has high cost and low efficiency; the rail inspection vehicle has the defects of high manufacturing cost, detection only in the empty window period of train operation and incapability of monitoring the real-time state.
Patent CN201910620162.8 discloses a track slab arching distributed monitoring system and a monitoring method, which detects an arching angle of a track slab in real time, and filters an angle value detected when a train passes through by processing the arching angle, so as to retain the angle value measured when the track slab is static; the track slab arching distributed monitoring system and the monitoring method have the following defects:
1) the existing algorithm can only realize the measurement of the angle value measured when the road plate is static, and the application scene is single;
2) when a train passes by, the data at the moment can be directly filtered by the existing algorithm, and the acquired original data is missing, so that the later data analysis is not facilitated;
3) the track slab arching distributed monitoring method and system adopt node output to process data, and have the disadvantages of high circuit power, high power consumption and reduced node service life.
Disclosure of Invention
Aiming at the technical problem, the invention provides a method for detecting the arch state on a track slab based on a deep learning technology. And converting the classification problem into a time series classification problem, namely inputting the time series classification problem into a series of track slab displacement data signal segments and outputting the time series classification problem into the state classification corresponding to the track slab.
The invention is realized by the following technical scheme:
a method for detecting the arch state on a track slab based on a deep learning technology comprises the following steps:
(1) the traditional vehicle-plate type track dynamics model is improved into a vehicle-plate type track dynamics model under the condition of considering the CA mortar void effect;
(2) constructing a database for deep learning of the classified neural network: outputting track slab displacement simulation data under different voiding degrees by setting different parameters by using the vehicle-slab track dynamics model under the voiding effect of the CA mortar obtained in the step (1), and after training for multiple times to obtain enough simulation data, adding labels to the simulation data according to disease types to form a database, wherein the database comprises a training set for inputting a classification neural network and a test set for inputting the classification neural network;
(3) training a classification network by using the database constructed in the step (2): the classified network structure comprises four layers of networks, wherein the first layer of network is an input layer, the second layer of network is a BilSTM, the third layer of network is a full-connection layer, and the fourth layer of network is a softmax layer;
(4) and (3) carrying out arch state classification on the track slab according to the classification network result: taking the track slab displacement data as the network input quantity of the trained classification network, and sorting the original track slab displacement data into a characteristic sequence to be input into a BilSTM network; and after circulation of the BilSTM layer, inputting the circular shape into a softmax layer to judge the type of the arch state on the track slab, and giving a final result of judgment of the arch state on the track slab.
Further, in the step (4), the original track slab displacement data is sorted into a feature sequence, and the specific method is as follows:
(1) preprocessing a track slab displacement signal: expressing track slab displacement signals by using various types of characteristics, normalizing time sequence signals to a (0,1) range, segmenting original track slab displacement signals, and extracting signal characteristics, wherein the signal characteristic extraction is to calculate data characteristics in small segments of data obtained by segmentation and then rearrange the data characteristics into a group of data to form characteristic data; after signal characteristics are extracted, connecting the characteristic data of each section in sequence to form a final characteristic data group;
the data features in the signal feature extraction process comprise a maximum value, a minimum value, an average value, a peak-to-peak value, a rectified mean value, a variance, a standard deviation, a kurtosis, a root mean square, a form factor, a peak factor, a kurtosis factor, a pulse factor and a margin factor;
(2) and calculating corresponding result relevance weight for each feature of the obtained feature data set by utilizing a Relief algorithm, and eliminating the features with lower weight according to a set weight threshold value to obtain a final trained feature sequence.
Further, step (4) is followed by step (5): and performing post-processing according to the obtained final result of the judgment of the arch-up state of the track slab, wherein the post-processing is to judge whether the track slab is arched up by combining other judgment programs, and the classification result can correct obvious identification errors, further optimize the identification effect and perform early warning work.
Further, in the step (1), the detailed form of the vibration type coordinate differential equation set of the vehicle-plate type track dynamics model under the consideration of CA mortar void action is as follows:
Figure RE-GDA0002596184080000041
in the formula: t isn(t) introducing a free beam orthogonal function system (X) after a Ritz method is adopted in a vertical vibration differential equation of the track slabnNMS is the mode order of the track plate, and NMS generalized coordinates T are selectedn(t);
Figure RE-GDA0002596184080000042
Is Tn(t) a second derivative; esIsBending rigidity of the track slab; m issβ being the unit length mass of the track slabnIs a constant; m is0Finger m0A plurality of discrete units; csqIs the distributed damping of the CA mortar at q; xpIs a free beam orthogonal function system { X }nMeaning of, p is similar to n, the range is the same as n, XnIs a free beam orthogonal function system, X, of the traditional track slabpThe free beam orthogonal function system of the improved track slab is provided; t isp(t) and XpIn the same way, TnP of (t) is like n, takenThe value range is the same as n, the band n represents the generalized coordinate of the traditional track slab, and the band p represents the generalized coordinate of the improved track slab; ksqThe distributed rigidity of the CA mortar at the position q; f corresponds to the concrete supporting layer, zf(x, t) is the vibration displacement (m), z of the concrete supporting layerf(xqAnd t) is the vibration displacement of the q CA mortar discrete unit concrete supporting layer;
Figure RE-GDA0002596184080000043
is zf(xqT) obtaining a derivation, and representing the vibration speed of the q CA mortar discrete unit concrete supporting layer; l issIs the length of a single track slab; n is0The number of fasteners of a steel rail on one rail plate; subscript s corresponds to the track slab; x is the number ofiThe method is characterized by comprising the following steps of (1) establishing a coordinate system of an original coordinate system on the basis of a steel rail model, wherein N is the number of fasteners in a plate type track length range L; x is the number ofqIs xiTime q ═ i; cpiIs in a coordinate system xiDamping of the under-rail pad layer when in position; NM is the modal order of the steel rail; kpiIs in a coordinate system xiThe stiffness of the under-rail pad in position; y isp(x(j-1)×n0+i)qp(t) expressing the vertical vibration displacement of the steel rail by applying the regular vibration mode function of the simply supported beam by a Ritz method
Figure RE-GDA0002596184080000051
In the formula, k is p and x has a subscript of (j-1) n0+ i wherein n0The number of fasteners of the steel rail on one rail plate is equal to the number of fasteners of the steel rail on the other rail plate.
Further, in the step (2), when constructing the database, selecting 104 meters, that is, 16 track slabs as the length of the vehicle-slab track dynamics model under the consideration of the CA mortar void effect, wherein the vehicle speed is 300 kilometers per hour, selecting the displacement of the track slabs as sample data when a train passes through one track slab, and when the sampling interval is 10-4When s is needed, different sample data are obtained by setting different emptying lengths, CA mortar emptying conditions with the longitudinal lengths of 0 meter, 0.325 meter, 0.65 meter, 1.3 meter and 1.95 meter are set to obtain samples with five emptying conditions, and then, for each sample, the data are acquiredThe number of samples obtained for each of the blanking conditions is sufficiently large to obtain 64 sets of samples.
The invention has the beneficial technical effects that:
(1) the method for detecting the arch state on the track slab based on the deep learning technology converts the classification problem into a time series classification problem, namely, the classification problem is input into a series of track slab displacement data signal segments and output into the state classification corresponding to the track slab.
(2) The invention provides a track slab arch state detection method based on a deep learning technology, aiming at the characteristics of a track slab arch state identification task, the identification process is divided into two stages: a feature extraction stage and a classification stage. Considering that the noise of the original signal is more, especially the signal including the track slab irregularity runs through all the time, the calculation resource needed by the signal training model after the preprocessing is directly used is very large, and the type of the learned information is single, so that the calculation amount of the classification process can be simplified by using various types of features to represent the electroencephalogram signal in two stages, and the calculation speed is improved.
(3) In the method provided by the invention, the data acquisition adopts track slab monitoring equipment laid along the railway, the track slab monitoring equipment has the problems of large environmental interference, irregular installation and the like, and the requirements on the aspects of data accuracy, stability and the like are more strict; in the method provided by the invention, the artificial intelligence technology is utilized, the interference signals such as track irregularity and the like are added during the training of the classification network, and the classification network can continuously learn according to the updating of the data set, so that the high accuracy and the high robustness are achieved in the aspect of identifying the abnormal arch state on the track slab in the displacement signals of the track slab with poor signal quality and large interference, the subsequent labor cost is greatly reduced, and the period for generating the final conclusion is shortened.
Drawings
FIG. 1 is a conventional model of high speed railway vehicle-slab track dynamics in an embodiment of the present invention;
FIG. 2 is a diagram illustrating a stress relationship between a steel rail and a conventional model according to an embodiment of the present invention;
FIG. 3 is a diagram illustrating a stress analysis of a track slab model according to a conventional model in an embodiment of the present invention;
FIG. 4 is a high-speed railway vehicle-plate type track dynamics improved model (i.e. vehicle-plate type track dynamics model under consideration of CA mortar void action) in the embodiment of the invention;
FIG. 5 is a force analysis of an improved track slab model according to an embodiment of the present invention;
FIG. 6 is a BiLSTM neural network structure according to an embodiment of the present invention;
FIG. 7 is a diagram of a classification neural network structure according to an embodiment of the present invention;
FIG. 8 is a feature weight arrangement after a Relief algorithm in an embodiment of the present invention;
FIG. 9 is a feature sequence with low weight removed according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
On the contrary, the invention is intended to cover alternatives, modifications, equivalents and alternatives which may be included within the spirit and scope of the invention as defined by the appended claims. Furthermore, in the following detailed description of the present invention, certain specific details are set forth in order to provide a better understanding of the present invention. It will be apparent to one skilled in the art that the present invention may be practiced without these specific details.
Aiming at the technical problems of the track slab arching distributed monitoring system and monitoring method in the prior art, the invention provides a track slab arching state detection method based on a deep learning technology, which comprises the following steps:
(1) the traditional vehicle-plate type track dynamics model is improved into a vehicle-plate type track dynamics model under the condition of considering the CA mortar void effect;
the traditional vehicle-plate track dynamics model of the high-speed railway is provided by Zhai Ming professor of southwest transportation university in vehicle-track coupling dynamics, and is called the traditional model for short. The conventional model is a vertical model of the vehicle-plate type track coupling system as shown in fig. 1, which is established by using the vehicle-track coupling dynamics principle and integrating a vehicle system and a plate type track system.
The conventional model has the following assumptions:
(1) as the CA mortar under the track slab mainly plays a role in supporting, and the vertical vibration of the ballastless track is mainly influenced by the void of the CA mortar, only the vibration in the vertical direction of the vehicle-slab track coupling system is considered.
(2) The vehicle system and the plate-type track system are symmetrical to the central line of the track, and only half of structural research of the vehicle track coupling system is considered in the calculation process.
(3) The vehicle system is a multi-rigid-body model with a two-series spring damping system and comprises a vehicle body, two bogies and four wheel pairs. A primary suspension system is arranged between the bogie and the wheel pair, and a secondary suspension system is arranged between the vehicle body and the bogie. Considering the ups and downs of the vehicle body and the bogie and the vibration of the point head, the ups and downs of the wheel pair, and the total ten degrees of freedom of the vehicle system.
(4) The steel rail is a simply supported beam with limited length on the basis of discrete elastic point support, and the elasticity and damping of the base plate under the rail and the fastener are respectively determined by the elastic coefficient KpAnd damping coefficient CpAnd (4) showing.
(5) The track slab is a finite-length free beam supported on continuously distributed linear springs and linear dampers, and the distributed rigidity and the distributed damping of the CA mortar respectively use elastic coefficients ksAnd damping coefficient csAnd (4) showing.
(6) The coupling relationship between the vehicle system and the plate-type track system is realized through wheel-rail interaction force, and a classic Hertz nonlinear elastic contact model is adopted.
The physical quantities represented by the symbols in the model are shown in Table 1
TABLE 1 physical quantities represented by symbols in the conventional model
Figure RE-GDA0002596184080000081
The model equation of motion includes:
(1) vehicle body sinking and floating movement
Figure RE-GDA0002596184080000082
zcIs used for sinking and floating of the vehicle body,
Figure RE-GDA0002596184080000083
is zcAfter derivation, the sinking and floating speed of the vehicle body is represented,
Figure RE-GDA0002596184080000084
representing the sinking and floating acceleration of the vehicle body;
Figure RE-GDA0002596184080000091
and zt1Similarly, where subscript t1 denotes the front truck, t2 denotes the rear truck.
(2) Vehicle body nodding motion
Figure RE-GDA0002596184080000092
(3) Sinking and floating movement of front bogie
Figure RE-GDA0002596184080000093
(4) Front bogie nodding motion
Figure RE-GDA0002596184080000094
(5) Sinking and floating movement of rear bogie
Figure RE-GDA0002596184080000095
(6) Rear bogie nodding motion
Figure RE-GDA0002596184080000096
(7) Sinking and floating movement of the first wheel pair
Figure RE-GDA0002596184080000097
(8) Sinking and floating movement of the second wheel pair
Figure RE-GDA0002596184080000098
(9) The third wheel moves up and down
Figure RE-GDA0002596184080000099
(10) Fourth wheel pair sinking and floating movement
Figure RE-GDA00025961840800000910
In the formula, piAnd (t) is the vertical wheel-rail force (i is 1-4) of the single-side wheel.
The steel rail in the traditional model is regarded as a finite-length simply supported beam on a discrete elastic point supporting basis, and the stress analysis model is shown in figure 2. Wherein p isiThe force is wheel-rail force and moves forwards at a speed v along with the vehicle; frsi(i is 1-N) is steel rail fulcrum counter-force, and N is the number of fasteners in the plate track length range L; ox is a fixed coordinate system fixedly connected to the steel rail; and o 'x' is a moving coordinate system connected to the vehicle. The relationship between these two coordinates is:
x=x′+x0+νt (11)
wherein x is0Fixed coordinates of the rear wheel at the starting moment; t is a time variable.
As can be seen from fig. 2, the differential equation of the vibration of the steel rail is:
Figure RE-GDA0002596184080000101
wherein
Figure RE-GDA0002596184080000102
In the formula, subscript "r" corresponds to steel rail, subscript "s" corresponds to track plate, and zr(xiT) and zs(xiT) respectively representing vertical displacement variables (m) of the steel rail and the track plate at the fastener;
Figure RE-GDA0002596184080000103
and
Figure RE-GDA0002596184080000104
respectively, the vertical velocity variables (m/s) of the rail and the rail plate at the fasteners.
Coordinates x of each wheelwj(j 1 to 4) are each:
Figure RE-GDA0002596184080000105
coordinates of each fastener
xi=ilp(15)
By adopting a Ritz method and applying a regular vibration mode function of the simply supported beam, the vertical vibration displacement of the steel rail can be expressed as:
Figure RE-GDA0002596184080000106
wherein, NM is the mode order of the steel rail, and is generally 0.5N; vibration type of rail
Figure RE-GDA0002596184080000107
Substituting formula (16) for formula (12) to obtain:
Figure RE-GDA0002596184080000111
both sides of the above formula are multiplied by Yh(x) (h ═ 1,2,3,. and, NM), integrate x from 0 to L, and note the orthogonality of the modes:
Figure RE-GDA0002596184080000112
is provided with
Figure RE-GDA0002596184080000113
Depending on the nature of the function, equation (20) can be organized as:
Figure RE-GDA0002596184080000114
because of the fact that
Figure RE-GDA0002596184080000115
Figure RE-GDA0002596184080000116
Therefore, equation (21) can be simplified as:
Figure RE-GDA0002596184080000117
the method is a basic form of a steel rail vibration type coordinate second-order ordinary differential equation set.
Further, formula (16) is substituted for formula (13)
Figure RE-GDA0002596184080000118
Then, the formula (24) becomes
Figure RE-GDA0002596184080000121
This is the detailed form of the differential equation set of the vibration mode coordinate of the steel rail.
The slab track model simplifies the CA mortar into springs and dampers that are continuously distributed along the track slab, which is seen as a finite free beam supported on continuously distributed linear springs and linear dampers as shown in fig. 3.
The differential equation of the vertical vibration of the track slab is as follows:
Figure RE-GDA0002596184080000122
in the formula, ks、csRespectively distributing rigidity (N/m/m) and damping (N.s/m/m) of a CA mortar layer under the track slab along the track direction; n is0The number of fasteners of the steel rail on one rail plate is increased.
Introducing a free beam orthogonal function system { X ] by a Ritz methodn1-NMS, selecting NMS generalized coordinates Tn(t) in which
Figure RE-GDA0002596184080000123
In the formula, Cn、βnIs a constant. Cn、βnLsThe values of (A) are shown in Table 2.
TABLE 2 free Beam function coefficients
Figure RE-GDA0002596184080000124
The vertical displacement of the track plate can be approximated as:
Figure RE-GDA0002596184080000125
wherein NMS is the mode order of the track plate, and NMS is 0.5n0Substituting formula (29) for formula (27) to obtain:
Figure RE-GDA0002596184080000131
both sides of the above formula are multiplied by Xp(x) (p 1-NMS), then L is in the whole length of the track platesIntegrate x and note the orthogonality of the modes
Figure RE-GDA0002596184080000132
Is provided with
Figure RE-GDA0002596184080000133
Depending on the nature of the function, equation (32) can be organized as:
Figure RE-GDA0002596184080000134
because of the fact that
Figure RE-GDA0002596184080000135
Figure RE-GDA0002596184080000136
Therefore, equation (33) can be simplified as:
Figure RE-GDA0002596184080000137
this is the dynamic equation of the track slab model.
Further, formula (29) is substituted for formula (25)
Figure RE-GDA0002596184080000141
Then, equation (36) becomes:
Figure RE-GDA0002596184080000142
this is the detailed form of the differential equation set of the track slab mode coordinate.
In the formula: in the formula: t isn(t) introducing a free beam orthogonal function system { X ] after a Ritz method is adopted in a track slab vertical vibration differential equation (27)n}(n=1-NMS), NMS is the mode order of the track plate, and NMS generalized coordinates T are selectedn(t);
Figure RE-GDA0002596184080000143
Is Tn(t) a second derivative; esIsBending rigidity of the track slab; m issβ being the unit length mass of the track slabnIs a constant; m is0Finger m0A plurality of discrete units; csqIs the distributed damping of the CA mortar at q; xpIs a free beam orthogonal function system { X }nMeaning of, p is similar to n, the range is the same as n, XnIs a free beam orthogonal function system, X, of the traditional track slabpThe free beam orthogonal function system of the improved track slab is provided; t isp(t) and XpIn the same way, TnP of (t) is similar to n, the value range is the same as n, the n represents the generalized coordinate of the traditional track slab, and the p represents the generalized coordinate of the improved track slab; ksqThe distributed rigidity of the CA mortar at the position q; f corresponds to the concrete supporting layer, zf(x, t) is the vibration displacement (m), z of the concrete supporting layerf(xqAnd t) is the vibration displacement of the q CA mortar discrete unit concrete supporting layer;
Figure RE-GDA0002596184080000144
is zf(xqT) obtaining a derivation, and representing the vibration speed of the q CA mortar discrete unit concrete supporting layer; l issIs the length of a single track slab; n is0The number of fasteners of a steel rail on one rail plate; subscript s corresponds to the track slab; x is the number ofiThe method is characterized by comprising the following steps of (1) establishing a coordinate system of an original coordinate system on the basis of a steel rail model, wherein N is the number of fasteners in a plate type track length range L; x is the number ofqIs xiTime q ═ i; cpiIs in a coordinate system xiDamping of the under-rail pad layer when in position; NM is the modal order of the steel rail; kpiIs in a coordinate system xiThe stiffness of the under-rail pad in position; y isp(x(j-1)×n0+i)qp(t) expressing the vertical vibration displacement of the steel rail by applying the regular vibration mode function of the simply supported beam by a Ritz methodIn the formula, k is p and x has a subscript of (j-1) n0+ i wherein n0The number of fasteners of a steel rail on one rail plate; xp(xi)Tp(t) is obtained by substituting formula (29) into (42), wherein formula (29) is obtained by introducing a free beam orthogonal function system { X ] by a Ritz methodn1-NMS, selecting NMS generalized coordinates Tn(t) in which
Figure RE-GDA0002596184080000152
In the formula, Cn、βnIs a constant. Cn、βnLsThe values of (A) are shown in Table 2. The vertical displacement of the track plate can be approximated as:
Figure RE-GDA0002596184080000153
in a ballastless track dynamics equation of a conventional vehicle-plate track dynamics model, CA mortar is simplified into springs and damping which are continuously distributed along a track plate, and acting force on the track plate is embodied in the form of distributed force in the dynamics equation, so that the void of a CA mortar layer with the length of the whole track plate can be simulated only by changing the rigidity damping value of the CA mortar under the bottom of the whole track plate, and an arbitrary range cannot be simulated. In the improved model provided by the invention, a CA mortar model is discretized in the longitudinal direction, distributed spring damping is replaced by the discrete spring damping, the supporting distribution force of the CA mortar on a track plate is converted into concentrated force, and a vehicle-plate type track dynamic model considering the CA mortar void effect is established, as shown in figure 4.
In the vehicle-plate type track dynamic model considering the CA mortar void effect, provided by the invention, the CA mortar corresponding to one track plate is averagely divided into m0The CA mortar of each discrete unit is simplified into a spring and a damper which are concentrated to one point, and the elastic coefficients K are respectively usedsAnd damping coefficient CsThe longitudinal length of each unit is expressed as ls=Ls/m0And under the condition that the CA mortar is not empty, the rigidity of the bearing surface is as follows:
Figure RE-GDA0002596184080000161
ECAis the modulus of elasticity of the CA mortar; h isCAThe thickness of the CA mortar layer (CRTS II type plate ballastless track is 0.03m), and the vertical stiffness K of each discrete spring and the damping of the CA mortar model under the working condition of no voidsAnd damping coefficient CsRespectively as follows:
Figure RE-GDA0002596184080000162
in the formula,
Figure RE-GDA0002596184080000163
bearing surface stiffness, b0The overall width of the track slab (2.55 m for a CRTS type II track slab); c. CsDistribution of damping for CA mortar,. lsA longitudinal length for each discrete unit;
after the CA mortar model is discretized, the distributed force to the track slab is converted into concentrated force in the track slab vibration differential equation, as shown in fig. 5, in the improved track slab model, the track slab is regarded as a finite-length free beam on the basis of the discrete elastic point support, and the differential equation is:
Figure RE-GDA0002596184080000164
wherein,
Figure RE-GDA0002596184080000165
in the formula, the subscript f corresponds to a concrete supporting layer; z is a radical off(x, t) is the vibrational displacement (m) of the concrete support layer; fsfq(t)(q=1~m0) The supporting force (N) of the q & ltth & gt CA mortar discrete unit to the track slab is shown; z is a radical ofs(xqT) is the vertical displacement variable of the track plate at q、KsqThe distributed rigidity of the CA mortar at q, CsqIs the distributed damping of the CA mortar at q;
wherein E issIsRepresenting the bending rigidity of the track slab; z is a radical ofs(x, t) vertical displacement variation of the track plate at the fastener; x, t are variables; m issThe mass of the track slab in unit length; m is0Is m0A plurality of discrete units; q is the q-th CA mortar discrete unit concrete supporting layer; (x-x)q) Is a differential variable; n is0The number of fasteners of a steel rail on one rail plate; frsj(t) the supporting force of the jth CA mortar discrete unit concrete supporting layer; (x-x)i) Is a differential variable.
Introducing a free beam orthogonal function system { X ] by a Ritz methodn1-NMS, selecting NMS generalized coordinates Tn(t),{XnThe values are given in table 2.
Substituting equation (29) into equation (41), and multiplying both sides of the equation by Xp(x) (p 1-NMS), and then in the whole length range L of the track slabsInner pair x is integrated:
Figure RE-GDA0002596184080000171
using the modal orthogonality and the nature of the function one can derive:
Figure RE-GDA0002596184080000172
substituting formula (34) and formula (35) for formula (44):
Figure RE-GDA0002596184080000173
this is the dynamic equation of the improved model of the track slab.
Further, formula (29) is substituted for formula (42):
Figure RE-GDA0002596184080000174
substituting formula (37) and formula (46) for formula (45):
Figure RE-GDA0002596184080000175
this is the detailed form of the improved track slab model mode shape coordinate differential equation set.
After the CA mortar is discretized, the distribution force of the CA mortar on the track plate is converted into a concentrated force, a train track coupling dynamic model considering the CA mortar void effect can be established, and the problem that the existing dynamic model cannot simulate CA mortar void under any working condition can be solved; because the modeling is based on matlab software modeling, the practical problems of low simulation efficiency, limited line length, short track irregularity wavelength and the like of a finite element model method can be solved.
(2) Constructing a database for deep learning of the classified neural network: outputting track slab displacement simulation data under different voiding degrees by setting different parameters by using the vehicle-slab track dynamics model under the voiding effect of the CA mortar obtained in the step (1), and after training for multiple times to obtain enough simulation data, adding labels to the simulation data according to disease types to form a database, wherein the database comprises a training set for inputting a classification neural network and a test set for inputting the classification neural network;
when constructing the database, in this embodiment, 104 meters, that is, 16 track slabs are selected as the length of the vehicle-slab track dynamics model under consideration of CA mortar void, the vehicle speed is 300 km/h, the displacement of the track slab is selected as sample data when a train passes through one track slab, and when the sampling interval is 10-4When s is carried out, different sample data are obtained by setting different void lengths, and CA mortar voids with longitudinal lengths of 0 meter, 0.325 meter, 0.65 meter, 1.3 meter and 1.95 meter are respectively arranged under the track slab to obtain samples under five conditions; then a sufficient number of samples are obtained for each of the conditions set for the void, obtaining 64 sets of samples.
(3) Training a classification network by using the database constructed in the step (2): the classified network structure comprises four layers of networks, wherein the first layer of network is an input layer, the second layer of network is a BilSTM, the third layer of network is a full-connection layer, and the fourth layer of network is a softmax layer; both first and second tier networks are commonly used to model context information in natural language processing tasks, where the full name of LSTM is Long Short-Term Memory, which is one of rnn (current Neural network). LSTM is well suited for modeling time series data, such as text data, due to its design features. BilSTM is an abbreviation of Bi-directional Long Short-Term Memory, and is formed by combining forward LSTM and backward LSTM, and the structure of the BilSTM is shown in figure 6;
(4) and (3) carrying out arch state classification on the track slab according to the classification network result: taking the track slab displacement data as the network input quantity of the trained classification network, and sorting the original track slab displacement data into a characteristic sequence to be input into a BilSTM network; and after circulation of the BilSTM layer, inputting the circular shape into a softmax layer to judge the type of the arch state on the track slab, and giving a final result of judgment of the arch state on the track slab.
The classification network structure is shown in fig. 7, and specifically, in this step, the rail plate displacement data as the network input quantity of the trained classification network may be collected by the method provided in patent CN201910620162.8, where the collection method specifically is to collect rail plate vibration displacement data by using the inclination angle sensing nodes installed on the rail plate, and then convert the rail plate vibration displacement data into rail plate inclination angle data for use in the present invention.
Specifically, in the step (4), the original track slab displacement data is sorted into a feature sequence, and the specific method is as follows:
(1) preprocessing a track slab displacement signal: the method comprises the steps of representing track slab displacement signals by using various types of characteristics, normalizing time sequence signals to a (0,1) range, segmenting original track slab displacement signals and extracting signal characteristics, wherein the original track slab displacement signals can be generally divided into three segments, if the signals are too long, the number of the segments can be increased, and the signal characteristic extraction is to calculate data characteristics in small segments of data obtained by segmentation and then rearrange the data characteristics into a group of data to form characteristic data; after signal characteristics are extracted, connecting the characteristic data of each section in sequence to form a final characteristic data group;
the data features in the signal feature extraction process comprise a maximum value, a minimum value, an average value, a peak-to-peak value, a rectified mean value, a variance, a standard deviation, a kurtosis, a root mean square, a form factor, a peak factor, a kurtosis factor, a pulse factor and a margin factor;
(2) and calculating a corresponding result correlation weight for each feature of the obtained feature data set by using a Relief algorithm, as shown in fig. 8, removing features with lower weights according to a set weight threshold, and then obtaining a final training feature sequence, as shown in fig. 9.
The Relief algorithm is a Feature weighting algorithm (Feature weighting algorithms), different weights are given to features according to the relevance of each Feature and category, and features with weights smaller than a certain threshold value are removed. The relevance of features and classes in the Relief algorithm is based on the discriminative power of features on close-range samples. Randomly selecting a sample R from a training set D by adopting a Relief algorithm, then searching a nearest neighbor sample H from samples similar to the R, wherein the nearest neighbor sample H is called NearHit, searching a nearest neighbor sample M from samples different from the R, wherein the nearest neighbor sample M is called NearMiss, and the samples NearHit and NearMiss are also samples in the training set D; the weight of each feature is then updated according to the following rules: if the distance between R and Near Hit on a feature is smaller than the distance between R and Near Miss, the feature is beneficial to distinguishing the nearest neighbors of the same class and different classes, and the weight of the feature is increased; conversely, if the distance between R and Near Hit is greater than the distance between R and Near miss, indicating that the feature has a negative effect on distinguishing between similar and dissimilar nearest neighbors, the weight of the feature is reduced. Repeating the above processes m times to obtain the average weight of each feature. The larger the weight of a feature is, the stronger the classification capability of the feature is, and conversely, the weaker the classification capability of the feature is. The running time of the Relief algorithm is increased linearly along with the increase of the sampling times m of the samples and the number N of the original features, so that the running efficiency is very high.
Further, step (4) is followed by step (5): and performing post-processing according to the obtained final result of the judgment of the arch-up state of the track slab, wherein the post-processing is to judge whether the track slab is arched up by combining other judgment programs, and the classification result can correct obvious identification errors, further optimize the identification effect and perform early warning work. The other judging program may be the judging program mentioned in patent CN201910620162.8, and judges whether the track slab is arched upward according to the comparison between the set inclination angle threshold and the collected inclination angle data.
Compared with the prior art, the method can realize that:
(1) events are converted into time series classification problems, and the application scenes are more diverse by utilizing new technology classification;
(2) the deep learning technology is adopted, so that data information loss is avoided, and the method can be conveniently expanded to other similar tasks;
(3) according to the difference of the collected data sets, more result types can be judged, so that the model has wider application range and more extensibility;
(4) and the module can be further more accurate and reliable by combining with other result judgment.

Claims (5)

1. A method for detecting the arch state on a track slab based on a deep learning technology is characterized by comprising the following steps:
(1) the traditional vehicle-plate type track dynamics model is improved into a vehicle-plate type track dynamics model under the condition of considering the CA mortar void effect;
(2) constructing a database for deep learning of the classified neural network: outputting track slab displacement simulation data under different voiding degrees by setting different parameters by using the vehicle-slab track dynamics model under the voiding effect of the CA mortar obtained in the step (1), and after training for multiple times to obtain enough simulation data, adding labels to the simulation data according to disease types to form a database, wherein the database comprises a training set for inputting a classification neural network and a test set for inputting the classification neural network;
(3) training a classification network by using the database constructed in the step (2): the classified network structure comprises four layers of networks, wherein the first layer of network is an input layer, the second layer of network is a BilSTM, the third layer of network is a full-connection layer, and the fourth layer of network is a softmax layer;
(4) and (3) carrying out arch state classification on the track slab according to the classification network result: taking the track slab displacement data as the network input quantity of the trained classification network, and sorting the original track slab displacement data into a characteristic sequence to be input into a BilSTM network; and after circulation of the BilSTM layer, inputting the circular shape into a softmax layer to judge the type of the arch state on the track slab, and giving a final result of judgment of the arch state on the track slab.
2. The method for detecting the camber state on the track slab based on the deep learning technology as claimed in claim 1, wherein in the step (4), the original track slab displacement data is sorted into a feature sequence, and the specific method comprises:
(1) preprocessing a track slab displacement signal: expressing track slab displacement signals by using various types of characteristics, normalizing time sequence signals to a (0,1) range, segmenting original track slab displacement signals, and extracting signal characteristics, wherein the signal characteristic extraction is to calculate data characteristics in small segments of data obtained by segmentation and then rearrange the data characteristics into a group of data to form characteristic data; after signal characteristics are extracted, connecting the characteristic data of each section in sequence to form a final characteristic data group;
the data features in the signal feature extraction process comprise a maximum value, a minimum value, an average value, a peak-to-peak value, a rectified mean value, a variance, a standard deviation, a kurtosis, a root mean square, a form factor, a peak factor, a kurtosis factor, a pulse factor and a margin factor;
(2) and calculating corresponding result relevance weight for each feature of the obtained feature data set by utilizing a Relief algorithm, and eliminating the features with lower weight according to a set weight threshold value to obtain a final trained feature sequence.
3. The method for detecting the arch state on the track slab based on the deep learning technology as claimed in claim 1, wherein the step (4) is followed by the step (5): and performing post-processing according to the obtained final result of the judgment of the arch-up state of the track slab, wherein the post-processing is to judge whether the track slab is arched up by combining other judgment programs, and the classification result can correct obvious identification errors, further optimize the identification effect and perform early warning work.
4. The method for detecting the arch state on the track slab based on the deep learning technology as claimed in claim 1, wherein in the step (1), the detailed form of the vibration type coordinate differential equation set of the vehicle-slab track dynamics model under consideration of CA mortar void is as follows:
Figure FDA0002444576820000031
in the formula: t isn(t) introducing a free beam orthogonal function system (X) after a Ritz method is adopted in a vertical vibration differential equation of the track slabnNMS is the mode order of the track plate, and NMS generalized coordinates T are selectedn(t);
Figure FDA0002444576820000032
Is Tn(t) a second derivative; esIsBending rigidity of the track slab; m issβ being the unit length mass of the track slabnIs a constant; m is0Finger m0A plurality of discrete units; csqIs the distributed damping of the CA mortar at q; xpIs a free beam orthogonal function system { X }nMeaning of, p is similar to n, the range is the same as n, XnIs a free beam orthogonal function system, X, of the traditional track slabpThe free beam orthogonal function system of the improved track slab is provided; t isp(t) and XpIn the same way, TnP of (t) is similar to n, the value range is the same as n, the n represents the generalized coordinate of the traditional track slab, and the p represents the generalized coordinate of the improved track slab; ksqThe distributed rigidity of the CA mortar at the position q; f corresponds to the concrete supporting layer, zf(x, t) is the vibration displacement (m), z of the concrete supporting layerf(xqAnd t) is the vibration displacement of the q CA mortar discrete unit concrete supporting layer;
Figure FDA0002444576820000033
is zf(xqT) obtaining a derivation, and representing the vibration speed of the q CA mortar discrete unit concrete supporting layer; l issIs the length of a single track slab; n is0The number of fasteners of a steel rail on one rail plate; subscript s corresponds to the track slab; x is the number ofiThe method is characterized by comprising the following steps of (1) establishing a coordinate system of an original coordinate system on the basis of a steel rail model, wherein N is the number of fasteners in a plate type track length range L; x is the number ofqIs xiTime q ═ i; cpiIs in a coordinate system xiDamping of the under-rail pad layer when in position; NM is the modal order of the steel rail; kpiIs in a coordinate system xiThe stiffness of the under-rail pad in position; y isp(x(j-1)×n0+i)qp(t) expressing the vertical vibration displacement of the steel rail by applying the regular vibration mode function of the simply supported beam by a Ritz method
Figure FDA0002444576820000041
In the formula, k is p and x has a subscript of (j-1) n0+ i wherein n0The number of fasteners of the steel rail on one rail plate is equal to the number of fasteners of the steel rail on the other rail plate.
5. The method according to claim 1, wherein in the step (2), when constructing the database, 104 meters, that is, 16 track slabs are selected as the length of the vehicle-slab track dynamics model considering CA mortar void, the vehicle speed is 300 km/h, the displacement of the track slab is selected as the sample data when the train passes through one track slab, and when the sampling interval is 10-4And s, obtaining different sample data by setting different emptying lengths, setting CA mortar emptying conditions with the longitudinal lengths of 0 meter, 0.325 meter, 0.65 meter, 1.3 meter and 1.95 meter to obtain samples of five emptying conditions, and then obtaining enough samples for each emptying setting condition to obtain 64 groups of samples.
CN202010275388.1A 2020-04-09 2020-04-09 Track slab arch state detection method based on deep learning technology Active CN111651840B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010275388.1A CN111651840B (en) 2020-04-09 2020-04-09 Track slab arch state detection method based on deep learning technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010275388.1A CN111651840B (en) 2020-04-09 2020-04-09 Track slab arch state detection method based on deep learning technology

Publications (2)

Publication Number Publication Date
CN111651840A true CN111651840A (en) 2020-09-11
CN111651840B CN111651840B (en) 2024-04-02

Family

ID=72352449

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010275388.1A Active CN111651840B (en) 2020-04-09 2020-04-09 Track slab arch state detection method based on deep learning technology

Country Status (1)

Country Link
CN (1) CN111651840B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113191051A (en) * 2021-04-30 2021-07-30 西南交通大学 Establishment method and application of nonlinear Hertz spring

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10111957A1 (en) * 2001-03-13 2002-09-19 Baumann Theodor High speed rail track in tubing tube matches tubing width to support spacing appropriate to track scale and support legs splayed out at set angle to bear centrally between tubing ring joins
CN101719183A (en) * 2009-12-10 2010-06-02 北京交通大学 Test simulation system for rail structures of high-speed railways and urban railway system
CN103678823A (en) * 2013-12-25 2014-03-26 西南交通大学 Plate-type ballastless track structure evaluation method based on CA mortar damage
CN105631124A (en) * 2015-12-28 2016-06-01 北京交通大学 Ballastless track damage analysis method performing combined solution based on definite element expansion and fatigue analysis
CN106650067A (en) * 2016-12-12 2017-05-10 贵州理工学院 Vehicle-track-bridge coupled vibrationfrequency domain model with middle and high frequency response output
CN107729649A (en) * 2017-10-16 2018-02-23 西南交通大学 A kind of method for calculating stock rail coupled system dynamics

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10111957A1 (en) * 2001-03-13 2002-09-19 Baumann Theodor High speed rail track in tubing tube matches tubing width to support spacing appropriate to track scale and support legs splayed out at set angle to bear centrally between tubing ring joins
CN101719183A (en) * 2009-12-10 2010-06-02 北京交通大学 Test simulation system for rail structures of high-speed railways and urban railway system
CN103678823A (en) * 2013-12-25 2014-03-26 西南交通大学 Plate-type ballastless track structure evaluation method based on CA mortar damage
CN105631124A (en) * 2015-12-28 2016-06-01 北京交通大学 Ballastless track damage analysis method performing combined solution based on definite element expansion and fatigue analysis
CN106650067A (en) * 2016-12-12 2017-05-10 贵州理工学院 Vehicle-track-bridge coupled vibrationfrequency domain model with middle and high frequency response output
CN107729649A (en) * 2017-10-16 2018-02-23 西南交通大学 A kind of method for calculating stock rail coupled system dynamics

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
向俊;赫丹;曾庆元;: "横向有限条与无砟轨道板段单元的车轨系统竖向振动分析法", 铁道学报, no. 04 *
朱胜阳: "高速铁路无砟轨道结构伤损行为及其对动态性能的影响", 中国优秀博士学位论文全文数据库信息科技辑 *
赵东田等: "板式无碴轨道CA砂浆的配制和性能", 天津大学学报, vol. 41, no. 7, pages 793 - 799 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113191051A (en) * 2021-04-30 2021-07-30 西南交通大学 Establishment method and application of nonlinear Hertz spring

Also Published As

Publication number Publication date
CN111651840B (en) 2024-04-02

Similar Documents

Publication Publication Date Title
CN111144039B (en) Train dynamic weighing system and weighing method based on deep learning
CN110704801B (en) Bridge cluster structure operation safety intelligent monitoring and rapid detection complete method
CN111256924B (en) Intelligent monitoring method for expansion joint of large-span high-speed railway bridge
Au et al. Parameter identification of vehicles moving on continuous bridges
JP5384166B2 (en) Structural deformation detection system
CN104598753B (en) Bridge moving vehicle load recognition method based on Brakhage V method
CN105923014B (en) A kind of track transition Amplitude Estimation method based on evidential reasoning rule
CN101266190A (en) Normal traffic flow stayd-cable bridge mode parametric measurement device and method
Yin et al. Bridge damage identification under the moving vehicle loads based on the method of physics-guided deep neural networks
CN110806324A (en) Wheel polygon abrasion detection method based on rail displacement and data acquisition equipment
CN116881819B (en) Stay cable working state monitoring method based on isolated forest
CN112329520B (en) Truck bearing fault identification method based on generation countermeasure learning
CN116052435B (en) Urban road congestion influence range definition and road influence calculation method
CN114882069A (en) Taxi track abnormity detection method based on LSTM network and attention mechanism
CN111651840A (en) Method for detecting arch state on track slab based on deep learning technology
Nguyen et al. Damage detection of structural based on indirect vibration measurement results combined with Artificial Neural Network
Zhang et al. Structural Damage Identification System Suitable for Old Arch Bridge in Rural Regions: Random Forest Approach.
Shen et al. A boom damage prediction framework of wheeled cranes combining hybrid features of acceleration and Gaussian process regression
CN117390551A (en) Operation monitoring level high-speed railway bridge influence line tracking identification method
CN112781551B (en) Road surface condition monitoring system and monitoring method thereof
CN114707205A (en) Track foundation differential settlement identification method based on GRU neural network
CN115455826A (en) Asphalt pavement track data expansion method based on radial basis function neural network
Cheema et al. Bridge indirect monitoring using Uniform Manifold Approximation and Projection (UMAP)
CN113420387B (en) Migration diagnosis method and system for rolling bearing of compacting machine
CN117709084B (en) Contact force calculation method considering harshness of pantograph slide plate

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