JP2018204952A - Displacement amount calculation system, program and recording medium - Google Patents

Abstract

To specify the displacement amount of an upper structure of a bridge by disposing an accelerometer in the upper structure of the bridge.SOLUTION: A receiving unit 121 of a displacement amount calculating device 12 continuously receives measurement values of an accelerometer attached to an upper structure of a bridge. A first calculation unit 1231 of a calculation unit 123 calculates displacement amounts by second-order integration of the measurement values of the accelerometer. From the displacement amount calculated by the first calculation unit 1231, the gravity component included in the measurement value of the accelerometer and a part of the integration error due to the measurement error are removed. A second calculation unit 1232 of the calculation unit 123 specifies a spline curve which has the control point as the bottom appearing before and after the timing at which the smoothed displacement amount obtained by smoothing the displacement amount calculated by the first calculation unit 1231 shows the peak. The second calculation unit 1232 subtracts the value indicated by the spline curve from the displacement amount calculated by the first calculation unit 1231, and removes the integration error remaining in the displacement amount calculated by the first calculation unit 1231.SELECTED DRAWING: Figure 3

Description

  The present invention relates to a technique for calculating a displacement amount of a bridge when passing through a vehicle.

  If the amount of displacement in the vertical direction when a vehicle (automobile, train, etc.) passes over the bridge (hereinafter referred to as “displacement amount” in this application means “displacement amount in the vertical direction”), the vehicle has passed. The relative weight of the vehicle can be estimated. In addition, if the weight of the vehicle that has passed over the bridge is known, and the change over time of the amount of displacement of the bridge when a vehicle having the same weight passes, the degree of deterioration of the bridge can be estimated.

  In this application, a bridge is a plurality of lower structures (including piers, supports, etc.) arranged in the direction of the bridge axis, and an upper structure (including main girders, floor slabs, etc.) supported by the plurality of lower structures. Means a structure comprising

  In order to directly measure the displacement amount of the bridge, it is necessary to arrange a rod, a piano wire or the like for applying a force corresponding to a change in the distance between the stationary reference point and the bridge to the displacement meter. However, for example, when the bridge is a bridge over water, it is often impossible to secure a stationary reference point. Even if a stationary reference point can be secured, it is not easy to place a rod, a piano wire or the like between the reference point and the bridge.

  There is a method of calculating the displacement amount by integrating the acceleration measured by the accelerometer installed in the superstructure of the bridge in the time domain or the frequency domain. According to this method, there is no need to arrange a rod, a piano wire or the like between the reference point and the bridge as in the case of directly measuring the amount of displacement. However, the measurement value of the accelerometer includes a measurement error (noise) that changes depending on the type of accelerometer, the identification of the solid, the installation direction, etc. The value obtained by simply integrating the measurement value of the accelerometer to the second floor is Does not represent the actual displacement correctly.

  In view of this, a technique has been proposed for removing an integration error caused by an acceleration measurement error that occurs when the displacement is calculated by second-order integration of the acceleration of the superstructure of the bridge.

  For example, Patent Document 1 discloses that wavelet transform is performed on a measurement value of an acceleration sensor attached to an object, and the wavelet coefficient obtained by the previous wavelet transform is 2 within a frequency range suitable for displacement analysis of the object. There has been proposed an apparatus that performs inverse wavelet transform on the result of the step integration and calculates the displacement of the object.

  Further, Patent Document 2 specifies a forced vibration section, a boundary condition of speed and a boundary condition of displacement based on measurement values of a plurality of acceleration sensors attached to an object, and a velocity waveform in a section including the forced vibration section. Is corrected to a corrected velocity waveform that satisfies the velocity boundary condition, and the displacement waveform in the section including the forced vibration section is corrected to the corrected displacement waveform that satisfies the displacement boundary condition.

