Measurement Method, Measurement Device, Measurement System, And Measurement Program

Abstract

A measurement method includes: generating first measurement data based on observation data of an observation point of a structure; generating second measurement data by performing filter processing on the first measurement data; calculating a first deflection amount of the structure; calculating a second deflection amount by performing filter processing on the first deflection amount; approximating the second measurement data with a linear function of the second deflection amount to calculate a first-order coefficient and a zero-order coefficient; calculating a third deflection amount based on the first-order coefficient, the zero-order coefficient, and the second deflection amount; calculating an offset based on the zero-order coefficient, the second deflection amount, and the third deflection amount; and calculating a static response by adding the offset and a product of the first-order coefficient and the first deflection amount.

Description

The present application is based on, and claims priority from JP Application Serial Number 2021-109015, filed Jun. 30, 2021, the disclosure of which is hereby incorporated by reference herein in its entirety.

BACKGROUND

1. Technical Field

The present disclosure relates to a measurement method, a measurement device, a measurement system, and a measurement program.

2. Related Art

JP-A-2018-31676 describes a railway bridge structural performance investigation method in which a theoretical analysis model of a dynamic response of a railway bridge during traveling of a train is formulated with the train as a moving load train and the bridge as a simple beam, an acceleration of the bridge during the traveling of the train is measured, an unknown parameter of the theoretical analysis model is identified by a vibration characteristic identification method by considering a residual acceleration waveform of the bridge after the train passes the bridge as a free vibration waveform, and the unknown parameter is estimated based on acceleration data by an inverse analysis method. More specifically, in the structural performance investigation method described in JP-A-2018-31676, a probability model is defined by introducing an error term into the theoretical analysis model, a simultaneous posterior probability density function of the unknown parameter when the acceleration data is given is defined by substituting, into an equation obtained by Bayes' theorem, a simultaneous occurrence probability at which the acceleration data is generated when the unknown parameter is given and a prior probability density function of the unknown parameter, and structural performance of the railway bridge is evaluated by reflecting the estimated parameter and uncertainty of the parameter.

Since an amount of data communication is enormous when the acceleration data acquired by an acceleration sensor installed on the bridge is transmitted to a host via a communication network, a method in which a measurement device installed near the acceleration sensor acquires the acceleration data, performs data processing, and transmits measurement data after the data processing to the host is preferable. With such a system configuration, it is possible to reduce the amount of data communication and reduce the costs of the entire system. However, in a method of estimating an unknown parameter of a theoretical analysis model based on acceleration data by an inverse analysis method as in the structural performance investigation method described in JP-A-2018-31676, a calculation amount is extremely large, and therefore, an expensive measurement device having high performance is required, and it is difficult to sufficiently reduce the costs of the entire system.

SUMMARY

A measurement method according to an aspect of the present disclosure includes: a first measurement data generation step of generating, based on observation data output from an observation device configured to observe an observation point of a structure, first measurement data based on a physical quantity which is a response to actions of a plurality of parts of a moving object moving on the structure on the observation point; a second measurement data generation step of generating second measurement data in which a vibration component is reduced by performing filter processing on the first measurement data; an observation information generation step of generating observation information including an entry time point and an exit time point of the moving object with respect to the structure; an average velocity calculation step of calculating an average velocity of the moving object based on the observation information and environment information which is created in advance and includes a dimension of the moving object and a dimension of the structure; a first deflection amount calculation step of calculating, based on an approximate equation of deflection of the structure, the observation information, the environment information, and the average velocity, a first deflection amount of the structure caused by the moving object; a second deflection amount calculation step of calculating a second deflection amount in which a vibration component is reduced by performing filter processing on the first deflection amount; a coefficient calculation step of approximating the second measurement data with a linear function of the second deflection amount to calculate a first-order coefficient and a zero-order coefficient of the linear function; a third deflection amount calculation step of calculating a third deflection amount based on the first-order coefficient, the zero-order coefficient, and the second deflection amount; an offset calculation step of calculating an offset based on the zero-order coefficient, the second deflection amount, and the third deflection amount; and a static response calculation step of calculating a static response by adding the offset and a product of the first-order coefficient and the first deflection amount.

A measurement device according to an aspect of the present disclosure includes: a first measurement data generation unit configured to generate, based on observation data output from an observation device configured to observe an observation point of a structure, first measurement data based on a physical quantity which is a response to actions of a plurality of parts of a moving object moving on the structure on the observation point; a second measurement data generation unit configured to generate second measurement data in which a vibration component is reduced by performing filter processing on the first measurement data; an observation information generation unit configured to generate observation information including an entry time point and an exit time point of the moving object with respect to the structure; an average velocity calculation unit configured to calculate an average velocity of the moving object based on the observation information and environment information which is created in advance and includes a dimension of the moving object and a dimension of the structure; a first deflection amount calculation unit configured to calculate, based on an approximate equation of deflection of the structure, the observation information, the environment information, and the average velocity, a first deflection amount of the structure caused by the moving object; a second deflection amount calculation unit configured to calculate a second deflection amount in which a vibration component is reduced by performing filter processing on the first deflection amount; a coefficient calculation unit configured to approximate the second measurement data with a linear function of the second deflection amount to calculate a first-order coefficient and a zero-order coefficient of the linear function; a third deflection amount calculation unit configured to calculate a third deflection amount based on the first-order coefficient, the zero-order coefficient, and the second deflection amount; an offset calculation unit configured to calculate an offset based on the zero-order coefficient, the second deflection amount, and the third deflection amount; and a static response calculation unit configured to calculate a static response by adding the offset and a product of the first-order coefficient and the first deflection amount.

A measurement system according to an aspect of the present disclosure includes: the measurement device according to the above aspect; and the observation device.

A non-transitory computer-readable storage medium according to an aspect of the present disclosure stores a measurement program, and the measurement program causes a computer to execute: a first measurement data generation step of generating, based on observation data output from an observation device configured to observe an observation point of a structure, first measurement data based on a physical quantity which is a response to actions of a plurality of parts of a moving object moving on the structure on the observation point; a second measurement data generation step of generating second measurement data in which a vibration component is reduced by performing filter processing on the first measurement data; an observation information generation step of generating observation information including an entry time point and an exit time point of the moving object with respect to the structure; an average velocity calculation step of calculating an average velocity of the moving object based on the observation information and environment information which is created in advance and includes a dimension of the moving object and a dimension of the structure; a first deflection amount calculation step of calculating, based on an approximate equation of deflection of the structure, the observation information, the environment information, and the average velocity, a first deflection amount of the structure caused by the moving object; a second deflection amount calculation step of calculating a second deflection amount in which a vibration component is reduced by performing filter processing on the first deflection amount; a coefficient calculation step of approximating the second measurement data with a linear function of the second deflection amount to calculate a first-order coefficient and a zero-order coefficient of the linear function; a third deflection amount calculation step of calculating a third deflection amount based on the first-order coefficient, the zero-order coefficient, and the second deflection amount; an offset calculation step of calculating an offset based on the zero-order coefficient, the second deflection amount, and the third deflection amount; and a static response calculation step of calculating a static response by adding the offset and a product of the first-order coefficient and the first deflection amount.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagram showing a configuration example of a measurement system.

FIG. 2 is a cross-sectional view of a superstructure of FIG. 1 taken along line A-A.

FIG. 3 is a diagram illustrating an acceleration detected by an acceleration sensor.

FIG. 4 is a diagram showing an example of measurement data u(t).

FIG. 5 is a diagram showing a power spectrum density of the measurement data u(t).

FIG. 6 is a diagram showing an example of measurement data u_lp(t).

FIG. 7 is a diagram showing an example of a relation between the measurement data u_lp(t), an entry time point t_i, and an exit time point t_o.

FIG. 8 is a diagram showing an example of a length L_c (C_m) of a vehicle and a distance La(a_w(C_m,n)) between axles.

FIG. 9 is a diagram illustrating a structural model of a superstructure of a bridge.

FIG. 10 is a diagram showing an example of a deflection amount w_std(a_w(C_m,n),t).

FIG. 11 is a diagram showing an example of a deflection amount C_std(C_m,t).

FIG. 12 is a diagram showing an example of a deflection amount T_std(t).

FIG. 13 is a diagram showing an example of a deflection amount T_{std_lp}(t).

FIG. 14 is a diagram showing the measurement data u_lp(t) and the deflection amount T_{std_lp}(t) in an overlapping manner.

FIG. 15 is a diagram showing an example of a deflection amount T_{Estd_lp}(t).

FIG. 16 is a diagram showing an example of a deflection amount T_Estd(t).

FIG. 17 is a diagram showing an example of a relation between the deflection amount T_{Estd_lp}(t) and the deflection amount T_{std_lp}(t) and a predetermined interval T_avg for calculating average values thereof.

FIG. 18 is a diagram showing an example of an offset T_{offset_std}(t).

FIG. 19 is a diagram showing an example of a deflection amount T_EOstd(t).

FIG. 20 is a diagram showing a relation between the measurement data u(t) and the deflection amount T_EOstd(t).

FIG. 21 is a diagram showing an example of the measurement data u(t) including a static response and a dynamic response.

FIG. 22 is a diagram showing an example of the measurement data u_lp(t).

FIG. 23 shows an example of a relation between the measurement data u_lp(t) and the entry time point t_i and the exit time point t_o.

FIG. 24 is a diagram showing an example of the deflection amount T_std(t).

FIG. 25 is a diagram showing the measurement data u_lp(t) and the deflection amount T_{std_lp}(t) in an overlapping manner.

FIG. 26 shows an example of the offset T_{offset_std}(t).

FIG. 27 is a diagram showing an example of the deflection amount T_EOstd(t).

FIG. 28 is a diagram showing a relation between the measurement data u(t) and the deflection amount T_EOstd(t).

FIG. 29 is a flowchart showing an example of a procedure of a measurement method according to a first embodiment.

FIG. 30 is a flowchart showing an example of a procedure of a first measurement data generation step.

FIG. 31 is a flowchart showing an example of a procedure of a second measurement data generation step.

FIG. 32 is a flowchart showing an example of a procedure of an observation information generation step.

FIG. 33 is a flowchart showing an example of a procedure of an average velocity calculation step.

FIG. 34 is a flowchart showing an example of a procedure of a first deflection amount calculation step.

FIG. 35 is a flowchart showing an example of a procedure of a second deflection amount calculation step.

FIG. 36 is a flowchart showing an example of a procedure of an offset calculation step.

FIG. 37 is a diagram showing a configuration example of a sensor, a measurement device, and a monitoring device.

FIG. 38 is a diagram showing an example of a load waveform P_std(a_w (C_m, n),t).

FIG. 39 is a diagram showing an example of a load waveform PT_std(t).

FIG. 40 is a diagram showing a relation between the deflection amount w_std(a_w(C_m,n),t) and a maximum amplitude max{w_std(a_w (C_m, n), t)}.

FIG. 41 is a flowchart showing an example of a procedure of a measurement method according to a second embodiment.

FIG. 42 is a flowchart showing an example of a procedure of a conversion ratio calculation step.

FIG. 43 is a flowchart showing an example of a procedure of a load waveform calculation step.

FIG. 44 is a diagram showing a configuration example of a measurement device according to the second embodiment.

FIG. 45 is a diagram showing another configuration example of the measurement system.

FIG. 46 is a diagram showing another configuration example of the measurement system.

FIG. 47 is a diagram showing another configuration example of the measurement system.

FIG. 48 is a cross-sectional view of a superstructure of FIG. 47 taken along line A-A.

DESCRIPTION OF EXEMPLARY EMBODIMENTS

Hereinafter, preferred embodiments of the present disclosure will be described in detail with reference to the drawings. The embodiments described below do not in any way limit contents of the present disclosure described in the claims. Not all configurations to be described below are necessarily essential components of the present disclosure.

1. First Embodiment

1-1. Configuration of Measurement System

A moving object passing through a superstructure of a bridge that is a structure according to a first embodiment is a vehicle, a railway vehicle, or the like that has a large weight and can be measured by BWIM. The BWIM is an abbreviation of bridge weigh in motion, and is a technology in which a bridge is regarded as a “scale”, deformation of the bridge is measured, and thereby a weight and the number of axles of the moving object passing through the bridge are measured. The superstructure of the bridge, which enables analysis of the weight of the moving object passing through the bridge, based on a response such as deformation or strain, is a structure in which the BWIM functions. The BWIM system, which applies a physical process between an action on the superstructure of the bridge and the response, enables the measurement of the weight of the moving object that travels on the bridge. Hereinafter, a measurement system for implementing a measurement method according to the present embodiment will be described by taking a case where the moving object is a railway vehicle as an example.

FIG. 1 is a diagram showing an example of the measurement system according to the present embodiment. As shown in FIG. 1, a measurement system 10 according to the present embodiment includes a measurement device 1, and at least one sensor 2 provided on a superstructure 7 of a bridge 5. The measurement system 10 may include a monitoring device 3.

The bridge 5 includes the superstructure 7 and a substructure 8. FIG. 2 is a cross-sectional view of the superstructure 7 taken along line A-A of FIG. 1. As shown in FIGS. 1 and 2, the superstructure 7 includes a bridge floor 7a, a support 7b, rails 7c, ties 7d, and a ballast 7e, and the bridge floor 7a includes a floor plate F, a main girder G, a cross girder (not shown), and the like. As shown in FIG. 1, the substructure 8 includes bridge piers 8a and bridge abutments 8b. The superstructure 7 is a structure across any one of the bridge abutment 8b and the bridge pier 8a adjacent to each other, two adjacent bridge abutments 8b, and two adjacent bridge piers 8a. Both end portions of the superstructure 7 are located at positions of the bridge abutment 8b and the bridge pier 8a adjacent to each other, at positions of the two adjacent bridge abutments 8b, or at positions of the two adjacent bridge piers 8a.

When a railway vehicle 6 enters the superstructure 7, the superstructure 7 is bent due to a load of the railway vehicle 6. Since the railway vehicle 6 includes a plurality of vehicles coupled to each other, a phenomenon occurs in which the bending of the superstructure 7 is periodically repeated as the vehicles pass through the superstructure 7. This phenomenon is called a static response. On the other hand, as a structure, the superstructure 7 has a natural vibration frequency, and therefore, natural vibration of the superstructure 7 may be excited when the railway vehicle 6 passes through the superstructure 7. When the natural vibration of the superstructure 7 is excited, a phenomenon occurs in which the bending of the superstructure 7 is periodically repeated. This phenomenon is called a dynamic response.

The measurement device 1 and the sensors 2 are coupled by, for example, a cable (not shown) and communicate with each other via a communication network such as a CAN. CAN is an abbreviation for controller area network. Alternatively, the measurement device 1 and the sensors 2 may communicate with each other via a wireless network.

Each sensor 2 outputs data used to calculate a static response when the railway vehicle 6, which is a moving object, moves on the superstructure 7, which is a structure. In the present embodiment, each sensor 2 is an acceleration sensor, and may be, for example, a crystal acceleration sensor or a MEMS acceleration sensor. MEMS is an abbreviation for micro electro mechanical systems.

In the present embodiment, each sensor 2 is installed at a central portion of the superstructure 7 in a longitudinal direction, specifically, at a central portion of the main girder G in the longitudinal direction. However, each sensor 2 is not limited to being installed at the central portion of the superstructure 7 as long as each sensor 2 can detect an acceleration for calculating the static response. When each sensor 2 is provided on the floor plate F of the superstructure 7, the sensor 2 may be damaged due to traveling of the railway vehicle 6, and the measurement accuracy may be affected by local deformation of the bridge floor 7a, so that in the example of FIGS. 1 and 2, each sensor 2 is provided at the main girder G of the superstructure 7.

The floor plate F, the main girder G, and the like of the superstructure 7 are bent in a vertical direction due to a load of the railway vehicle 6 passing through the superstructure 7. Each sensor 2 detects an acceleration of the bending of the floor plate F or the main girder G caused by the load of the railway vehicle 6 passing through the superstructure 7.

The measurement device 1 calculates the static response when the railway vehicle 6 passes through the superstructure 7 based on acceleration data output from each sensor 2. The measurement device 1 is installed on, for example, the bridge abutment 8b.

The measurement device 1 and the monitoring device can communicate with each other via, for example, a wireless network of a mobile phone and a communication network 4 such as the Internet. The measurement device 1 transmits measurement data including the static response when the railway vehicle 6 passes through the superstructure 7 to the monitoring device 3. The monitoring device 3 may store the information in a storage device (not shown), and may perform, for example, processing such as monitoring of the railway vehicle 6 and abnormality determination of the superstructure 7 based on the information.

In the present embodiment, the bridge 5 is a railway bridge, and is, for example, a steel bridge, a girder bridge, or an RC bridge. RC is an abbreviation for reinforced-concrete.

As shown in FIG. 2, in the present embodiment, an observation point R is set in association with the sensor 2. In the example of FIG. 2, the observation point R is set at a position on a surface of the superstructure 7 located vertically above the sensor 2 provided at the main girder G. That is, the sensor 2 is an observation device that observes the observation point R. The sensor 2 detects a physical quantity which is a response to actions of a plurality of parts of the railway vehicle 6 moving on the superstructure 7, which is a structure, on the observation point R, and outputs data including the detected physical quantity. For example, each of the plurality of parts of the railway vehicle 6 is an axle or a wheel, and is hereinafter assumed to be an axle. In the present embodiment, each sensor 2 is an acceleration sensor and detects an acceleration as the physical quantity. The sensor 2 may be provided at a position where the acceleration generated at the observation point R due to the traveling of the railway vehicle 6 can be detected, but the sensor 2 is preferably provided at a position close to the observation point R in the vertical direction.

The number and installation positions of the sensors 2 are not limited to the example shown in FIGS. 1 and 2, and various modifications can be made.

The measurement device 1 acquires, based on the acceleration data output from the sensor 2, an acceleration in a direction intersecting a surface of the superstructure 7 on which the railway vehicle 6 moves. The surface of the superstructure 7 on which the railway vehicle 6 moves is defined by a direction in which the railway vehicle 6 moves, that is, an X direction which is the longitudinal direction of the superstructure 7, and a direction orthogonal to the direction in which the railway vehicle 6 moves, that is, a Y direction which is a width direction of the superstructure 7. Since the observation point R is bent in a direction orthogonal to the X direction and the Y direction due to the traveling of the railway vehicle 6, the measurement device 1 preferably acquires an acceleration in the direction orthogonal to the X direction and the Y direction, that is, a Z direction which is a normal direction of the floor plate F, in order to accurately calculate a magnitude of the acceleration of the bending.

