CN101408951A  Method for obtaining equivalent load spectrum and estimating weariness residual longevity of bridge crane based on neural network  Google Patents
Method for obtaining equivalent load spectrum and estimating weariness residual longevity of bridge crane based on neural network Download PDFInfo
 Publication number
 CN101408951A CN101408951A CNA2008101806221A CN200810180622A CN101408951A CN 101408951 A CN101408951 A CN 101408951A CN A2008101806221 A CNA2008101806221 A CN A2008101806221A CN 200810180622 A CN200810180622 A CN 200810180622A CN 101408951 A CN101408951 A CN 101408951A
 Authority
 CN
 China
 Prior art keywords
 crane
 neural network
 fatigue
 load spectrum
 load
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Granted
Links
Images
Abstract
The invention relates to a method for acquiring valent weight load spectrum and estimating fatigue surplus service life of a bridge crane based on the nerve network, and the method is applied to acquire the valent weight load spectrum of a general bridge crane and estimate the fatigue surplus service life of a boxshaped main girder thereof. The invention aims at solving the technical difficult problem in the prior art that a great deal of actual load spectrum actual measurement is difficult to actualize and the random variation tendency of the valent weight load spectrum of the bridge crane is difficult to express by using a determinate mathematic analytic relational expression due to the limitation of field experimental condition. In order to solve the technical difficult problem, the adopted technical proposal includes that the method comprises the following steps: (1) a nerve network model is established; (2) a leaning sample is obtained; (3) an LMBP nerve network is trained; (4) the valent weight load spectrum of the bridge crane is obtained; and at last, a fatigue stress spectrum and the like are calculated by adopting the valent weight load spectrum of the bridge crane obtained by the nerve network, so the fatigue surplus service life of the main girder of the crane can be estimated.
Description
Technical field
The present invention relates to a kind of overhead travelling crane equivalent load spectrum based on neural network obtains and the fatigue surplus life evaluation method, it is applicable to obtaining of generalpurpose overhead crane equivalent load spectrum and the estimation of box girder fatigue surplus life thereof, and it is obtaining and the fatigue surplus life estimation applicable to other type crane equivalent load spectrum also.
Background technology
Hoisting machinery plays indispensable vital role in economic construction, be the special equipment in the major technologies and equipment industry, in case have an accident, economic loss heaviness, social influence are abominable.Especially in recent years, sudden fracture of hoisting machinery structural system and failure event take place in domestic and international many countries and regions in succession, are that fatigue break destroys and cause the one of the main reasons of this class accident.Therefore, national governments and inspection body unprecedentedly pay attention to the research that the hoisting machinery structural system carries out deciding the longevity, lengthening the life.
Solve that crane is decided the longevity, the condition precedent of lengthening the life is to obtain representative typical loadthe time history that can reflect the true operating position of crane metal construction, i.e. loading spectrum.At present the Crane Load spectrum obtain the two kinds of methods that mainly contain: the one, field measurement, promptly obtain the load of crane by direct testtime history (loading spectrum), this not only needs to consume lot of manpower and material resources, and because the repeated load of variation is born by the crane structure system, each results measured is all inequality, this randomness and uncertainty cause measured result being directly applied to theoretical analysis and engineering practice; The 2nd, use traditional regression method, traditional regression method can find the general relationship of lifted load and working cycle number of times basically, but it is to be based upon a certain amount of mathematical analysis to concern relatively accurate statement on the basis that tradition returns (match) mode, and the concrete actual condition randomness of crane is very strong, is difficult to represent with a definite mathematical analysis relational expression.
Summary of the invention
The objective of the invention is to solve that existing method exists limits, is difficult to carry out a large amount of real loads spectrums because of the field experiment condition and survey and be difficult to represent the technical barrier of the randomness variation tendency of crane equivalent load spectrum with a definite mathematical analysis relational expression, and provide a kind of overhead travelling crane equivalent load spectrum of conveniently obtaining equivalent load spectrum, effectively estimating the crane girder fatigue surplus life to obtain and the fatigue surplus life evaluation method based on neural network.
For solving abovementioned technical barrier, the present invention proposes a kind of overhead travelling crane equivalent load spectrum acquisition methods based on neural network, this method is utilized the Nonlinear Mapping function of neural network, select the feedforward neural network of the single output of dual input, adopt the specified lifted load of crane (lifting capacity), lifted load and working cycle time logarithmic data is as learning sample, using the error back propagation learning algorithm of LM (LevenbergMarquardts) BP trains network, obtain optimum weights, realization is obtained the overhead travelling crane equivalent load spectrum, and it comprises the following steps:
(1) sets up neural network model
The specified lifted load and the lifted load of getting overhead travelling crane are input quantity, and the working cycle number of times of different lifted loads in a working hour of corresponding such crane one specified lifted load is output quantity, make up the nerve network system of the single output of dual input;
(2) obtain learning sample
Specified lifted load and lifted load with overhead travelling crane are the input of learning sample, the working cycle number of times of different lifted loads in a working hour of corresponding such crane one specified lifted load is the hope output of learning sample, obtain learning sample, and specified lifted load value is divided into 8 parts at least learning sample carried out classification;
(3) training LMBP neural network
On the basis of the learning sample that neural network model that step (1) obtains and step (2) obtain, training LMBP neural network, and obtain optimum model parameter;
(4) the crane equivalent load spectrum obtains
The different lifted loads with at least 8 grades of specified lifted load of needs being estimated the overhead travelling crane of fatigue surplus life are input to the neural network that step (3) trains, and can obtain the equivalent load spectrum of this overhead travelling crane.
Described neural network structure is three layers, i.e. input layer, hidden layer and output layer, and the input layer number is 2, and output layer node number is 1, and the number of hidden nodes is between 13～17.
Sample data in the described step (2) is that overhead travelling crane is in lifted load and the corresponding work cycle index thereof under the normal operating conditions, makes learning sample meet the actual applying working condition of crane.
A kind of method that adopts the equivalent load spectrum estimation girder fatigue surplus life of the overhead travelling crane that obtains by neural network, it comprises the following steps:
(1) adopts the equivalent load spectrum of the overhead travelling crane that obtains by neural network, calculate the Simulating of Fatigue Stress Spectra of this crane girder span centre dangerouse crosssection;
(2) from Simulating of Fatigue Stress Spectra, extract the cycle index and the global cycle number of times of stress amplitudes at different levels, stress amplitudes at different levels again with rain flow method;
(3) calculate the equivalent stress width of cloth σ of this Simulating of Fatigue Stress Spectra according to Miner fatigue damage stress amplitude equivalence formula
_{Re}
(4) with equivalent stress amplitude σ
_{Re}The computing formula of hardware fatigue surplus life is calculated in substitution: the fatigue surplus life that can estimate crane girder.
Because the present invention has adopted technique scheme, utilize the LMBP neural network that trains, can obtain the equivalent load spectrum of overhead travelling crane and the fatigue surplus life of estimation overhead travelling crane key structure spare quickly and efficiently.Therefore, compare with background technology, the present invention has following substantial characteristics:
(1) based on the neural network that has trained, set up the securing mechanism of the equivalent load spectrum of the type crane, and needn't seek to set up specified lifted load and lifted load and the analytic relationship formula that complexity between corresponding working cycle number of times embarrasses, have realization simply, characteristics effectively efficiently;
(2) utilize the network model that trains, can obtain the equivalent load spectrum of the type overhead travelling crane fast and estimate the fatigue surplus life of its girder by computing machine, thereby save the loaded down with trivial details process and a large amount of input of crane field measurement greatly, it is convenient and can accurately estimate the purpose of crane girder fatigue surplus life to realize obtaining equivalent load spectrum.
Description of drawings
Fig. 1 is a schematic flow sheet of the present invention;
Fig. 2 is a LMBP neural network topology structure synoptic diagram of the present invention;
Fig. 3 is the equivalent stress spectral curve of 75 tons of generalpurpose overhead cranes and the measured stress line chart of setting a song to music.
Embodiment
The present invention is described in further detail below in conjunction with drawings and Examples.
As depicted in figs. 1 and 2, the overhead travelling crane equivalent load spectrum acquisition methods in the present embodiment based on neural network, it comprises the following steps:
(1) sets up neural network model
Must obtain the equivalent load spectrum of crane earlier for realizing the fatigue surplus life estimation of crane key structure spare (girder), the specified lifted load (lifting capacity) and the lifted load of getting overhead travelling crane are input quantity, the working cycle number of times of different lifted loads in a working hour of corresponding such crane one specified lifted load (lifting capacity) is output quantity, makes up the single output nerve network system of dual input.Described neural network structure is three layers, i.e. input layer, hidden layer and output layer, and selecting the input layer number is 2, output layer node number is 1, the number of hidden nodes optional 15.
(2) obtain learning sample
The primary work that solves practical problems with artificial neural network is to collect the learning sample data.For making the neural network after the study have good performance, collected sample data should be to begin to lift by crane a weight with overhead travelling crane to rise, end to can begin to lift by crane next article the time, a complete procedure that comprises overhead travelling crane operation and rest normally be in lifted load and corresponding work cycle index thereof under the normal operating conditions, make learning sample meet the actual applying working condition of crane.
Present embodiment is the input of learning sample with the specified lifted load (lifting capacity) and the lifted load of overhead travelling crane, the working cycle number of times of different lifted loads in a working hour of a corresponding specified lifted load (lifting capacity) is the hope output of learning sample, obtain learning sample, and specified lifted load value is divided into 8 parts learning sample carried out classification (the sample classification is The more the better), be about to 0 load value and be divided into 8 parts at least, by equal score value lifted load and corresponding working cycle number of times thereof are carried out classification then to specified lifted load interval; Learning sample data in the present embodiment realize by data acquisition system (DAS) and data recording: the support level place that at first signal wire of pulling force sensor is installed in lifting mechanism top sheave assembly axis, carry out the lifted load that the record pulling force sensor shows in the process of operate as normal at crane then, and in a working hour working cycle number of times of corresponding different lifted loads.Table 1 be depicted as the specified lifted load of obtaining be 100 tons, 120 tons, 140 tons, 160 tons with the lifted load of 180 tons of generalpurpose overhead cranes and the working cycle of corresponding different lifted loads in a working hour time fractional part sample data bunch:
Table 1 overhead travelling crane lifted load working cycle time fractional part sample data bunch
(3) training LMBP neural network
Have only the neural network after the training could realize its function.The learning sample normalization that step (2) obtains, the LMBP neural network that training is determined by step (1) obtains its model parameter.The neural metwork training process is the weight coefficient w1 by continuous adjustment input layer and hidden layer, and the weight coefficient w2 of hidden layer and output layer makes neural network to the output of nonlearning sample and the very little requirement up to satisfied application of error of desired value.
(4) the crane equivalent load spectrum obtains
The neural network that step (3) trains has promptly been set up the different lifted loads of generalpurpose overhead crane of different specified lifted loads and the randomness variation tendency mapping relations of working cycle number of timesand be the equivalent load spectrum of general class overhead travelling crane.The specified lifted load of needs being estimated the generalpurpose overhead crane of fatigue surplus life is input to the neural network that step (3) trains with 8 grades of (same wellbehaved progression) different lifted loads, can obtain the equivalent load spectrum of this crane.As being that the specified lifted load of the generalpurpose overhead crane of 28.5 meters of 75 tons, span is input to the neural network that step (3) trains with 8 grades of different lifted loads with specified lifted load, can obtain the equivalent load spectrum of this crane, its result is as follows, can see that their working cycle number of times variation tendency is consistent with the input sample.
Specified lifted load is the equivalent load spectrum of 75 tons generalpurpose overhead crane
(5) fatigue surplus life estimation
The equivalent load spectrum of the generalpurpose overhead crane that employing is obtained by neural network, calculate the Simulating of Fatigue Stress Spectra of this crane girder span centre dangerouse crosssection, utilize rain flow method from Simulating of Fatigue Stress Spectra, to extract the cycle index and the global cycle number of times of stress amplitudes at different levels, stress amplitudes at different levels, and then calculate the equivalent stress width of cloth σ of this Simulating of Fatigue Stress Spectra according to Miner fatigue damage stress amplitude equivalence formula by 8 grades (same wellbehaved progression)
_{Re}, with equivalent stress amplitude σ
_{Re}The computing formula of hardware fatigue surplus life is calculated in substitution.Can estimate the fatigue surplus life of this crane girder.
As the specified lifted load that step (4) is obtained is the Simulating of Fatigue Stress Spectra that the equivalent load spectrum of the generalpurpose overhead crane of 28.5 meters of 75 tons, span changes into the girder dangerouse crosssection, and the equivalent stress width of cloth, and substitution calculates the computing formula of crane fatigue surplus life, can calculate the fatigue surplus life of this crane girder.As shown in Figure 3, be the Simulating of Fatigue Stress Spectra of the dangerouse crosssection that changes into by equivalent load spectrum and the correlation curve of measured stress spectrum, can see that they are very approaching.
Claims (4)
1, a kind of overhead travelling crane equivalent load spectrum acquisition methods based on neural network is characterized in that comprising the following steps:
(1) sets up neural network model
The specified lifted load and the lifted load of getting overhead travelling crane are input quantity, and the working cycle number of times of different lifted loads in a working hour of corresponding such crane one specified lifted load is output quantity, make up the nerve network system of the single output of dual input;
(2) obtain learning sample
Specified lifted load and lifted load with overhead travelling crane are the input of learning sample, the working cycle number of times of different lifted loads in a working hour of corresponding such crane one specified lifted load is the hope output of learning sample, obtain learning sample, and specified lifted load value is divided into 8 parts at least learning sample carried out classification;
(3) training LMBP neural network
On the basis of the learning sample that neural network model that step (1) obtains and step (2) obtain, training LMBP neural network, and obtain optimum model parameter;
(4) the crane equivalent load spectrum obtains
The different lifted loads with at least 8 grades of specified lifted load of needs being estimated the overhead travelling crane of fatigue surplus life are input to the neural network that step (3) trains, and can obtain the equivalent load spectrum of this overhead travelling crane.
2, the overhead travelling crane equivalent load spectrum acquisition methods based on neural network according to claim 1, it is characterized in that: described neural network structure is three layers, be input layer, hidden layer and output layer, the input layer number is 2, output layer node number is 1, and the number of hidden nodes is between 13～17.
3, the overhead travelling crane equivalent load spectrum acquisition methods based on neural network according to claim 1, it is characterized in that: the sample data in the described step (2) is that overhead travelling crane is in lifted load and the corresponding work cycle index thereof under the normal operating conditions, makes learning sample meet the actual applying working condition of crane.
4, a kind of bridging crane main beam fatigue surplus life evaluation method based on neural network is characterized in that comprising the following steps:
(1) adopts the equivalent load spectrum of the overhead travelling crane that obtains by neural network, calculate the Simulating of Fatigue Stress Spectra of this crane girder span centre dangerouse crosssection;
(2) from Simulating of Fatigue Stress Spectra, extract the cycle index and the global cycle number of times of stress amplitudes at different levels, stress amplitudes at different levels again with rain flow method;
(3) calculate the equivalent stress width of cloth σ of this Simulating of Fatigue Stress Spectra according to Miner fatigue damage stress amplitude equivalence formula
_{Re}
(4) with equivalent stress amplitude σ
_{Re}The computing formula of hardware fatigue surplus life is calculated in substitution, can estimate the fatigue surplus life of crane girder.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN 200810180622 CN101408951B (en)  20081118  20081118  Method for obtaining equivalent load spectrum and estimating weariness residual longevity of bridge crane based on neural network 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN 200810180622 CN101408951B (en)  20081118  20081118  Method for obtaining equivalent load spectrum and estimating weariness residual longevity of bridge crane based on neural network 
Publications (2)
Publication Number  Publication Date 