Japanese Patent Laid-Open No. 2006-148549 JP 2017-003556 A

  In Patent Document 1 described above, it is necessary to appropriately set a frequency range for performing wavelet inverse transform for each object.

  An object of the present invention is to provide a technique for specifying the amount of displacement of a bridge superstructure by installing an accelerometer in the superstructure of the bridge.

  In order to solve the above-described problem, the present invention provides the predetermined value calculated using a measured value of acceleration continuously measured by the accelerometer at the predetermined position in a period including a period in which the vehicle passes through the predetermined position of the bridge. An acquisition unit that acquires displacement amount data indicating a temporal change in the displacement amount of the bridge at a position; a smoothing unit that smoothes the displacement amount data to generate smoothed displacement amount data; and the bridge loads a vehicle load. When the direction of movement received is the positive direction, the displacement amount at the start timing of the period in which the rate of change of the displacement amount indicated by the smoothed displacement amount data in a predetermined time length is maintained at a predetermined threshold value or more is the first amount. The end timing of the period that is specified as the control point of the period, and is a period subsequent to the period, in which the rate of change of the displacement amount indicated by the smoothed displacement amount data is maintained at a predetermined threshold value or less. A control point specifying unit for specifying a unit amount as a second control point, a spline curve having the first control point and the second control point as control points, and a part of the smoothed displacement amount data A displacement amount calculation system comprising: a baseline identifying unit that identifies, as baseline data, data replaced with data represented by the spline curve; and a correction unit that generates corrected displacement amount data by subtracting the baseline data from the displacement amount data. provide.

  In the displacement amount calculation system described above, the control point specifying unit specifies a timing at which the displacement amount indicated by the smoothed displacement amount data has a maximum value as a peak timing, and the control point specifying unit includes the smoothing displacement amount. The timing at which the displacement amount indicated by the data indicates the minimum value in the period before the peak timing is specified as the timing of the first control point, and the control point specifying unit determines that the displacement amount indicated by the smoothed displacement amount data is A timing indicating a minimum value in a period after the peak timing may be specified as the timing of the second control point.

  Further, in the displacement amount calculation system described above, the control point specifying unit specifies a period including a timing at which the displacement amount indicated by the smoothed displacement amount data indicates a maximum value as a peak period, and the control point specifying unit includes: The timing at which the displacement amount indicated by the smoothed displacement amount data indicates the minimum value in the period before the peak period is specified as the timing of the first control point, and the control point specifying unit is configured to output the smoothed displacement amount data. May be specified as the timing of the second control point. The timing at which the amount of displacement indicated by indicates a minimum value in a period after the peak period.

  In the displacement amount calculation system described above, the control point specifying unit uses the data obtained by smoothing data obtained by correcting the displacement amount data to determine the first control point and the second control point. You may specify.

  In the displacement amount calculation system described above, the baseline identifying unit may identify a primary spline curve as the spline curve having the first control point and the second control point as control points.

  In the displacement amount calculation system described above, the base line specifying unit may use the first control point and the second control point as the spline curve having the first control point and the second control point as control points. A spline curve that does not pass through at least one of the control points may be specified.

  Further, the present invention provides the computer with the predetermined position at the predetermined position calculated using the acceleration measurement value continuously measured by the accelerometer at the predetermined position during a period including a period in which the vehicle passes through the predetermined position of the bridge. A process of acquiring displacement data indicating a change in displacement of the bridge over time, a process of smoothing the displacement data to generate smoothed displacement data, and a direction in which the bridge moves in response to a vehicle load. In the case of the positive direction, the displacement amount at the start timing of the period in which the change rate in the predetermined time length of the displacement amount indicated by the smoothed displacement amount data is maintained at a predetermined threshold value or more is specified as the first control point. A displacement amount at the end timing of a period that is a period subsequent to the period and in which the rate of change of the displacement amount indicated by the smoothed displacement amount data in a predetermined time length is maintained below a predetermined threshold value is a second value. A process for specifying as a control point, a spline curve having the first control point and the second control point as control points are specified, and a part of the smoothed displacement data is replaced with data indicated by the spline curve A program for executing a process for identifying the data as baseline data and a process for generating corrected displacement data by subtracting the baseline data from the displacement data is provided.

  Further, the present invention provides the computer with the predetermined position at the predetermined position calculated using the acceleration measurement value continuously measured by the accelerometer at the predetermined position during a period including a period in which the vehicle passes through the predetermined position of the bridge. A process of acquiring displacement data indicating a change in displacement of the bridge over time, a process of smoothing the displacement data to generate smoothed displacement data, and a direction in which the bridge moves in response to a vehicle load. In the case of the positive direction, the displacement amount at the start timing of the period in which the change rate in the predetermined time length of the displacement amount indicated by the smoothed displacement amount data is maintained at a predetermined threshold value or more is specified as the first control point. A displacement amount at the end timing of a period that is a period subsequent to the period and in which the rate of change of the displacement amount indicated by the smoothed displacement amount data in a predetermined time length is maintained below a predetermined threshold value is a second value. A process for specifying as a control point, a spline curve having the first control point and the second control point as control points are specified, and a part of the smoothed displacement data is replaced with data indicated by the spline curve A computer-readable recording medium for continuously recording a program for executing the process of specifying the acquired data as baseline data and the process of generating the corrected displacement data by subtracting the baseline data from the displacement data provide.

  According to the present invention, the amount of displacement of the bridge superstructure is specified by installing the accelerometer in the superstructure of the bridge.