FIG. 3 is a diagram showing the acceleration detected by the sensor 2. The sensor 2 is an acceleration sensor that detects accelerations generated in three axes orthogonal to one another.

In order to detect the acceleration of the bending at the observation point R caused by the traveling of the railway vehicle 6, the sensor 2 is installed such that one of three detection axes thereof, which are the x axis, the y axis, and the z axis, is in a direction intersecting the X direction and the Y direction. In FIGS. 1 and 2, the sensor 2 is installed such that one axis thereof is in a direction intersecting the X direction and the Y direction. The observation point R is bent in the direction orthogonal to the X direction and the Y direction. Therefore, in order to accurately detect the acceleration of the bending, ideally, the sensor 2 is installed such that one axis thereof is aligned with the direction orthogonal to the X direction and the Y direction, that is, the normal direction of the floor plate F.

However, when the sensor 2 is installed on the superstructure 7, an installation location may be inclined. In the measurement device 1, even if one of the three detection axes of the sensor 2 is not installed in alignment with the normal direction of the floor plate F, since the one axis is substantially oriented in the normal direction, an error is small and thus can be ignored. The measurement device 1 can correct a detection error due to inclination of the sensor 2 based on a three-axis combined acceleration obtained by combining the accelerations in the x axis, the y axis, and the z axis even if one of the three detection axes of the sensor 2 is not installed in alignment with the normal direction of the floor plate F. Further, the sensor 2 may be a one-axis acceleration sensor that detects an acceleration generated at least in a direction substantially parallel to the vertical direction or an acceleration in the normal direction of the floor plate F.

Hereinafter, details of the measurement method according to the present embodiment executed by the measurement device 1 will be described.

1-2. Details of Measurement Method

First, the measurement device 1 integrates acceleration data a(k) output from the sensor 2, which is an acceleration sensor, to generate velocity data v(k) as in Equation (1), and further integrates the velocity data v(k) to generate measurement data u(k) as in Equation (2). The acceleration data a(k) is data of an acceleration change excluding an acceleration bias unnecessary for calculating a displacement change when the railway vehicle 6 passes through the bridge 5. For example, the acceleration directly before the railway vehicle 6 passes through the bridge 5 may be set to 0, and the subsequent acceleration change may be set as the acceleration data a(k). In Equation (1) and Equation (2), k is a sample number, and ΔT is a time interval of samples. The measurement data u(k) is data of the displacement of the observation point R due to the traveling of the railway vehicle 6.

v(k)=a(k)ΔT+v(k−1)  (1)

u(k)=v(k)ΔT+u(k−1)  (2)

The measurement data u(k) having the sample number k as a variable is converted into measurement data u(t) having the time point t as a variable at the time point t=kΔT. FIG. 4 shows an example of the measurement data u(t). Since the measurement data u(t) is generated based on the acceleration data a(t) output from the sensor 2 that observes the observation point R, the measurement data u(t) is data based on the acceleration that is a response to the actions of a plurality of parts of the railway vehicle 6 moving on the superstructure 7 on the observation point R.

Next, the measurement device 1 generates measurement data u_lp(t) obtained by performing filter processing on the measurement data u(t) in order to reduce a vibration component having a fundamental frequency F_f included in the measurement data u(t) and a harmonic of the vibration component. The filter processing may be, for example, low-pass filter processing or band-pass filter processing.

Specifically, first, the measurement device 1 calculates a power spectrum density by performing fast Fourier transform processing on the measurement data u(t), and calculates a peak of the power spectrum density as the fundamental frequency Ff. FIG. 5 shows the power spectrum density obtained by performing fast Fourier transform processing on the measurement data u(t) of FIG. 4. In the example of FIG. 5, the fundamental frequency F_f is calculated as about 3 Hz. Then, the measurement device 1 calculates a basic cycle T_f based on the fundamental frequency F_f according to Equation (3), and calculates a moving average interval k_mf adjusted to a time resolution of the data by dividing the basic cycle T_f by ΔT as in Equation (4). The basic cycle T_f is a cycle corresponding to the fundamental frequency F_f, and T_f>2ΔT.

$\begin{array}{cc}{T}_f=\frac{1}{F_f}& \left(3\right)\end{array}$ $\begin{array}{cc}{k}_{\mathrm{mf}}=2\lfloor \frac{T_f}{2\Delta T}\rfloor +1& \left(4\right)\end{array}$

Then, the measurement device 1 performs, as the filter processing, moving average processing on the measurement data u(t) in the basic cycle T_f according to Equation (5) to generate the measurement data u_lp(t) in which the vibration component included in the measurement data u(t) is reduced. In the moving average processing, not only the necessary calculation amount is small, but also an attenuation amount of a signal component of the fundamental frequency F_f and a harmonic component of the signal component is very large, so that the measurement data u_lp(t) in which the vibration component is effectively reduced is obtained. FIG. 6 shows an example of the measurement data u_lp(t). As shown in FIG. 6, the measurement data u_lp(t) from which the vibration component included in the measurement data u(t) is almost removed is obtained.

$\begin{array}{cc}{u}_{\mathrm{lp}}\left(k\right)=\frac{1}{k_{\mathrm{mf}}}\underset{n=k-\frac{k_{\mathrm{mf}}-1}{2}}{\sum ^{k+\frac{k_{\mathrm{mf}}-1}{2}}}u\left(n\right)& \left(5\right)\end{array}$

The measurement device 1 may generate the measurement data u_lp(t) by performing, as the filter processing, FIR filter processing for attenuating a signal component having a frequency equal to or higher than the fundamental frequency F_f on the measurement data u(t). FIR is an abbreviation of finite impulse response. In the FIR filter processing, although a calculation amount is larger than that of the moving average processing, all signal components having a frequency equal to or higher than the fundamental frequency F_f can be attenuated.

Next, the measurement device 1 calculates two time points at which an amplitude of the measurement data u_lp(t) matches a threshold C_Lu_a which is a product of a predetermined coefficient C_L and an amplitude u_a calculated based on the measurement data u_lp(t), or two times at which the amplitude of the measurement data u_lp(t) exceeds the threshold C_Lu_a, as an entry time point t₁ and an exit time point t_o of the railway vehicle 6 with respect to the superstructure 7. However, 0<C_L<1, and the amplitude u_a is calculated as, for example, an average value in an interval from a time point t₁ to a time point t₂ in which the amplitude of the measurement data u_lp(t) is shifted, according to Equation (6).

$\begin{array}{cc}{u}_a=\frac{1}{t_2-t_1}\sum _{t=t_1}^{t_2}{u}_{\mathrm{lp}}\left(t\right)& \left(6\right)\end{array}$

The entry time point t_i is a time point at which a leading axle of the plurality of axles of the railway vehicle 6 passes through an entry end of the superstructure 7. The exit time point t_o is a time point at which a rearmost axle of the plurality of axles of the railway vehicle 6 passes through an exit end of the superstructure 7. FIG. 7 shows an example of a relation between the measurement data u_lp(t), the entry time point t_i, and the exit time point t_o.

Next, the measurement device 1 calculates a difference between the exit time point t_o and the entry time point t_i, as a passing time t_s during which the railway vehicle 6 passes through the superstructure 7 of the bridge 5 according to Equation (7).

t _s =t ₀ −t _i  (7)

The measurement device 1 calculates, as the number of vehicles C_T of the railway vehicle 6, a maximum integer less than or equal to a number obtained by subtracting 1 from a product of the passing time t_s and the fundamental frequency F_f, according to Equation (8).

C _T =└t _s F _f−1┘=floor(t_sF_f−1)={t_sF_f−1}  (8)

The measurement device 1 stores observation information including the entry time point t_i, the exit time point t_o, the passing time t_s, and the number of vehicles CT in a storage unit (not shown). In the example of FIG. 7, the entry time point t_i is at 7.155 seconds, the exit time point t_o is at 12.845 seconds, the passing time t_s is 5.69 seconds, and the number of vehicles C_T is 16.

Then, the measurement device 1 performs the following processing based on the observation information and environment information which is created in advance and includes dimensions of the railway vehicle 6 and dimensions of the superstructure 7.

The environment information includes, for example, a length L_B of the superstructure 7 and a position L_x of the observation point R as the dimensions of the superstructure 7. The length L_B of the superstructure 7 is a distance between the entry end and the exit end of the superstructure 7. The position L_x of the observation point R is a distance from the entry end of the superstructure 7 to the observation point R. The environment information includes, for example, a length L_c(C_m) of each vehicle of the railway vehicle 6, the number of axles a_T(C_m) of each vehicle, and a distance La(a_w(C_m,n)) between axles of each vehicle, as the dimensions of the railway vehicle 6. C_m is a vehicle number, and the length L_c(C_m) of each vehicle is a distance between two ends of a C_m-th vehicle. The number of axles a_T(C_m) of each vehicle is the number of axles of the C_m-th vehicle, n is an axle number of each vehicle, and 1≤n≤a_T(C_m). The distance La (a_w (C_m, n)) between axles of each vehicle is a distance between a front end of the C_m-th vehicle and a first axle when n=1, and is a distance between the (n−1)-th axle and the n-th axle when n 2. FIG. 8 shows an example of the length L_c(Cm) of the C_m-th vehicle of the railway vehicle 6 and the distance La(a_w(C_m,n)) between the axles. The dimensions of the railway vehicle 6 and the dimensions of the superstructure 7 can be measured by a known method. A database of the dimensions of the railway vehicle 6 passing through the bridge 5 may be created in advance, and the dimensions of a corresponding vehicle may be referred to according to the passing time point.

When it is assumed that the railway vehicle 6 in which any number of vehicles having the same dimensions are coupled to each other travels on the superstructure 7 of the bridge 5, the environment information may include the length L_C(C_m) of the vehicle, the number of axles a_T(C_m) of the vehicle, and the distance La(a_w(C_m,n)) between the axles, which are related to one vehicle.

A total number of axles Ta_T of the railway vehicle 6 is calculated according to Equation (9) using the number of vehicles C_T included in the observation information and the number of axles a_T(C_m) of each vehicle included in the environment information.

$\begin{array}{cc}{\mathrm{Ta}}_T=\sum _{C_m=1}^{C_T}{a}_T\left(C_m\right)& \left(9\right)\end{array}$

A distance D_wa(a_w(C_m,n)) from the leading axle to the n-th axle of the C_m-th vehicle of the railway vehicle 6 is calculated according to Equation (10) using the length L_C(C_m) of each vehicle, the number of axles a_T(C_m) of each vehicle, and the distance La(a_w(C_m,n)) between axles of each vehicle included in the environment information. In Equation (10), it is assumed that L_C(C_m)=L_C(1).

$\begin{array}{cc}{D}_{\mathrm{wa}}\left(a_w\left(C_m,n\right)\right)=\sum _{y=1}^{C_m}{L}_C\left(y\right)+\sum _{x=1}^n\mathrm{La}\left(a_w\left(C_m,x\right)\right)-\left\{L_C\left(1\right)+\mathrm{La}\left(a_w\left(1,1\right)\right)\right\}& \left(10\right)\end{array}$

The measurement device 1 calculates a distance D_wa(a_w(C_T,a_T(C_T))) from the leading axle to the rearmost axle of the rearmost vehicle of the railway vehicle 6 according to Equation (11) obtained by substituting C_m=C_T and n=a_T(C_T) into Equation (10).

$\begin{array}{cc}{D}_{\mathrm{wa}}\left(a_w\left(C_T,a_T\left(C_T\right)\right)\right)=\sum _{y=1}^{C_T}{L}_C\left(y\right)+\sum _{x=1}^{a_T\left(C_T\right)}\mathrm{La}\left(a_w\left(C_T,x\right)\right)-\left\{L_C\left(1\right)+\mathrm{La}\left(a_w\left(1,1\right)\right)\right\}& \left(11\right)\end{array}$

An average velocity v_a of the railway vehicle 6 is calculated according to Equation (12) using the length L_B of the superstructure 7 included in the environment information, the passing time t_s included in the observation information, and the calculated distance D_wa(a_w(C_T,a_T(C_T))).

$\begin{array}{cc}{v}_a=\frac{L_B}{t_s}+\frac{D_{\mathrm{wa}}\left(a_w\left(C_T,a_T\left(C_T\right)\right)\right)}{t_s}& \left(12\right)\end{array}$

The measurement device 1 calculates the average velocity v_a of the railway vehicle 6 according to Equation (13) obtained by substituting Equation (11) into Equation (12).

$\begin{array}{cc}{v}_a=\frac{L_B}{t_s}+\frac{1}{t_s}\left[\sum _{y=1}^{C_T}{L}_C\left(y\right)+\sum _{x=1}^{a_T\left(C_T\right)}\mathrm{La}\left(a_w\left(C_T,x\right)\right)-\left\{L_C\left(1\right)+\mathrm{La}\left(a_w\left(1,1\right)\right)\right\}\right]& \left(13\right)\end{array}$

Next, the measurement device 1 calculates a deflection amount of the superstructure 7 caused by the traveling of the railway vehicle 6 in the following manner.

In the present embodiment, considering that the superstructure 7 of the bridge 5 has a configuration in which one or a plurality of bridge floors 7a including the floor plate F, the main girder G, and the like are continuously arranged, the measurement device 1 calculates a displacement of one bridge floor 7a as a displacement at the central portion in the longitudinal direction. The load applied to the superstructure 7 moves from one end to the other end of the superstructure 7. In this case, the deflection amount, which is the displacement of the central portion of the superstructure 7, can be expressed by a position of the load on the superstructure 7 and an amount of the load. In the present embodiment, in order to express the deflection deformation when the axles of the railway vehicle 6 moves on the superstructure 7 as a trajectory of the deflection amount caused by the movement on the bridge under one-point load, a structural model shown in FIG. 9 is considered, and the deflection amount at the central portion is calculated in the structural model. In FIG. 9, P is a load, a is a load position from the entry end of the superstructure 7 on a side where the railway vehicle 6 enters, and b is a load position from the exit end of the superstructure 7 on a side where the railway vehicle 6 exits. L_B is the length of the superstructure 7, that is, the distance between two ends of the superstructure 7. The structural model shown in FIG. 9 is a simple beam in which two ends are supported with the two ends as fulcrums.

In the structural model shown in FIG. 9, when the position of the entry end of the superstructure 7 is zero and an observation position of the deflection amount is x, a bending moment M of the simple beam is expressed by Equation (14).

$\begin{array}{cc}M=\frac{b}{L_B}\mathrm{Px}-{\mathrm{PH}}_a\left(x-a\right)& \left(14\right)\end{array}$

In Equation (14), a function H_a is defined as Equation (15).