CN101408951A true CN101408951A (en)  20090415 
CN101408951B CN101408951B (en)  20130403 
Family
ID=40571957
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN 200810180622 Expired  Fee Related CN101408951B (en)  20081118  20081118  Method for obtaining equivalent load spectrum and estimating weariness residual longevity of bridge crane based on neural network 
Country Status (1)
Country  Link 

CN (1)  CN101408951B (en) 
Cited By (17)
Publication number  Priority date  Publication date  Assignee  Title 

CN101830397A (en) *  20100402  20100915  太原科技大学  System and method for acquiring loading spectrum parameter data of bridge crane 
CN102261988A (en) *  20110425  20111130  吴峰崎  Method for testing residual life of crane 
CN102324040A (en) *  20110909  20120118  太原科技大学  Methods for acquiring equivalent load spectrum and estimating residual fatigue life of bridge crane 
CN102567632A (en) *  20111222  20120711  上海交通大学  Shore bridge structure wind vibration fatigue life forecasting method based on accumulated damage of probability 
CN102767133A (en) *  20110506  20121107  东南大学  Arranging method of sensor for detecting fatigue stress of steel box girder bridge 
CN103389220A (en) *  20130801  20131113  太原科技大学  Test detection method of load shedding equivalent weight of large bridge crane 
CN104809317A (en) *  20150519  20150729  武汉理工大学  Entire load spectrum compiling method for bridge type crane based on clustering analysis 
CN106096262A (en) *  20160608  20161109  南京航空航天大学  A kind of aeroengine loading spectrum Calculation of correlation factor method based on rainflow counting circulation 
CN107356445A (en) *  20170622  20171117  长安大学  Working device of loader fatigue test spectrum method for sorting, loading method and device 
CN108760371A (en) *  20180704  20181106  南京泽楚科技有限公司  A kind of inservice hoisting machinery structure residual life computational methods 
CN109408998A (en) *  20181108  20190301  太原科技大学  Estimating method for fatigue life is carried out based on sample incremental quick obtaining stress spectra 
CN109766637A (en) *  20190112  20190517  太原科技大学  Bridge crane structural optimization based on reliability method based on Krigng agent model 
CN109882998A (en) *  20190214  20190614  中国联合网络通信集团有限公司  A kind of prediction technique and system of airconditioning refrigerant 
CN110032555A (en) *  20190416  20190719  上海建科工程咨询有限公司  A kind of neural network tower crane Risk Forecast Method and system 
CN111131237A (en) *  20191223  20200508  深圳供电局有限公司  Microgrid attack identification method based on BP neural network and gridconnected interface device 
CN111591893A (en) *  20200527  20200828  太原科技大学  Method for measuring hoisting load of automobile crane based on neural network 
CN113138005A (en) *  20210513  20210720  华东理工大学  Vehicle random load pattern recognition system based on feature extraction and neural network 

