CN106895900B - Method based on the bridge dynamic strain identification bicycle number of axle - Google Patents
Method based on the bridge dynamic strain identification bicycle number of axle Download PDFInfo
- Publication number
- CN106895900B CN106895900B CN201710089866.8A CN201710089866A CN106895900B CN 106895900 B CN106895900 B CN 106895900B CN 201710089866 A CN201710089866 A CN 201710089866A CN 106895900 B CN106895900 B CN 106895900B
- Authority
- CN
- China
- Prior art keywords
- bridge
- strain
- axle
- dynamic strain
- vehicle
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000012360 testing method Methods 0.000 claims abstract description 22
- 238000009795 derivation Methods 0.000 claims abstract description 11
- 238000001514 detection method Methods 0.000 claims abstract description 7
- 230000000694 effects Effects 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims description 20
- 238000005452 bending Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 4
- 230000007935 neutral effect Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000000638 stimulation Effects 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 230000036541 health Effects 0.000 abstract description 4
- 230000009467 reduction Effects 0.000 abstract description 3
- 238000009434 installation Methods 0.000 description 8
- 238000005259 measurement Methods 0.000 description 7
- 230000004044 response Effects 0.000 description 5
- 238000001914 filtration Methods 0.000 description 4
- 238000012423 maintenance Methods 0.000 description 3
- 230000035882 stress Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000002955 isolation Methods 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000005303 weighing Methods 0.000 description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01G—WEIGHING
- G01G19/00—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups
- G01G19/02—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups for weighing wheeled or rolling bodies, e.g. vehicles
- G01G19/03—Weighing apparatus or methods adapted for special purposes not provided for in the preceding groups for weighing wheeled or rolling bodies, e.g. vehicles for weighing during motion
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
- Vehicle Body Suspensions (AREA)
Abstract
The invention discloses the methods based on the bridge dynamic strain identification bicycle number of axle, belong to bridge testing detection, health monitoring and overload vehicle and administer field.In beam bridge span centre or maximum strain reaction cross-section arrangement high-resolution strain detection testing device (sensor), test the bridge moving strain time history curve under vehicle effect, noise reduction or fitting are filtered to dynamic strain time-history curves, and carry out single order and second order derivation, in conjunction with finite element analog result, curve derivation as a result, calculating the axle number for obtaining vehicle.This method is only needed to arrange 1 strain testing section on bridge, be installed simple, low in cost；It is not damaged to existing highway bridge road without destroying bridge floor or pavement structure；Recognizer is simple, accuracy is high.
Description
Technical field
The present invention relates to bridge testing detection, health monitoring and overload vehicles to administer field, is based particularly on bridge dynamic strain
The method for identifying the bicycle number of axle.
Background technique
In recent years, there are multiple great Bridge Accidents successively in China.The accident that these occur is related with several factors, but
It is a lack of effective monitoring measure and necessary maintenance, maintenance measure is important one of reason.These startling accidents
So that people are to the quality of modern bridge and service life, also gradually concern is got up.Testing inspection, health monitoring are carried out to bridge structure
With communications and transportation overload control vehicle, it has also become the hot spot that domestic and international academia, engineering circles are studied.Traditional bridge machinery is very
The experience of manager and technical staff are depended in big degree, often the situation of bridge especially large bridge are lacked comprehensive
It holds and understands, information cannot timely feedback.If the disease to bridge is underestimated, just it is likely to lose the best of maintenance
Opportunity accelerates the process of bridge damage, shortens the service life of bridge.If will be caused to the disease overestimate of bridge
Unnecessary fund waste, prevent the bearing capacity of bridge is from giving full play to.
Currently, vehicle overload occurs again and again the phenomenon that traveling, and the load of bridge is that have certain limit, when these are overweight
Vehicle by bridge, certain damage can be undoubtedly generated to bridge；And over time, the aging of bridge, holds
Loading capability is also changing, therefore monitors the weight for passing through bridge vehicle, goes forward side by side line number according to statistics, for understanding the healthy shape of bridge
Condition is significant.But generally conventional weighbridge weight measuring equipment is expensive, bicycle road it is generally the least expensive be also required to 100,000 yuan or more, and
It needs to destroy former pavement structure when installation；Routine weighing weighbridge needs artificial or differentiates vehicle and vehicle by picture pick-up device when use
The number of axle；And conventional weighbridge is mounted on the road surface other than bridge main body structure, multilane bridge is difficult to pair using conventional weighbridge
Influence of the driving alongside vehicle to bridge counts on bridge.I.e. using conventional weighing weighbridge and the not applicable vehicle on bridge
Axle number, axle away from or axle weight measurement, so, be badly in need of in society a kind of total suitable for measuring the vehicle travelled in bridge
Weight, axle number, axle away from or axle weight cost effective method or equipment.
The open method based on bridge dynamic strain identification vehicle weight of Chinese patent 201210249735.9, including following step
It is rapid: dynamic strain measure device is arranged on the longitudinally asymmetric section of bridge；In the bridge runway for installing dynamic strain measure device
Upper progress sport car experiment, while data are acquired, and send data to computer；It is answered using finite element stimulation theory is dynamic
Become, extract the dynamic strain time-history curves that dynamic strain result is formed, and be compared with measured value and correlation analysis, obtains driving speed
The relational expression T=f (x, ε) or curve of degree, dynamic strain peak value and vehicle weight；According to each section dynamic strain reach peak value when
Between difference and cross-sectional distance running speed calculated using software automatically, be averaged with each measuring point dynamic strain peak value in section, bring relationship into
Formula T=f (x, ε) or curve are calculated or are compared the car weight acquired through the section.Here driving speed is found by test
The relational expression of degree, dynamic strain peak value and vehicle weight, has that heavy workload, adaptability be not high, especially each bridge
Architectural characteristic be all different, cause to install every time in this way require bridge in use is needed to carry out closure test,
And test number (TN) is more, installation effectiveness is low, it is at high cost to install and use, rate of its popularization is not high, can not cope with Chinese bridge at present
Measurement demand.
Chinese patent 201610114464.4 discloses a kind of axle recognition methods and system for bridge, this method packet
It includes step: altogether including two section groups in 4~6 sections along bridge longitudinally label, form two virtual simply supported beams；Acquisition
The bridge time-histories data at two section groups is measured, two groups of isolation responses are calculated according to time-histories data, two groups of isolation are rung
Two response time-history curves should be constructed respectively；Local peaking's quantity on response time-history curves is counted, axle for vehicle is obtained
Or the quantity N of axis group.Here it by establishing modeling scheme, and realizes by response time-history curves the judgement of axle number, but needs
Multiple groups sensor is installed to form cross-referenced purpose, there is installation complexity, deficiency at high cost.
Summary of the invention
Goal of the invention of the invention is, in view of the above-mentioned problems, the method based on the bridge dynamic strain identification bicycle number of axle is provided,
By provide it is a kind of with simple installation, the recognition methods of at low cost, high-efficient axle number, by installation one dynamic stress
Sensor can meet the measurement demand to axle number.
In order to achieve the above objectives, the technical scheme adopted by the invention is that:
Based on the method for the bridge dynamic strain identification bicycle number of axle, the axle number of the vehicle is surveyed when bicycle passes through bridge
Measure step are as follows: arrange that dynamic strain sensor, the dynamic strain sensor are being arranged at 1/2 across footpath of bridge or maximum strain reaction
Being longitudinally arranged on section and along bridge；The dynamic strain sensor is sequentially connected high speed acquisition device and processing by shielded cable
Unit；The high speed acquisition device is acquired the strain signal of dynamic strain sensor and is sent in processing unit；The place
Reason unit is filtered to strain signal and obtains dynamic strain time-history curves, and extraction obtains the continuous of dynamic strain time-history curves
Section occurs for normal strain, and the dynamic strain time-history curves progress second order derivation in section occurs to normal strain and handles to obtain load system
Number curve；The axle number is the line segment quantity for being not less than amplitude threshold in load factor curve in ordinate negative direction；It is described
Amplitude threshold is to correspond to minimum in load factor curve when meeting the minimum single-point load of bridge survey demand by the bridge
The amplitude size in the corresponding curve section of single-point load.
This programme focuses on obtaining load factor curve to obtain the bicycle number of axle, senses particular by dynamic strain
Device obtains strain signal data, and filtering processing obtains dynamic strain time-history curves, establishes corresponding two-dimensional coordinate system, confirms and mention
Section occurs for the continuous normal strain passed through of picking up the car, and guarantees data accuracy, and reduce operation, improves working efficiency；
Second order derivation, which is carried out, by dynamic strain time-history curves normal strain occurred in section obtains load factor curve.Load factor is bent
When line is that Vehicle Axles pass through dynamic strain sensor test point, it is related to correspond to corresponding axle load size when strain variation value mutation
Coefficient curve, it is unique when caused strain when each axle passes through test point.Processing unit is according to load factor curve and width
Degree threshold value determines generation corresponding with the axle of vehicle or caused curve section on load factor curve, and calculates the curve section
Amplitude size, the amplitude size be maximum curve section amplitude size.Amplitude threshold is the minimum for meeting bridge survey demand
When single-point load passes through the bridge, the amplitude size in the corresponding curve section of minimum single-point load in load factor curve is corresponded to.
Here amplitude threshold guarantees to choose the correctness in Vehicle Axles homologous thread section, specifically screens out due in vehicle travel process
Vibration or bridge vibration caused by dynamic strain influence, reduce error influence.According in load factor curve, each axle of vehicle
Corresponding curve section be all it is unique, i.e., amplitude size is all unique, it can when obtaining bicycle and passing through bridge: axle number
For the line segment quantity for being not less than amplitude threshold in ordinate negative direction in load factor curve, or meet the curve area of condition
Between quantity.In contrast, the prior art obtains Vehicle Axles number and generally uses " Peak Intensity Method ", i.e. acquisition dynamic strain time-history curves,
By the number of peaks of software number dynamic strain time-history curves thus the Vehicle Axles number to vehicle；But there are following for this method
Insufficient: vehicle or bridge vibration can bring the variation of the wave crest quantity of dynamic strain time-history curves, and software statistics quantity, which can exist, to be missed
Difference.
Preferably, the bridge moving strain value size under vehicle effect and fluctuation amplitude are influenced by many, only work as vehicle
When gross weight reaches certain numerical value, under the conditions of existing technical level and instrument and equipment, just using dynamic strain identification one bicycle axle weight
Can there are certain applicability and accuracy.According to a large amount of verification experimental verification, when using high-resolution strain transducer (resolution ratio
No more than 0.1 μ ε) carry out dynamic strain acquisition, the bridge maximum strain under gross combination weight effectIt needs to meet(MmaxFor the maximal bending moment of the bridge strain testing section under vehicle effect, ymaxFor Edge Distance neutral axis
Maximum height, IE are the bending stiffness in strain testing section) when, this method has universal acceptable accuracy.
Preferably, due to bridge force-bearing structure type multiplicity, if any by curved beam bridge, the arcuately bridge being pressurized, skewed
Bridge and suspension bridge are drawn, is compared through test and force analysis, due to beam bridge wide applicability, structure stress is simple, power transmission
The features such as path is clear, active force and structural response linear relationship are good is had preferably again using dynamic strain identification bicycle axle
Mechanical foundation and higher applicability and accuracy of identification, experiment proves that and compare analysis, before beam bridge calculates across footpath L and vehicle
Maximum axle spacing l meets relational expression afterwardsWhen, recognition accuracy can meet bridge testing detection and health to one bicycle axle again
The needs of monitoring field.
Preferably, the dynamic strain time-history curves of actual measurement include biggish ambient noise and system noise, without noise reduction process
It cannot function as the initial data of identification, the influence factor and feature of comprehensive strain time history curve, the used filtering drop of this method
Method for de-noising is clipping Glitch Filter.
Preferably, the dynamic strain time-history curves pass through the dynamic strain of finite element stimulation theory by processing unit, and mention
Dynamic strain result is taken to be formed.
Due to the adoption of the above technical scheme, the invention has the following advantages:
1. the present invention is suitable for the installation detection of new and old bridge, bridge or pavement of road itself do not injured, is not necessarily to
Breaking off a friendship to lead to can complete to install.
2. the present invention only span centre or strain maximum cross-section arrange 1 strain monitoring section, so integral installation cost
It is low and easy for installation.
3. the present invention needs a dynamic strain time-history curves only to identify bicycle axle number, thus analyze identification more rapidly,
Fast, energy consumption is lower, increases the timeliness of data processing.
Detailed description of the invention
Fig. 1 is the strain time history curve after the simply supported beam noise reduction filtering under single Concentrated load of the invention.
Fig. 2 is the single order derived function figure of Fig. 1.
Fig. 3 is the second order derived function figure of Fig. 1.
Fig. 4 is single-point load of the present invention and multi-load strain time history curve.
Fig. 5 is multi-load strain time history curve single order derived function figure in Fig. 4.
Fig. 6 is multi-load strain time history curve second order derived function figure in Fig. 4.
Fig. 7 is actual measurement strain time history curve graph of the embodiment of the present invention.
Fig. 8 is the area that dynamic strain time-history curves envelope in section occurs for normal strain of the embodiment of the present invention.
Fig. 9 is vehicle wheelbase of embodiment of the present invention identification schematic diagram.
Figure 10 is strain time history curve derived function figure of the embodiment of the present invention.
Specific embodiment
It is further illustrated below in conjunction with specific implementation of the attached drawing to invention.
Based on the method for the bridge dynamic strain identification bicycle number of axle, the axle number of the vehicle is surveyed when bicycle passes through bridge
Measure step are as follows: arrange that dynamic strain sensor, the dynamic strain sensor are being arranged at 1/2 across footpath of bridge or maximum strain reaction
Being longitudinally arranged on section and along bridge；The dynamic strain sensor is sequentially connected high speed acquisition device and processing by shielded cable
Unit；The high speed acquisition device is acquired the strain signal of dynamic strain sensor and is sent in processing unit；The place
Reason unit is filtered to strain signal and obtains dynamic strain time-history curves, and extraction obtains the continuous of dynamic strain time-history curves
Section occurs for normal strain, and the dynamic strain time-history curves progress second order derivation in section occurs to normal strain and handles to obtain load system
Number curve；The axle number is the line segment quantity for being not less than amplitude threshold in load factor curve in ordinate negative direction；It is described
Amplitude threshold is to correspond to minimum in load factor curve when meeting the minimum single-point load of bridge survey demand by the bridge
The amplitude size in the corresponding curve section of single-point load.
Wherein, in order to improve the accuracy of above-mentioned measurement method, the above method should be applied preferentially is in bridge structure form
On the bridge of beam bridge, and the bridge should meet: the bridge maximum strain under gross combination weight effect MmaxFor the maximal bending moment of the bridge strain testing section under vehicle effect, ymaxFor Edge Distance neutral axis
Maximum height, IE are the bending stiffness in strain testing section.Beam bridge calculates across footpath L and maximum axle spacing l before and after vehicle and needs to meet
Relational expressionDynamic strain sensor is high-resolution strain detection testing device, and resolution ratio is less than 0.1 μ ε.Wherein,
Processing unit carries out the processing of clipping Glitch Filter to strain signal, obtains dynamic strain time-history curves.
In following further explanations, dynamic strain sensor, dynamic strain sensor are arranged specifically at 1/2 across footpath of bridge
It is sequentially connected high speed acquisition device and processing unit by shielded cable, processing unit is the computer to bridge monitoring.It is related to
Bridge parameter is all made ofThat is εmax=5；The bridge being related to is all satisfied with vehicle: before span of bridge L and vehicle
The proportionate relationship of maximum axle spacing l afterwards are as follows:Dynamic strain sensor is that resolution ratio is 0.01 μ ε；Principle and act are carried out with this
Example explanation.
In order to further appreciate that above scheme working principle, said here by simply supported beam strain time history curvilinear function
It is bright.
The strain time history curve after simply supported beam noise reduction filtering under single Concentrated load may be expressed as:
Wherein, t is the time, and L is that simply supported beam calculates across footpath, and v is speed, and measuring point sensor is at a distance from the cross-sectional neutral axis of place
Y, E are simply supported beam elasticity modulus, and I is simply supported beam cross sectional moment of inertia.
Above formula derivation is obtained:
Formula (2) is normal function discontinuous function, continues to have its derivation:
Make respectively formula (1), formula (2), formula (3) functional arrangement is as shown in Figure 1-3, the peak point pair of strain curve as shown in Figure 1
When what is answered is the section where load P passes through dynamic strain Sensor C, caused strain size, while where measuring point C
Section is also the separation that strain time history curve switchs to decline by rising, and the derivative value in corresponding diagram 2 becomes negative value from positive value.This
The reason of illustrating that the strain variation value at measuring point is mutated, and strain is caused to be mutated is that load successively acts on survey
The left and right sides in section where point, according to certain resilient relationship corresponding with power is strained, the size of strain variation value is in certain journey
To corresponding load it is directly proportional on degree, therefore is one by one for the size of second derivative values shown in Fig. 3 and load P
It is corresponding.
For the expression formula of its single order of function of strain and second order derived function under the effect of multiple unit concentrated loads, such as following formula
(4), formula (5), shown in formula (6).
Wherein, t is the time, and N is positive integer, i.e. unit load quantity.
Wherein, t is the time, and N is positive integer, i.e. unit load quantity.
Wherein, t is the time, and N is positive integer, i.e. unit load quantity.
It is assumed that N=4, i.e. hypothesis vehicle has 4 axle loads, and axle load equation is followed successively by ε1(x)、ε2(x)、ε3(x)
And ε4(x), then ε is made according to formula (1) and formula (4) respectively1(x)、ε2(x)、ε3(x) and ε4(x) single-point load functional arrangement and overloading
Lotus functional arrangement, as shown in part on Fig. 4.
Corresponding one complete strain of each axle load P as seen from Figure 4 influences line, with ε1(x) for, its table
Show load P1It is into bridge at origin x=0, in time x=d3Locate bridge out, corresponding d3Across footpath L is calculated for practical bridge.At this time
It is (d for the strain time history curve total length surveyed under multiaxis load1+d2+d3), the total size of actual measurement strain stress (x) is (ε1
(x)+ε2(x)+ε3(x)+ε4(x)).Time change x=vt is done to x, v is speed, just obtains total strain time history curve graph as schemed
Shown in 4 lower parts.Corresponding axis when just obtaining strain variation value mutation to the actual measurement secondary derivation of overall strain time-history curves function again
The load factor curve of load, as shown in Figure 6.It is positive here it is worth noting that leading value for strain time history curve second order
When be that measuring point strain value starts to become larger at span centre or maximum strain reaction caused when axle load is loaded into bridge or out bridge, negative value is corresponding
Be then axle load by measuring point when caused strain start to become smaller.
The measuring point of time-history curves second-order derivative value when passing through to(for) unit axle load is constant for one in extremely short time Δ t
Constant, that is, form the curve section determined in load factor curve, which is unit axle load by surveying
Time-history curves second-order derivative value when point.
When measuring based on the above principles to bridge, dynamic strain sensor is arranged at 1/2 across footpath of bridge.High speed is adopted
Storage is acquired the strain signal of dynamic strain sensor and is sent in processing unit.Processing unit carries out strain signal
Clipping Glitch Filter handles and obtains dynamic strain time-history curves, and extraction obtains the continuous normal strain generating region of dynamic strain time-history curves
Between, and the dynamic strain time-history curves progress second order derivation in section is occurred to normal strain and handles to obtain load factor curve.Processing
Unit is determining according to load factor curve and amplitude threshold and extracts curve corresponding with the axle of vehicle on load factor curve
Section calculates the amplitude size in the curve section.Amplitude threshold is to meet the minimum single-point load of bridge survey demand described in
When bridge, the amplitude size in the corresponding curve section of minimum single-point load in load factor curve is corresponded to.Curve section meets it
Amplitude size is not less than amplitude threshold size.
Wherein, Vehicle Axles number are as follows: be not less than the line segment of amplitude threshold in load factor curve in ordinate negative direction
Quantity or the number in curve section.Amplitude threshold is when meeting the minimum single-point load of bridge survey demand by the bridge,
The amplitude size in the corresponding curve section of minimum single-point load in corresponding load factor curve isParameter physics
Meaning is as hereinbefore.Here minimum single-point load should meet in MmaxIt is the magnitude of load that bridge is subject to when being minimized.
Weight size or the Vehicle Axles weight of each axle load of the vehicle can be obtained by formula (7):
It is wherein N is positive integer, i.e., unit load quantity, GVW are gross combination weight.
It can be obtained by load factor curve, axle weight is=(all axles of the corresponding amplitude size/vehicle of the axle of vehicle
The sum of corresponding amplitude size) × gross combination weight；Gross combination weight is the interior dynamic strain time-history curves envelope in normal strain generating region
Area × constant coefficient.Constant coefficient=standard tonnage vehicle car weight/(the normal strain that the vehicle of standard tonnage is obtained by the bridge
The area of dynamic strain time-history curves envelope in section occurs).
The vehicle speed can be obtained by formula (8):
L is that bridge calculates across footpath, t1It is axle from bridge is entered to the time for reaching dynamic strain Sensor.
It can be obtained by load factor curve, speed is first curve area of the span of bridge/normal strain generating region negative direction
Between occur to normal strain twice of spacing of section starting point
The vehicle wheelbase can be obtained by formula (9):
Δ d=v (tn-tn-1)=v Δ t formula (9)；
V is speed, and Δ t is the time interval that axle passes through dynamic strain Sensor.
It can be obtained by load factor curve, the wheelbase=[across footpath of the bridge/(positive direction first in normal strain generation section
Twice of spacing of section starting point occurs to normal strain for a curve section)] spacing between × two adjacent curve sections.
It is further illustrated below with actual tests.
Same car weight is taken first, here by taking 20 tons of model (4.6kg) vehicles as an example, 20T -20Hz simply supported beam span centre is taken to strain
AI01 (left side) and AI02 (the right), two o'clock carry out the identification of speed, identify and are divided into 3 grades, 20km/h (0.29m/s),
40km/h (0.56m/s) and 60km/h (0.83m/s).
In Fig. 7, horizontal axis indicates the time, and sample frequency 20Hz is 0.05 second per small lattice time interval, and wherein AI01 line is being just
Strain value is 127 lattice, and the time of passing a bridge is 127x0.05=6.35 seconds；AI02 line normal strain be 128 lattice, pass a bridge the time be 128 ×
0.05=6.4 seconds.
Relative error
As shown in figure 8, acquiring 20t car weight (ε-vt) to surround area using related software is 77.46m × ε, constant coefficient is
77.46/4.6=16.84.
The identification of car weight now divides four grades in terms of the area of strain time history curve and reference axis envelope with car weight identification,
20t (4.6kg), 40t (8.85kg), 60t (13.42kg) and 80t (17.9kg), wherein 20t is that standard is heavily loaded, 40t, 60t,
80t is identification heavy duty.
| Car weight | Area m × ε | Identify car weight kg | Opposite accidentally % |
| 20t(4.6kg) | 77.46 | 4.6 | 0 |
| 40t(8.85kg) | 148.65 | 8.82 | 0.3 |
| 60t(13.42kg) | 234.72 | 13.94 | 3.8 |
| 80t(17.9kg) | 303.65 | 18.03 | 0.7 |
As seen from the above table, the available vehicle of area surrounded by constant coefficient and strain time history curve and reference axis is total
Weight, error can be within the specified scope.
As shown in figure 9, same car weight, different speeds carry out wheelbase identification with same wheelbase.Now with 20t-20Hz, vehicle
Two axle distance of axis is 18cm, and speed is respectively 20km/h (0.29m/s), 40km/h (0.56m/s) and 60km/h (0.83m/
S) identification of wheelbase is carried out respectively.
| 20t-20Hz | Recognition speed m/s | Recognition time s | Identify wheelbase m | Relative error % |
| 20km/h | 0.32 | 0.5 | 0.16 | 11 |
| 40km/h | 0.606 | 0.25 | 0.152 | 15 |
| 60km/h | 0.869 | 0.2 | 0.174 | 3.4 |
With same wheelbase 18cm, same car weight 20t (4.6kg), for same speed 20km/h (0.29m/s).Figure 10 is
Trolley is with 20km/h, front axle weight 1.38kg rear axle weight 3.22kg, sample frequency extracted strain time history curve graph when being 20Hz,
If being 0.05s (i.e. former sample frequency) to time-history curves derivation function according to step-length, due to each in original timeamplitude map
There is the case where fluctuation up and down in point, resulting result is very big by curve influence of noise, as Figure 10 is marked；Therefore it will cause secondary lead
Also fluctuation, recognition result are difficult to achieve the desired results function up and down therewith.Here it can be filtered out and be met by setting amplitude threshold
Recognition result, amplitude threshold size is that can measure to bear in the load factor curve that minimum tonnage is determined by the bridge
Minimum amplitude value on direction.As shown in Figure 10, meet curve section definition only there are two, mark 1, mark 2, this and practical vehicle
The number of axle is consistent.
It is obtained according to formula (7):
Relative error,
To the error of axle weight 7.6%, receiving in error range.As it can be seen that vehicle can be measured to obtain under above scheme
Axle weight, meets the needs of bridge survey.
Above description is the detailed description for the present invention preferably possible embodiments, but embodiment is not limited to this hair
Bright patent claim, it is all the present invention suggested by technical spirit under completed same changes or modifications change, should all belong to
In the covered the scope of the patents of the present invention.
Claims (5)
1. the method based on the bridge dynamic strain identification bicycle number of axle, which is characterized in that when bicycle passes through bridge to the vehicle
Axle number measuring process are as follows: dynamic strain sensor is arranged at 1/2 across footpath of bridge or maximum strain reaction, the dynamic strain passes
Sensor setting section on and being longitudinally arranged along bridge；The dynamic strain sensor is sequentially connected high speed by shielded cable and adopts
Storage and processing unit；The high speed acquisition device is acquired the strain signal of dynamic strain sensor and is sent to processing unit
In；The processing unit carries out the processing of clipping Glitch Filter to strain signal and obtains dynamic strain time-history curves, and extraction is moved
Section occurs for the continuous normal strain of strain time history curve, and the dynamic strain time-history curves in section occur to normal strain and carry out second order
Derivation handles to obtain load factor curve；The axle number is to be not less than amplitude threshold in load factor curve in ordinate negative direction
The line segment quantity of value；The amplitude threshold is when meeting the minimum single-point load of bridge survey demand by the bridge, to correspond to
The amplitude size in the corresponding curve section of minimum single-point load in load factor curve.
2. the method according to claim 1 based on the bridge dynamic strain identification bicycle number of axle, it is characterised in that: described move is answered
Become time-history curves by processing unit by the dynamic strain of finite element stimulation theory, and extracts dynamic strain result and formed.
3. the method according to claim 1 based on the bridge dynamic strain identification bicycle number of axle, it is characterised in that: the vehicle
Bridge maximum strain under gross weight effectIt needs to meetMmaxFor the bridge under vehicle effect
The maximal bending moment in strain testing section, ymaxFor the maximum height of Edge Distance neutral axis, IE is that the bending resistance in strain testing section is rigid
Degree.
4. the method according to claim 1 based on the bridge dynamic strain identification bicycle number of axle, it is characterised in that: the bridge
Structure type is beam bridge, and maximum axle spacing l need to meet relational expression before and after beam bridge calculates across footpath L and vehicle
5. the method according to claim 1 based on the bridge dynamic strain identification bicycle number of axle, it is characterised in that: described move is answered
Change sensor is high-resolution strain detection testing device, and resolution ratio is less than 0.1 μ ε.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710089866.8A CN106895900B (en) | 2017-02-20 | 2017-02-20 | Method based on the bridge dynamic strain identification bicycle number of axle |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710089866.8A CN106895900B (en) | 2017-02-20 | 2017-02-20 | Method based on the bridge dynamic strain identification bicycle number of axle |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106895900A CN106895900A (en) | 2017-06-27 |
| CN106895900B true CN106895900B (en) | 2019-08-23 |
Family
ID=59184577
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710089866.8A Active CN106895900B (en) | 2017-02-20 | 2017-02-20 | Method based on the bridge dynamic strain identification bicycle number of axle |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106895900B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106895900B (en) * | 2017-02-20 | 2019-08-23 | 广西交通科学研究院有限公司 | Method based on the bridge dynamic strain identification bicycle number of axle |
| CN107389168A (en) * | 2017-07-18 | 2017-11-24 | 重庆交通大学 | A kind of vehicle for bridge moves the recognition methods of load |
| CN107393312B (en) * | 2017-08-22 | 2023-04-07 | 湖南大学 | Method and system for identifying axle and speed of vehicle running on bridge |
| CN110874932B (en) * | 2018-08-29 | 2021-06-15 | 北京万集科技股份有限公司 | Method, device and system for identifying number of axles |
| CN112132227B (en) * | 2020-09-30 | 2024-04-05 | 石家庄铁道大学 | Bridge train load action time course extraction method and device and terminal equipment |
| CN113280747B (en) * | 2021-04-22 | 2023-05-12 | 杭州申昊科技股份有限公司 | System and method for judging rail safety based on fiber bragg grating sensor |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN2888442Y (en) * | 2005-08-23 | 2007-04-11 | 北京华力兴科技发展有限责任公司 | Safety check-up system for roll on/roll off transportation vehicle |
| CN102735320A (en) * | 2012-07-19 | 2012-10-17 | 广西交通科学研究院 | Method for identifying weights of cars based on dynamic strain of bridges |
| JP2014228480A (en) * | 2013-05-24 | 2014-12-08 | 国立大学法人福井大学 | Device and method for passing vehicle weight analysis processing |
| CN105651338A (en) * | 2016-03-01 | 2016-06-08 | 湖南大学 | Axle recognition method and system for axle beam |
| CN106895900A (en) * | 2017-02-20 | 2017-06-27 | 广西交通科学研究院有限公司 | The method that the bicycle number of axle is recognized based on bridge dynamic strain |
-
2017
- 2017-02-20 CN CN201710089866.8A patent/CN106895900B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN2888442Y (en) * | 2005-08-23 | 2007-04-11 | 北京华力兴科技发展有限责任公司 | Safety check-up system for roll on/roll off transportation vehicle |
| CN102735320A (en) * | 2012-07-19 | 2012-10-17 | 广西交通科学研究院 | Method for identifying weights of cars based on dynamic strain of bridges |
| JP2014228480A (en) * | 2013-05-24 | 2014-12-08 | 国立大学法人福井大学 | Device and method for passing vehicle weight analysis processing |
| CN105651338A (en) * | 2016-03-01 | 2016-06-08 | 湖南大学 | Axle recognition method and system for axle beam |
| CN106895900A (en) * | 2017-02-20 | 2017-06-27 | 广西交通科学研究院有限公司 | The method that the bicycle number of axle is recognized based on bridge dynamic strain |
Non-Patent Citations (1)
| Title |
|---|
| 应变式汽车轴重动态测量系统性能增长研究;周志峰;《中国博士学位论文全文数据库 工程科技II辑》;20080630(第6期);第4、85-102页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106895900A (en) | 2017-06-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106840337B (en) | Method based on bridge dynamic strain identification one bicycle axle weight | |
| CN106895900B (en) | Method based on the bridge dynamic strain identification bicycle number of axle | |
| CN106710242B (en) | Method based on bridge dynamic strain identification fleet vehicle number | |
| CN102735320B (en) | Method for identifying weights of cars based on dynamic strain of bridges | |
| CN110689723B (en) | Truck overload identification method based on power distribution and self-learning | |
| CN108515984B (en) | Wheel damage detection method and device | |
| CN104792937B (en) | Bridge head bump detection evaluation method based on vehicle-mounted gravitational acceleration sensor | |
| CN105651338B (en) | The recognition methods of axletree quantity, wheelbase recognition methods and system for bridge | |
| CN102901550A (en) | Method for implementing vehicle-mounted dynamic weighing | |
| CN104309435A (en) | System and method for online recognition of road unevenness | |
| CN105258770B (en) | Road vehicle dynamic weighing method and equipment | |
| CN106872005B (en) | Method based on bridge dynamic strain identification fleet's bicycle car weight | |
| WO2019153876A1 (en) | Method and device for evaluating severity of head injury incurred by cyclist after impact with road surface, and testing method | |
| CN106885551B (en) | Based on bridge dynamic strain identification one bicycle axle away from method | |
| CN110926735A (en) | Bridge structure rapid diagnosis method based on multidimensional dynamic parameters | |
| CN106871847B (en) | Method based on bridge dynamic strain identification fleet's bicycle spacing | |
| CN109059840A (en) | A kind of city rail vehicle wheel out of round is along detection method | |
| CN107393312A (en) | Vehicle Axles and the method and system of speed are travelled on a kind of identification bridge | |
| CN112697249B (en) | Dynamic vehicle overrun determination method and determination system | |
| CN110987499A (en) | Bridge dynamic load test method | |
| CN205388511U (en) | Dynamic vehicle weighing system | |
| CN109357823A (en) | A method of actual measurement road surface structare layer axis carries active position, the maximum strain of layer bottom and axis and carries speed | |
| CN112781720B (en) | Road condition evaluation system and evaluation method thereof | |
| CN111174885A (en) | Method for classifying and sectionally acquiring signals of vehicle dynamic weighing sensor | |
| CN103884408A (en) | Axle type detecting method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP03 | Change of name, title or address |
Address after: No.6, Gaoxin 2nd Road, XiXiangTang District, Nanning City, Guangxi Zhuang Autonomous Region Patentee after: Guangxi Jiaoke Group Co.,Ltd. Country or region after: China Address before: No. 6 Gaoxin 2nd Road, Nanning City, Guangxi Zhuang Autonomous Region Patentee before: GUANGXI TRANSPORTATION RESEARCH & CONSULTING Co.,Ltd. Country or region before: China |
|
| CP03 | Change of name, title or address |