CN106908128A - The recognition methods of the bridge floor multiaxis traveling load based on Landweber iterative methods - Google Patents
The recognition methods of the bridge floor multiaxis traveling load based on Landweber iterative methods Download PDFInfo
- Publication number
- CN106908128A CN106908128A CN201710061242.5A CN201710061242A CN106908128A CN 106908128 A CN106908128 A CN 106908128A CN 201710061242 A CN201710061242 A CN 201710061242A CN 106908128 A CN106908128 A CN 106908128A
- Authority
- CN
- China
- Prior art keywords
- bridge
- delta
- equation
- vehicle
- 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.)
- Pending
Links
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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a kind of recognition methods of the bridge floor multiaxis traveling load of Landweber iterative methods, comprise the following steps:1), in bridge bottom surface correspondence position x1,x2,…xmM displacement transducer is pasted at place respectively, measures bridge floor multiaxis mobile vehicle load fkT () displacement of t at x position is v (x, t), k=1,2,3 ..., be the vehicle number of axle;2) oscillatory differential equation, is set up;3), equation (1) is solved;4) bridge, is set up under k axle Vehicle Loads, and multiaxis traveling load system equation is recognized by dynamic respond;5) exact value of multiaxis traveling load, is tried to achieve using Landweber iterative methods.The present invention need to only measure bridge displacement response and can recognize that bridge floor multiaxis traveling load, and recognition methods is simple and precision is higher, with good feasibility, can be widely applied to the moving load identification of all kinds bridge.
Description
Technical field
Bridge floor multiaxis is recognized by bridge displacement the invention belongs to bridge floor moving load identification technical field, more particularly to one kind
The method of traveling load.
Background technology
Bridges in Our Country present situation is " rebuilding light supporting ", and from 1999 to 2013, domestic media disclosed report China because various
The bridge that reason collapses up to more than 110 seats, wherein the bridge collapse not caused including Wenchuan earthquake still.Cause bridge damnification and break
Bad reason can be summarized as external factor and internal factor, wherein in external factor due to automobile overload cause bridge fatigue damage and
Durability reduction occupies leading position, and internal factor is then mainly the reduction of bridge itself bearing capacity and the strength of materials is degenerated.
With the explosive growth of highway in China traffic, the actual vehicle flowrate for bearing of many bridges increases compared with Earlier designs value
A lot, the increase of speed and car weight can have a negative impact to bridge, the multiple-axle vehicle and large-scale multiple-axle vehicle especially overloads
Appearance substantially exacerbate the risk of bridge collapse.
Highway in China transfinites to stand and make many work in control vehicle overload method, but measuring method is to use ground mostly at present
Pound technology, i.e., weigh to realize the measurement of gross combination weight by parking.Under the trend of development rapid transit, how in vehicle row
Accurate vehicle load has important practical meaning in engineering during sailing, especially to the accurate measurement of each axle load of multiaxis lorry
Security and durability to protecting bridge all have very great help.
It is identified for conventional two-axle car more than existing moving load identification technology, it is impossible to which multiple-axle vehicle load is entered
Row identification, therefore it is badly in need of a kind of method that can be identified to bridge floor multiaxis mobile vehicle load.
The content of the invention
Only need to measure the identification bridge floor multiaxis by bridge displacement is responded rapidly and efficiently it is an object of the invention to provide a kind of
Mobile vehicle load, accuracy of identification is high and does not influence bridge floor vehicle normal pass.
To reach above-mentioned purpose, the technical solution adopted by the present invention is:It is a kind of based on Landweber iterative methods
The recognition methods of bridge floor multiaxis traveling load, comprises the following steps:
1), in bridge bottom surface correspondence position x1,x2,…xmM displacement transducer is pasted at place respectively, measures the shifting of bridge floor multiaxis
Dynamic vehicular load fkT () displacement of t at x position is v (x, t), k=1,2,3 ..., be the vehicle number of axle;
2) vehicle-bridge system oscillatory differential equation, is set up:Bridge length is taken for L, bending rigidity is EI, bridge unit length matter
It is ρ to measure, it is considered to which viscous damping simultaneously takes damped coefficient for C, ignores detrusion and the rotary inertia of bridge, bridge floor multiaxis locomotive
Load fk (t) is moved right with speed c from beam left end supporting, then the oscillatory differential equation of vehicle-bridge system is:
Wherein δ (x-ct) is Dirac function;
The boundary condition of equation (1) is:
V (0, t)=0, v (L, t)=0,
3), equation (1) is solved;
4) bridge, is set up under k axle Vehicle Loads, and multiaxis traveling load system equation is recognized by dynamic respond:
v(m×1)=S(m×k)·f(k×1) (2)
v(m×1)It is traveling load fkT () is in x1,x2,…xmThe actual displacement at place, and m >=k;S(m×k)It is known system square
Battle array;f(k×1)It is required k axle traveling loads;
The discrete form of formula (2) is expressed as:
Wherein
5) exact value of multiaxis traveling load, is tried to achieve using Landweber iterative methods;
When bridge displacement response is measured, made an uproar comprising less in the dynamic respond generally measured due to measurement noise interference
Sound, actually measured dynamic respond is vδ, i.e., equation Sf=v right-hand vectors be typically with disturbance data vδ, and meet
||v-vδ||≤δ (4)
Wherein δ is noise threshold limit value, and this noise threshold limit value is determined by environmental factor, and span is 0 to 1
Between;
Landweber Iterations are:
Wherein S*It is the adjoint matrix of vehicle bridge moving load identification sytem matrix S,
Stopping criterion for iteration is:
Wherein C is the influence coefficient of noise threshold limit value, is a constant, and this value is manually set, according to identification essence
Degree requires different freely to set;
Taking multiaxis traveling load initial value isIteration can be written as:
Wherein l is unit matrix, if initial value is taken asIteration can be written as:
WhereinThe multiaxis traveling load for as being recognized based on Landweber alternative manners.
Described step 3) in equation (1) is solved comprise the following steps that it is described:
Based on modal superposition principle, it is assumed that the n-th order Mode Shape function of bridge isThen equation (1)
Solution be expressed as:
Matrix form is:
Here n is mode number, qnT () (n=1,2 ... ∞) are n-th order modal displacements, equation (12) is substituted into equation (1),
And x is integrated in [0, L], using boundary condition and Dirac function characteristic, vehicle-bridge system oscillatory differential equation qn
T () is expressed as:
HereIt is qnThe second dervative of (t),It is qnThe first derivative of (t), Respectively circular frequency, damping rate and bridge floor mobile vehicle load mod table
Up to formula;
As vehicle has a k axletree, and k-th axletree to first distance of axletree be xk(x1=0), then equation (14)
It is written as:
Then corresponding to the modal displacement at m measuring point can be expressed as by equation (13):
X on bridge1,x2,…xmThe speed at place is tried to achieve by the once differentiation of displacement:
Further, x on bridge1,x2,…xmThe acceleration at place is tried to achieve by the second differential of displacement:
Similarly, x on beam1,x2,…xmThe moment of flexure at place can utilize relational expressionTry to achieve:
If f1,f2,…,fkIt is each axle correspondence load of known k axles vehicle, ignores the influence of damping, then the solution of equation (1) can
It is expressed as:
Wherein
The present invention can respond identification multiaxis traveling load, the method letter of measurement bridge displacement response by measuring bridge displacement
List and precision is higher, therefore identification bridge floor traveling load is responded by bridge displacement there is good feasibility and accuracy of identification energy
Guarantee is accessed, dynamic respond need to be only obtained using method proposed by the present invention and be can recognize that bridge floor multiaxis traveling load, therefore
Recognition methods proposed by the present invention has good feasibility, can be widely applied to the moving load identification of all kinds bridge.
But in actual application, the dynamic respond of measurement is the influence of noise, Landweber iteration with small distracter
Method can accurately recognize multiaxis traveling load by limiting stopping criterion for iteration by the dynamic respond comprising noise, in the complicated ring of measurement
There is larger advantage when border or more noise jamming.
Brief description of the drawings
Fig. 1 is flow chart of the method for the present invention.
Specific embodiment
As shown in figure 1, the invention discloses a kind of identification of the bridge floor multiaxis traveling load based on Landweber iterative methods
Method, comprises the following steps:
1), in bridge bottom surface correspondence position x1,x2,…xmM displacement transducer is pasted at place respectively, measures the shifting of bridge floor multiaxis
Dynamic vehicular load fk(t) at x position t displacement for v (x, t), k=1,2,3 ... be the vehicle number of axle;
2) vehicle-bridge system oscillatory differential equation, is set up:Bridge length is taken for L, bending rigidity is EI, bridge unit length matter
It is ρ to measure, it is considered to which viscous damping simultaneously takes damped coefficient for C, ignores detrusion and the rotary inertia of bridge, bridge floor multiaxis locomotive
Load fk (t) is moved right with speed c from beam left end supporting, then the oscillatory differential equation of vehicle-bridge system is:
Wherein δ (x-ct) is Dirac function;
The boundary condition of equation (1) is:
V (0, t)=0, v (L, t)=0,
3), equation (1) is solved;
31), based on modal superposition principle, it is assumed that the n-th order Mode Shape function of beam isThen equation
(1) solution is represented by:
Matrix form is:
Here n is mode number, qnT () (n=1,2 ... ∞) are n-th order modal displacements, equation (12) is substituted into equation (1),
And x is integrated in [0, L], using boundary condition and Dirac function characteristic, vehicle-bridge system oscillatory differential equation qn
T () is expressed as:
HereIt is qnThe second dervative of (t),It is qnThe first derivative of (t), Respectively circular frequency, damping rate and bridge floor mobile vehicle load mod table
Up to formula.
As vehicle has a k axletree, and k-th axletree to first distance of axletree be xk(x1=0), then equation (14)
It is written as:
Then corresponding to the modal displacement at m measuring point can be expressed as by equation (13):
X on bridge1,x2,…xmThe speed at place is tried to achieve by the once differentiation of displacement:
Further, x on bridge1,x2,…xmThe acceleration at place is tried to achieve by the second differential of displacement:
Similarly, x on beam1,x2,…xmThe moment of flexure at place can utilize relational expressionTry to achieve:
If f1,f2,…,fkIt is each axle correspondence load of known k axles vehicle, ignores the influence of damping, then the solution of equation (1) can
It is expressed as:
Wherein
4) bridge, is set up under k axle Vehicle Loads, and multiaxis traveling load system equation is recognized by dynamic respond:
v(m×1)=S(m×k)·f(k×1) (2)
v(m×1)It is traveling load fkT () is in x1,x2,…xmThe actual displacement at place (is exactly measured position in step (1)
Move), and m >=k;S(m×k)It is known sytem matrix;f(k×1)It is required k axle traveling loads;
The discrete form of formula (2) is expressed as
Wherein
5) exact value of multiaxis traveling load, is tried to achieve using Landweber iterative methods;
When bridge displacement response is measured, made an uproar comprising less in the dynamic respond generally measured due to measurement noise interference
Sound, actually measured dynamic respond is vδ, i.e., equation Sf=v right-hand vectors be typically with disturbance data vδ, and meet
||v-vδ||≤δ (4)
Wherein δ is noise threshold limit value, and this noise threshold limit value is determined by environmental factor, and span is 0 to 1
Between;
Landweber Iterations are:
Wherein S*It is the adjoint matrix of vehicle bridge moving load identification sytem matrix S,
Stopping criterion for iteration is:
Wherein C is the influence coefficient of noise threshold limit value, is a constant, and this value is manually set, according to identification essence
Degree requires different freely to set;
Taking multiaxis traveling load initial value isIteration can be written as:
Wherein l is unit matrix, if initial value is taken asIteration can be written as:
WhereinThe multiaxis traveling load for as being recognized based on Landweber alternative manners.
Claims (2)
1. a kind of recognition methods of the bridge floor multiaxis traveling load based on Landweber iterative methods, it is characterised in that:Including following
Step:
1), in bridge bottom surface correspondence position x1,x2,…xmM displacement transducer is pasted at place respectively, measures bridge floor multiaxis locomotive
Load fkT () displacement of t at x position is v (x, t), k=1,2,3 ..., be the vehicle number of axle;
2) vehicle-bridge system oscillatory differential equation, is set up:Bridge length is taken for L, bending rigidity is EI, and bridge linear mass is
ρ, it is considered to which viscous damping simultaneously takes damped coefficient for C, ignores detrusion and the rotary inertia of bridge, bridge floor multiaxis mobile vehicle lotus
Carry fkT () is moved right with speed c from beam left end supporting, then the oscillatory differential equation of vehicle-bridge system is:
Wherein δ (x-ct) is Dirac function;
The boundary condition of equation (1) is:
V (0, t)=0, v (L, t)=0,V (x, 0)=0,
3), equation (1) is solved;
4) bridge, is set up under k axle Vehicle Loads, and multiaxis traveling load system equation is recognized by dynamic respond:
v(m×1)=S(m×k)·f(k×1) (2)
v(m×1)It is traveling load fkT () is in x1,x2,…xmThe actual displacement at place, and m >=k;S(m×k)It is known sytem matrix;
f(k×1)It is required k axle traveling loads;
The discrete form of formula (2) is expressed as:
Wherein
5) exact value of multiaxis traveling load, is tried to achieve using Landweber iterative methods;
When bridge displacement response is measured, comprising less noise in the dynamic respond generally measured due to measurement noise interference,
Actually measured dynamic respond is v6, i.e., equation Sf=v right-hand vectors be typically with disturbance data v6, and meet
||v-vδ||≤δ (4)
Wherein δ be noise threshold limit value, this noise threshold limit value is determined by environmental factor, span for 0 to 1 it
Between;
Landweber Iterations are:
Wherein S*It is the adjoint matrix of vehicle bridge moving load identification sytem matrix S,
Stopping criterion for iteration is:
Wherein C is the influence coefficient of noise threshold limit value, is a constant, and this value is manually set, will according to accuracy of identification
Ask different freely to set;
Taking multiaxis traveling load initial value isIteration can be written as:
Wherein l is unit matrix, if initial value is taken asIteration can be written as:
WhereinThe multiaxis traveling load for as being recognized based on Landweber alternative manners.
2. the recognition methods of the bridge floor multiaxis traveling load based on Landweber iterative methods as claimed in claim 1, its feature
It is:Described step 3) in equation (1) is solved comprise the following steps that it is described:
Based on modal superposition principle, it is assumed that the n-th order Mode Shape function of bridge isThe then solution of equation (1)
It is expressed as:
Matrix form is:
Here n is mode number, qnT () (n=1,2 ... ∞) are n-th order modal displacements, by equation (12) substitution equation (1), and
X is integrated in [0, L], using boundary condition and Dirac function characteristic, vehicle-bridge system oscillatory differential equation qn(t) table
It is shown as:
HereIt is qnThe second dervative of (t),It is qnThe first derivative of (t), Respectively circular frequency, damping rate and bridge floor mobile vehicle load mod table
Up to formula;
As vehicle has a k axletree, and k-th axletree to first distance of axletree be xk(x1=0), then equation (14) is written as:
Then corresponding to the modal displacement at m measuring point can be expressed as by equation (13):
X on bridge1,x2,…xmThe speed at place is tried to achieve by the once differentiation of displacement:
Further, x on bridge1,x2,…xmThe acceleration at place is tried to achieve by the second differential of displacement:
Similarly, x on beam1,x2,…xmThe moment of flexure at place can utilize relational expressionTry to achieve:
If f1,f2,…,fkIt is each axle correspondence load of known k axles vehicle, ignores the influence of damping, then the solution of equation (1) can be represented
For:
Wherein
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710061242.5A CN106908128A (en) | 2017-01-25 | 2017-01-25 | The recognition methods of the bridge floor multiaxis traveling load based on Landweber iterative methods |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710061242.5A CN106908128A (en) | 2017-01-25 | 2017-01-25 | The recognition methods of the bridge floor multiaxis traveling load based on Landweber iterative methods |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106908128A true CN106908128A (en) | 2017-06-30 |
Family
ID=59207627
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710061242.5A Pending CN106908128A (en) | 2017-01-25 | 2017-01-25 | The recognition methods of the bridge floor multiaxis traveling load based on Landweber iterative methods |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106908128A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108520099A (en) * | 2018-03-12 | 2018-09-11 | 榆林学院 | The measurement method of dynamic load(loading) suffered by a kind of ramp Truck |
CN108897958A (en) * | 2018-06-30 | 2018-11-27 | 中国石油大学(华东) | Oscillator acts on cased vibration analysis and describes method |
CN109948254A (en) * | 2019-03-21 | 2019-06-28 | 华东交通大学 | One kind bridge displacement extreme value response analysis method under traveling load excitation |
CN110940405A (en) * | 2018-09-25 | 2020-03-31 | 比亚迪股份有限公司 | Method and device for calculating total mass of vehicle |
CN115900906A (en) * | 2022-06-15 | 2023-04-04 | 东南大学 | Bridge dynamic weighing method based on mid-span boundary beam measuring point strain |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104598753A (en) * | 2015-02-06 | 2015-05-06 | 华北水利水电大学 | Bridge moving vehicle load recognition method based on Brakhage V method |
CN104598931A (en) * | 2015-02-06 | 2015-05-06 | 华北水利水电大学 | Bridge moving vehicle load identifying method based on Kaczmarz algebraic iterative reconstruction method |
CN104615888A (en) * | 2015-02-06 | 2015-05-13 | 华北水利水电大学 | Method for recognizing moving vehicle load on bridge based on generalized minimal residual algorithm |
-
2017
- 2017-01-25 CN CN201710061242.5A patent/CN106908128A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104598753A (en) * | 2015-02-06 | 2015-05-06 | 华北水利水电大学 | Bridge moving vehicle load recognition method based on Brakhage V method |
CN104598931A (en) * | 2015-02-06 | 2015-05-06 | 华北水利水电大学 | Bridge moving vehicle load identifying method based on Kaczmarz algebraic iterative reconstruction method |
CN104615888A (en) * | 2015-02-06 | 2015-05-13 | 华北水利水电大学 | Method for recognizing moving vehicle load on bridge based on generalized minimal residual algorithm |
Non-Patent Citations (2)
Title |
---|
T.H.T.CHAN ET AL: "AN INTERPRETIVE METHOD FOR MOVING FORCE INDENTIFICATION", 《JOURNAL OF SOUND AND VIBRATION》 * |
常晓通 等: "Landweber迭代正则化方法在动态载荷识别中的应用", 《应用数学和力学》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108520099A (en) * | 2018-03-12 | 2018-09-11 | 榆林学院 | The measurement method of dynamic load(loading) suffered by a kind of ramp Truck |
CN108897958A (en) * | 2018-06-30 | 2018-11-27 | 中国石油大学(华东) | Oscillator acts on cased vibration analysis and describes method |
CN108897958B (en) * | 2018-06-30 | 2022-09-13 | 中国石油大学(华东) | Vibration analysis description method of casing under action of oscillator |
CN110940405A (en) * | 2018-09-25 | 2020-03-31 | 比亚迪股份有限公司 | Method and device for calculating total mass of vehicle |
CN110940405B (en) * | 2018-09-25 | 2022-07-15 | 比亚迪股份有限公司 | Method and device for calculating total mass of vehicle |
CN109948254A (en) * | 2019-03-21 | 2019-06-28 | 华东交通大学 | One kind bridge displacement extreme value response analysis method under traveling load excitation |
CN115900906A (en) * | 2022-06-15 | 2023-04-04 | 东南大学 | Bridge dynamic weighing method based on mid-span boundary beam measuring point strain |
CN115900906B (en) * | 2022-06-15 | 2023-07-25 | 东南大学 | Bridge dynamic weighing method based on strain of measuring point of midspan boundary beam |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106908128A (en) | The recognition methods of the bridge floor multiaxis traveling load based on Landweber iterative methods | |
CN101281117B (en) | Wide span rail traffic bridge damnification recognition method | |
CN104615888B (en) | Method for recognizing moving vehicle load on bridge based on generalized minimal residual algorithm | |
CN105910692B (en) | Vehicle load measurement method based on acceleration transducer | |
CN104598931B (en) | Bridge based on Kaczmarz algebraically iterative reconstruction approach moves vehicle Load Identification Methods | |
CN104077459B (en) | A kind of automobile chassis handling stability energy analysis model method for building up based on suspension KC characteristics | |
CN102874277A (en) | Estimation of wheel rail interaction forces | |
CN106840336A (en) | The recognition methods of the bridge floor multiaxis traveling load based on QMR algorithms | |
Calvo et al. | Influence of a shock absorber model on vehicle dynamic simulation | |
CN106768234A (en) | Make the recognition methods of the bridge floor multiaxis traveling load of usage based on classical row | |
CN105389458A (en) | Highway random traffic flow wind-vehicle-bridge coupling computation method | |
Velinsky et al. | Vehicle energy dissipation due to road roughness | |
Wei et al. | On-board measurement of aerodynamic loads for high-speed trains negotiating transitions in windbreak walls | |
Zou et al. | Analysis of the effects of wind barrier on driving safety and comfort of vehicles on long-span bridges under crosswinds | |
Li et al. | Effects of parameters on dynamic responses for a heavy vehicle-pavement-foundation coupled system | |
CN106844982A (en) | Recognition methods based on the bridge floor multiaxis traveling load for solving maximum entropy | |
Mohanta et al. | A parametric analysis on the performance of vehicle tires | |
Brach | An impact moment coefficient for vehicle collision analysis | |
CN106595932A (en) | Truncated total least squares-based bridge floor multiple axle moving load identifying method | |
CN109145466A (en) | A kind of 1/4 car model modeling method based on McPherson suspension | |
CN106909776A (en) | The recognition methods of the bridge floor multiaxis traveling load based on principal element weighting iterative method | |
CN112364426A (en) | Railway pier damage assessment method and system based on driving safety and dynamic response, terminal device and readable storage medium | |
Giavotto et al. | Vehicle dynamics and crash dynamics with minicomputer | |
CN106874584A (en) | The recognition methods of the bridge floor multiaxis traveling load based on composition Average Iteration algorithm | |
CN113449378B (en) | Transverse load transfer rate calculation method based on deformation of vehicle suspension |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170630 |
|
RJ01 | Rejection of invention patent application after publication |