2008
 20081118 CN CN 200810180622 patent/CN101408951B/en not_active Expired  Fee Related
Cited By (26)
Publication number  Priority date  Publication date  Assignee  Title 

CN101830397A (en) *  20100402  20100915  太原科技大学  System and method for acquiring loading spectrum parameter data of bridge crane 
CN102261988A (en) *  20110425  20111130  吴峰崎  Method for testing residual life of crane 
CN102767133A (en) *  20110506  20121107  东南大学  Arranging method of sensor for detecting fatigue stress of steel box girder bridge 
CN102324040A (en) *  20110909  20120118  太原科技大学  Methods for acquiring equivalent load spectrum and estimating residual fatigue life of bridge crane 
CN102324040B (en) *  20110909  20130605  太原科技大学  Methods for acquiring equivalent load spectrum and estimating residual fatigue life of bridge crane 
CN102567632A (en) *  20111222  20120711  上海交通大学  Shore bridge structure wind vibration fatigue life forecasting method based on accumulated damage of probability 
CN102567632B (en) *  20111222  20141210  上海交通大学  Shore bridge structure wind vibration fatigue life forecasting method based on accumulated damage of probability 
CN103389220B (en) *  20130801  20150916  太原科技大学  Large bridge crane load shedding equivalent method for testing and detecting 
CN103389220A (en) *  20130801  20131113  太原科技大学  Test detection method of load shedding equivalent weight of large bridge crane 
CN104809317B (en) *  20150519  20180123  武汉理工大学  Based on cluster analysis crane span structure class crane complete machine loading spectrum preparation method 
CN104809317A (en) *  20150519  20150729  武汉理工大学  Entire load spectrum compiling method for bridge type crane based on clustering analysis 
CN106096262B (en) *  20160608  20181123  南京航空航天大学  Aeroengine load spectral correlative coefficient calculation method based on rainflow counting circulation 
CN106096262A (en) *  20160608  20161109  南京航空航天大学  A kind of aeroengine loading spectrum Calculation of correlation factor method based on rainflow counting circulation 
CN107356445A (en) *  20170622  20171117  长安大学  Working device of loader fatigue test spectrum method for sorting, loading method and device 
CN107356445B (en) *  20170622  20190507  长安大学  Method for sorting, loading method and device are composed in working device of loader fatigue test 
CN108760371A (en) *  20180704  20181106  南京泽楚科技有限公司  A kind of inservice hoisting machinery structure residual life computational methods 
CN109408998A (en) *  20181108  20190301  太原科技大学  Estimating method for fatigue life is carried out based on sample incremental quick obtaining stress spectra 
CN109408998B (en) *  20181108  20230103  太原科技大学  Fatigue life evaluation method by rapidly acquiring stress spectrum based on sample continuous increment 
CN109766637B (en) *  20190112  20220610  太原科技大学  Bridge crane structure reliability optimization method based on Krigng agent model 
CN109766637A (en) *  20190112  20190517  太原科技大学  Bridge crane structural optimization based on reliability method based on Krigng agent model 
CN109882998A (en) *  20190214  20190614  中国联合网络通信集团有限公司  A kind of prediction technique and system of airconditioning refrigerant 
CN109882998B (en) *  20190214  20210827  中国联合网络通信集团有限公司  Prediction method and system for air conditioner refrigerant 
CN110032555A (en) *  20190416  20190719  上海建科工程咨询有限公司  A kind of neural network tower crane Risk Forecast Method and system 
CN111131237A (en) *  20191223  20200508  深圳供电局有限公司  Microgrid attack identification method based on BP neural network and gridconnected interface device 
CN111591893A (en) *  20200527  20200828  太原科技大学  Method for measuring hoisting load of automobile crane based on neural network 
CN113138005A (en) *  20210513  20210720  华东理工大学  Vehicle random load pattern recognition system based on feature extraction and neural network 
Also Published As
Publication number  Publication date 