The figure which showed the displacement amount calculation system which concerns on one Embodiment. The figure which showed the basic composition of the computer used for implement | achieving the displacement amount calculation apparatus which concerns on one Embodiment. The figure which showed the structure of the displacement amount calculation apparatus which concerns on one Embodiment. The graph which shows the gravity correction acceleration calculated by the 1st smoothing part which concerns on one Embodiment. The graph which shows the 1st smoothing acceleration computed by the 1st smoothing part concerning one embodiment. The graph which shows the passage period which the passage period specific | specification part which concerns on one Embodiment specified. The graph which shows the analysis period which the analysis period specific | specification part which concerns on one Embodiment specified. The graph which shows the speed which the 1st integration part concerning one embodiment computed. The graph which shows the 1st base line which the 1st base line specific part concerning one embodiment specified. The graph which shows the correction speed which the 1st correction part concerning one embodiment generated. The graph which shows the displacement amount which the 2nd integration part concerning one embodiment computed. The graph which shows the 2nd base line which the 2nd base line specific part concerning one embodiment specified. The graph which shows the 1st correction | amendment displacement amount which the 2nd correction | amendment part which concerns on one Embodiment produced | generated. The graph which shows the 3rd base line which the 3rd base line specific part concerning one embodiment specified. The graph which shows the 2nd correction | amendment displacement amount which the 3rd correction | amendment part which concerns on one Embodiment produced | generated. The graph which shows the smoothing displacement amount which the 2nd smoothing part which concerns on one Embodiment computed. The graph which shows the 1st control point and the 2nd control point which the control point specific | specification part which concerns on one Embodiment specified. The graph which shows the 4th base line which the 4th base line specific part concerning one embodiment specified. The graph which shows the 3rd correction | amendment displacement amount which the 4th correction | amendment part which concerns on one Embodiment produced | generated. The graph which shows the 1st control point and the 2nd control point which the smoothing displacement amount which the 2nd smoothing part which concerns on one modification calculated, and the control point specific | specification part specified.

  A displacement amount calculation system 1 according to an embodiment of the present invention will be described below. FIG. 1 is a diagram showing a displacement amount calculation system 1. The displacement amount calculation system 1 receives the acceleration data indicating the acceleration value in the vertical direction that the accelerometer 11 continuously measures and the accelerometer 11 attached to the bridge 9 at a predetermined position, and the received acceleration data. And a displacement amount calculating device 12 for calculating a displacement amount in a vertical direction at a predetermined position of the bridge 9.

  The accelerometer 11 and the displacement amount calculation device 12 may directly transmit / receive acceleration data, or may transmit / receive acceleration data via a relay device.

  The bridge 9 includes a plurality of lower structures 91 and an upper structure 92 arranged so as to span the plurality of lower structures 91. Various types and weights of vehicles C pass over the superstructure 92 of the bridge 9.

  The predetermined position where the accelerometer 11 is attached to the upper structure 92 is a position substantially at the center of two adjacent lower structures 91 in the bridge axis direction. The position where the accelerometer 11 is attached may be any position as long as it is a position where displacement occurs as the vehicle passes. In addition to the basic function of measuring acceleration, the accelerometer 11 includes a transmission unit that transmits acceleration data indicating acceleration measured by wire or wireless to the displacement amount calculation device 12.

  The displacement amount calculation device 12 may be configured as a so-called dedicated device, or may be realized by a computer executing data processing according to a program. Hereinafter, a case where the displacement amount calculation device 12 is realized by a computer will be described as an example. However, each of the configurations of the displacement amount calculation device 12 described below may be realized by hardware.

  FIG. 2 is a diagram showing a basic configuration of the computer 10 used for realizing the displacement amount calculation device 12. The computer 10 includes a memory 101 that stores various data, a processor 102 that performs various data processing in accordance with a program stored in the memory 101, and a communication interface 103 that performs data communication with an external device.

  FIG. 3 is a diagram showing the configuration of the displacement amount calculation device 12. In the present embodiment, each component shown in FIG. 3 is realized by the computer 10 performing data processing according to a program.

  The receiving unit 121 receives acceleration data transmitted from the accelerometer 11. The receiving unit 121 is mainly realized by the communication interface 103. The storage unit 122 stores various data. The storage unit 122 is mainly realized by the memory 101. The data stored in the storage unit 122 includes acceleration data received by the reception unit 121 and various types of data indicating results calculated by the calculation unit 123 described below.

  The calculation unit 123 calculates the vertical displacement amount at the position where the accelerometer 11 of the bridge 9 is attached, using the acceleration data received from the accelerometer 11 by the reception unit 121. The arithmetic unit 123 is mainly realized by the processor 102.

  The calculation unit 123 is large and includes a first calculation unit 1231 and a second calculation unit 1232. The first calculation unit 1231 includes the following components.

A first acquisition unit 301 that acquires acceleration data indicating acceleration continuously measured by the accelerometer 11.
A first smoothing unit 302 that generates first smoothed acceleration data obtained by smoothing an acceleration measurement value indicated by acceleration data acquired by the first acquisition unit 301 (data indicating a temporal change in the acceleration measurement value).
A passing period specifying unit 303 that specifies a passing period that is a period in which the vehicle has passed in the acceleration data based on the first smoothed acceleration data.
An analysis period specifying unit 304 that specifies an analysis period that is a period to be analyzed in the acceleration data based on the passage period.

A first integration unit 305 that calculates the velocity by integrating the acceleration indicated by the acceleration data in the analysis period in the time domain.
A first baseline specifying unit 306 that specifies a linear approximate curve of the velocity calculated by the first integrating unit 305 as a first baseline.
A first correction unit 307 that generates a correction speed that is a speed correction value by subtracting the value indicated by the first baseline from the speed calculated by the first integration unit 305.

A second integration unit 308 that integrates the correction speed generated by the first correction unit 307 in the time domain and calculates a displacement amount.
A second baseline identifying unit 309 that identifies a linear approximate curve of the displacement calculated by the second integrating unit 308 as a second baseline.
A second correction unit 310 that subtracts the value indicated by the second baseline from the displacement amount calculated by the second integration unit 308 to generate a first correction displacement amount that is a first correction value of the displacement amount.

A third baseline specifying unit 311 that specifies a curve obtained by smoothing the first correction displacement amount generated by the second correction unit 310 as a third baseline.
By subtracting the value indicated by the third baseline from the first correction displacement amount generated by the second correction unit 310, second correction displacement amount data indicating the second correction displacement amount, which is the second correction value of the displacement amount, is generated. A third correction unit 312.

  Moreover, the 2nd calculating part 1232 is provided with the following structure parts.

A second acquisition unit 313 that acquires second correction displacement amount data generated by the third correction unit 312.
A second smoothing unit 314 that smoothes the second corrected displacement amount data to generate smoothed displacement amount data.

When the direction in which the bridge 9 moves under the load of the vehicle C is a positive direction, the timing at which the displacement amount indicated by the smoothed displacement amount data indicates the maximum value is specified as the peak timing, and the period before the peak timing A control point specifying unit 315 that specifies a point indicating the minimum value as the first control point and specifies a point indicating the minimum value in the period after the peak timing as the second control point.
A spline curve having the first control point and the second control point specified by the control point specifying unit 315 as control points is specified, and a part of the smoothed displacement data is replaced with data indicated by the specified spline curve. A fourth baseline identifying unit 316 that identifies the obtained data as baseline data.
A fourth correction unit 317 that generates third correction displacement amount data by subtracting the baseline data from the second correction displacement amount data;

  Note that the second acquisition unit 313 uses the acceleration measured continuously by the accelerometer at the predetermined position during a period including the period during which the vehicle passes the predetermined position of the bridge. It is an example of the acquisition part which acquires the displacement amount data which shows the time-dependent change of the displacement amount.

  The third corrected displacement data generated by the fourth correction unit 317 is stored in the storage unit 122, and the relative weight of the vehicle C that has passed through the bridge 9 is estimated and the degree of deterioration of the bridge 9 is estimated. Used for etc.

  Below, the specific example of the process which each structure part of the calculating part 123 mentioned above performs is demonstrated. The 1st acquisition part 301 acquires the acceleration data which show the measured value of the acceleration regarding the period before and the period before the vehicle C passes the position where the accelerometer 11 of the bridge 9 was attached.