$\begin{array}{cc}{H}_a=\{\begin{array}{c}0\left(\mathrm{if}x\le a\right)\\ 1\left(\mathrm{if}x>a\right)\end{array}& \left(15\right)\end{array}$

Equation (16) is obtained by transforming Equation (14).

$\begin{array}{cc}-\frac{{\mathrm{ML}}_B}{P}=-\mathrm{bx}+H_a{L}_B\left(x-a\right)& \left(16\right)\end{array}$

Meanwhile, the bending moment M is expressed by Equation (17). In Equation (17), 0 is an angle, I is a secondary moment, and E is a Young's modulus.

$\begin{array}{cc}-M=\mathrm{EI}\frac{d\theta }{\mathrm{dx}}& \left(17\right)\end{array}$

Equation (18) is obtained by substituting Equation (17) into Equation (16).

$\begin{array}{cc}\frac{{\mathrm{EIL}}_B}{P}\frac{d\theta }{\mathrm{dx}}=-\mathrm{bx}+H_a{L}_B\left(x-a\right)& \left(18\right)\end{array}$

Equation (20) is obtained by calculating Equation (19) which is obtained by integrating Equation (18) with respect to the observation position x. In Equation (20), C₁ is an integral constant.

$\begin{array}{cc}\int \frac{{\mathrm{EIL}}_B}{P}\frac{d\theta }{\mathrm{dx}}\mathrm{dx}=\int \left(-\mathrm{bx}+H_a{L}_B\left(x-a\right)\right)\mathrm{dx}& \left(19\right)\end{array}$ $\begin{array}{cc}\frac{{\mathrm{EIL}}_B}{P}\theta =-\frac{{\mathrm{bx}}^2}{2}+H_a\frac{{L_B\left(x-a\right)}^2}{2}+C_1& \left(20\right)\end{array}$

Equation (22) is obtained by calculating Equation (21) which is obtained by integrating Equation (20) with respect to the observation position x. In Equation (22), C₂ is an integral constant.

$\begin{array}{cc}\int \frac{{\mathrm{EIL}}_B}{P}\theta \mathrm{dx}=\int \left\{-\frac{{\mathrm{bx}}^2}{2}+H_a\frac{{L_B\left(x-a\right)}^2}{2}+C_1\right\}\mathrm{dx}& \left(21\right)\end{array}$ $\begin{array}{cc}\frac{{\mathrm{EIL}}_B}{P}\theta =-\frac{{\mathrm{bx}}^2}{2}+H_a\frac{{L_B\left(x-a\right)}^2}{2}+C_1& \left(22\right)\end{array}$

In Equation (22), Ox represents the deflection amount, and Equation (23) is obtained by replacing Ox with a deflection amount w.

$\begin{array}{cc}\frac{{\mathrm{EIL}}_B}{P}w=-\frac{{\mathrm{bx}}^3}{6}+H_a\frac{{L_B\left(x-a\right)}^3}{6}+C_{1}x+C_2& \left(23\right)\end{array}$

As shown in FIG. 9, since b=L_B−a, Equation (23) is transformed as in Equation (24).

$\begin{array}{cc}\frac{{\mathrm{EIL}}_B}{P}w=-\frac{\left(L_B-a\right)x^3}{6}+H_a\frac{\left(L_B-a\right)x^3}{6}+C_{1}x+C_2& \left(24\right)\end{array}$

When x=0 and the deflection amount w=0, H_a=0 as x≤a, and therefore, when x=w=H_a=0 is substituted into Equation (24), Equation (25) is obtained.

C ₂=0  (25)

When x=L_B and the deflection amount w=0, H_a=1 as x>a, and therefore, when x=L_B, w=0, and H_a=1 are substituted into Equation (24), Equation (26) is obtained.

$\begin{array}{cc}{C}_1=\frac{a\left(L_B-a\right)\left(a+2\left(L_B-a\right)\right)}{6}& \left(26\right)\end{array}$

Equation (27) is obtained by substituting b=L_B−a into Equation (26).

$\begin{array}{cc}{C}_1=\frac{\mathrm{ab}\left(a+2b\right)}{6}& \left(27\right)\end{array}$

Equation (28) is obtained by substituting an integral constant C₁ of Equation (25) and an integral constant C₂ of Equation (26) into Equation (23).

$\begin{array}{cc}\frac{{\mathrm{EIL}}_B}{P}w=-\frac{{\mathrm{bx}}^3}{6}+H_a\frac{{L_B\left(x-a\right)}^3}{6}+\frac{\mathrm{ab}\left(a+2b\right)}{6}x& \left(28\right)\end{array}$

Equation (28) is transformed, and the deflection amount w at the observation position x when the load P is applied to the position a is expressed by Equation (29).

$\begin{array}{cc}w=\frac{P}{6{\mathrm{EIL}}_R}\left\{-{\mathrm{bx}}^3+H_a{L_B\left(x-a\right)}^3+\mathrm{ab}\left(a+2b\right)x\right\}& \left(29\right)\end{array}$

A deflection amount w_0.5LB at the observation position x at the center when the load P is at the center of the superstructure 7 is expressed by Equation (30), wherein x=0.5LB, a=b=0.5LB, and H_a=0. The deflection amount w_0.5LB is a maximum amplitude of the deflection amount w.

$\begin{array}{cc}{w}_{0.5L_B}=\frac{P}{48\mathrm{EI}}{L}_B^3& \left(30\right)\end{array}$

The deflection amount w at any observation position x is normalized by the deflection amount w_0.5LB. When the position a of the load P is on an entry end side of the observation position x, as x>a, Equation (31) is obtained by substituting H_a=1 is into Equation (30).

$\begin{array}{cc}w=\frac{P}{6{\mathrm{EIL}}_B}\left\{-{\mathrm{bx}}^3+{L_B\left(x-a\right)}^3+\mathrm{ab}\left(a+2b\right)x\right\}& \left(31\right)\end{array}$

When the position a of the load P is represented by a=L_Br, and a=L_Br, b=L_B(1−r) is substituted into Equation (31), Equation (32) is obtained, and a deflection amount w_std in which the deflection amount w is normalized is obtained according to Equation (32). r represents a ratio of the position a of the load P to the length L_B of the superstructure 7.

$\begin{array}{cc}{w}_{\mathrm{std}}=\frac{8}{L_B}\left\{{\mathrm{xr}}^3+\left(\frac{x^3}{L_B^2}+2x\right)r\right\}-\frac{8}{L_B}\left(L_B{r}^3+\frac{3x^2}{L_B}r\right)& \left(32\right)\end{array}$

Similarly, when the position a of the load P is on an exit end side of the observation position x, as x≤a, Equation (33) is obtained by substituting H_a=0 into Equation (30).

$\begin{array}{cc}w=\frac{P}{6{\mathrm{EIL}}_B}\left\{-{\mathrm{bx}}^3+\mathrm{ab}\left(L_B+b\right)x\right\}& \left(33\right)\end{array}$

When the position a of the load P is represented by a=L_Br, and a=L_Br and b=L_B(1−r) are substituted into Equation (33), Equation (34) is obtained, and the deflection amount w_std in which the deflection amount w is normalized is obtained according to Equation (34).

$\begin{array}{cc}{w}_{\mathrm{std}}=\frac{8}{L_B}\left\{{\mathrm{xr}}^3+\left(\frac{x^3}{L_B^2}+2x\right)r\right\}-\frac{8}{L_B}\left(3{xr}^3+\frac{x^3}{L_B^2}\right)& \left(34\right)\end{array}$

Equation (32) and Equation (34) are combined, and a deflection amount w_std(r) at any observation position x=L_x is expressed by Equation (35). In Equation (35), a function R(r) is expressed by Equation (36). Equation (35) is an approximate equation of the deflection of the superstructure 7 which is a structure, and is an equation based on the structural model of the superstructure 7. Specifically, Equation (35) is an approximate equation normalized by the maximum amplitude of the deflection at the center position between the entry end and the exit end of the superstructure 7.

$\begin{array}{cc}{w}_{\mathrm{std}}\left(r\right)=\frac{8}{L_B}\left\{L_x{r}^3+\left(\frac{L_x^3}{L_B^2}+2L_x\right)r-R\left(r\right)\right\}& \left(35\right)\end{array}$ $\begin{array}{cc}R\left(r\right)=\{\begin{array}{cc}{L}_B{r}^3+\frac{3L_x^2}{L_B}r& \left(\mathrm{if}{L}_x>L_{B}r\right)\\ 3L_x{r}^2+\frac{L_x^3}{L_B^2}& \left(\mathrm{if}{L}_x\le L_{B}r\right)\end{array}& \left(36\right)\end{array}$

In the present embodiment, the load P is a load of any axle of the railway vehicle 6. A time t_xn required for a certain axle of the railway vehicle 6 to reach the position L_x of the observation point R from the entry end of the superstructure 7 is calculated according to Equation (37) using the average velocity v_a calculated according to Equation (12).

$\begin{array}{cc}{t}_{\mathrm{xn}}=\frac{L_x}{v_a}& \left(37\right)\end{array}$

A time tin required for a certain axle of the railway vehicle 6 to pass through the superstructure 7 having the length L_B is calculated according to Equation (38).

$\begin{array}{cc}{t}_{\ln }=\frac{L_B}{v_a}& \left(38\right)\end{array}$

A time point t₀(C_m,n) at which the n-th axle of the C_m-th vehicle of the railway vehicle 6 reaches the entry end of the superstructure 7 is calculated according to Equation (39) using the entry time point t_i included in the observation information, the distance D_wa(a_w(C_m,n)) calculated according to Equation (10), and the average velocity v_a calculated according to Equation (12).

$\begin{array}{cc}{t}_0\left(C_m,n\right)=t_i+\frac{1}{v_a}{D}_{\mathrm{wa}}\left(a_w\left(C_m,n\right)\right)& \left(39\right)\end{array}$

Using Equation (37), Equation (38), and Equation (39), the measurement device 1 calculates, according to Equation (40), a deflection amount w_std(a_w(C_m,n),t) obtained by replacing, with time, the deflection amount w_std(r) caused by the n-th axle of the C_m-th vehicle and represented by Equation (35). In Equation (40), a function R(t) is expressed by Equation (41). FIG. 10 shows an example of the deflection amount w_std(a_w(C_m,n),t).

$\begin{array}{cc}{w}_{\mathrm{std}}\left(a_w\left(C_m,n\right),t\right)=\{\begin{array}{c}0\\ \frac{8}{t_{\ln }}\left\{{t_{\mathrm{xn}}\left(\frac{t-t_0\left(C_m,n\right)}{t_{\ln }}\right)}^3+\left(\frac{t_{\mathrm{xn}}^3}{t_{\ln }^2}+2t_{\mathrm{xn}}\right)\left(\frac{t-t_0\left(C_m,n\right)}{t_{\ln }}\right)-R\left(t\right)\right\}\begin{array}{c}\mathrm{if}\left(t<t_0\left(C_m,n\right)\right)\\ \mathrm{if}\left(t_0\left(C_m,n\right)<t<t_0\left(C_m,n\right)+t_{\ln }\right)\\ \mathrm{if}\left(t_0\left(C_m,n\right)+t_{\ln }<t\right)\end{array}\\ 0\end{array}& \left(40\right)\end{array}$ $\begin{array}{cc}R\left(t\right)=\{\begin{array}{cc}0& \mathrm{if}\left(t<t_0\left(C_m,n\right)\right)\\ {t_{\ln }\left(\frac{t-t_0\left(C_m,n\right)}{t_{\ln }}\right)}^3+\frac{3t_{\mathrm{xn}}^2}{t_{\ln }}\left(\frac{t-t_0\left(C_m,n\right)}{t_{\ln }}\right)& \mathrm{if}\left(t_0\left(C_m,n\right)\le t\le t_0\left(C_m,n\right)+t_{\ln }\bigcap t_{\mathrm{xn}}>t-t_0\left(C_m,n\right)\right)\\ 3{t_{\mathrm{xn}}\left(\frac{t-t_0\left(C_m,n\right)}{t_{\ln }}\right)}^2+\frac{t_{\mathrm{xn}}^3}{t_{\ln }^2}& \mathrm{if}\left(t_0\left(C_m,n\right)\le t\le t_0\left(C_m,n\right)+t_{\ln }\bigcap t_{\mathrm{xn}}\le t-t_0\left(C_m,n\right)\right)\\ 0& \mathrm{if}\left(t_0\left(C_m,n\right)+t_{\ln }<t\right)\end{array}& \left(41\right)\end{array}$

The measurement device 1 calculates a deflection amount C_std(C_m,t) caused by the C_m-th vehicle according to Equation (42). FIG. 11 shows an example of the deflection amount C_std(C_m,t) caused by the C_m-th vehicle with the number of axles n=4.

$\begin{array}{cc}{C}_{\mathrm{std}}\left(C_m,t\right)=\sum _{n=1}^{a_T\left(C_m\right)}{w}_{\mathrm{std}}\left(a_w\left(C_m,n\right),t\right)& \left(42\right)\end{array}$

The measurement device 1 further calculates a deflection amount T_std(t) caused by the railway vehicle 6 according to Equation (43). FIG. 12 shows an example of the deflection amount T_std(t) caused by the railway vehicle 6 with the number of vehicles C_T=16. In FIG. 12, the broken line indicates 16 deflection amounts C_std(1,t) to C_std(16,t).

$\begin{array}{cc}{T}_{\mathrm{std}}\left(t\right)=\sum _{C_m=1}^{C_T}{C}_{\mathrm{std}}\left(C_m,t\right)& \left(43\right)\end{array}$

Next, the measurement device 1 generates a deflection amount T_{std_lp}(t), obtained by performing filter processing on the deflection amount T_std(t), in order to reduce the vibration component having a fundamental frequency F_M included in the deflection amount T_std(t) and a harmonic of the vibration component. The filter processing may be, for example, low-pass filter processing or band-pass filter processing.

Specifically, first, the measurement device 1 calculates a power spectrum density by performing fast Fourier transform processing on the deflection amount T_std(t) and calculates a peak of the power spectrum density as the fundamental frequency F_M. Then, the measurement device 1 calculates a basic cycle T_M based on the fundamental frequency F_M according to Equation (44), and calculates a moving average interval k_nM, adjusted to a time resolution of the data, by dividing the basic cycle T_M by ΔT as in Equation (45). The basic cycle T_M is a cycle corresponding to the fundamental frequency F_M, and T_M>2ΔT.

$\begin{array}{cc}{T}_M=\frac{1}{f_M}& \left(44\right)\end{array}$ $\begin{array}{cc}{k}_{\mathrm{mM}}=2\lfloor \frac{T_M}{2\Delta T}\rfloor +1& \left(45\right)\end{array}$

Then, the measurement device 1 performs, as the filter processing, moving average processing on the deflection amount T_std(t) in the basic cycle T_M according to Equation (46) to calculate the deflection amount T_{std_lp}(t) in which the vibration component included in the deflection amount T_std(t) is reduced. In the moving average processing, not only a necessary calculation amount is small, but also an attenuation amount of the signal component of the fundamental frequency F_M and the harmonic component of the signal component is very large, so that the deflection amount T_{std_lp}(t) in which the vibration component is effectively reduced is obtained. FIG. 13 shows an example of the deflection amount T_{std_lp}(t). As shown in FIG. 13, the deflection amount T_{std_lp}(t) from which the vibration component included in the deflection amount T_std(t) is almost removed is obtained.

$\begin{array}{cc}{T}_{\mathrm{std}\_\mathrm{lp}}\left(k\right)=\frac{1}{k_{\mathrm{mM}}}\sum _{n=k-\frac{k_{\mathrm{mM}}-1}{2}}^{k+\frac{k_{\mathrm{mM}}-1}{2}}{T}_{\mathrm{std}}\left(n\right)& \left(46\right)\end{array}$

The measurement device 1 may generate the deflection amount T_{std_lp}(t) by performing, as the filter processing, FIR filter processing for attenuating a signal component having a frequency equal to or higher than the fundamental frequency F_M, on the deflection amount T_std(t). In the FIR filter processing, although the calculation amount is larger than that of the moving average processing, all signal components having a frequency equal to or higher than the fundamental frequency F_f can be attenuated.

FIG. 14 shows the measurement data u_lp(t) shown in FIG. 6 and the deflection amount T_{std_lp}(t) shown in FIG. 13 in an overlapping manner. The deflection amount T_{std_lp}(t) is considered to be a deflection amount proportional to the load of the railway vehicle 6 passing through the superstructure 7, and it is assumed that a linear function of the deflection amount T_{std_lp}(t) is substantially equal to the measurement data u_lp(t). That is, the measurement device 1 approximates the measurement data u_lp(t) by the linear function of the deflection amount T_{std_lp}(t) as in Equation (47). An approximate time interval is a time interval between the entry time point t_i and the exit time point t_o or a time interval in which the amplitude of the deflection amount T_{std_lp}(t) is 0.

u _lp(t)≅c₁T_{std_lp}(t)+c₀  (47)

Then, the measurement device 1 calculates a first-order coefficient c₁ and a zero-order coefficient c₀ of the linear function represented by Equation (47). For example, the measurement device 1 calculates, using a least-squares method, the first-order coefficient c₁ and the zero-order coefficient c₀ at which an error e(t) represented by Equation (48), that is, a difference between the measurement data u_lp(t) and the linear function of Equation (47) is minimized.

e(t)=u_lp(t)c₁−T_{std_lp}(t)+c₀

t _i ≤t≤t _o  (48)

The first-order coefficient c₁ and the zero-order coefficient c₀ are calculated according to Equation (49) and Equation (50), respectively. A data section corresponding to the approximate time interval is set as k_a≤k≤k_b.

$\begin{array}{cc}{c}_1=\left\{n\sum _{k=k_a}^{k_b}{u}_{\mathrm{lp}}\left(k\right)T_{{\mathrm{std}}_{\mathrm{lp}}}\left(k\right)-\sum _{k=k_a}^{k_b}{T}_{{\mathrm{std}}_{\mathrm{lp}}}\left(k\right)\sum _{k=k_a}^{k_b}{u}_{\mathrm{lp}}\left(k\right)\right\}/\left\{n\sum _{k=k_a}^{k_b}{T_{{\mathrm{std}}_{\mathrm{lp}}}\left(k\right)}^2-\sum _{k=k_a}^{k_b}{T_{{\mathrm{std}}_{\mathrm{lp}}}\left(k\right)}^2\right\}& \left(49\right)\end{array}$ $n=\sum _{k=k_a}^{k_b}1$ $\begin{array}{cc}\begin{array}{cc}{c}_0=\left\{\sum _{k=k_a}^{k_b}{u}_{\mathrm{lp}}\left(k\right)-c_1\sum _{k=k_a}^{k_b}{T}_{\mathrm{std}\_\mathrm{lp}}\left(k\right)\right\}/n& \left(50\right)\end{array}& \end{array}$ $n=\sum _{k=k_a}^{k_b}1$

Then, the measurement device 1 calculates a deflection amount T_{Estd_lp}(t) in which the deflection amount T_{std_lp}(t) is adjusted using the first-order coefficient c₁ and the zero-order coefficient c₀, as in Equation (51). As shown in Equation (51), the deflection amount T_{Estd_lp}(t) basically corresponds to a right side of Equation (47), and the zero-order coefficient c₀ is set to 0 in an interval before the entry time point t_i and an interval after the exit time point t_o. FIG. 15 shows an example of the deflection amount T_{Estd_lp}(t).

$\begin{array}{cc}{T}_{\mathrm{Estd}\_\mathrm{lp}}\left(t\right)=\{\begin{array}{cc}t<t_i& c_1{T}_{\mathrm{std}\_\mathrm{lp}}\left(t\right)\\ t_i\le t\le t_o& c_1{T}_{\mathrm{std}\_\mathrm{lp}}\left(t\right)+c_0\\ t_o<t& c_1{T}_{\mathrm{std}\_\mathrm{lp}}\left(t\right)\end{array}& \left(51\right)\end{array}$

As in Equation (52), it is assumed that a linear function of the deflection amount T_std(t) using the first-order coefficient c₁ calculated according to Equation (49) and the zero-order coefficient c₀ calculated according to Equation (50) is substantially equal to the measurement data u(t).

u(t)≅C₁T_std(t)+C₀

t _i ≤t≤t _o  (52)

A deflection amount T_Estd(t) obtained by adjusting the deflection amount T_std(t) using the first-order coefficient c₁ and the zero-order coefficient c₀ is calculated according to Equation (53). A right side of Equation (53) is obtained by replacing T_{std_lp}(t) on a right side of Equation (51) with T_std(t). FIG. 16 shows an example of the deflection amount T_Estd(t).

$\begin{array}{cc}{T}_{\mathrm{Estd}}\left(t\right)=\{\begin{array}{cc}t<t_i& c_1{T}_{\mathrm{std}}\left(t\right)\\ t_i\le t\le t_o& c_1{T}_{\mathrm{std}}\left(t\right)+c_0\\ t_o<t& c_1{T}_{\mathrm{std}}\left(t\right)\end{array}& \left(53\right)\end{array}$

Next, the measurement device 1 calculates an amplitude ratio R_T between the deflection amount T_{Estd_lp}(t) and the deflection amount T_{std_lp}(t) in a predetermined interval according to Equation (54) with t=kΔT. In Equation (54), a numerator is an average value of n+1 samples of the deflection amount T_{Estd_lp}(t) included in a predetermined interval which is a part of an interval in which the waveform of the deflection amount T_{Estd_lp}(t) and the waveform of the deflection amount T_{std_lp}(t) are shifted, and a denominator is an average value of n+1 samples of the deflection amount T_{std_lp}(t) included in the predetermined interval. FIG. 17 shows an example of a relation between the deflection amount T_{Estd_lp}(t) and the deflection amount T_{std_lp}(t) and a predetermined interval T_avg for calculating the average values thereof.

$\begin{array}{cc}{R}_T=\left(\frac{1}{n+1}\sum _{k=k_0}^{k_0+n}{T}_{\mathrm{Estd}\_\mathrm{lp}}\left(k\right)\right)/\left(\frac{1}{n+1}\sum _{k=k_0}^{k_0+n}{T}_{\mathrm{std}\_\mathrm{lp}}\left(k\right)\right)& \left(54\right)\end{array}$

Next, the measurement device 1 compares a product R_TT_{std_lp}(t) of the amplitude ratio R_T and the deflection amount T_{std_lp}(t) with the zero-order coefficient c₀ to calculate an offset T_{offset_std} (t). Specifically, the measurement device 1 calculates the offset T_{offset_std}(t) by replacing, with the zero-order coefficient c₀, an interval of the product R_TT_{std_lp}(t) in which an absolute value of the product R_TT_{std_lp}(t) of the amplitude ratio R_T and the deflection amount T_{std_lp}(t) is bigger than an absolute value of the zero-order coefficient c₀, as in Equation (55). FIG. shows an example of the offset T_{offset_std}(t). In the example of FIG. 18, since the amplitude of the deflection amount T_{std_lp}(t) is 0 or negative, the measurement device 1 calculates the offset T_{offset_std}(t) by replacing, with the zero-order coefficient c₀, an interval in which the product R_TT_{std_lp}(t) is smaller than the zero-order coefficient c₀.

$\begin{array}{cc}{T}_{\mathrm{offset}\_\mathrm{std}}\left(t\right)=\{\begin{array}{cc}{R}_T{T}_{\mathrm{std}\_\mathrm{lp}}\left(t\right)\ge c_0& R_T{T}_{\mathrm{std}\_\mathrm{lp}}\left(t\right)\\ R_T{T}_{\mathrm{std}\_\mathrm{lp}}\left(t\right)<c_0& c_0\end{array}& \left(55\right)\end{array}$

Then, the measurement device 1 adds a product c₁T_std(t) of the first-order coefficient c₁ and the deflection amount T_std(t) and the offset T_offsetstd(t) to calculate a deflection amount T_EOstd(t) as the static response, as in Equation (56). The deflection amount T_EOstd(t) corresponds to the static response when the railway vehicle 6 passes through the superstructure 7. FIG. 19 shows an example of the deflection amount T_EOstd(t). FIG. 20 shows a relation between the measurement data u(t) and the deflection amount T_EOstd(t).

T _EOstd(t)=C₁T_std(t)+T_{offset_std}(t)  (56)

When a vibration frequency of the static response caused by the railway vehicle 6 passing through the superstructure 7 is greatly separated from the natural vibration frequency of the superstructure 7, the natural vibration of the superstructure 7 is less likely to be excited even when the railway vehicle 6 passes through the superstructure 7, so that the dynamic response hardly occurs, and the measurement data u(t) has a waveform close to the static response. The measurement data u(t) shown in FIG. 20 has a waveform close to the deflection amount T_EOstd(t) as the static response, and it can be seen that the measurement data u (t) includes almost no dynamic response.

On the other hand, when the vibration frequency of the static response is close to the natural vibration frequency of the superstructure 7, the natural vibration of the superstructure 7 is excited, so that not only the static response but also the dynamic response is generated in a superimposed manner. Therefore, the measurement data u(t) includes not only the static response but also the dynamic response. FIG. 21 shows an example of the measurement data u(t) including the static response and the dynamic response.

FIG. 22 shows an example of the measurement data u_lp(t) calculated according to Equation (5) based on the measurement data u(t) of FIG. 21. FIG. 23 shows an example of a relation between the measurement data u_lp(t) of FIG. 22, the entry time point t_i, and the exit time point t_o. FIG. 24 shows an example of the deflection amount T_std(t) calculated according to Equation (43) using the entry time point t_i and the exit time point t_o of FIG. 23 and the number of vehicles C_T calculated according to Equation (8). FIG. 25 shows the measurement data u_lp(t) shown in FIG. 22 and the deflection amount T_{std_lp}(t) calculated according to Equation (46) based on the deflection amount T_std(t) of FIG. 23 in an overlapping manner. FIG. 26 shows an example of the offset T_{offset_std}(t) calculated according to Equation (55) based on the deflection amount T_{std_lp}(t) of FIG. 25. FIG. shows an example of the deflection amount T_EOstd(t) calculated according to Equation (56) based on the deflection amount T_std(t) of FIG. 24 and the offset T_{offset_std} (t) in FIG. 26. FIG. 28 shows a relation between the measurement data u(t) shown in FIG. 21 and the deflection amount T_EOstd(t) shown in FIG. 27.

As shown in FIG. 28, it can be seen that, with respect to the measurement data u(t) including not only the static response but also the dynamic response, the deflection amount T_EOstd(t) does not include the vibration of the dynamic response and is the static response. As described above, with the measurement method according to the present embodiment, the deflection amount T_Eo_std(t) which is the static response can be calculated based on the measurement data u(t) which has a large dynamic response and an unclear static response.

1-3. Procedure of Measurement Method

FIG. 29 is a flowchart showing an example of a procedure of a measurement method according to the first embodiment. In the present embodiment, the measurement device 1 executes the procedure shown in FIG. 29.

As shown in FIG. 29, first, in an observation data acquisition step S10, the measurement device 1 acquires acceleration data a(k) which is observation data output from the sensor 2 which is an observation device.

Next, in a first measurement data generation step S20, the measurement device 1 generates, based on the acceleration data a(k) which is the observation data acquired in step S10, the measurement data u(t) which is first measurement data based on the acceleration as the physical quantity, which is the response to the actions of the plurality of axles of the railway vehicle 6 moving on the superstructure 7 on the observation point R. An example of a procedure of the first measurement data generation step S20 will be described later.

Next, in a second measurement data generation step S30, the measurement device 1 generates the measurement data u_lp(t), which is second measurement data in which a vibration component is reduced by the measurement device 1 performing filter processing on the measurement data u(t) generated in step S20. For example, the measurement device 1 performs, as the filter processing, low-pass filter processing for attenuating the vibration component having a frequency equal to or higher than the fundamental frequency F_f of the measurement data u(t). An example of a procedure of the second measurement data generation step S30 will be described later.

Next, in an observation information generation step S40, the measurement device 1 generates the observation information including the entry time point t_i and the exit time point t_o of the railway vehicle 6 with respect to the superstructure 7. The entry time point t_i is the time point at which the leading axle of the plurality of axles of the railway vehicle 6 passes through the entry end of the superstructure 7, and the exit time point t_o is the time point at which the rearmost axle of the plurality of axles of the railway vehicle 6 passes through the exit end of the superstructure 7. In the present embodiment, the measurement device 1 calculates the entry time point t_i and the exit time point t_o based on the measurement data u_lp(t) generated in step S30. Further, the measurement device 1 generates the number of vehicles C_T. An example of a procedure of the observation information generation step S40 will be described later.

Next, in an average velocity calculation step S50, the measurement device 1 calculates the average velocity v_a of the railway vehicle 6 based on the observation information generated in step S40 and the environment information which is created in advance and includes the dimension of the railway vehicle 6 and the dimension of the superstructure 7. The environment information includes the length L_B of the superstructure 7, the position L_x of the observation point R, the length L_c(C_m) of each vehicle of the railway vehicle 6, the number of axles a_T(C_m) of each vehicle, and the distance L_a(a_w(C_m,n)) between the axles corresponding to the position of each of the plurality of axles of the railway vehicle 6. An example of a procedure of the average velocity calculation step S50 will be described later.

Next, in a first deflection amount calculation step S60, the measurement device 1 calculates the deflection amount T_std(t), which is a first deflection amount of the superstructure 7 caused by the railway vehicle 6, based on the approximate equation of the deflection of the superstructure 7, which is Equation (35), the observation information generated in step S40, the environment information, and the average velocity v_a of the railway vehicle 6 calculated in step S50. Specifically, the measurement device 1 calculates the deflection amount w_std(a_w(C_m,n),t) of the superstructure 7 caused by each of the plurality of axles based on the approximate equation of the deflection of the superstructure 7, the observation information, the environment information, and the average velocity v_a, and calculates the deflection amount T_std(t) by adding the deflection amount w_std(a_w(C_m,n),t) of the superstructure 7 caused by each of the plurality of axles. An example of a procedure of the first deflection amount calculation step S60 will be described later.

Next, in a second deflection amount calculation step S70, the measurement device 1 calculates the deflection amount T_{std_lp}(t), which is a second deflection amount in which a vibration component is reduced by the measurement device 1 performing filter processing on the deflection amount T_std(t) calculated in step S60. For example, the measurement device 1 performs, as the filter processing, low-pass filter processing for attenuating the vibration component having a frequency equal to or higher than the fundamental frequency F_M of the deflection amount T_std(t). An example of a procedure of the second deflection amount calculation step S70 will be described later.

Next, in a coefficient calculation step S80, the measurement device 1 approximates the measurement data u_lp(t) generated in step S30 with the linear function of the deflection amount T_{std_lp}(t) calculated in step S70, and calculates the first-order coefficient c₁ and the zero-order coefficient c₀ of the linear function. Specifically, the measurement device 1 approximates the measurement data u_lp(t) with the linear function of the deflection amount T_{std_lp}(t) as in Equation (47), and calculates the first-order coefficient c₁ and the zero-order coefficient c₀ according to Equation (49) and Equation (50) using the least-squares method.

Next, in a third deflection amount calculation step S90, the measurement device 1 calculates the deflection amount T_{Estd_lp}(t), which is a third deflection amount, based on the first-order coefficient c₁ and the zero-order coefficient c₀ calculated in step S80 and the deflection amount T_{std_lp}(t) calculated in step S70. Specifically, the measurement device 1 calculates the deflection amount T_{Estd_lp}(t), which is a product c₁T_{std_lp}(t) of the first-order coefficient c₁ and the deflection amount T_{std_lp}(t) in the interval before the entry time point t_i and the interval after the exit time point t_o, and is a sum of the product c₁T_{std_lp}(t) and the zero-order coefficient c₀ in an interval from the entry time point t_i to the exit time point t_o, as in Equation (51).

Next, in an offset calculation step S100, the measurement device 1 calculates the offset T_{offset_std}(t) based on the zero-order coefficient c₀ calculated in step S80, the deflection amount T_{std_lp}(t) calculated in step S70, and the deflection amount T_{Estd_lp}(t) calculated in step S90. An example of a procedure of the offset calculation step S100 will be described later.

Next, in a static response calculation step S110, the measurement device 1 calculates the deflection amount T_EOstd(t) as the static response by adding the product c₁T_std(t) of the first-order coefficient c₁ calculated in step S80 and the deflection amount T_std(t) calculated in step S60 and the offset T_{offset_std}(t) calculated in step S100, as in Equation (56).

Next, in a measurement data output step S120, the measurement device 1 outputs the measurement data including the deflection amount T_EOstd(t) as the static response calculated in step S110 to the monitoring device 3. Specifically, the measurement device 1 transmits the measurement data to the monitoring device 3 via the communication network 4. The measurement data may include the measurement data u(t) and u_lp(t), the deflection amounts T_std(t), T_{std_lp}(t), and T_{Estd_lp}(t), and the like, in addition to the deflection amount T_EOstd(t).

Then, the measurement device 1 repeats the processing of steps S10 to S120 until the measurement is completed in step S130.

FIG. 30 is a flowchart showing the example of the procedure of the first measurement data generation step S20 of FIG. 29.

As shown in FIG. 30, in step S201, the measurement device 1 integrates the acceleration data a(t) output from the sensor 2 to generate velocity data v(t) as in Equation (1).

Then, in step S202, the measurement device 1 integrates the velocity data v(t) generated in step S201 to generate the measurement data u(t), as in Equation (2).

As described above, in the present embodiment, the measurement data u(t) is data of the displacement of the superstructure 7 caused by the railway vehicle 6 which is a moving object moving on the superstructure 7 which is a structure, and is data obtained by integrating twice the acceleration in the direction intersecting the surface of the superstructure 7 on which the railway vehicle 6 moves. Therefore, the measurement data u(t) includes data having a waveform projecting in a positive direction or a negative direction, specifically, a rectangular waveform, a trapezoidal waveform, or a sine half-wave waveform. The rectangular waveform includes not only an accurate rectangular waveform but also a waveform approximate to the rectangular waveform. Similarly, the trapezoidal waveform includes not only an accurate trapezoidal waveform but also a waveform approximate to the trapezoidal waveform. Similarly, the sine half-wave waveform includes not only an accurate sine half-wave waveform but also a waveform approximate to the sine half-wave waveform.

FIG. 31 is a flowchart showing the example of the procedure of the second measurement data generation step S30 of FIG. 29.

As shown in FIG. 31, in step S301, the measurement device 1 calculates the power spectrum density by performing fast Fourier transform processing on the measurement data u(t) calculated in step S202 of FIG. 30, and calculates the peak of the power spectrum density as the fundamental frequency F_f.

Then, in step S302, the measurement device 1 generates the measurement data u_lp(t) by performing low-pass filter processing for attenuating the vibration component having a frequency equal to or higher than the fundamental frequency F_f of the measurement data u(t). The measurement device 1 may generate the measurement data u_lp(t) by performing, as the low-pass filter processing, moving average processing on the measurement data u(t) in the basic cycle T_f corresponding to the fundamental frequency F_f, as in Equation (5). Alternatively, the measurement device 1 may generate the measurement data u_lp(t) by performing, as the low-pass filter processing, FIR filter processing for attenuating the signal component having a frequency equal to or higher than the fundamental frequency F_f on the measurement data u(t).

FIG. 32 is a flowchart showing an example of the procedure of the observation information generation step S40 of FIG. 29.

As shown in FIG. 32, first, in step S401, the measurement device 1 calculates, as the amplitude u_a, the average value in the interval from a time point t₁ to a time point t₂ in which the amplitude of the measurement data u_lp(t) generated in step S302 of FIG. 31 is shifted, according to Equation (6).

Then, in step S402, the measurement device 1 calculates, as the entry time point t_i, a first time point at which the amplitude of the measurement data u_lp(t) matches or exceeds the threshold C_Lu_a which is the product of the predetermined coefficient C_L and the amplitude u_a calculated in step S401.

In step S403, the measurement device 1 calculates, as the exit time point t_o, a second time point after the first time point at which the amplitude of the measurement data u_lp(t) matches or exceeds the threshold C_Lu_a.

In step S404, the measurement device 1 calculates the difference between the exit time point t_o and the entry time point t_i as the passing time t_s as in Equation (7).

Then, in step S405, the measurement device 1 calculates, as the number of vehicles C_T of the railway vehicle 6, the maximum integer less than or equal to the number obtained by subtracting 1 from the product t_sF_f of the passing time t_s calculated in step S404 and the fundamental frequency F_f calculated in step S301 of FIG. 31, as in Equation (8).

Then, in step S406, the measurement device 1 generates the observation information including the entry time point t_i calculated in step S402, the exit time point t_o calculated in step S403, the passing time t_s calculated in step S404, and the number of vehicles C_T calculated in step S405.

FIG. 33 is a flowchart showing an example of the procedure of the average velocity calculation step S50 of FIG. 29.

In step S501, the measurement device 1 calculates, based on the environment information, the distance D_wa(a_w(C_T,a_T(C_T))) from the leading axle to the rearmost axle of the railway vehicle 6 according to Equation (11).

In step S502, the measurement device 1 calculates, based on the environment information, the distance from the entry end to the exit end of the superstructure 7. In the present embodiment, the distance from the entry end to the exit end of the superstructure 7 is the length L_B of the superstructure 7 included in the environment information.

Then, in step S503, the measurement device 1 calculates the average velocity v_a of the railway vehicle 6 according to Equation (12) based on the entry time point t_i and the exit time point t_o included in the observation information generated in step S406 of FIG. 32, the distance D_wa(a_w(C_T,a_T(C_T))) from the leading axle to the rearmost axle of the railway vehicle 6 calculated in step S501, and the length L_B of the superstructure 7 which is the distance from the entry end to the exit end of the superstructure 7 calculated in step S502.

FIG. 34 is a flowchart showing an example of the procedure of the first deflection amount calculation step S60 of FIG. 29.

First, in step S601, the measurement device 1 calculates, based on the environment information, each distance D_wa(a_w(C_m,n)) from the leading axle to the n-th axle of the C_m-th vehicle of the railway vehicle 6, according to Equation (10).

Next, in step S602, the measurement device 1 calculates the time t_xn required for a certain axle of the railway vehicle 6 to reach the position L_x of the observation point R from the entry end of the superstructure 7 using the position L_x of the observation point R included in the environment information and the average velocity v_a calculated in step S503 of FIG. 33, according to Equation (37).

In step S603, the measurement device 1 calculates the time t_ln required for a certain axle of the railway vehicle 6 to pass through the superstructure 7 using the length L_B of the superstructure 7, which is the distance from the entry end to the exit end of the superstructure 7 calculated in step S502 of FIG. 33, and the average velocity v_a, according to Equation (38).

In step S604, the measurement device 1 calculates the time point t₀(C_m,n) at which the n-th axle of the C_m-th vehicle of the railway vehicle 6 reaches the entry end of the superstructure 7 using the entry time point t_i included in the observation information generated in step S406 of FIG. 32, the distance D_wa(a_w(C_m,n)) calculated in step S601, and the average velocity v_a, according to Equation (39).

Next, in step S605, the measurement device 1 calculates the deflection amount w_std (a_w(C_m,n), t) of the superstructure 7 caused by the n-th axle of the C_m-th vehicle using the approximate equation of the deflection of the superstructure 7, which is Equation (35), the time t_xn calculated in step S602, the time t_ln calculated in step S603, and the time point t₀(C_m,n) calculated in step S604, according to Equation (40).

Next, in step S606, the measurement device 1 calculates, according to Equation (42), the deflection amount C_std(C_m,t) of the superstructure 7 caused by each vehicle by adding the deflection amount w_std(a_w(C_m,n),t) of the superstructure 7 caused by each axle of each vehicle calculated in step S605.

Then, in step S607, the measurement device 1 calculates, according to Equation (43), the deflection amount T_std(t) of the superstructure 7 caused by the railway vehicle 6 by adding the deflection amount C_std(C_m, t) of the superstructure 7 caused by each vehicle, which is calculated in step S606.

FIG. 35 is a flowchart showing the example of the procedure of the second deflection amount calculation step S70 of FIG. 29.

As shown in FIG. 35, in step S701, the measurement device 1 calculates the power spectrum density by performing fast Fourier transform processing on the deflection amount T_std(t) calculated in step S607 of FIG. 34, and calculates the peak of the power spectrum density as the fundamental frequency F_M.

Then, in step S702, the measurement device 1 calculates the deflection amount T_{std_lp}(t) by performing low-pass filter processing for attenuating the vibration component having a frequency equal to or higher than the fundamental frequency F_M of the deflection amount T_std(t). The measurement device 1 may calculate the deflection amount T_{std_lp}(t) by performing, as the low-pass filter processing, moving average processing on the deflection amount T_std(t) in the basic cycle T_M corresponding to the fundamental frequency F_M, according to Equation (46). Alternatively, the measurement device 1 may calculate the deflection amount T_{std_lp}(t) by performing, as the low-pass filter processing, FIR filter processing for attenuating the signal component having a frequency equal to or higher than the fundamental frequency F_M on the deflection amount T_std(t).

FIG. 36 is a flowchart showing the example of the procedure of the offset calculation step S100 of FIG. 29.

As shown in FIG. 36, in step S1001, the measurement device 1 calculates the amplitude ratio R_T in a predetermined interval between the deflection amount T_{Estd_lp}(t) calculated in step S90 of FIG. 29 and the deflection amount T_{std_lp}(t) calculated in step S702 of FIG. 35, according to Equation (54).

Then, in step S1002, the measurement device 1 calculates the offset T_{offset_std}(t) by replacing, with the zero-order coefficient c₀, the interval of the product R_TT_{std_lp}(t) in which the absolute value of the product R_TT_{std_lp}(t) of the amplitude ratio R_T calculated in step S1001 and the deflection amount T_{std_lp}(t) is bigger than the absolute value of the zero-order coefficient c₀ calculated in step S80 of FIG. 29, as in Equation (55).

1-4. Configuration of Observation Device, Measurement Device, and Monitoring Device

FIG. 37 is a diagram showing a configuration example of the sensor 2 which is the observation device, the measurement device 1, and the monitoring device 3.

As shown in FIG. 37, the sensor 2 includes a communication unit 21, an acceleration sensor 22, a processor 23, and a storage unit 24.

The storage unit 24 is a memory that stores various programs, data, and the like for the processor 23 to perform calculation processing and control processing. Further, the storage unit 24 stores programs, data, and the like for the processor 23 to implement predetermined application functions.

The acceleration sensor 22 detects an acceleration generated in each axial direction of the three axes.

The processor 23 controls the acceleration sensor 22 by executing an observation program 241 stored in the storage unit 24, generates observation data 242 based on the acceleration detected by the acceleration sensor 22, and stores the generated observation data 242 in the storage unit 24. In the present embodiment, the observation data 242 is the acceleration data a(k).

The communication unit 21 transmits the observation data 242 stored in the storage unit 24 to the measurement device 1 under the control of the processor 23.

As shown in FIG. 37, the measurement device 1 includes a first communication unit 11, a second communication unit 12, a storage unit 13, and a processor 14.

The first communication unit 11 receives the observation data 242 from the sensor 2, and outputs the received observation data 242 to the processor 14. As described above, the observation data 242 is the acceleration data a(k).

The storage unit 13 is a memory that stores programs, data, and the like for the processor 14 to perform calculation processing and control processing. The storage unit 13 stores programs, data, and the like for the processor 14 to implement predetermined application functions. The processor 14 may receive various programs, data, and the like via the communication network 4 and store the programs, data, and the like in the storage unit 13.

The processor 14 generates measurement data 135 based on the observation data 242 received by the first communication unit 11 and environment information 132 stored in advance in the storage unit 13, and stores the generated measurement data 135 in the storage unit 13.

In the present embodiment, the processor 14 functions as an observation data acquisition unit 141, a first measurement data generation unit 142, a second measurement data generation unit 143, an observation information generation unit 144, an average velocity calculation unit 145, a first deflection amount calculation unit 146, a second deflection amount calculation unit 147, a coefficient calculation unit 148, a third deflection amount calculation unit 149, an offset calculation unit 150, a static response calculation unit 151, and a measurement data output unit 152 by executing a measurement program 131 stored in the storage unit 13. That is, the processor 14 includes the observation data acquisition unit 141, the first measurement data generation unit 142, the second measurement data generation unit 143, the observation information generation unit 144, the average velocity calculation unit 145, the first deflection amount calculation unit 146, the second deflection amount calculation unit 147, the coefficient calculation unit 148, the third deflection amount calculation unit 149, the offset calculation unit 150, the static response calculation unit 151, and the measurement data output unit 152.

The observation data acquisition unit 141 acquires the observation data 242 received by the first communication unit 11, and stores the observation data 242 in the storage unit 13 as observation data 133. That is, the observation data acquisition unit 141 performs the processing of the observation data acquisition step S10 in FIG. 29.

The first measurement data generation unit 142 reads the observation data 133 stored in the storage unit 13, and generates, based on the acceleration data a(t) which is the observation data 133, the measurement data u(t) which is the first measurement data based on the acceleration as the physical quantity, which is the response to the actions of the plurality of axles of the railway vehicle 6 moving on the superstructure 7 on the observation point R. Specifically, the first measurement data generation unit 142 integrates the acceleration data a(t), which is the observation data 133, to generate the velocity data v(t) as in Equation (1), and further integrates the velocity data v(t) to generate the measurement data u(t) as in Equation (2). That is, the first measurement data generation unit 142 performs the processing of the first measurement data generation step S20 in FIG. 29, specifically, the processing of steps S201 and S202 in FIG. 30.

The second measurement data generation unit 143 generates the measurement data u_lp(t), which is the second measurement data in which the vibration component is reduced by performing filter processing on the measurement data u(t) generated by the first measurement data generation unit 142. For example, the second measurement data generation unit 143 performs, as the filter processing, low-pass filter processing for attenuating the vibration component having a frequency equal to or higher than the fundamental frequency F_f of the measurement data u(t). Specifically, the second measurement data generation unit 143 calculates the power spectrum density by performing fast Fourier transform processing on the measurement data u(t), calculates the peak of the power spectrum density as the fundamental frequency F_f, and generates the measurement data u_lp(t) by performing low-pass filter processing for attenuating the vibration component having a frequency equal to or higher than the fundamental frequency F_f of the measurement data u(t). The second measurement data generation unit 143 may generate the measurement data u_lp(t) by performing, as the low-pass filter processing, moving average processing on the measurement data u(t) in the basic cycle T_f corresponding to the fundamental frequency F_f, as in Equation (5). Alternatively, the second measurement data generation unit 143 may generate the measurement data u_lp(t) by performing, as the low-pass filter processing, FIR filter processing for attenuating the signal component having a frequency equal to or higher than the fundamental frequency F_f on the measurement data u(t). That is, the second measurement data generation unit 143 performs the processing of the second measurement data generation step S30 in FIG. 29, specifically, the processing of steps S301 and S302 in FIG. 31.

Based on the measurement data u_lp(t) generated by the second measurement data generation unit 143, the observation information generation unit 144 generates observation information 134 including the entry time point t_i and the exit time point t_o of the railway vehicle 6 with respect to the superstructure 7, and stores the observation information 134 in the storage unit 13. Specifically, first, the observation information generation unit 144 calculates, as the amplitude u_a, the average value of the interval from the time point t₁ to the time point t₂ in which the amplitude of the measurement data u_lp(t) is shifted, according to Equation (6). Next, the observation information generation unit 144 calculates, as the entry time point t_i, the first time point at which the amplitude of the measurement data u_lp(t) matches or exceeds the threshold C_Lu_a which is the product of the predetermined coefficient C_L and the amplitude u_a. The observation information generation unit 144 calculates, as the exit time point t_o, the second time point after the first time point at which the amplitude of the measurement data u_lp(t) matches or exceeds the threshold C_Lu_a. The observation information generation unit 144 calculates the difference between the exit time point t_o and the entry time point t_i as the passing time t_s as in Equation (7). Next, the observation information generation unit 144 calculates, as the number of vehicles C_T of the railway vehicle 6, a maximum integer less than or equal to the number obtained by subtracting 1 from the product t_sF_f of the passing time t_s and the fundamental frequency F_f, as in Equation (8). Then, the observation information generation unit 144 generates the observation information 134 including the entry time point t_i, the exit time point t_o, the passing time t_s, and the number of vehicles C_T. That is, the observation information generation unit 144 performs the processing of the observation information generation step S40 in FIG. 29, specifically, the processing of steps S401 to S406 in FIG. 32.

The average velocity calculation unit 145 calculates the average velocity v_a of the railway vehicle 6 based on the observation information 134 stored in the storage unit 13 and the environment information 132 which is created in advance and stored in the storage unit 13 and includes the dimensions of the railway vehicle 6 and the dimensions of the superstructure 7. Specifically, the average velocity calculation unit 145 calculates, based on the environment information 132, the distance D_wa(a_w(C_T,a_T(C_T))) from the leading axle to the rearmost axle of the railway vehicle 6 according to Equation (11). The average velocity calculation unit 145 calculates the length L_B of the superstructure 7, which is the distance from the entry end to the exit end of the superstructure 7, based on the environment information 132. Then, the average velocity calculation unit 145 calculates the average velocity v_a of the railway vehicle 6 according to Equation (12), based on the entry time point t_i and the exit time point t_o included in the observation information 134, the distance D_wa(a_w(C_T,a_T(C_T))), and the length L_B of the superstructure 7. That is, the average velocity calculation unit 145 performs the processing of the average velocity calculation step S50 in FIG. 29, specifically, the processing of steps S501, S502, and S503 in FIG. 33.

The first deflection amount calculation unit 146 calculates the deflection amount T_std(t), which is the first deflection amount of the superstructure 7 caused by the railway vehicle 6, based on the approximate equation of the deflection of the superstructure 7, which is Equation (35), the observation information 134 stored in the storage unit 13, the environment information 132 stored in the storage unit 13, and the average velocity v_a of the railway vehicle 6 calculated by the average velocity calculation unit 145. Specifically, first, the first deflection amount calculation unit 146 calculates, based on the environment information 132, the distance D_wa(a_w(C_m,n)) from the leading axle to the n-th axle of the C_m-th vehicle of the railway vehicle 6, according to Equation (10). Next, the first deflection amount calculation unit 146 calculates the time t_xn required for a certain axle of the railway vehicle 6 to reach the position L_x of the observation point R from the entry end of the superstructure 7 using the position L_x of the observation point R included in the environment information 132 and the average velocity v_a, according to Equation (37). The first deflection amount calculation unit 146 calculates the time t_ln required for a certain axle of the railway vehicle 6 to pass through the superstructure 7 using the length L_B of the superstructure 7, which is the distance from the entry end to the exit end of the superstructure 7, and the average velocity v_a, according to Equation (38). Further, the first deflection amount calculation unit 146 calculates the time point t₀(C_m,n) at which the n-th axle of the C_m-th vehicle of the railway vehicle 6 reaches the entry end of the superstructure 7 using the entry time point t_i included in the observation information 134, the distance D_wa(a_w(C_m,n)), and the average velocity v_a, according to Equation (39). Next, the first deflection amount calculation unit 146 calculates the deflection amount w_std(a_w(C_m,n),t) of the superstructure 7 caused by the n-th axle of the C_m-th vehicle using the approximate equation of the deflection of the superstructure 7, which is Equation (35), the time t_xn, the time t_ln, and the time point t₀(C_m,n), according to Equation (40). Next, the first deflection amount calculation unit 146 calculates the deflection amount C_std(C_m,t) of the superstructure 7 caused by the C_m-th vehicle using the deflection amount w_std(a_w(C_m,n),t), according to Equation (42). Then, the first deflection amount calculation unit 146 calculates the deflection amount T_std(t) of the superstructure 7 caused by the railway vehicle 6 using the deflection amount C_std(C_m,t), according to Equation (43). That is, the first deflection amount calculation unit 146 performs the processing of the first deflection amount calculation step S60 in FIG. 29, specifically, the processing of steps S601 to S607 in FIG. 34.

The second deflection amount calculation unit 147 calculates the deflection amount T_{std_lp}(t), which is the second deflection amount in which the vibration component is reduced by performing filter processing on the deflection amount T_std(t) calculated by the first deflection amount calculation unit 146. For example, the second deflection amount calculation unit 147 performs, as the filter processing, low-pass filter processing for attenuating the vibration component having a frequency equal to or higher than the fundamental frequency F_M of the deflection amount T_std(t). Specifically, the second deflection amount calculation unit 147 calculates the deflection amount T_std(t) as the fundamental frequency F_M by performing fast Fourier transform processing on the deflection amount T_std(t), and calculates the deflection amount T_{std_lp}(t) by performing low-pass filter processing for attenuating the vibration component having a frequency equal to or higher than the fundamental frequency F_M of the deflection amount T_std(t). The second deflection amount calculation unit 147 may calculate the deflection amount T_{std_lp}(t) by performing, as the low-pass filter processing, moving average processing on the deflection amount T_std(t) in the basic cycle T_M corresponding to the fundamental frequency F_M, according to Equation (46). Alternatively, the second deflection amount calculation unit 147 may calculate the deflection amount T_{std_lp}(t) by performing, as the low-pass filter processing, FIR filter processing for attenuating the signal component having a frequency equal to or higher than the fundamental frequency F_M on the deflection amount T_std(t). That is, the second deflection amount calculation unit 147 performs the processing of the second deflection amount calculation step S70 in FIG. 29, specifically, the processing of steps S701 and S702 in FIG. 35.

The coefficient calculation unit 148 approximates the measurement data u_lp(t) generated by the second measurement data generation unit 143 with the linear function of the deflection amount T_{std_lp}(t) calculated by the second deflection amount calculation unit 147, and calculates the first-order coefficient c₁ and the zero-order coefficient c₀ of the linear function. Specifically, the coefficient calculation unit 148 approximates the measurement data u_lp(t) with the linear function of the deflection amount T_{std_lp}(t) as in Equation (47), and calculates the first-order coefficient c₁ and the zero-order coefficient c₀ according to Equation (49) and Equation (50) using the least-squares method. That is, the coefficient calculation unit 148 performs the processing of the coefficient calculation step S80 in FIG. 29.

The third deflection amount calculation unit 149 calculates the deflection amount T_{Estd_lp}(t), which is the third deflection amount, based on the first-order coefficient c₁ and the zero-order coefficient c₀ calculated by the coefficient calculation unit 148 and the deflection amount T_{std_lp}(t) calculated by the second deflection amount calculation unit 147. Specifically, the third deflection amount calculation unit 149 calculates the deflection amount T_{Estd_lp}(t), which is the product c₁T_{std_lp}(t) of the first-order coefficient c₁ and the deflection amount T_{std_lp}(t) in the interval before the entry time point t_i and the interval after the exit time point t_o, and is the sum of the product c₁T_{std_lp}(t) and the zero-order coefficient c₀ in the interval from the entry time point t_i to the exit time point t_o, as in Equation (51). That is, the third deflection amount calculation unit 149 performs the processing of the third deflection amount calculation step S90 in FIG. 29.

The offset calculation unit 150 calculates the offset T_{offset_std}(t) based on the zero-order coefficient c₀ calculated by the coefficient calculation unit 148, the deflection amount T_{std_lp}(t) calculated by the second deflection amount calculation unit 147, and the deflection amount T_{Estd_lp}(t) calculated by the third deflection amount calculation unit 149. Specifically, the offset calculation unit 150 calculates the amplitude ratio R_T between the deflection amount T_{Estd_lp}(t) and the deflection amount T_{std_lp}(t) in a predetermined interval, according to Equation (54). Then, the offset calculation unit 150 calculates the offset T_{offset_std}(t) by replacing, with the zero-order coefficient c₀, the interval in which the product R_TT_{std_lp}(t) the amplitude ratio R_T and the deflection amount T_{std_lp}(t) is smaller than the zero-order coefficient c₀, as in Equation (55). That is, the offset calculation unit 150 performs the processing of the offset calculation step S100 in FIG. 29, specifically, the processing of steps S1001 and S1002 in FIG. 36.

The static response calculation unit 151 calculates the deflection amount T_EOstd(t) as the static response by adding the product c₁T_std(t) of the first-order coefficient c₁ calculated by the coefficient calculation unit 148 and the deflection amount T_std(t) calculated by the first deflection amount calculation unit 146, and the offset T_{offset_std}(t) calculated by the offset calculation unit 150, as in Equation (56). That is, the static response calculation unit 151 performs the processing of the static response calculation step S110 in FIG. 29.

The deflection amount T_EOstd(t) as the static response is stored in the storage unit 13 as at least a part of the measurement data 135. The measurement data 135 may include the measurement data u(t) and u_lp(t), the deflection amounts T_std (t), T_{std_lp} (t), and T_{Estd_lp} (t), and the like, in addition to the deflection amount T_EOstd(t).

The measurement data output unit 152 reads the measurement data 135 stored in the storage unit 13 and outputs the measurement data 135 to the monitoring device 3. Specifically, the second communication unit 12 transmits the measurement data 135 stored in the storage unit 13 to the monitoring device 3 via the communication network 4 under the control of the measurement data output unit 152. That is, the measurement data output unit 152 performs the processing of the measurement data output step S120 in FIG. 29.

As described above, the measurement program 131 is a program that causes the measurement device 1, which is a computer, to execute each procedure of the flowchart shown in FIG. 29.

As shown in FIG. 37, the monitoring device 3 includes a communication unit 31, a processor 32, a display unit 33, an operation unit 34, and a storage unit 35.

The communication unit 31 receives the measurement data 135 from the measurement device 1 and outputs the received measurement data 135 to the processor 32.

The display unit 33 displays various types of information under the control of the processor 32. The display unit 33 may be, for example, a liquid crystal display or an organic EL display. EL is an abbreviation for electro luminescence.

The operation unit 34 outputs operation data corresponding to an operation of a user to the processor 32. The operation unit 34 may be, for example, an input device such as a mouse, a keyboard, or a microphone.

The storage unit 35 is a memory that stores various programs, data, and the like for the processor 32 to perform calculation processing and control processing. The storage unit 35 stores programs, data, and the like for the processor 32 to implement predetermined application functions.

The processor 32 acquires the measurement data 135 received by the communication unit 31, generates evaluation information by evaluating a temporal change in the displacement of the superstructure 7 based on the acquired measurement data 135, and displays the generated evaluation information on the display unit 33.

In the present embodiment, the processor 32 functions as a measurement data acquisition unit 321 and a monitoring unit 322 by executing a monitoring program 351 stored in the storage unit 35. That is, the processor 32 includes the measurement data acquisition unit 321 and the monitoring unit 322.

The measurement data acquisition unit 321 acquires the measurement data 135 received by the communication unit 31, and adds the acquired measurement data 135 to a measurement data sequence 352 stored in the storage unit 35.

The monitoring unit 322 statistically evaluates a temporal change in the deflection amount of the superstructure 7 based on the measurement data sequence 352 stored in the storage unit 35. Then, the monitoring unit 322 generates evaluation information indicating the evaluation result, and displays the generated evaluation information on the display unit 33. The user can monitor a state of the superstructure 7 based on the evaluation information displayed on the display unit 33.

The monitoring unit 322 may perform processing such as monitoring of the railway vehicle 6 and abnormality determination of the superstructure 7 based on the measurement data sequence 352 stored in the storage unit 35.

The processor 32 transmits, based on the operation data output from the operation unit 34, information for adjusting operation states of the measurement device 1 and the sensor 2 to the measurement device 1 via the communication unit 31. The operation state of the measurement device 1 is adjusted according to the information received via the second communication unit 12. The measurement device 1 transmits information for adjusting the operation state of the sensor 2 received via the second communication unit 12 to the sensor 2 via the first communication unit 11. The operation state of the sensor 2 is adjusted according to the information received via the communication unit 21.

In the processors 14, 23, and 32, for example, the functions of the respective units may be implemented by individual hardware, or the functions of the respective units may be implemented by integrated hardware. For example, the processors 14, 23, and 32 include hardware, and the hardware may include at least one of a circuit that processes a digital signal and a circuit that processes an analog signal. The processors 14, 23, and 32 may be a CPU, a GPU, a DSP, or the like. CPU is an abbreviation for central processing unit, GPU is an abbreviation for graphics processing unit, and DSP is an abbreviation for digital signal processor. The processors 14, 23, and 32 may be configured as custom ICs such as ASICs so as to implement the functions of the respective units, or may implement the functions of the respective units by a CPU and an ASIC. ASIC is an abbreviation for application specific integrated circuit, and IC is an abbreviation for integrated circuit.

The storage units 13, 24, and 35 are configured by, for example, various IC memories such as a ROM, a flash ROM, and a RAM, and a recording medium such as a hard disk, a memory card, and the like. ROM is an abbreviation for read only memory, RAM is an abbreviation for random access memory, and IC is an abbreviation for integrated circuit. The storage units 13, 24, and 35 include a non-volatile information storage device that is a computer-readable device or a medium, and various programs, data, and the like may be stored in the information storage device. The information storage device may be an optical disk such as an optical disk DVD or a CD, a hard disk drive, or various memories such as a card type memory or a ROM.

Although only one sensor 2 is shown in FIG. 37, each of a plurality of sensors 2 may generate the observation data 242 and transmit the observation data 242 to the measurement device 1. In this case, the measurement device 1 receives a plurality of pieces of the observation data 242 transmitted from the plurality of sensors 2, generates a plurality of pieces of measurement data 135, and transmits the plurality of pieces of measurement data 135 to the monitoring device 3. The monitoring device 3 receives the plurality of pieces of measurement data 135 transmitted from the measurement device 1, and monitors a plurality of states of the superstructures 7 based on the plurality of pieces of received measurement data 135.

1-5. Function and Effect

In the measurement method of the first embodiment described above, the measurement device 1 generates the measurement data u(t) based on the acceleration data a(t) output from the sensor 2, and calculates the deflection amount T_std(t) of the superstructure 7 caused by the railway vehicle 6, based on the measurement data u(t) and Equation (35) which is an approximate equation of the deflection based on the structural model reflecting the structure of the superstructure 7 of the bridge 5. The measurement device 1 calculates the static response when the railway vehicle 6 moves on the superstructure 7 by relatively simple processing using the measurement data u(t) and the deflection amount T_std(t). Therefore, according to the measurement method of the first embodiment, the measurement device 1 can calculate the static response by performing processing with a relatively small calculation amount instead of performing processing with a very large calculation amount such as estimating an unknown parameter of a theoretical analysis model based on the acceleration data a(t) by an inverse analysis method.

According to the measurement method of the first embodiment, since the velocity of the railway vehicle 6 actually changes slightly but hardly changes, the measurement device 1 calculates the deflection amount T_std(t) based on the average velocity v_a assuming that the railway vehicle 6 travels at a constant average velocity v_a, and thus it is possible to significantly reduce the calculation amount while maintaining the calculation accuracy of the deflection amount T_std(t).

According to the measurement method of the first embodiment, the measurement device 1 can calculate the average velocity v_a of the railway vehicle 6 by simple calculation according to Equation (13) based on the acceleration data a(t) output from the sensor 2 instead of directly measuring the average velocity v_a of the railway vehicle 6.

According to the measurement method of the first embodiment, by approximating the measurement data u_lp(t), in which the vibration component is reduced by performing filter processing on the measurement data u(t), with the linear function of the deflection amount T_{std_lp}(t) in which the vibration component is reduced by performing filter processing on the deflection amount T_std(t), the measurement device 1 can calculate the static response even when the dynamic response is included in the measurement data u(t) together with the static response.

According to the measurement method of the first embodiment, since the product c₁T_std(t) of the first-order coefficient c₁, which is a first-order term of the linear function for approximating the measurement data u_lp(t), and the deflection amount T_std(t) corresponds to the displacement of the superstructure 7 that is proportional to the load of the railway vehicle 6, and the offset T_{offset_std}(t) corresponds to the displacement of the superstructure 7 that is not proportional to the load of the railway vehicle 6, such as play or floating, the measurement device 1 can accurately calculate the static response by adding the product c₁T_std(t) and the offset T_{offset_std}(t).

According to the measurement method of the first embodiment, the measurement device 1 approximates the measurement data u_lp(t), in which the vibration component having a frequency equal to or higher than the fundamental frequency F_f included in the measurement data u(t) is attenuated, by the linear function of the deflection amount T_{std_lp}(t), and thus the calculation accuracy of the first-order coefficient c₁ and the zero-order coefficient c₀ of the linear function is improved, so that the static response can be accurately calculated.

According to the measurement method of the first embodiment, the measurement device 1 calculates, according to Equation (55), the offset T_{offset_std}(t) which reflects that in an interval where the railway vehicle 6 passes through the superstructure 7, a displacement of the superstructure 7 such as play and floating that are not proportional to the load of the railway vehicle 6 occur, and that the displacement of the superstructure 7 does not occur in other intervals, and thus the static response can be accurately calculated.

According to the measurement method of the first embodiment, since the measurement device 1 can calculate the number of vehicles C_T of the railway vehicle 6 based on the entry time point t_i of the railway vehicle 6 to the superstructure 7 and the exit time point t_o of the railway vehicle 6 from the superstructure 7 according to Equation (8), the static response when the railway vehicle 6 of which the number of vehicles C_T is unknown moves on the superstructure 7 can be accurately calculated.

According to the measurement method of the first embodiment, since the measurement device 1 can accurately calculate the entry time point t_i of the railway vehicle 6 to the superstructure 7 and the exit time point t_o of the railway vehicle 6 from the superstructure 7 based on the measurement data u_lp(t) in which the vibration component is reduced, the static response can be accurately calculated.

2. Second Embodiment

Hereinafter, in a second embodiment, the same components as those in the first embodiment will be denoted by the same reference numerals, repetitive description as those in the first embodiment will be omitted or simplified, and contents different from those in the first embodiment will be mainly described.

In the second embodiment, the measurement device 1 calculates a load waveform PT_std(t) of the railway vehicle 6 traveling on the superstructure 7.

Specifically, first, the measurement device 1 multiplies the deflection amount w_std (a_w (C_m, n), t) of the superstructure 7 caused by each of the plurality of axles of the railway vehicle 6, which is calculated according to Equation (40), by a predetermined conversion ratio U_rn to calculate a load waveform P_std(a_w(C_m,n),t) of each of the plurality of axles, as in Equation (57). FIG. 38 shows an example of the load waveform P_std(a_w(C_m,n),t).

P _std(a_w(C_m,n),t)=U_mW_std(a_w(C_m,n),t)  (57)

Next, the measurement device 1 calculates a load waveform PC_std(C_m,t) of the C_m-th vehicle by adding the load waveform P_std(a_w(C_m,n),t) of each axle of the C_m-th vehicle, as in Equation (58).

$\begin{array}{cc}{\mathrm{PC}}_{\mathrm{std}}\left(C_m,t\right)=\sum _{n=1}^{a_T\left(C_m\right)}{P}_{\mathrm{std}}\left(a_w\left(C_m,n\right),t\right)& \left(58\right)\end{array}$

Finally, the measurement device 1 calculates a load waveform PT_std(t) of the railway vehicle 6 by adding the load waveform PC_std(C_m,t) of each vehicle, as in Equation (59). FIG. 39 shows an example of the load waveform PT_std(t).

$\begin{array}{cc}{\mathrm{PT}}_{\mathrm{std}}\left(t\right)=\sum _{C_m=1}^{C_T}{\mathrm{PC}}_{\mathrm{std}}\left(C_m,t\right)& \left(59\right)\end{array}$

The measurement device 1 may calculate the load waveform PT_std(t) of the railway vehicle 6 by multiplying the deflection amount T_std(t) caused by the railway vehicle 6, which is calculated according to Equation (43), by the conversion ratio U_m as in Equation (60). The load waveform PT_std(t) calculated according to Equation (59) is equivalent to the load waveform PT_std(t) calculated according to Equation (60).

$\begin{array}{cc}{\mathrm{PT}}_{\mathrm{std}}\left(t\right)=U_m{T}_{\mathrm{std}}\left(t\right)=U_m\sum _{C_m=1}^{C_T}{C}_{\mathrm{std}}\left(C_m,t\right)& \left(60\right)\end{array}$

The conversion ratio U_m for converting the deflection amount into the load waveform is calculated as follows prior to the calculation of the load waveform PT_std(t) of the railway vehicle 6.

As in Equation (61), the conversion ratio U_m is calculated as a ratio between the load P of each axle of a known railway vehicle 6′ and a maximum amplitude r_m of the displacement of the superstructure 7 caused by each axle when the known railway vehicle 6′ passes through the superstructure 7.

$\begin{array}{cc}{U}_m=\frac{P}{r_m}& \left(61\right)\end{array}$

For example, the measurement device 1 calculates the deflection amount w_std(a_w(C_m,n),t) of each axle using the acceleration data a(t) output from the sensor 2 when the known railway vehicle 6′ passes through the superstructure 7, according to Equation (40). The measurement device 1 calculates a maximum amplitude max{w_std(a_w(C_m,n),t)} of the deflection amount w_std(a_w(C_m,n),t). FIG. 40 shows a relation between the deflection amount w_std(a_w(C_m,n),t) and the maximum amplitude max{w_std(a_w(C_m,n),t)}.

Further, the measurement device 1 calculates the first-order coefficient c₁ according to Equation (49). Then, the measurement device 1 multiplies the maximum amplitude max{w_std(a_w (C_m, n)} of the deflection amount w_std(a_w(C_m,n),t) by the first-order coefficient c₁ to calculate a maximum amplitude r_m of the displacement of the superstructure 7 when the known railway vehicle 6′ passes through the superstructure 7, as in Equation (62).

r _m =c ₁ max{w_std(a_w(C_m,n),t)}  (62)

The conversion ratio U_m is calculated according to Equation (63) obtained by substituting Equation (62) into Equation (61).

$\begin{array}{cc}{U}_m=\frac{P}{c_1\max \left\{w_{\mathrm{std}}\left(a_w\left(C_m,n\right),t\right)\right\}}& \left(63\right)\end{array}$

For example, when a weight PT of each vehicle of the known railway vehicle 6′ is 30 tons and the number of axles of each vehicle is 4, if it is assumed that the weight PT is equally distributed to each axle, the load P of each axle is 7.5 tons. When the maximum amplitude max{w_std(a_w(C_m,n),t)} of the deflection amount w_std(a_w(C_m,n),t) is −0.245 and the first-order coefficient c₁ is 2.2631, the conversion ratio U_m is calculated as in Equation (64) and is −13.527.

$\begin{array}{cc}{U}_m=\frac{7.5}{2.2631\times \left(-0.245\right)}=-13.527& \left(64\right)\end{array}$

FIG. 41 is a flowchart showing an example of a procedure of a measurement method according to the second embodiment. In FIG. 41, the same reference numerals are given to the steps of performing processing the same as the steps of FIG. 29. In the present embodiment, the measurement device 1 executes the procedure shown in FIG. 41.

As shown in FIG. 41, first, in a conversion ratio calculation step S2, the measurement device 1 calculates the conversion ratio U_m based on the acceleration data a(t) output from the sensor 2 when the known railway vehicle 6′ passes through the superstructure 7.

Next, as in the first embodiment, the measurement device 1 performs the processing of steps S10 to S110.

Next, in a load waveform calculation step S112, the measurement device 1 calculates the load waveform PT_std(t) of the railway vehicle 6 based on the deflection amount w_std(a_w(C_m,n),t) of the superstructure 7 caused by each of the plurality of axles of the railway vehicle 6, which is calculated in the first deflection amount calculation step S60, and the conversion ratio U_m calculated in step S2.

Next, as in the first embodiment, the measurement device 1 performs the processing of the measurement data output step S120. The measurement data output by the measurement device 1 may include the load waveform PT_std(t) calculated in step S112.

Then, the measurement device 1 repeats the processing of steps S10 to S120 until the measurement is completed in step S130.

FIG. 42 is a flowchart showing an example of the procedure of the conversion ratio calculation step S2 in FIG. 41.

As shown in FIG. 42, first, in step S21, the measurement device 1 generates, based on the acceleration data a(t) output from the sensor 2, the measurement data u(t) which is acceleration-based third measurement data which is a response to actions of the plurality of axles of the known railway vehicle 6′ moving on the superstructure 7 on the observation point R. Specifically, the measurement device 1 generates the measurement data u(t) by performing processing the same as that in the first measurement data generation step S20 of FIG. 29, that is, processing the same as that in the steps shown in FIG. 30.

Next, in step S22, the measurement device 1 generates the measurement data u_lp(t), which is fourth measurement data, in which a vibration component is reduced by performing filter processing on the measurement data u(t) generated in step S21. Specifically, the measurement device generates the measurement data u_lp(t) by performing processing the same as that in the second measurement data generation step S30 of FIG. 29, that is, processing the same as that in the steps shown in FIG. 31.

Next, in step S23, the measurement device 1 calculates the deflection amount w_std(a_w(C_m,n),t) of the superstructure 7 caused by each of the plurality of axles of the known railway vehicle 6′ based on the measurement data u_lp(t) generated in step S22, the approximate equation of the deflection of the superstructure 7, which is Equation (35), and the environment information, and adds the deflection amounts w_std(a_w(C_m,n),t) to calculate the deflection amount T_std(t), which is a fourth deflection amount of the superstructure 7 caused by the known railway vehicle 6′. Specifically, the measurement device 1 calculates the deflection amount T_std(t) by performing processing the same as that in the observation information generation step S40, the average velocity calculation step S50, and the first deflection amount calculation step S60 in FIG. 29, that is, processing the same as that in the steps shown in FIG. 32, the steps shown in FIG. 33, and the steps shown in FIG. 34.

Next, in step S24, the measurement device 1 calculates the deflection amount T_{std_lp}(t), which is a fifth deflection amount in which a vibration component is reduced by performing filter processing on the deflection amount T_std(t) calculated in step S23. Specifically, the measurement device 1 generates the deflection amount T_{std_lp}(t) by performing processing the same as that in the second deflection amount calculation step S70 in FIG. 29, that is, processing the same as that in the steps shown in FIG. 35.

Next, in step S25, the measurement device 1 approximates the measurement data u_lp(t) generated in step S22 with the linear function of the deflection amount T_{std_lp}(t) calculated in step S24, and calculates the first-order coefficient c₁ of the linear function. Specifically, the measurement device 1 approximates the measurement data u_lp(t) with the linear function of the deflection amount T_{std_lp}(t) as in Equation (47), and calculates the first-order coefficient c₁ according to Equation (49) using the least-squares method.

Finally, in step S26, the measurement device 1 calculates, as in Equation (63), the conversion ratio U_m by dividing the load P of each of the plurality of axles of the known railway vehicle 6′ by a product c₁ max{w_std(a_w(C_m,n),t)} of the first-order coefficient c₁ calculated in step S25 and the maximum amplitude max{w_std (a_w(C_m, n), t)} of the deflection amount w_std(a_w(C_m,n), t) of the superstructure 7 caused by each of the plurality of axles, which is calculated in step S23.

FIG. 43 is a flowchart showing an example of a procedure of the load waveform calculation step S112 in FIG. 41.

As shown in FIG. 43, in step S1121, the measurement device 1 multiplies the deflection amount w_std(a_w(C_m,n), t) of the superstructure 7 caused by each of the plurality of axles of the known railway vehicle 6′ by the predetermined conversion ratio U_m calculated in the conversion ratio calculation step S2 of FIG. 41 to calculate the load waveform P_std(a_w(C_m,n),t) of each of the plurality of axles, as in Equation (57).

Then, in step S1122, the measurement device 1 calculates the load waveform PT_std(t) of the railway vehicle 6 by adding the load waveform P_std(a_w(C_m,n), t) of each of the plurality of axles, as in Equation (58) and Equation (59).

FIG. 44 is a diagram showing a configuration example of the measurement device 1 according to the second embodiment. As shown in FIG. 44, the measurement device 1 according to the second embodiment includes the first communication unit 11, the second communication unit 12, the storage unit 13, and the processor 14, similarly to the first embodiment. Since the functions of the first communication unit 11, the second communication unit 12, and the storage unit 13 are similar to those in the first embodiment, description thereof will be omitted.

In the present embodiment, the processor 14 functions as the observation data acquisition unit 141, the first measurement data generation unit 142, the second measurement data generation unit 143, the observation information generation unit 144, the average velocity calculation unit 145, the first deflection amount calculation unit 146, the second deflection amount calculation unit 147, the coefficient calculation unit 148, the third deflection amount calculation unit 149, the offset calculation unit 150, the static response calculation unit 151, the measurement data output unit 152, a conversion ratio calculation unit 153, and a load waveform calculation unit 154 by executing the measurement program 131 stored in the storage unit 13. That is, the processor 14 includes the observation data acquisition unit 141, the first measurement data generation unit 142, the second measurement data generation unit 143, the observation information generation unit 144, the average velocity calculation unit 145, the first deflection amount calculation unit 146, the second deflection amount calculation unit 147, the coefficient calculation unit 148, the third deflection amount calculation unit 149, the offset calculation unit 150, the static response calculation unit 151, the measurement data output unit 152, the conversion ratio calculation unit 153, and the load waveform calculation unit 154.

Since the functions of the observation data acquisition unit 141, the first measurement data generation unit 142, the second measurement data generation unit 143, the observation information generation unit 144, the average velocity calculation unit 145, the first deflection amount calculation unit 146, the second deflection amount calculation unit 147, the coefficient calculation unit 148, the third deflection amount calculation unit 149, the offset calculation unit 150, the static response calculation unit 151, and the measurement data output unit 152 are similar to those in the first embodiment, description thereof will be omitted. The observation data acquisition unit 141 performs the processing of the observation data acquisition step S10 of FIG. 41. The first measurement data generation unit 142 performs the processing of the first measurement data generation step S20 of FIG. 41. The second measurement data generation unit 143 performs the processing of the second measurement data generation step S30 of FIG. 41. The observation information generation unit 144 performs the processing of the observation information generation step S40 of FIG. 41. The average velocity calculation unit 145 performs the processing of the average velocity calculation step S50 of FIG. 41. The first deflection amount calculation unit 146 performs the processing of the first deflection amount calculation step S60 of FIG. 41. The second deflection amount calculation unit 147 performs the processing of the second deflection amount calculation step S70 of FIG. 41. The coefficient calculation unit 148 performs the processing of the coefficient calculation step S80 of FIG. 41. The third deflection amount calculation unit 149 performs the processing of the third deflection amount calculation step S90 of FIG. 41. The offset calculation unit 150 performs the processing of the offset calculation step S100 of FIG. 41. The static response calculation unit 151 performs the processing of the static response calculation step S110 of FIG. 41. The measurement data output unit 152 performs the processing of the measurement data output step S120 of FIG. 41.

The conversion ratio calculation unit 153 calculates the conversion ratio U_m based on the acceleration data a(t) output from the sensor 2 when the known railway vehicle 6′ passes through the superstructure 7. The acceleration data a(t) is acquired by the observation data acquisition unit 141 and is stored in the storage unit 13 as the observation data 133, and the conversion ratio calculation unit 153 reads the observation data 133 from the storage unit 13 and calculates the conversion ratio U_m.

Specifically, first, the conversion ratio calculation unit 153 generates, based on the observation data 133, the measurement data u(t), which is the acceleration-based third measurement data, which is the response to the actions of the plurality of axles of the known railway vehicle 6′ moving on the superstructure 7 on the observation point R. Next, the conversion ratio calculation unit 153 generates the measurement data u_lp(t), which is the fourth measurement data in which the vibration component is reduced by performing filter processing on the measurement data u(t). Next, the conversion ratio calculation unit 153 calculates the deflection amount w_std(a_w(C_m,n),t) of the superstructure 7 caused by each of the plurality of axles of the known railway vehicle 6′ based on the measurement data u_lp(t), the approximate equation of the deflection of the superstructure 7, which is Equation (35), and the environment information 132 stored in the storage unit 13, and adds the deflection amounts w_std(a_w(C_m,n),t) to calculate the deflection amount T_std(t), which is the fourth deflection amount of the superstructure caused by the known railway vehicle 6′. Next, the conversion ratio calculation unit 153 calculates the deflection amount T_{std_lp}(t), which is the fifth deflection amount in which the vibration component is reduced by performing filter processing on the deflection amount T_std (t). Next, the conversion ratio calculation unit 153 approximates the measurement data u_lp(t) with the linear function of the deflection amount T_{std_lp}(t), and calculates the first-order coefficient c₁ of the linear function. Finally, the conversion ratio calculation unit 153 calculates the conversion ratio U_m by dividing the load P of each of the plurality of axles of the known railway vehicle 6′ by the product c₁ max{w_std (a_w (C_m, n), t)} of the first-order coefficient c₁ and the maximum amplitude max{w_std (a_w (C_m, n), t)} of the deflection amount w_std(a_w(C_m,n),t), as in Equation (63). Then, the conversion ratio calculation unit 153 stores the calculated conversion ratio U_m in the storage unit 13 as a conversion ratio 136.

As described above, the conversion ratio calculation unit 153 performs the processing of the conversion ratio calculation step S2 in FIG. 41, specifically, the processing of steps S21 to S26 in FIG. 42.

The load waveform calculation unit 154 performs the processing of calculating the load waveform PT_std(t) of the railway vehicle 6 based on the deflection amount w_std(a_w(C_m,n),t) of the superstructure 7 caused by each of the plurality of axles of the railway vehicle 6, which is calculated by the first deflection amount calculation unit 146, and the conversion ratio 136 stored in the storage unit 13. Specifically, the load waveform calculation unit 154 multiplies the deflection amount w_std(a_w(C_m,n),t) of the superstructure 7 caused by each of the plurality of axles of the known railway vehicle 6′ by the conversion ratio 136, that is, the conversion ratio U_m, to calculate the load waveform P_std(a_w(C_m,n),t) of each of the plurality of axles, as in Equation (57). Then, the load waveform calculation unit 154 calculates the load waveform PT_std(t) of the railway vehicle 6 by adding the load waveform P_std(a_w(C_m,n),t) of each of the plurality of axles, as in Equation (58) and Equation (59). The load waveform PT_std(t) calculated by the load waveform calculation unit 154 may be stored in the storage unit 13 as at least a part of the measurement data 135.

As described above, the load waveform calculation unit 154 performs the processing of the load waveform calculation step S112 in FIG. 41, specifically, the processing of steps S1121 and S1122 in FIG. 43.

As described above, the measurement program 131 is a program that causes the measurement device 1, which is a computer, to execute each procedure of the flowchart shown in FIG. 41.

The calculation of the conversion ratio U_m may be performed by a device different from the measurement device 1. In this case, the processor 14 may not function as the conversion ratio calculation unit 153.

In the measurement method according to the second embodiment described above, the first-order term of the linear function of the deflection amount T_std(t) for approximating the measurement data u_lp(t) corresponds to the displacement of the superstructure 7 that is proportional to the load of the known railway vehicle 6′. Therefore, according to the measurement method of the second embodiment, the measurement device 1 can accurately calculate the conversion ratio U_m between the displacement and the load by dividing the load of each of the plurality of axles of the known railway vehicle 6′ by the product c₁ max{w_std(a_w(C_m,n),t)} of the first-order coefficient c₁ of the linear function and the maximum amplitude max{w_std (a_w(C_m,t)} of the deflection amount w_std(a_w(C_m,n),t) of the superstructure 7 caused by each of the plurality of axles.

According to the measurement method of the second embodiment, the measurement device 1 can use the conversion ratio U_m with high accuracy to accurately calculate the load waveform PT_std(t) of the railway vehicle 6 by processing with a relatively small calculation amount.

According to the measurement method of the second embodiment, it is possible to achieve the same effects as those of the measurement method according to the first embodiment.

3. Modification

The present disclosure is not limited to the above embodiments, and various modifications can be made within the scope of the gist of the present disclosure.

In the embodiments described above, the sensor 2, which is an observation device, is an acceleration sensor that outputs the acceleration data a(k), but the observation device is not limited to the acceleration sensor. For example, the observation device may be an impact sensor, a pressure-sensitive sensor, a strain gauge, an image measuring device, a load cell, or a displacement meter.

The impact sensor detects an impact acceleration as a response to an action of each axle of the railway vehicle 6 on the observation point R. The pressure-sensitive sensor, the strain gauge, and the load cell detect a stress change as a response to an action of each axle of the railway vehicle 6 on the observation point R. The image measuring device detects, by image processing, a displacement as a response to an action of each axle of the railway vehicle 6 on the observation point R. The displacement gauge is, for example, a contact-type displacement meter, a ring-type displacement meter, a laser displacement meter, a pressure-sensitive sensor, or a displacement measurement device using an optical fiber, and detects a displacement as a response to an action of each axle of the railway vehicle 6 on the observation point R.

As an example, FIG. 45 shows a configuration example of the measurement system 10 using a ring-type displacement meter as the observation device. FIG. 46 shows a configuration example of the measurement system 10 using an image measuring device as the observation device. In FIGS. 45 and 46, the same components as those in FIG. 1 are denoted by the same reference numerals, and description thereof will be omitted. In the measurement system 10 shown in FIG. 45, a piano wire 41 is fixed between an upper surface of a ring-type displacement meter 40 and a lower surface of the main girder G immediately above the ring-type displacement meter 40, and the ring-type displacement meter 40 measures a displacement of the piano wire 41 caused by bending of the superstructure 7 and transmits the measured displacement data to the measurement device 1. The measurement device 1 generates the measurement data 135 based on the displacement data transmitted from the ring-type displacement meter 40. In the measurement system 10 shown in FIG. 46, a camera 50 transmits, to the measurement device 1, an image obtained by imaging a target 51 provided on a side surface of the main girder G. The measurement device 1 processes the image transmitted from the camera 50, calculates a displacement of the target 51 caused by bending of the superstructure 7 to generate displacement data, and generates the measurement data 135 based on the generated displacement data. In the example of FIG. 46, the measurement device 1 generates the displacement data as an image measuring device, but the displacement data may be generated by an image measuring device (not shown) different from the measurement device 1 by image processing.

In the embodiments described above, the bridge 5 is a railway bridge, and the moving object moving on the bridge 5 is the railway vehicle 6, but the bridge 5 may be a road bridge, and the moving object moving on the bridge 5 may be a vehicle such as an automobile, a road train, a truck, or a construction vehicle. FIG. 47 shows a configuration example of the measurement system 10 when the bridge 5 is a road bridge and a vehicle 6a moves on the bridge 5. In FIG. 47, the same components as those in FIG. 1 are denoted by the same reference numerals. As shown in FIG. 47, the bridge 5, which is a road bridge, includes the superstructure 7 and the substructure 8, similarly to the railway bridge. FIG. 48 is a cross-sectional view of the superstructure 7 taken along line A-A of FIG. 47. As shown in FIGS. 47 and 48, the superstructure 7 includes the bridge floor 7a and the support 7b, and the bridge floor 7a includes the floor plate F, the main girder G, and a cross girder (not shown). As shown in FIG. 47, the substructure 8 includes bridge piers 8a and bridge abutments 8b. The superstructure 7 is a structure across any one of the bridge abutment 8b and the bridge pier 8a adjacent to each other, two adjacent bridge abutments 8b, and two adjacent bridge piers 8a. Both end portions of the superstructure 7 are located at positions of the bridge abutment 8b and the bridge pier 8a adjacent to each other, at positions of the two adjacent bridge abutments 8b, or at positions of the two adjacent bridge piers 8a. The bridge 5 is, for example, a steel bridge, a girder bridge, or an RC bridge.

Each sensor 2 is installed at a central portion of the superstructure 7 in a longitudinal direction, specifically, at a central portion of the main girder G in the longitudinal direction. However, each sensor 2 is not limited to being installed at the central portion of the superstructure 7 as long as each sensor 2 can detect an acceleration for calculating the displacement of the superstructure 7. When each sensor 2 is provided on the floor plate F of the superstructure 7, the sensor 2 may be damaged due to traveling of the vehicle 6a, and the measurement accuracy may be affected by local deformation of the bridge floor 7a, so that in the example of FIGS. 47 and 48, each sensor 2 is provided at the main girder G of the superstructure 7.

As shown in FIG. 48, the superstructure 7 has two lanes L1 and L2 on which the vehicle 6a as a moving object can move and three main girders G. In the example of FIGS. 47 and 48, in the central portion of the superstructure 7 in the longitudinal direction, the sensors 2 are respectively provided at two main girders at two ends, an observation point R1 is provided at a position of a surface of the lane L1 vertically above one of the sensors 2, and an observation point R2 is provided at a position of a surface of the lane L2 vertically above the other of the sensors 2. That is, the two sensors 2 are observation devices for observing the observation points R1 and R2, respectively. The two sensors 2 for respectively observing the observation points R1 and R2 may be provided at positions where accelerations generated at the observation points R1 and R2 due to the traveling of the vehicle 6a can be detected, and are preferably provided at positions close to the observation points R1 and R2. The number and installation position of the sensor 2, and the number of lanes are not limited to the example shown in FIGS. 47 and 48, and various modifications can be made.

The measurement device 1 calculates displacements of bending of the lanes L1 and L2 caused by the traveling of the vehicle 6a based on the acceleration data output from the sensors 2, and transmits information on the displacements of the lanes L1 and L2 to the monitoring device 3 via the communication network 4. The monitoring device 3 may store the information in a storage device (not shown), and may perform processing such as monitoring of the vehicle 6a and abnormality determination of the superstructure 7 based on the information, for example.

In the embodiments described above, each sensor 2 is provided at the main girder G of the superstructure 7, but the sensor 2 may be provided on the surface of or inside the superstructure 7, at the lower surface of the floor plate F, at the bridge pier 8a, or the like. In the embodiments described above, the superstructure of the bridge is described as an example of the structure, but the present disclosure is not limited thereto, and any structure may be used as long as the structure is deformed due to the movement of the moving object.

In the embodiments described above, the measurement device 1 calculates the entry time point t_i based on the observation data output from the observation device that observes the observation point R, but the measurement device 1 may calculate the entry time point t_i based on observation data output from another observation device that observes the entry end of the superstructure 7. Similarly, in the embodiments described above, the measurement device 1 calculates the exit time point t_o based on the observation data output from the observation device that observes the observation point R, but the measurement device 1 may calculate the exit time point t_o based on observation data output from another observation device that observes the exit end of the superstructure 7.

The embodiments and the modifications described above are merely examples, and the present disclosure is not limited thereto. For example, the embodiments and modifications can be appropriately combined.

The present disclosure includes a configuration substantially the same as the configuration described in the embodiment, for example, a configuration having the same function, method, and result, or a configuration having the same object and effect. The present disclosure includes a configuration in which a non-essential portion of the configuration described in the embodiment is replaced. The present disclosure includes a configuration having the same function and effect as the configuration described in the embodiment, or a configuration capable of achieving the same object. Further, the present disclosure includes a configuration in which a known technique is added to the configuration described in the embodiment.

The following contents are derived from the embodiments and modifications described above.

A measurement method according to an aspect includes: a first measurement data generation step of generating, based on observation data output from an observation device configured to observe an observation point of a structure, first measurement data based on a physical quantity which is a response to actions of a plurality of parts of a moving object moving on the structure on the observation point; a second measurement data generation step of generating second measurement data in which a vibration component is reduced by performing filter processing on the first measurement data; an observation information generation step of generating observation information including an entry time point and an exit time point of the moving object with respect to the structure; an average velocity calculation step of calculating an average velocity of the moving object based on the observation information and environment information which is created in advance and includes a dimension of the moving object and a dimension of the structure; a first deflection amount calculation step of calculating, based on an approximate equation of deflection of the structure, the observation information, the environment information, and the average velocity, a first deflection amount of the structure caused by the moving object; a second deflection amount calculation step of calculating a second deflection amount in which a vibration component is reduced by performing filter processing on the first deflection amount; a coefficient calculation step of approximating the second measurement data with a linear function of the second deflection amount to calculate a first-order coefficient and a zero-order coefficient of the linear function; a third deflection amount calculation step of calculating a third deflection amount based on the first-order coefficient, the zero-order coefficient, and the second deflection amount; an offset calculation step of calculating an offset based on the zero-order coefficient, the second deflection amount, and the third deflection amount; and a static response calculation step of calculating a static response by adding the offset and a product of the first-order coefficient and the first deflection amount.

In the measurement method, the static response when the moving object moves on the structure is calculated by relatively simple processing using the first measurement data generated based on the observation data and the first deflection amount generated based on the approximate equation of the deflection of the structure. Therefore, according to the measurement method, it is possible to calculate the static response by performing processing with a relatively small calculation amount instead of performing processing with a very large calculation amount such as estimating an unknown parameter of a theoretical analysis model based on the acceleration data by an inverse analysis method.

According to the measurement method, since the velocity of the moving object actually changes slightly but hardly changes, it is possible to calculate the first deflection amount based on the average velocity assuming that the moving object moves at a constant average velocity, and thus it is possible to significantly reduce the calculation amount while maintaining the calculation accuracy of the first deflection amount.

According to the measurement method, the second measurement data in which the vibration component is reduced by performing filter processing on the first measurement data is approximated with the linear function of the second deflection amount in which the vibration component is reduced by performing filter processing on the first deflection amount, and thus the static response can be calculated even when the dynamic response is included in the first measurement data together with the static response.

According to the measurement method, since the product of the first-order coefficient, which is the first-order term of the linear function approximating the first deflection amount, and the first deflection amount corresponds to the displacement of the structure that is proportional to the load of the moving object, and the offset corresponds to the displacement of the structure that is not proportional to the load of the moving object, such as play or floating, it is possible to accurately calculate the static response by adding the offset and the product of the first-order coefficient and the first deflection amount.

In the measurement method according to the above aspect, in the offset calculation step, a product of an amplitude ratio of the third deflection amount and the second deflection amount in a predetermined interval and the second deflection amount may be calculated, and the offset may be calculated by replacing, with the zero-order coefficient, an interval of the product in which an absolute value of the product of the amplitude ratio and the second deflection amount is bigger than an absolute value of the zero-order coefficient.

According to the measurement method, the offset which reflects that in an interval where the moving object moves on the structure, the displacement of the structure such as play and floating that are not proportional to the load of the moving object occur, and that the displacement of the structure does not occur in other intervals is calculated, and thus the static response can be accurately calculated.

In the measurement method according to the above aspect, the entry time point may be a time point at which a leading part of the plurality of parts of the moving object passes through an entry end of the structure, the exit time point may be a time point at which a rearmost part of the plurality of parts of the moving object passes through an exit end of the structure, and in the average velocity calculation step, a distance from the leading part to the rearmost part and a distance from the entry end to the exit end may be calculated based on the environment information, and the average velocity may be calculated based on the entry time point, the exit time point, the distance from the leading part to the rearmost part, and the distance from the entry end to the exit end.

According to the measurement method, the average velocity of the moving object can be calculated by simple calculation based on observation data output from an observation device that observes the observation point instead of directly measuring the average velocity of the moving object.

In the measurement method according to the above aspect, the environment information may include a length of the structure, a position of the observation point, a length of each vehicle of the moving object, and a position of each of the plurality of parts of the moving object.

In the measurement method according to the above aspect, in the second measurement data generation step, low-pass filter processing for attenuating the vibration component having a frequency equal to or higher than a fundamental frequency of the first measurement data may be performed as the filter processing.

According to the measurement method, in the coefficient calculation step, since the second measurement data in which the vibration component having a frequency equal to or higher than the fundamental frequency included in the first measurement data is attenuated is approximated by the linear function of the second deflection amount, the calculation accuracy of the first-order coefficient and the zero-order coefficient of the linear function is improved, and thus the static response can be accurately calculated.

In the measurement method according to the above aspect, the observation information may include the number of vehicles of the moving object, and in the observation information generation step, a maximum integer less than or equal to a number obtained by subtracting 1 from a product of the fundamental frequency and a passing time, which is a difference between the exit time point and the entry time point, may be calculated as the number of vehicles.

According to the measurement method, since the number of vehicles of the moving object can be calculated based on the entry time point of the moving object to the structure and the exit time point of the moving object from the structure, it is possible to accurately calculate the static response when the moving object of which the number of vehicles is unknown moves on the structure.

In the measurement method according to the above aspect, in the observation information generation step, a first time point at which an amplitude of the second measurement data matches or exceeds a threshold may be calculated as the entry time point, and a second time point after the first time point at which the amplitude of the second measurement data matches or exceeds the threshold may be calculated as the exit time point.

According to the measurement method, since the entry time point of the moving object to the structure and the exit time point of the moving object from the structure can be accurately calculated based on the second measurement data in which the vibration component is reduced, it is possible to accurately calculate the static response.

In the measurement method according to the above aspect, the structure may be a superstructure of a bridge.

According to the measurement method, it is possible to calculate the static response when the moving object moves on the superstructure of the bridge by processing with a relatively small calculation amount.

In the measurement method according to the above aspect, the moving object may be a vehicle or a railway vehicle, and each of the plurality of parts may be an axle or a wheel.

According to the measurement method, it is possible to calculate the static response when the vehicle or the railway vehicle moves on the structure by processing with a relatively small calculation amount.

In the measurement method according to the above aspect, the approximate equation of the deflection of the structure may be an equation based on a structural model of the structure.

According to the measurement method, it is possible to calculate the first deflection amount reflecting the configuration of the structure on which the moving object moves, and it is possible to accurately calculate the static response.

In the measurement method according to the above aspect, the structural model may be a simple beam whose both ends are supported.

According to the measurement method, it is possible to accurately calculate the static response when the moving object moves on a structure having a configuration similar to the simple beam.

The measurement method according to the above aspect may further include a load waveform calculation step of calculating a load waveform of the moving object, in the first deflection amount calculation step, deflection amounts of the structure respectively caused by the plurality of parts may be calculated based on the approximate equation of the deflection of the structure, the observation information, the environment information, and the average velocity, and the first deflection amount may be calculated by adding the deflection amounts of the structure respectively caused by the plurality of parts, and in the load waveform calculation step, load waveforms of the plurality of parts may be respectively calculated by multiplying the deflection amounts of the structure respectively caused by the plurality of parts by a predetermined conversion ratio, and a load waveform of the moving object may be calculated by adding the load waveforms of the plurality of parts.

According to the measurement method, it is possible to accurately calculate the load waveform of the moving object by processing with a relatively small calculation amount.

The measurement method according to the above aspect may further include a conversion ratio calculation step of calculating the conversion ratio, and the conversion ratio calculation step may include: generating, based on the observation data output from the observation device, third measurement data based on a physical quantity which is a response to actions of a plurality of parts of a known moving object moving on the structure on the observation point, generating fourth measurement data in which a vibration component is reduced by performing filter processing on the third measurement data, calculating, based on the fourth measurement data, the approximate equation of the deflection of the structure, and the environment information, the deflection amounts of the structure respectively caused by the plurality of parts of the known moving object, and calculating a fourth deflection amount of the structure caused by the known moving object by adding the deflection amounts of the structure respectively caused by the plurality of parts, calculating a fifth deflection amount in which a vibration component is reduced by performing filter processing on the fourth deflection amount, approximating the fourth measurement data with a linear function of the fifth deflection amount to calculate a first-order coefficient of the linear function, and calculating the conversion ratio by dividing a load of each of the plurality of parts of the known moving object by a product of the first-order coefficient and a maximum amplitude of the deflection amount of the structure caused by each of the plurality of parts.

According to the measurement method, since the first-order term of the linear function approximating the fourth measurement data corresponds to the displacement of the structure proportional to the load of the known moving object, it is possible to accurately calculate the conversion ratio between the displacement and the load by dividing the load of each of the plurality of parts of the known moving object by the product of the first-order coefficient of the linear function and the maximum amplitude of the deflection amount of the structure caused by each of the plurality of parts.

In the measurement method according to the above aspect, the observation device may be an acceleration sensor, an impact sensor, a pressure-sensitive sensor, a strain gauge, an image measuring device, a load cell, or a displacement meter.

According to the measurement method, it is possible to accurately measure the static response using data of acceleration, a stress change, or a displacement.

In the measurement method according to the above aspect, the structure may be a structure in which bridge weigh in motion (BWIM) functions.

A measurement device according to an aspect includes: a first measurement data generation unit configured to generate, based on observation data output from an observation device configured to observe an observation point of a structure, first measurement data based on a physical quantity which is a response to actions of a plurality of parts of a moving object moving on the structure on the observation point; a second measurement data generation unit configured to generate second measurement data in which a vibration component is reduced by performing filter processing on the first measurement data; an observation information generation unit configured to generate observation information including an entry time point and an exit time point of the moving object with respect to the structure; an average velocity calculation unit configured to calculate an average velocity of the moving object based on the observation information and environment information which is created in advance and includes a dimension of the moving object and a dimension of the structure; a first deflection amount calculation unit configured to calculate, based on an approximate equation of deflection of the structure, the observation information, the environment information, and the average velocity, a first deflection amount of the structure caused by the moving object; a second deflection amount calculation unit configured to calculate a second deflection amount in which a vibration component is reduced by performing filter processing on the first deflection amount; a coefficient calculation unit configured to approximate the second measurement data with a linear function of the second deflection amount to calculate a first-order coefficient and a zero-order coefficient of the linear function; a third deflection amount calculation unit configured to calculate a third deflection amount based on the first-order coefficient, the zero-order coefficient, and the second deflection amount; an offset calculation unit configured to calculate an offset based on the zero-order coefficient, the second deflection amount, and the third deflection amount; and a static response calculation unit configured to calculate a static response by adding the offset and a product of the first-order coefficient and the first deflection amount.

The measurement device calculates the static response when the moving object moves on the structure by relatively simple processing using the first measurement data generated based on the observation data and the first deflection amount generated based on the approximate equation of the deflection of the structure. Therefore, according to the measurement device, it is possible to calculate the static response by performing processing with a relatively small calculation amount instead of performing processing with a very large calculation amount such as estimating an unknown parameter of a theoretical analysis model based on the acceleration data by an inverse analysis method.

According to the measurement device, since the velocity of the moving object actually changes slightly but hardly changes, it is possible to calculate the first deflection amount based on the average velocity assuming that the moving object moves at a constant average velocity, and thus it is possible to significantly reduce the calculation amount while maintaining the calculation accuracy of the first deflection amount.

According to the measurement device, the second measurement data in which the vibration component is reduced by performing filter processing on the first measurement data is approximated with the linear function of the second deflection amount in which the vibration component is reduced by performing filter processing on the first deflection amount, and thus the static response can be calculated even when the dynamic response is included in the first measurement data together with the static response.

According to the measurement device, since the product of the first-order coefficient, which is the first-order term of the linear function approximating the first deflection amount, and the first deflection amount corresponds to the displacement of the structure that is proportional to the load of the moving object, and the offset corresponds to the displacement of the structure that is not proportional to the load of the moving object, such as play or floating, it is possible to accurately calculate the static response by adding the offset and the product of the first-order coefficient and the first deflection amount.

A measurement system according to an aspect includes: the measurement device according to the above aspect; and the observation device.

A non-transitory computer-readable storage medium according to an aspect of the present disclosure stores a measurement program, and the measurement program causes a computer to execute: a first measurement data generation step of generating, based on observation data output from an observation device configured to observe an observation point of a structure, first measurement data based on a physical quantity which is a response to actions of a plurality of parts of a moving object moving on the structure on the observation point; a second measurement data generation step of generating second measurement data in which a vibration component is reduced by performing filter processing on the first measurement data; an observation information generation step of generating observation information including an entry time point and an exit time point of the moving object with respect to the structure; an average velocity calculation step of calculating an average velocity of the moving object based on the observation information and environment information which is created in advance and includes a dimension of the moving object and a dimension of the structure; a first deflection amount calculation step of calculating, based on an approximate equation of deflection of the structure, the observation information, the environment information, and the average velocity, a first deflection amount of the structure caused by the moving object; a second deflection amount calculation step of calculating a second deflection amount in which a vibration component is reduced by performing filter processing on the first deflection amount; a coefficient calculation step of approximating the second measurement data with a linear function of the second deflection amount to calculate a first-order coefficient and a zero-order coefficient of the linear function; a third deflection amount calculation step of calculating a third deflection amount based on the first-order coefficient, the zero-order coefficient, and the second deflection amount; an offset calculation step of calculating an offset based on the zero-order coefficient, the second deflection amount, and the third deflection amount; and a static response calculation step of calculating a static response by adding the offset and a product of the first-order coefficient and the first deflection amount.

In the measurement program, the static response when the moving object moves on the structure is calculated by relatively simple processing using the first measurement data generated based on the observation data and the first deflection amount generated based on the approximate equation of the deflection of the structure. Therefore, according to the measurement program, it is possible to calculate the static response by performing processing with a relatively small calculation amount instead of performing processing with a very large calculation amount such as estimating an unknown parameter of a theoretical analysis model based on the acceleration data by an inverse analysis method.

According to the measurement program, since the velocity of the moving object actually changes slightly but hardly changes, it is possible to calculate the first deflection amount based on the average velocity assuming that the moving object moves at a constant average velocity, and thus it is possible to significantly reduce the calculation amount while maintaining the calculation accuracy of the first deflection amount.

According to the measurement program, the second measurement data in which the vibration component is reduced by performing filter processing on the first measurement data is approximated with the linear function of the second deflection amount in which the vibration component is reduced by performing filter processing on the first deflection amount, and thus the static response can be calculated even when the dynamic response is included in the first measurement data together with the static response.

According to the measurement program, since the product of the first-order coefficient, which is the first-order term of the linear function approximating the first deflection amount, and the first deflection amount corresponds to the displacement of the structure that is proportional to the load of the moving object, and the offset corresponds to the displacement of the structure that is not proportional to the load of the moving object, such as play or floating, it is possible to accurately calculate the static response by adding the offset and the product of the first-order coefficient and the first deflection amount.

Claims

1. A measurement method, comprising: a first measurement data generation step of generating, based on observation data output from an observation device configured to observe an observation point of a structure, first measurement data based on a physical quantity which is a response to actions of a plurality of parts of a moving object moving on the structure on the observation point;a second measurement data generation step of generating second measurement data in which a vibration component is reduced by performing filter processing on the first measurement data;an observation information generation step of generating observation information including an entry time point and an exit time point of the moving object with respect to the structure;an average velocity calculation step of calculating an average velocity of the moving object based on the observation information and environment information which is created in advance and includes a dimension of the moving object and a dimension of the structure;a first deflection amount calculation step of calculating, based on an approximate equation of deflection of the structure, the observation information, the environment information, and the average velocity, a first deflection amount of the structure caused by the moving object;a second deflection amount calculation step of calculating a second deflection amount in which a vibration component is reduced by performing filter processing on the first deflection amount;a coefficient calculation step of approximating the second measurement data with a linear function of the second deflection amount to calculate a first-order coefficient and a zero-order coefficient of the linear function;a third deflection amount calculation step of calculating a third deflection amount based on the first-order coefficient, the zero-order coefficient, and the second deflection amount;an offset calculation step of calculating an offset based on the zero-order coefficient, the second deflection amount, and the third deflection amount; anda static response calculation step of calculating a static response by adding the offset and a product of the first-order coefficient and the first deflection amount.

2. The measurement method according to claim 1, wherein in the offset calculation step,a product of an amplitude ratio of the third deflection amount and the second deflection amount in a predetermined interval and the second deflection amount is calculated, and the offset is calculated by replacing, with the zero-order coefficient, an interval of the product in which an absolute value of the product of the amplitude ratio and the second deflection amount is bigger than an absolute value of the zero-order coefficient.

3. The measurement method according to claim 1, wherein the entry time point is a time point at which a leading part of the plurality of parts of the moving object passes through an entry end of the structure,the exit time point is a time point at which a rearmost part of the plurality of parts of the moving object passes through an exit end of the structure, andin the average velocity calculation step,a distance from the leading part to the rearmost part and a distance from the entry end to the exit end are calculated based on the environment information, andthe average velocity is calculated based on the entry time point, the exit time point, the distance from the leading part to the rearmost part, and the distance from the entry end to the exit end.

4. The measurement method according to claim 1, wherein the environment information includes a length of the structure, a position of the observation point, a length of each vehicle of the moving object, and a position of each of the plurality of parts of the moving object.

5. The measurement method according to claim 1, wherein in the second measurement data generation step,low-pass filter processing for attenuating the vibration component having a frequency equal to or higher than a fundamental frequency of the first measurement data is performed as the filter processing.

6. The measurement method according to claim 5, wherein the observation information includes the number of vehicles of the moving object, andin the observation information generation step,a maximum integer less than or equal to a number obtained by subtracting 1 from a product of the fundamental frequency and a passing time, which is a difference between the exit time point and the entry time point, is calculated as the number of vehicles.

7. The measurement method according to claim 1, wherein in the observation information generation step,a first time point at which an amplitude of the second measurement data matches or exceeds a threshold is calculated as the entry time point, anda second time point after the first time point at which the amplitude of the second measurement data matches or exceeds the threshold is calculated as the exit time point.

8. The measurement method according to claim 1, wherein the structure is a superstructure of a bridge.

9. The measurement method according to claim 1, wherein the moving object is a vehicle or a railway vehicle, andeach of the plurality of parts is an axle or a wheel.

10. The measurement method according to claim 1, wherein the approximate equation of the deflection of the structure is an equation based on a structural model of the structure.

11. The measurement method according to claim 10, wherein the structural model is a simple beam whose both ends are supported.

12. The measurement method according to claim 1, further comprising: a load waveform calculation step of calculating a load waveform of the moving object, whereinin the first deflection amount calculation step,deflection amounts of the structure respectively caused by the plurality of parts are calculated based on the approximate equation of the deflection of the structure, the observation information, the environment information, and the average velocity, and the first deflection amount is calculated by adding the deflection amounts of the structure respectively caused by the plurality of parts, andin the load waveform calculation step,load waveforms of the plurality of parts are respectively calculated by multiplying the deflection amounts of the structure respectively caused by the plurality of parts by a predetermined conversion ratio, and a load waveform of the moving object is calculated by adding the load waveforms of the plurality of parts.

13. The measurement method according to claim 12, further comprising: a conversion ratio calculation step of calculating the conversion ratio, whereinthe conversion ratio calculation step includes:generating, based on the observation data output from the observation device, third measurement data based on a physical quantity which is a response to actions of a plurality of parts of a known moving object moving on the structure on the observation point,generating fourth measurement data in which a vibration component is reduced by performing filter processing on the third measurement data,calculating, based on the fourth measurement data, the approximate equation of the deflection of the structure, and the environment information, the deflection amounts of the structure respectively caused by the plurality of parts of the known moving object, and calculating a fourth deflection amount of the structure caused by the known moving object by adding the deflection amounts of the structure respectively caused by the plurality of parts,calculating a fifth deflection amount in which a vibration component is reduced by performing filter processing on the fourth deflection amount,approximating the fourth measurement data with a linear function of the fifth deflection amount to calculate a first-order coefficient of the linear function, andcalculating the conversion ratio by dividing a load of each of the plurality of parts of the known moving object by a product of the first-order coefficient and a maximum amplitude of the deflection amount of the structure caused by each of the plurality of parts.

14. The measurement method according to claim 1, wherein the observation device is an acceleration sensor, an impact sensor, a pressure-sensitive sensor, a strain gauge, an image measuring device, a load cell, or a displacement meter.

15. The measurement method according to claim 1, wherein the structure is a structure in which bridge weigh in motion (BWIM) functions.

16. A measurement device, comprising: a first measurement data generation unit configured to generate, based on observation data output from an observation device configured to observe an observation point of a structure, first measurement data based on a physical quantity which is a response to actions of a plurality of parts of a moving object moving on the structure on the observation point;a second measurement data generation unit configured to generate second measurement data in which a vibration component is reduced by performing filter processing on the first measurement data;an observation information generation unit configured to generate observation information including an entry time point and an exit time point of the moving object with respect to the structure;an average velocity calculation unit configured to calculate an average velocity of the moving object based on the observation information and environment information which is created in advance and includes a dimension of the moving object and a dimension of the structure;a first deflection amount calculation unit configured to calculate, based on an approximate equation of deflection of the structure, the observation information, the environment information, and the average velocity, a first deflection amount of the structure caused by the moving object;a second deflection amount calculation unit configured to calculate a second deflection amount in which a vibration component is reduced by performing filter processing on the first deflection amount;a coefficient calculation unit configured to approximate the second measurement data with a linear function of the second deflection amount to calculate a first-order coefficient and a zero-order coefficient of the linear function;a third deflection amount calculation unit configured to calculate a third deflection amount based on the first-order coefficient, the zero-order coefficient, and the second deflection amount;an offset calculation unit configured to calculate an offset based on the zero-order coefficient, the second deflection amount, and the third deflection amount; anda static response calculation unit configured to calculate a static response by adding the offset and a product of the first-order coefficient and the first deflection amount.

17. A measurement system, comprising: the measurement device according to claim 16; andthe observation device.

18. A non-transitory computer-readable storage medium storing a measurement program, the measurement program causing a computer to execute: a first measurement data generation step of generating, based on observation data output from an observation device configured to observe an observation point of a structure, first measurement data based on a physical quantity which is a response to actions of a plurality of parts of a moving object moving on the structure on the observation point;a second measurement data generation step of generating second measurement data in which a vibration component is reduced by performing filter processing on the first measurement data;an observation information generation step of generating observation information including an entry time point and an exit time point of the moving object with respect to the structure;an average velocity calculation step of calculating an average velocity of the moving object based on the observation information and environment information which is created in advance and includes a dimension of the moving object and a dimension of the structure;a first deflection amount calculation step of calculating, based on an approximate equation of deflection of the structure, the observation information, the environment information, and the average velocity, a first deflection amount of the structure caused by the moving object;a second deflection amount calculation step of calculating a second deflection amount in which a vibration component is reduced by performing filter processing on the first deflection amount;a coefficient calculation step of approximating the second measurement data with a linear function of the second deflection amount to calculate a first-order coefficient and a zero-order coefficient of the linear function;a third deflection amount calculation step of calculating a third deflection amount based on the first-order coefficient, the zero-order coefficient, and the second deflection amount;an offset calculation step of calculating an offset based on the zero-order coefficient, the second deflection amount, and the third deflection amount; anda static response calculation step of calculating a static response by adding the offset and a product of the first-order coefficient and the first deflection amount.