CN101408951B (en)  20130403 
Similar Documents
Publication  Publication Date  Title 

CN101408951B (en)  Method for obtaining equivalent load spectrum and estimating weariness residual longevity of bridge crane based on neural network  
CN105043631B (en)  Vibratory drilling method Cable force measuring method based on linear model  
CN102324040B (en)  Methods for acquiring equivalent load spectrum and estimating residual fatigue life of bridge crane  
CN105241660B (en)  Highspeed rail large bridge performance test methods based on health monitoring data  
CN107016235A (en)  The equipment running status health degree appraisal procedure adaptively merged based on multiple features  
Deng et al.  Fatigue reliability assessment for bridge welded details using longterm monitoring data  
CN108226399A (en)  A kind of beamstring structure damage combined recognising method based on modal parameter  
CN103489035A (en)  Power grid electric energy quality comprehensive evaluation method based on gray weighting correlation analysis algorithm  
CN106932135B (en)  Flexible inhaul cable force testing method for identifying vibration frequency based on weighted narrowband peak searching method  
CN104063612A (en)  Tunnel engineering risk situation fuzzy evaluation method and system  
CN105300692A (en)  Bearing fault diagnosis and prediction method based on extended Kalman filtering algorithm  
CN106570259A (en)  Gross error elimination method for dam displacement data  
CN102486833A (en)  Method of efficacy anticipation and failure examination for apparatus  
CN101667214B (en)  Method for arranging displacement measuring points  
CN104933254A (en)  Vibrationmethod cable force measurement method of improved linear model on the basis of frequency ratio offset coefficient  
CN103838931A (en)  Method for evaluating remanufacturing access period of engineering mechanical arm rest class structure  
CN108319767A (en)  The method of suspension bridge sunpender stress assessment based on traveling load  
CN108549847B (en)  Beam type structure crack damage identification method under nonreference data condition  
CN107421672B (en)  Weighted search force calculation method based on global peak searching of vibration frequency  
CN106446384B (en)  A kind of damnification recognition method of bridging crane main beam structure  
CN109978374B (en)  Risk assessment method for oil and gas pipeline system  
CN110796293A (en)  Power load prediction method  
CN105488734A (en)  Full life circle based quality assessment method and system for electric energy metering device  
CN108573119A (en)  Inservice hoisting machinery overall structure remaining life realtime evaluation system  
CN104807661A (en)  Highlevel and highrise structure dynamic detection bearing capacity evaluating method 
Legal Events
Date  Code  Title  Description 

C06  Publication  
PB01  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
C14  Grant of patent or utility model  
GR01  Patent grant  
CF01  Termination of patent right due to nonpayment of annual fee 
Granted publication date: 20130403 Termination date: 20141118 

EXPY  Termination of patent right or utility model 