The first smoothing unit 302 first performs a process of removing the gravity component included in the measured acceleration value indicated by the acceleration data. Specifically, the first smoothing unit 302 specifies the median value (or average value) of the measurement values of acceleration indicated by the acceleration data as a gravity component, subtracts the specified gravity component from the measurement value of acceleration, Calculate gravity correction acceleration. FIG. 4 is a graph showing the gravity correction acceleration A ₁ calculated by the first smoothing unit 302.

  4 to 19, the direction in which the bridge 9 moves under the load of the vehicle C is a positive direction. Therefore, when the direction in which the bridge 9 moves under the load of the vehicle C is a negative direction, in the processing described below, data indicated by a graph in which the positive and negative of the vertical axis in FIGS. It will be. In that case, the peak in the following description is replaced with the bottom, the bottom is replaced with the peak, the maximum value is replaced with the minimum value, and the minimum value is replaced with the maximum value. It should be noted that whether the direction in which the bridge moves under the load of the vehicle is the positive direction or the negative direction is merely a difference in expression and does not cause a technical difference.

The first smoothing unit 302 calculates a first smoothed acceleration obtained by smoothing a negative value of the gravity correction acceleration A ₁ (FIG. 4), and shows a first smoothed acceleration. Generate acceleration data. FIG. 5 is a graph showing the first smoothing acceleration A ₂ calculated by the first smoothing unit 302.

Specifically, the first smoothing unit 302 calculates the square root of each squared moving average value of the gravity correction acceleration A ₁ (FIG. 4) as the first smoothing acceleration A ₂ . In addition, the time length which the 1st smoothing part 302 uses for calculation of a moving average value is about 1-2 seconds, for example.

The passing period specifying unit 303 specifies the peak of the first smoothed acceleration A ₂ (FIG. 5) calculated by the first smoothing unit 302, and the vehicle C attaches the accelerometer 11 based on the specified peak value. The start timing and end timing of the period (passing period) passing through the specified position are specified. FIG. 6 is a graph showing the passing period T ₁ specified by the passing period specifying unit 303.

Specifically, the passage period specifying unit 303 specifies a point where the first smoothing acceleration A ₂ (FIG. 5) shows the maximum value as the peak P ₁ . Then, passing period specifying unit 303, in a period before the timing of a peak P _1, a predetermined ratio to the value U ₁ of the peak P ₁ (e.g., 10%) first smoothed acceleration to a threshold L ₁ multiplied by The point at which A ₂ decays is identified as bottom B ₁ . In addition, the passage period specifying unit 303 specifies a point where the first smoothing acceleration A ₂ is attenuated to the threshold L ₁ as the bottom B _{2 in} a period after the timing indicating the peak P ₁ . The passage period specifying unit 303 specifies a period having the timing of the bottom B ₁ as the start timing and the timing of the bottom B ₂ as the end timing as the passage period T ₁ .

The analysis period specifying unit 304 specifies an analysis period, which is a period of acceleration data used for the following analysis, based on the passing period T ₁ specified by the passing period specifying unit 303. FIG. 7 is a graph showing the analysis period T ₂ specified by the analysis period specifying unit 304.

Specifically, the analysis period specifying unit 304, before and after the passage period T _1, a predetermined ratio of the time length of the passing time T ₁ (e.g., time of 100% of the length of the passage time period) was added to the time period length identifying the period as the analysis period T _2. The process described below is performed for the analysis period T ₂ specified by the analysis period specifying unit 304.

The first integration unit 305 integrates the gravity correction acceleration A ₁ (FIG. 4), which is a measured value of the acceleration from which the gravity component has been removed by the first smoothing unit 302, in the time domain to calculate the velocity. FIG. 8 is a graph showing the speed V ₁ calculated by the first integration unit 305.

The first baseline identifying unit 306 identifies the linear approximate curve of the velocity V ₁ (FIG. 8) calculated by the first integrating unit 305 as the first baseline. FIG. 9 is a graph showing the first baseline E ₁ identified by the first baseline identifying unit 306. The method by which the first baseline identifying unit 306 identifies the linear approximation curve is, for example, the least square method.

The first correction unit 307 is indicated by the first base line E ₁ (FIG. 9) specified by the first base line specifying unit 306 from each timing value of the speed V ₁ (FIG. 8) calculated by the first integration unit 305. A correction speed is calculated by subtracting the value of each timing, and correction speed data indicating the correction speed is generated. FIG. 10 is a graph showing the correction speed V ₂ generated by the first correction unit 307.

The second integration unit 308 integrates the correction speed V ₂ (FIG. 10) generated by the first correction unit 307 in the time domain to calculate the displacement amount, and generates displacement amount data indicating the calculated displacement amount. FIG. 11 is a graph showing the displacement amount D ₁ calculated by the second integration unit 308.

The second baseline identifying unit 309 identifies the linear approximate curve of the displacement D ₁ (FIG. 11) generated by the second integrating unit 308 as the second baseline. FIG. 12 is a graph showing the second baseline E ₂ identified by the second baseline identifying unit 309. The method by which the second baseline identifying unit 309 identifies the linear approximation curve is, for example, the least square method.

