CN109341989A - A kind of Bridge Influence Line recognition methods that can reject vehicle power effect - Google Patents
A kind of Bridge Influence Line recognition methods that can reject vehicle power effect Download PDFInfo
- Publication number
- CN109341989A CN109341989A CN201811016637.4A CN201811016637A CN109341989A CN 109341989 A CN109341989 A CN 109341989A CN 201811016637 A CN201811016637 A CN 201811016637A CN 109341989 A CN109341989 A CN 109341989A
- Authority
- CN
- China
- Prior art keywords
- bridge
- imf
- response
- res
- quasi
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
- G01M5/0008—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings of bridges
Abstract
The invention belongs to safety of structure detection technique fields, disclose a kind of Bridge Influence Line recognition methods that can reject vehicle power effect, and method includes: that (1) pre-processes the bridge response data that sensor measures using empirical mode decomposition;(2) it is identified to by the response of pretreated bridge using Quasi-static Method.The power superposition item that bridge can be responded by being pre-processed using empirical mode decomposition to bridge response is rejected, to obtain the quasi_static response of bridge, is further useful for influencing line identification.This method can handle the structural response data measured on various bridge types by sensor, have good future in engineering applications without understanding any bridge structure information in advance simultaneously.
Description
Technical field
The invention belongs to safety of structure detection technique fields, and in particular to one kind can reject the bridge of vehicle power effect
Beam influences line recognition methods.
Technical background
In recent years, the promotion with the continuous development and people of high-speed rail technology to trip efficiency requirements, in the world
Built high-speed railway total mileage is unprecedented to be improved.Since most operating mileages of high-speed railway are all in height
On bridge formation, the safe condition of high-speed rail bridge is also had been to be concerned by more and more people.High-speed rail bridge is arranged due to constantly bearing high speed
The power of vehicle loads, and the load effect being subject to is more complicated than common highway bridge, how to monitor high metal bridge in real time
The military service performance of beam becomes a significant problem concerning passenger survival property safety.
It is also worth noting that occurring in recent years a kind of based on the damnification recognition method for influencing line variation.Bridge
The response curve that line is bridge a certain specific position when unit load passes through bridge is influenced, as a kind of statical feature of bridge,
It influences line and contains a large amount of structural information.As the index of a structural damage, the rigid of bridge can directly be reacted by influencing line
Degree and flexibility information.In addition to this, the location information that line further comprises bridge is influenced, it can be by influencing the variation of line to damage
It is positioned.Since load running track is single, being particularly suitable for influences line to it identifies railroad bridge.But bullet train mistake
It is biggish that generation is compared in quasi-static bridge response when dynamic effect caused by when bridge can make bridge response pass a bridge with low speed train
Deviation, how from dynamic effect identification obtain influence line become high-speed railway bridge influence line identification a difficult point ask
Topic.
For the semi-static load situation that low speed load train is passed a bridge, influencing line knowledge method for distinguishing has had many researchs
Achievement.OBerin proposes a kind of matrix method at first, which carrys out the identification from measured data, influences line.This method establish it is a kind of and
The relevant matrix of train information, by identifying Bridge Influence Line to matrix inversion.Sio-Song Leng is estimated using maximum likelihood
Meter method is deduced a kind of influence line identification formula, and this method can consider the effect of multiple load vehicles simultaneously, thus
Give a kind of more accurate influence line estimation.Chen uses the influence line that regularization method obtains identification first and is not inconsistent
The fluctuation for closing physical significance is eliminated.Recognize that bridge response can be expressed as influencing line and train information letter
Several convolution influences line by the method identification of fast Fourier transform deconvolution.This method can greatly improve influence line
The efficiency of identification, and also have certain guarantee on accuracy of identification.In the case where considering dynamic effect, Wang Ningbo et al. is logical
It crosses multinomial and SIN function to be respectively fitted static(al) effect and dynamic effect, identification has obtained certain specific bridge types
Influence line.And for other more common situations, up for further studying.
Summary of the invention
The purpose of the present invention is one kind can reject the Bridge Influence Line recognition methods of vehicle power effect.
Technical solution of the present invention:
A kind of Bridge Influence Line recognition methods that can reject vehicle power effect, steps are as follows:
Step 1 pre-processes the bridge response message that sensor measures using empirical mode decomposition
(1) local maximum and local minimum of bridge response signal R (t) are found out;
(2) minimum interpolation and maximum interpolating function are calculated separately using cubic spline interpolation method, and is corresponded to
Signal envelope emin(t) and emax(t);
(3) local mean value M (t)=(e of signal maximum and minimum envelope is calculatedmin(t)+emax(t))/2。
(4) original input signal is subtracted into local mean value and obtains oscillator signal H (t)=R (t)-M (t);
(5) when H (t) meets the condition of intrinsic mode functions, H (t) becomes an imf (t);Otherwise, with H (t) replacement step
(1) R (t) in is simultaneously repeated the above process since step (1);One intrinsic mode functions must satisfy 2 conditions: 1) function exists
In entire time range, the number of Local Extremum and zero crossing is equal, or at most poor one;2) at any point in time, office
The envelope of portion's maximum value and the envelope average value of local minimum are necessary for 0;
(6) imf is enabled1(t)=H (t), then imf1It (t) is first imf (t), corresponding surplus Res1(t)=R (t)-
imf1(t), by Res1(t) above-mentioned all steps are repeated as new signal, obtains second imf (t) component, and so on, it obtains
To Res1(t)-imf2(t)=Res2(t) ... Resn-1(t)-imfn(t)=Resn(t);When intrinsic mode functions can not be further continued for mentioning
When taking, stop screening process;
(7) pass through decomposable process above, R (t) is finally broken down into n imf (t) component and a surplus Resn(t);
Original signal R (t) is expressed asSurplus Res after treatmentnIt (t) is exactly that bridge standard is quiet
State response;
Step 2 carries out influence line identification using Quasi-static Method to by the response of pretreated bridge
After the pretreatment of empirical mode decomposition, bridge dynamic response is converted to the response of the bridge under quasi-static situation,
It is identified using classical quasi-static influence line identification model.
Beneficial effects of the present invention:
(1) Bridge Influence Line recognition methods of the invention is compared with original influence line recognition methods, can be arranged high speed
Bridge structure response when vehicle loads is identified, is influenced line identification for high-speed rail bridge and is provided theoretical basis;
(2) Bridge Influence Line recognition methods of the invention has stringent theoretical basis, is rung based on bridge load and bridge
The synchronous acquisition information answered, and combine advanced Statistical error algorithm, it is ensured that have by the influence line of the system identification higher
Precision;
(3) Bridge Influence Line recognition methods of the invention uses simple, can be to acquisition without predicting bridge structure information
Obtained bridge dynamic response data is directly handled, to obtain influencing line.
Detailed description of the invention
Fig. 1 for the used algorithm of the present invention implementation flow chart;
Fig. 2 is the model of vibration that the load vehicle simulated in embodiment of the present invention method is passed a bridge;
Fig. 3 is that bridge dynamic response caused by vehicle is loaded in embodiment of the present invention method;
Fig. 4 is the quasi_static response obtained after being handled in embodiment of the present invention method by empirical mode decomposition;
Fig. 5 is the influence line identified in embodiment of the present invention method to quasi_static response after processing;
Specific embodiment
The present invention is described in further detail with a numerical example with reference to the accompanying drawing.
Influence line recognition methods of the invention point " the bridge response message that sensor is measured using empirical mode decomposition into
Row pretreatment " and " being identified to by the response of pretreated bridge using Quasi-static Method " two steps, specific embodiment
It has been given above, the application method and feature of invention is next illustrated in conjunction with an example.
Implement to calculate: 72km/h Biaxial stress vehicle, which is passed a bridge, influences line identification
In this numerical example, we simulate the for-wheel vehicle that a speed is 72km/h and pass through simply supported beam, pass through foundation
Vehicle bridge coupling vibration model responds to analyze the amount of deflection of span centre.The wheelbase of for-wheel vehicle is 4m, and freely-supported beam length is 16m.By vehicle
Body is simplified to three mass blocks, and quality is respectively: m1=524kg;m2=297kg;m3=6451kg.The concrete condition of model is shown in
Fig. 2.
In this example, the differential equation of motion of four-degree-of-freedom vehicle and distributed mass system simply supported beam is established, is passed through
The time-histories that ordinary differential system obtains the response of simply supported beam span centre is solved, as shown in Figure 3.Deflection of bridge span when bullet train is passed a bridge
Response is segmented into quasi_static response and power superposition item two parts, can be divided two parts by empirical mode decomposition
Solution, decomposition result are shown in Fig. 4.
After being rejected to dynamic effect, it can be identified with quasi-static influence line recognition methods to line is influenced, in this calculation
In example, we are using regularization LSQR iterative method.Identify that obtained influence line and the comparison of true impact line are shown in Fig. 5.It can be with
Find out, the influence line identified and true value coincide substantially.This method is one kind of Bridge Influence Line when identifying dynamically load
Effective ways.
Claims (1)
1. the Bridge Influence Line recognition methods that one kind can reject vehicle power effect, which is characterized in that steps are as follows:
Step 1 pre-processes the bridge response message that sensor measures using empirical mode decomposition
(1) local maximum and local minimum of bridge response signal R (t) are found out;
(2) minimum interpolation and maximum interpolating function are calculated separately using cubic spline interpolation method, and obtained to induction signal
Envelope emin(t) and emax(t);
(3) local mean value M (t)=(e of signal maximum and minimum envelope is calculatedmin(t)+emax(t))/2;
(4) original input signal is subtracted into local mean value and obtains oscillator signal H (t)=R (t)-M (t);
(5) when H (t) meets the condition of intrinsic mode functions, H (t) becomes an imf (t);Otherwise, in H (t) replacement step (1)
R (t) and repeated the above process since step (1);One intrinsic mode functions must satisfy 2 conditions: 1) function is entire
In time range, the number of Local Extremum and zero crossing is equal, or at most poor one;2) at any point in time, part is most
The envelope being worth greatly and the envelope average value of local minimum are necessary for 0;
(6) imf is enabled1(t)=H (t), then imf1It (t) is first imf (t), corresponding surplus Res1(t)=R (t)-imf1(t),
By Res1(t) above-mentioned all steps are repeated as new signal, obtains second imf (t) component, and so on, obtain Res1
(t)-imf2(t)=Res2(t) ... Resn-1(t)-imfn(t)=Resn(t);When intrinsic mode functions can not be further continued for extracting,
Stop screening process;
(7) pass through decomposable process above, R (t) is finally broken down into n imf (t) component and a surplus Resn(t);It is original
Signal R (t) is expressed asSurplus Res after treatmentnIt (t) is exactly the quasi-static sound of bridge
It answers;
Step 2 carries out influence line identification using Quasi-static Method to by the response of pretreated bridge
After the pretreatment of empirical mode decomposition, bridge dynamic response is converted to the response of the bridge under quasi-static situation, uses
Classical quasi-static influence line identification model is identified.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811016637.4A CN109341989B (en) | 2018-09-03 | 2018-09-03 | Bridge influence line identification method capable of eliminating vehicle power effect |
PCT/CN2019/073826 WO2020048090A1 (en) | 2018-09-03 | 2019-01-30 | Bridge influence line identification method capable of eliminating vehicle power effect |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811016637.4A CN109341989B (en) | 2018-09-03 | 2018-09-03 | Bridge influence line identification method capable of eliminating vehicle power effect |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109341989A true CN109341989A (en) | 2019-02-15 |
CN109341989B CN109341989B (en) | 2020-06-12 |
Family
ID=65292177
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811016637.4A Active CN109341989B (en) | 2018-09-03 | 2018-09-03 | Bridge influence line identification method capable of eliminating vehicle power effect |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN109341989B (en) |
WO (1) | WO2020048090A1 (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110261305A (en) * | 2019-06-18 | 2019-09-20 | 南京东南建筑机电抗震研究院有限公司 | Based on across footpaths continuous bridge damnification recognition methods such as the multispan for influencing line |
CN112485030A (en) * | 2020-11-09 | 2021-03-12 | 深圳市桥博设计研究院有限公司 | Bridge structure dynamic monitoring method, system and equipment based on frequency coupling |
CN114414122A (en) * | 2022-01-07 | 2022-04-29 | 大连理工大学 | Sensing device for measuring multi-component force of wheel |
Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11326108A (en) * | 1998-05-07 | 1999-11-26 | Shimadzu Corp | Helium leak detector |
JP2003233640A (en) * | 2001-12-03 | 2003-08-22 | Tateishi Kozo Sekkei:Kk | Structural calculation device for building, computer program, record medium and building |
CN101532917A (en) * | 2009-04-15 | 2009-09-16 | 江苏省交通科学研究院股份有限公司 | Quick load test method for bridge carrying capacity |
CN102509005A (en) * | 2011-10-27 | 2012-06-20 | 招商局重庆交通科研设计院有限公司 | Bridge bearing capacity evaluation method based on field tested influence line |
JP2014010572A (en) * | 2012-06-28 | 2014-01-20 | Strauss Technology Co Ltd | Vector image drawing device and vector image drawing method |
CN104389275A (en) * | 2014-12-10 | 2015-03-04 | 邢兵 | Comprehensive structure prestress reinforcing method on basis of influence line theory |
CN104819813A (en) * | 2015-04-29 | 2015-08-05 | 中南大学 | Bridge influence line dynamic test method |
CN104897180A (en) * | 2015-05-26 | 2015-09-09 | 广州大学 | Method for pretreating bridge monitoring signals |
CN104964837A (en) * | 2015-06-12 | 2015-10-07 | 广东电网有限责任公司电力科学研究院 | Method and system for monitoring rigidity damage of structure based on EMD |
KR20160044528A (en) * | 2013-08-22 | 2016-04-25 | 엠파이어 테크놀로지 디벨롭먼트 엘엘씨 | Influence of line of sight for driver safety |
CN105628381A (en) * | 2016-01-13 | 2016-06-01 | 东北石油大学 | Reciprocating compressor bearing fault diagnosis method based on improved local mean value decomposition |
CN105973619A (en) * | 2016-04-27 | 2016-09-28 | 厦门大学 | Bridge local damage identification method based on influence line under structure health monitoring system |
CN105973627A (en) * | 2016-05-26 | 2016-09-28 | 东南大学 | Long-gauge-length-strain-influence-envelope-based bridge damage identification method |
CN106482917A (en) * | 2016-10-11 | 2017-03-08 | 广州大学 | A kind of detection method of cable-stayed bridge main-beam dynamic deflection |
CN107132011A (en) * | 2017-05-31 | 2017-09-05 | 中南大学 | A kind of bridge quick determination method based on influence line |
CN107300452A (en) * | 2017-08-10 | 2017-10-27 | 天津市德力电子仪器有限公司 | A kind of Test on Bridge Loading rapid detection system |
WO2017200380A1 (en) * | 2016-05-18 | 2017-11-23 | Heijmans N.V. | Method for determining the structural integrity of an infrastructural element |
CN107389168A (en) * | 2017-07-18 | 2017-11-24 | 重庆交通大学 | A kind of vehicle for bridge moves the recognition methods of load |
CN107588915A (en) * | 2017-10-18 | 2018-01-16 | 厦门大学 | A kind of Bridge Influence Line recognition methods and system |
CN207395997U (en) * | 2017-10-18 | 2018-05-22 | 厦门大学 | A kind of Bridge Influence Line identifying system |
-
2018
- 2018-09-03 CN CN201811016637.4A patent/CN109341989B/en active Active
-
2019
- 2019-01-30 WO PCT/CN2019/073826 patent/WO2020048090A1/en active Application Filing
Patent Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11326108A (en) * | 1998-05-07 | 1999-11-26 | Shimadzu Corp | Helium leak detector |
JP2003233640A (en) * | 2001-12-03 | 2003-08-22 | Tateishi Kozo Sekkei:Kk | Structural calculation device for building, computer program, record medium and building |
CN101532917A (en) * | 2009-04-15 | 2009-09-16 | 江苏省交通科学研究院股份有限公司 | Quick load test method for bridge carrying capacity |
CN102509005A (en) * | 2011-10-27 | 2012-06-20 | 招商局重庆交通科研设计院有限公司 | Bridge bearing capacity evaluation method based on field tested influence line |
JP2014010572A (en) * | 2012-06-28 | 2014-01-20 | Strauss Technology Co Ltd | Vector image drawing device and vector image drawing method |
KR20160044528A (en) * | 2013-08-22 | 2016-04-25 | 엠파이어 테크놀로지 디벨롭먼트 엘엘씨 | Influence of line of sight for driver safety |
CN104389275A (en) * | 2014-12-10 | 2015-03-04 | 邢兵 | Comprehensive structure prestress reinforcing method on basis of influence line theory |
CN104819813A (en) * | 2015-04-29 | 2015-08-05 | 中南大学 | Bridge influence line dynamic test method |
CN104897180A (en) * | 2015-05-26 | 2015-09-09 | 广州大学 | Method for pretreating bridge monitoring signals |
CN104964837A (en) * | 2015-06-12 | 2015-10-07 | 广东电网有限责任公司电力科学研究院 | Method and system for monitoring rigidity damage of structure based on EMD |
CN105628381A (en) * | 2016-01-13 | 2016-06-01 | 东北石油大学 | Reciprocating compressor bearing fault diagnosis method based on improved local mean value decomposition |
CN105973619A (en) * | 2016-04-27 | 2016-09-28 | 厦门大学 | Bridge local damage identification method based on influence line under structure health monitoring system |
WO2017200380A1 (en) * | 2016-05-18 | 2017-11-23 | Heijmans N.V. | Method for determining the structural integrity of an infrastructural element |
CN105973627A (en) * | 2016-05-26 | 2016-09-28 | 东南大学 | Long-gauge-length-strain-influence-envelope-based bridge damage identification method |
CN106482917A (en) * | 2016-10-11 | 2017-03-08 | 广州大学 | A kind of detection method of cable-stayed bridge main-beam dynamic deflection |
CN107132011A (en) * | 2017-05-31 | 2017-09-05 | 中南大学 | A kind of bridge quick determination method based on influence line |
CN107389168A (en) * | 2017-07-18 | 2017-11-24 | 重庆交通大学 | A kind of vehicle for bridge moves the recognition methods of load |
CN107300452A (en) * | 2017-08-10 | 2017-10-27 | 天津市德力电子仪器有限公司 | A kind of Test on Bridge Loading rapid detection system |
CN107588915A (en) * | 2017-10-18 | 2018-01-16 | 厦门大学 | A kind of Bridge Influence Line recognition methods and system |
CN207395997U (en) * | 2017-10-18 | 2018-05-22 | 厦门大学 | A kind of Bridge Influence Line identifying system |
Non-Patent Citations (4)
Title |
---|
CHENG, CHIA-CHI等: "Evaluation of Bridge Span by Recovered Stiffness Data Obtained with Moving Vehicle Loadings", 《HEALTH MONITORING OF STRUCTURAL AND BIOLOGICAL SYSTEMS 2012》 * |
WU BITAO等: "Damage identification and bearing capacity evaluation of bridges based on distributed long-gauge strain envelope line under moving vehicle loads", 《JOURNAL OF INTELLIGENT MATERIAL SYSTEMS AND STRUCTURES》 * |
吴必涛: "基于分布式长标距应变传感的桥梁动态性能监测和损伤识别方法研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
王宁波等: "基于影响线的桥梁移动荷载识别", 《振动与冲击》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110261305A (en) * | 2019-06-18 | 2019-09-20 | 南京东南建筑机电抗震研究院有限公司 | Based on across footpaths continuous bridge damnification recognition methods such as the multispan for influencing line |
CN112485030A (en) * | 2020-11-09 | 2021-03-12 | 深圳市桥博设计研究院有限公司 | Bridge structure dynamic monitoring method, system and equipment based on frequency coupling |
CN112485030B (en) * | 2020-11-09 | 2023-03-14 | 深圳市桥博设计研究院有限公司 | Bridge structure dynamic monitoring method, system and equipment based on frequency coupling |
CN114414122A (en) * | 2022-01-07 | 2022-04-29 | 大连理工大学 | Sensing device for measuring multi-component force of wheel |
CN114414122B (en) * | 2022-01-07 | 2022-09-27 | 大连理工大学 | Sensing device for measuring multi-component force of wheel |
Also Published As
Publication number | Publication date |
---|---|
CN109341989B (en) | 2020-06-12 |
WO2020048090A1 (en) | 2020-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shao et al. | Modified stacked autoencoder using adaptive Morlet wavelet for intelligent fault diagnosis of rotating machinery | |
CN113221277B (en) | Bearing performance degradation evaluation method and system based on digital twin model | |
Huang et al. | A multisource dense adaptation adversarial network for fault diagnosis of machinery | |
Li et al. | Joint attention feature transfer network for gearbox fault diagnosis with imbalanced data | |
CN109341989A (en) | A kind of Bridge Influence Line recognition methods that can reject vehicle power effect | |
CN101916241B (en) | Method for identifying time-varying structure modal frequency based on time frequency distribution map | |
CN109115879B (en) | Structural damage identification method based on modal shape and convolutional neural network | |
Liang et al. | Unsupervised fault diagnosis of wind turbine bearing via a deep residual deformable convolution network based on subdomain adaptation under time-varying speeds | |
Di et al. | Ensemble deep transfer learning driven by multisensor signals for the fault diagnosis of bevel-gear cross-operation conditions | |
Chen et al. | Application of nonlinear output frequency response functions and deep learning to RV reducer fault diagnosis | |
Fan et al. | Fan fault diagnosis based on lightweight multiscale multiattention feature fusion network | |
Zhang et al. | A multi-module generative adversarial network augmented with adaptive decoupling strategy for intelligent fault diagnosis of machines with small sample | |
Yao et al. | Implementation of a novel algorithm of wheelset and axle box concurrent fault identification based on an efficient neural network with the attention mechanism | |
CN114218778A (en) | Method and device for analyzing sonic boom test data | |
Qin et al. | Cross-domain fault diagnosis based on improved multi-scale fuzzy measure entropy and enhanced joint distribution adaptation | |
Bi et al. | A variety of engine faults detection based on optimized variational mode decomposition-robust independent component analysis and fuzzy C-mean clustering | |
CN115659782A (en) | Bridge damage identification method based on PGDNN under action of mobile vehicle load | |
Lv et al. | Gear fault feature extraction based on fuzzy function and improved Hu invariant moments | |
Jin et al. | Structural damage recognition based on filtered feature selection and convolutional neural network | |
Zhao et al. | Research on main bearing life prediction of direct-drive wind turbine based on digital twin technology | |
Tan et al. | Sensible multiscale symbol dynamic entropy for fault diagnosis of bearing | |
Shuai et al. | Adaptive morphological feature extraction and support vector regressive classification for bearing fault diagnosis | |
Liu et al. | A novel wind turbine health condition monitoring method based on common features distribution adaptation | |
Dalian et al. | Dual-rotor misalignment fault quantitative identification based on DBN and improved DS evidence theory | |
Qi et al. | A new deep fusion network for automatic mechanical fault feature learning |
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 |