The second correction unit 310 is indicated by the second baseline E ₂ (FIG. 12) identified by the second baseline identifying unit 309 from the timing value of the displacement D ₁ (FIG. 11) calculated by the second integration unit 308. The first correction displacement amount is calculated by subtracting the timing values to be generated, and first correction displacement amount data indicating the calculated first correction displacement amount is generated. FIG. 13 is a graph showing the first correction displacement amount D ₂ generated by the second correction unit 310.

The third baseline identifying unit 311 identifies the third baseline obtained by smoothing the first correction displacement amount D ₂ (FIG. 13) generated by the second correction unit 310. The third baseline mainly indicates a long-cycle integration error included in the displacement calculated by second-order integration of the measurement value of the accelerometer 11. FIG. 14 is a graph showing the third baseline E ₃ identified by the third baseline identifying unit 311.

Specifically, the third baseline specifying unit 311 calculates, for example, the moving average value of the first corrected displacement amount D ₂ as the value of each timing of the third baseline E ₃ . The time length used by the third baseline specifying unit 311 for calculating the moving average value is, for example, the same time length as the time length of the passage period T ₁ .

The third correction unit 312 determines the third baseline E ₃ (FIG. 14) identified by the third baseline identification unit 311 from each timing value of the first correction displacement amount D ₂ (FIG. 13) calculated by the second correction unit 310. The second corrected displacement amount is calculated by subtracting the value of each timing indicated by (), and second corrected displacement amount data indicating the calculated second corrected displacement amount is generated. FIG. 15 is a graph showing the second corrected displacement amount D ₃ generated by the third correction unit 312.

  The first calculator 1231 generates the second corrected displacement amount data as described above. The displacement amount indicated by the second corrected displacement amount data generated by the first calculation unit 1231 is obtained by removing the gravitational component from the measurement value of the accelerometer 11, removing the integration error after the first-order integration, This value is obtained by removing the integration error (long-cycle integration error) after integration. Therefore, the measured value of the accelerometer 11 is simply compared with the amount of displacement calculated by second-order integration, and a highly accurate value is shown.

  However, the second calculation unit 1232 included in the displacement amount calculation device 12 further removes the integration error remaining in the displacement amount indicated by the second corrected displacement amount data generated by the first calculation unit 1231, thereby achieving more accurate. Calculate high displacement.

  The second acquisition unit 313 of the second calculation unit 1232 acquires the second corrected displacement amount data (FIG. 15) from the first calculation unit 1231.

The second smoothing unit 314, the second correction displacement amount D ₃ indicated by the second correction displacement data to calculate a smoothed smooth displacement, to produce the calculated smoothed displacement data. FIG. 16 is a graph showing the smoothed displacement amount D ₄ calculated by the second smoothing unit 314.

Specifically, the second smoothing unit 314 calculates, for example, a moving average value of the second corrected displacement amount D ₃ (FIG. 15) as the smoothed displacement amount D ₄ . The time length used by the second smoothing unit 314 for calculating the moving average value is, for example, about 0.5 seconds.

The control point specifying unit 315 performs processing for specifying two control points used by the fourth baseline specifying unit 316 to specify a spline curve. FIG. 17 is a graph showing the first control point (bottom B ₃ ) and the second control point (bottom B ₄ ) specified by the control point specifying unit 315.

In order to specify the first control point and the second control point, the control point specifying unit 315 first superimposes the moving window having a predetermined time length on the smoothed displacement D ₄ while shifting the moving window within the moving window. A range (period) in which a value equal to or higher than the upper limit threshold calculated according to a predetermined rule or a value equal to or lower than the lower limit threshold calculated according to a predetermined rule is specified as the bottom appearance period T ₃ . The time length of the moving window used by the control point specifying unit 315 is, for example, about half the time length of the passage period T ₁ . Moreover, the upper limit threshold value used by the control point specifying unit 315 is, for example, a predetermined number of times (for example, + 1.5σ) the standard deviation σ of the smoothing displacement amount D ₄ . Moreover, the lower limit threshold value used by the control point specifying unit 315 is, for example, a predetermined number of times (for example, −1.5σ) the standard deviation σ of the smoothing displacement amount D ₄ .

In the bottom appearance period T ₃ specified as described above, the control point specifying unit 315 specifies a point where the smoothing displacement amount D ₄ has the maximum value as the peak P ₂ . Subsequently, the control point specifying unit 315 determines a point (bottom B ₃ ) in which the smoothed displacement amount D ₄ has the minimum value in the period before the peak P ₂ timing (peak timing) in the bottom appearance period T _3. Specify as the first control point. In addition, the control point specifying unit 315 sets a point (bottom B ₄ ) at which the smoothed displacement amount D ₄ shows the minimum value in the period after the timing of peak P ₂ (peak timing) in the bottom appearance period T ₃ . 2 is specified as a control point.

The fourth baseline identification unit 316 includes a period before the timing of the first control point (bottom B ₃ ) and a period after the timing of the second control point (bottom B ₄ ) of the smoothed displacement D ₄ . Based on the value, a spline curve having the first control point and the second control point as control points is specified, and a part of the smoothed displacement D ₄ (FIG. 16) (for example, from the timing of the first control point) A curve obtained by replacing the value of the period until the timing of the second control point) with the value indicated by the specified spline curve is specified as the fourth baseline, and baseline data indicating the value of the specified fourth baseline is generated.

The spline curve specified by the fourth baseline specifying unit 316 is, for example, a straight line (primary spline curve) that passes through the first control point (bottom B ₃ ) and the second control point (bottom B ₄ ). FIG. 18 is a graph showing the fourth base line E ₄ specified by the fourth base line specifying unit 316 using a primary spline curve.

  Note that the spline curve specified by the fourth baseline specifying unit 316 is not limited to the primary spline curve. For example, the nth order spline curve passing through the first control point and the second control point (n is an integer of 2 or more). Alternatively, it may be a B-spline curve that passes through the vicinity of the first control point and the second control point but does not pass at least one of them.

The fourth correction unit 317 uses the fourth base line E ₄ specified by the fourth base line specifying unit 316 from the timing values (values indicated by the second correction displacement amount data) of the second correction displacement amount D ₃ (FIG. 15). The third correction displacement amount is calculated by subtracting the timing values (values indicated by the baseline data) indicated by (FIG. 18), and third correction displacement amount data indicating the calculated third correction displacement amount is generated. FIG. 19 is a graph showing the third correction displacement amount D ₅ generated by the fourth correction unit 317.

  The second calculator 1232 generates the third corrected displacement amount data as described above. The displacement amount indicated by the third correction displacement amount data generated by the second calculation unit 1232 is obtained by removing the integration error remaining in the displacement amount indicated by the second correction displacement amount data generated by the first calculation unit 1231. Is. Accordingly, the third corrected displacement amount data indicates a displacement amount with higher accuracy than the second corrected displacement amount data.

  As described above, according to the displacement amount calculation system 1, the displacement amount of the upper structure 92 of the bridge 9 is specified by installing the accelerometer 11 in the upper structure 92 of the bridge 9.

[Modification]
The above-described embodiments can be variously modified within the scope of the technical idea of the present invention. Examples of these modifications are shown below. It should be noted that two or more of the following modified examples may be employed in combination.

(1) In the above-described embodiment, the displacement amount calculation device 12 is configured by one device. However, the present invention is not limited to this, and the displacement amount calculation device 12 is configured by two or more devices that operate in conjunction with each other. May be. For example, the displacement amount calculation device 12 includes two or more server devices that perform data communication via a network such as the Internet, and each of the server devices includes one or more of the components included in the calculation unit 123. May be adopted.

(2) In the embodiment described above, the first calculation unit 1231 calculates the amount of displacement by second-order integration of acceleration in the time domain, but the present invention is not limited to this, and the first calculation unit 1231 calculates the acceleration. The displacement may be calculated by second-order integration in the frequency domain.

(3) In the above-described embodiment, the configuration unit included in the first calculation unit 1231 is an example, and as long as the first calculation unit 1231 generates displacement amount data using acceleration data, Such a component may be provided. For example, the first calculation unit 1231 does not include the second baseline specifying unit 309, the second correction unit 310, the third baseline specifying unit 311, and the third correction unit 312, and the first calculation unit 1231 is added to the second calculation unit 1232. On the other hand, a configuration in which displacement amount data indicating the displacement amount D ₁ (FIG. 11) is delivered may be employed.

(4) In the above-described embodiment, the control point specifying unit 315 specifies the timing of the point (peak P ₂ ) at which the smoothing displacement amount D ₄ (FIG. 17) shows the maximum value as the peak timing, and before the peak timing. The point (bottom B ₃ ) where the smoothing displacement amount D ₄ shows the minimum value in the period is the first control point, and the point where the smoothing displacement amount D ₄ shows the minimum value (bottom B) in the period after the peak timing. ₄ ) was specified as the second control point. In the present invention, the method for specifying the first control point and the second control point is not limited to this.

  That is, when the direction in which the bridge 9 moves under the load of the vehicle C is a positive direction, the control point specifying unit 315 has a predetermined rate of change in the displacement amount indicated by the smoothed displacement amount data. The amount of change at the start timing of the period maintained above the threshold is specified as the first control point, and the rate of change in the predetermined time length of the amount of displacement indicated by the smoothed displacement amount data is a period following the period. The first control point and the second control point may be specified by any method as long as the amount of displacement at the end timing of the period during which is maintained below the predetermined threshold is specified as the second control point.

For example, the control point specifying unit 315, instead smoothing displacement D ₄ is to be specified as the peak timing to the timing of the point indicating the maximum value (peak P _2), the smoothed displacement D ₄ is a timing indicating the maximum value the period including identifying a peak period, the timing of the minimum value smoothed displacement D ₄ in the period before the peak period as the timing of the first control point, also, smooth displacement in the period after the peak period The timing at which the amount D ₄ shows the minimum value may be specified as the timing of the second control point. As a method for controlling point specifying unit 315 specifies the peak periods, for example, smoothing displacement D ₄ is the peak period indicating a value greater than a predetermined multiple of the standard deviation σ of the smoothed displacement D ₄ (e.g., 2 [sigma]) An example is a method of specifying the period.

Further, the control point specifying unit 315 specifies the first control point and the second control point without specifying the timing of the peak P ₂ (peak timing) or the period including the peak timing (peak period). Also good. An example of a method in which the control point specifying unit 315 specifies the first control point and the second control point without specifying the peak timing or peak period will be described below.

In the configuration exemplified in the modification (3), the first calculation unit 1231 does not include the second baseline specifying unit 309, the second correction unit 310, the third baseline specifying unit 311, and the third correction unit 312. The second acquisition unit 313 acquires displacement amount data indicating the displacement amount D ₁ (FIG. 11). In this case, a graph of the smoothing displacement amount D ₆ generated by the second smoothing unit 314 is shown in FIG.

First, the control point specifying unit 315 has a change rate in a predetermined time length (for example, a time length of about 0.5 to 1 second) of the smoothing displacement amount D _{6 with} a predetermined threshold (for example, about 0.1 mm / second). ) identify the time period T ₄ to be maintained at least to identify the F ₁ point indicating the displacement of the start timing of the period T ₄ as the first control point.

Subsequently, the control point specifying unit 315 is a period subsequent to the period T _4, and the rate of change of the smoothed displacement amount D _{6 in} a predetermined time length (for example, a time length of about 0.5 to 1 second) is changed. predetermined threshold value (e.g., -0.1 mm / degree s) to identify the time period T ₅ is maintained below identifies the F ₂ points indicating the amount of displacement of the end timing of the period T ₅ as a second control point.

(5) In the above-described embodiment, the control point specifying unit 315 uses the data obtained by smoothing the data obtained by correcting the displacement amount data acquired by the second acquisition unit 313 and uses the first control point and the second control point. Control points may be specified.

  As an example of a configuration that realizes this modified example, for example, the first calculation unit 1231 does not include the second baseline specifying unit 309, the second correction unit 310, the third baseline specifying unit 311, and the third correction unit 312. The point specifying unit 315 includes a second baseline specifying unit 309, a second correcting unit 310, a third baseline specifying unit 311 and a third correcting unit 312 and performs the same processing as the second smoothing unit 314. The structure provided with a conversion part is mentioned.

In this case, the second acquisition unit 313 acquires displacement amount data indicating the displacement amount D ₁ (FIG. 11). The second baseline identifying unit 309 included in the control point identifying unit 315 identifies the linear approximation curve of the displacement D ₁ (FIG. 11) acquired by the second acquiring unit 313 as the second baseline E ₂ (FIG. 12). The second correction unit 310 included in the control point specifying unit 315 _first subtracts each timing value indicated by the second baseline E ₂ (FIG. 12) from each timing value of the displacement D ₁ (FIG. 11). A corrected displacement amount D ₂ (FIG. 13) is calculated.

The third baseline identifying unit 311 included in the control point identifying unit 315 generates a third baseline E ₃ (FIG. 14) obtained by smoothing the first correction displacement amount D ₂ (FIG. 13). The third correction unit 312 included in the control point specifying unit 315 subtracts the value of each timing indicated by the third baseline E ₃ (FIG. 14) from the value of each timing of the first correction displacement amount D ₂ (FIG. 13). The second corrected displacement amount D ₃ (FIG. 15) is calculated.

The third smoothing unit included in the control point specifying unit 315 generates a smoothed displacement amount D ₄ (FIG. 16) obtained by smoothing the second corrected displacement amount D ₃ (FIG. 15). The control point specifying unit 315 specifies the peak P _{2 at} which the smoothing displacement amount D ₄ (FIG. 16) has the maximum value, and the bottom B at which the smoothing displacement amount D ₄ has the minimum value before and after the peak P _2. The timings ₃ and B ₄ are specified as the first timing and the second timing. Subsequently, the control point specifying unit 315 determines the smoothed displacement amount at the first timing from the smoothed displacement amount D ₆ (FIG. 20) generated by the second smoothing unit 314 as the first control point F _1. identified as to identify the smooth displacement of the second timing as a second control point F _2.

(6) In the above-described embodiment, the moving average is used as the smoothing method. However, low-pass filter processing other than the moving average may be employed instead of the moving average. In the above-described embodiment, the numerical values such as the time length specifically shown are examples, and other numerical values may be adopted.

(7) In the above-described embodiment, some processes may be omitted. For example, when the first base line identification unit 306 and the first correction unit 307 perform correction by the linear approximation curve of the first-order integral value on the final displacement amount is small, the first baseline identification unit 306 and the first baseline identification unit 306 The process by the 1 correction | amendment part 307 may be abbreviate | omitted.

(8) In the above-described embodiment, the order of some processes may be changed or integrated. For example, the order of the processes performed by the second baseline identifying unit 309 and the second correcting unit 310 and the processes performed by the third baseline identifying unit 311 and the third correcting unit 312 may be reversed. That is, after the displacement calculated by the second integration unit 308 is corrected by the third baseline identification unit 311 and the third correction unit 312, the correction is performed by the second baseline identification unit 309 and the second correction unit 310. May be performed. Also, a baseline identifying unit that identifies the baseline obtained by integrating the second baseline and the third baseline by integrating the second baseline identifying unit 309 and the third baseline identifying unit 311 may be employed. In this case, the second correction unit 310 and the third correction unit 312 are also integrated.

(9) In the above-described embodiment, the processes of the passage period specifying unit 303 and the analysis period specifying unit 304 may be replaced by other methods. For example, a sensor that detects passage of the vehicle C at a predetermined position may be attached to the bridge 9, and the passage period and the analysis period may be specified based on the detection result of the sensor.

(10) In the above-described embodiment, some processes may be simplified to reduce the calculation load. For example, the fourth base line specifying unit 316, when performing a particular spline curve, instead of using all of the smoothed displacement D _4, for example, those thinned portions of the sample of the smoothed displacement D ₄ (e.g. , Thinned to ¼) may be used.

(11) In the above-described embodiment, the program executed by the computer 10 to implement the displacement amount calculation device 12 may be downloaded to the computer 10 via a network such as the Internet, or may be a computer-readable record. The recording medium may be continuously recorded and distributed, and read by the computer 10 from the recording medium.

  DESCRIPTION OF SYMBOLS 1 ... Displacement amount calculation system, 9 ... Bridge, 10 ... Computer, 11 ... Accelerometer, 12 ... Displacement amount calculation apparatus, 91 ... Lower structure, 92 ... Upper structure, 101 ... Memory, 102 ... Processor, 103 ... Communication interface, DESCRIPTION OF SYMBOLS 121 ... Reception part, 122 ... Memory | storage part, 123 ... Operation part, 301 ... 1st acquisition part, 302 ... 1st smoothing part, 303 ... Passing period specific | specification part, 304 ... Analysis period specific part, 305 ... 1st integration part , 306 ... 1st baseline specifying part, 307 ... 1st correcting part, 308 ... 2nd integrating part, 309 ... 2nd baseline identifying part, 310 ... 2nd correcting part, 311 ... 3rd baseline identifying part, 312 ... 3rd Correction unit, 313 ... second acquisition unit, 314 ... second smoothing unit, 315 ... control point specifying unit, 316 ... fourth baseline specifying unit, 317 ... fourth correction unit, 1231 ... first calculation unit, 1232 ... first 2 operation part

Claims (8)

The time-dependent change of the displacement amount of the bridge at the predetermined position calculated by using the acceleration measurement value continuously measured by the accelerometer at the predetermined position in a period including a period in which the vehicle passes through the predetermined position of the bridge. An acquisition unit for acquiring displacement amount data;
A smoothing unit that smoothes the displacement amount data to generate smoothed displacement amount data;
Start of a period in which the rate of change of the displacement amount indicated by the smoothed displacement amount data is maintained at a predetermined threshold value or more when the direction in which the bridge moves under the load of the vehicle is a positive direction The displacement amount of the timing is specified as the first control point, and the rate of change of the displacement amount indicated by the smoothed displacement amount data in a predetermined time length is maintained below a predetermined threshold value in a period subsequent to the period. A control point specifying unit that specifies the amount of displacement at the end timing of the period as a second control point;
A spline curve having the first control point and the second control point as control points is specified, and data obtained by replacing a part of the smoothed displacement amount data with data indicated by the spline curve is specified as baseline data. A baseline identification part;
A displacement amount calculation system comprising: a correction unit that subtracts the baseline data from the displacement amount data to generate corrected displacement amount data. The control point specifying unit specifies a timing at which the displacement amount indicated by the smoothed displacement amount data indicates a maximum value as a peak timing,
The control point specifying unit specifies a timing at which the displacement amount indicated by the smoothed displacement amount data indicates a minimum value in a period before the peak timing as the timing of the first control point,
The displacement according to claim 1, wherein the control point specifying unit specifies a timing at which a displacement amount indicated by the smoothed displacement amount data indicates a minimum value in a period after the peak timing as a timing of the second control point. Quantity calculation system. The control point specifying unit specifies a period including a timing at which the displacement amount indicated by the smoothed displacement amount data indicates a maximum value as a peak period,
The control point specifying unit specifies a timing at which the displacement amount indicated by the smoothed displacement amount data indicates a minimum value in a period before the peak period as the timing of the first control point,
2. The displacement according to claim 1, wherein the control point specifying unit specifies a timing at which a displacement amount indicated by the smoothed displacement amount data indicates a minimum value in a period after the peak period as a timing of the second control point. Quantity calculation system. 2. The displacement amount according to claim 1, wherein the control point identifying unit identifies the first control point and the second control point using data obtained by smoothing data obtained by correcting the displacement amount data. Calculation system. The displacement according to any one of claims 1 to 4, wherein the base line specifying unit specifies a primary spline curve as the spline curve having the first control point and the second control point as control points. Quantity calculation system. The base line specifying unit is a spline curve that does not pass through at least one of the first control point and the second control point as the spline curve having the first control point and the second control point as control points. The displacement amount calculation system according to any one of claims 1 to 4, wherein: On the computer,
The time-dependent change of the displacement amount of the bridge at the predetermined position calculated by using the acceleration measurement value continuously measured by the accelerometer at the predetermined position in a period including a period in which the vehicle passes through the predetermined position of the bridge. Processing to obtain displacement data;
A process of smoothing the displacement amount data to generate smoothed displacement amount data;
Start of a period in which the rate of change of the displacement amount indicated by the smoothed displacement amount data is maintained at a predetermined threshold value or more when the direction in which the bridge moves under the load of the vehicle is a positive direction The displacement amount of the timing is specified as the first control point, and the rate of change of the displacement amount indicated by the smoothed displacement amount data in a predetermined time length is maintained below a predetermined threshold value in a period subsequent to the period. A process for identifying the amount of displacement at the end timing of the period as a second control point;
A spline curve having the first control point and the second control point as control points is specified, and data obtained by replacing a part of the smoothed displacement amount data with data indicated by the spline curve is specified as baseline data. Processing,
And a process for generating corrected displacement data by subtracting the baseline data from the displacement data. On the computer,
The time-dependent change of the displacement amount of the bridge at the predetermined position calculated by using the acceleration measurement value continuously measured by the accelerometer at the predetermined position in a period including a period in which the vehicle passes through the predetermined position of the bridge. Processing to obtain displacement data;
A process of smoothing the displacement amount data to generate smoothed displacement amount data;
Start of a period in which the rate of change of the displacement amount indicated by the smoothed displacement amount data is maintained at a predetermined threshold value or more when the direction in which the bridge moves under the load of the vehicle is a positive direction The displacement amount of the timing is specified as the first control point, and the rate of change of the displacement amount indicated by the smoothed displacement amount data in a predetermined time length is maintained below a predetermined threshold value in a period subsequent to the period. A process for identifying the amount of displacement at the end timing of the period as a second control point;
A spline curve having the first control point and the second control point as control points is specified, and data obtained by replacing a part of the smoothed displacement amount data with data indicated by the spline curve is specified as baseline data. Processing,
A computer-readable recording medium for continuously recording a program for executing a process of generating corrected displacement data by subtracting the baseline data from the displacement data.

Images