Derivation Method, Derivation Device, Derivation System, And Program

Abstract

A derivation method includes: an acquisition step of acquiring time-series data including a physical quantity generated at a predetermined observation point in a structure as a response caused by a movement of a formation moving object formed with one or more moving objects on the structure; an environment information acquisition step of acquiring, as environment information, information on a structure length that is a length of the structure, a moving object length that is a length of the moving object, and an installation position of a contact portion of the moving object with the structure; a fundamental frequency derivation step of deriving a fundamental frequency of the time-series data based on the time-series data; a passing period derivation step of deriving a passing period during which the formation moving object passes through the structure based on the time-series data; and a number derivation step of deriving the number of the moving objects included in the formation moving object based on the environment information, the fundamental frequency, and the passing period.

Claims

1. A derivation method, comprising: an acquisition step of acquiring time-series data including a physical quantity generated at a predetermined observation point in a structure as a response caused by a movement of a formation moving object formed with one or more moving objects on the structure; an environment information acquisition step of acquiring, as environment information, information on a structure length that is a length of the structure, a moving object length that is a length of the moving object, and an installation position of a contact portion of the moving object with the structure; a fundamental frequency derivation step of deriving a fundamental frequency of the time-series data based on the time-series data; a passing period derivation step of deriving a passing period during which the formation moving object passes through the structure based on the time-series data; and a number derivation step of deriving the number of the moving objects included in the formation moving object based on the environment information, the fundamental frequency, and the passing period.

2. The derivation method according to claim 1, wherein in the passing period derivation step, an entry time point and an exit time point are derived by performing, for each of the entry time point and the exit time point of the formation moving object with respect to the structure, processing of acquiring a time point within a period after one time point and before another time point among two time points corresponding to two consecutive pieces of data between which a predetermined threshold is, the two consecutive pieces of data being included in the time-series data on which low-pass filter processing of attenuating a vibration component of a frequency equal to or higher than a fundamental frequency of the time-series data is performed, and a period from the derived entry time point to the derived exit time point is derived as the passing period.

3. The derivation method according to claim 1, wherein the physical quantity is a displacement of the structure, the method further comprises a moving object passing period derivation step of deriving a moving object passing period which is a period during which the moving object passes through the structure, and in the number derivation step, a value obtained by adding one to a product of a value obtained by subtracting the moving object passing period from the passing period and the fundamental frequency is derived as the number of moving objects included in the formation moving object.

4. The derivation method according to claim 3, wherein the contact portion is a wheel installed on an axle of the moving object, and in the moving object passing period derivation step, a value obtained by dividing a sum of the structure length and a distance between a foremost axle and a rearmost axle in the moving object by a product of the fundamental frequency and the moving object length is derived as the moving object passing period.

5. The derivation method according to claim 1, wherein in the fundamental frequency derivation step, a frequency corresponding to a peak having a lowest corresponding frequency among peaks indicated by a result of Fourier transform of the time-series data is derived as the fundamental frequency.

6. The derivation method according to claim 1, wherein a deflection model of the structure is an equation based on an architecture of the structure.

7. The derivation method according to claim 1, wherein the structure has a simple beam shape with both ends supported.

8. The derivation method according to claim 1, wherein the structure is a bridge.

9. The derivation method according to claim 1, wherein Bridge Weigh in Motion (BWIM) is applicable to the structure.

10. The derivation method according to claim 1, wherein the moving object is a railway vehicle that moves on the structure via a wheel.

11. The derivation method according to claim 1, wherein the time-series data is data based on data detected via at least one of an acceleration sensor, an impact sensor, a pressure-sensitive sensor, a strain gauge, an image measuring device, a load cell, and a displacement gauge.

12. A derivation device, comprising: an acquisition unit configured to acquire time-series data including a physical quantity generated at a predetermined observation point in a structure as a response caused by a movement of a formation moving object formed with one or more moving objects on the structure; an environment information acquisition unit configured to acquire, as environment information, information on a structure length that is a length of the structure, a moving object length that is a length of the moving object, and an installation position of a contact portion of the moving object with the structure; a fundamental frequency derivation unit configured to derive a fundamental frequency of the time-series data based on the time-series data; a passing period derivation unit configured to derive a passing period during which the formation moving object passes through the structure based on the time-series data; and a number derivation unit configured to derive the number of the moving objects included in the formation moving object based on the environment information, the fundamental frequency, and the passing period.

13. A derivation system comprising: a derivation device; and a sensor, wherein the derivation device includes: an acquisition unit configured to acquire time-series data including a physical quantity that is generated at a predetermined observation point in a structure as a response caused by a movement of a formation moving object formed with one or more moving objects on the structure and that is measured via the sensor; an environment information acquisition unit configured to acquire, as environment information, information on a structure length that is a length of the structure, a moving object length that is a length of the moving object, and an installation position of a contact portion of the moving object with the structure; a fundamental frequency derivation unit configured to derive a fundamental frequency of the time-series data based on the time-series data; a passing period derivation unit configured to derive a passing period during which the formation moving object passes through the structure based on the time-series data; and a number derivation unit configured to derive the number of the moving objects included in the formation moving object based on the environment information, the fundamental frequency, and the passing period.

14. A non-transitory computer-readable storage medium storing a program, the program for causing a computer to execute: an acquisition step of acquiring time-series data including a physical quantity generated at a predetermined observation point in a structure as a response caused by a movement of a formation moving object formed with one or more moving objects on the structure; an environment information acquisition step of acquiring, as environment information, information on a structure length that is a length of the structure, a moving object length that is a length of the moving object, and an installation position of a contact portion of the moving object with the structure; a fundamental frequency derivation step of deriving a fundamental frequency of the time-series data based on the time-series data; a passing period derivation step of deriving a passing period during which the formation moving object passes through the structure based on the time-series data; and a number derivation step of deriving the number of the moving objects included in the formation moving object based on the environment information, the fundamental frequency, and the passing period.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0011] FIG. 1 is a block diagram showing a configuration of a derivation system.

[0012] FIG. 2 is a diagram showing a cross section of a bridge.

[0013] FIG. 3 is a diagram showing dimensions of a unit bridge girder.

[0014] FIG. 4 is a diagram showing dimensions of a railway vehicle.

[0015] FIG. 5 is a diagram showing an outline of the unit bridge girder.

[0016] FIG. 6 is a diagram showing a bending moment at the unit bridge girder.

[0017] FIG. 7 is a diagram showing an outline of deflection of the unit bridge girder caused by a wheel.

[0018] FIG. 8 is a diagram showing an outline of deflection of the unit bridge girder caused by a railway vehicle.

[0019] FIG. 9 is a diagram showing an outline of deflection of the unit bridge girder caused by a railway train.

[0020] FIG. 10 is a diagram showing deflection of the unit bridge girder caused by each railway vehicle.

[0021] FIG. 11 is a diagram showing deflection of the unit bridge girder caused by each railway vehicle and the railway train.

[0022] FIG. 12 is a diagram showing time-series data of displacement of the unit bridge girder.

[0023] FIG. 13 is a diagram showing an FFT result of the time-series data.

[0024] FIG. 14 is a diagram showing derivation processing of an entry time point and an exit time point.

[0025] FIG. 15 is a diagram showing the derivation processing of the entry time point and the exit time point.

[0026] FIG. 16 is a diagram showing the time-series data of the displacement of the unit bridge girder.

[0027] FIG. 17 is a diagram showing the FFT result of the time-series data.

[0028] FIG. 18 is a diagram showing time-series data subjected to low-pass filter processing.

[0029] FIG. 19 is a diagram showing details of elements of the derivation system.

[0030] FIG. 20 is a flowchart showing the derivation processing.

DESCRIPTION OF EXEMPLARY EMBODIMENTS

[0031] Embodiments of the present disclosure will be described in the following order.

(1) Configuration of Derivation System

(1-1) Outline of Derivation System

(1-2) Deflection Model

(1-3) Verification Experiment

(1-4) Details of Elements

(2) Derivation Process

(3) Other Embodiments

(1) Configuration of Derivation System

(1-1) Outline of Derivation System

[0032] FIG. 1 is a block diagram showing an example of a configuration of a derivation system 10 according to the present embodiment. The derivation system 10 is a system that derives the number of railway vehicles included in a railway train 6 based on time-series data including a physical quantity at a predetermined observation point on a bridge 5 on which the railway train 6 formed with one or more railway vehicles moves. The railway train 6 is an example of a formation moving object. Each of the railway vehicles included in the railway train 6 is an example of a moving object. The bridge 5 is an example of a structure on which the moving object moves. Each railway vehicle of the railway train 6 moves on the bridge 5 via wheels provided on an axle. The wheel is an example of a contact portion between the railway vehicle and the bridge. In the present embodiment, each of the railway vehicles formed in the railway train 6 is a railway vehicle having the same structure. As shown in FIG. 1, the derivation system 10 includes a measurement device 1, at least one sensor device 2 provided in a superstructure 7 of the bridge 5, and a server device 3.

[0033] The measurement device 1 calculates deflection, that is, a displacement of the superstructure 7 caused by traveling of the railway train 6 based on acceleration data output from each sensor device 2. The measurement device 1 is installed on, for example, a bridge abutment 8b. The measurement device 1 and the server device 3 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 information on the displacement of the superstructure 7 caused by the traveling of the railway train 6 to the server device 3. The server device 3 derives the number of railway vehicles formed in the railway train 6 based on the transmitted displacement data.

[0034] In the present embodiment, the bridge 5 is a railroad bridge, and is, for example, a steel bridge, a girder bridge, or an RC bridge. The RC is an abbreviation for reinforced-concrete. In the present embodiment, the bridge 5 is a structure to which Bridge Weigh In Motion (BWIM) is applicable. The BWIM 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 a moving object passing through the bridge is measured. The bridge, which enables analysis of the weight of the moving object traveling on the bridge, based on a response such as deformation or strain of the bridge, is considered to be a structure to which BWIM is applicable. Therefore, the BWIM system, which applies a physical process between an action on the bridge and the response, enables the measurement of the weight of the moving object that travels on the bridge. The weight of the moving object is measured by measuring a correlation coefficient between the displacement and a load in advance, and using the correlation coefficient to derive the load of the moving object passing through from the measurement result of the displacement of the bridge when the moving object passes through.

[0035] The bridge 5 includes the superstructure 7 that is a portion where the moving object moves, and a substructure 8 that supports the superstructure 7. FIG. 2 is a cross-sectional view of the superstructure 7 taken along a 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 which is not shown. As shown in FIG. 1, the substructure 8 includes bridge piers 8a and the bridge abutments 8b. The superstructure 7 is a structure across the bridge abutment 8b and the bridge pier 8a adjacent to each other, two adjacent bridge abutments 8b, or two adjacent bridge piers 8a. Hereinafter, the bridge abutment 8b and the bridge pier 8a are collectively referred to as a support portion. In the present embodiment, a set of support portions and a portion of the bridge girder of the superstructure 7 across the set of support portions are collectively referred to as one bridge girder. That is, a simple beam-like structure in which both ends are supported by two support portions is defined as one bridge girder. Therefore, the bridge 5 shown in FIG. 1 includes two bridge girders. Hereinafter, each bridge girder included in the bridge 5 is referred to as a unit bridge girder.

[0036] The measurement device 1 and the sensor device 2 are coupled to each other, for example, in a wired or wireless manner, and communicate with one another via a communication network such as a controller area network (CAN).

[0037] The sensor device 2 is used to measure a predetermined physical quantity used to derive a displacement (deflection) at an observation point set on the superstructure 7. In the present embodiment, the predetermined physical quantity is an acceleration. In the present embodiment, the sensor device 2 is installed at the observation point. The sensor device 2 includes an acceleration sensor such as a quartz acceleration sensor or a micro electro mechanical systems (MEMS) acceleration sensor. The sensor device 2 outputs acceleration data for deriving the displacement of the superstructure 7 caused by a movement of the railway train 6 which is a moving object at the observation point.

[0038] In the present embodiment, the sensor device 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. The sensor device 2 is not limited to being installed at the central portion of the superstructure 7 as long as the sensor device 2 can detect an acceleration for calculating the displacement of the superstructure 7. When the sensor device 2 is provided on the floor plate F of the superstructure 7, the sensor device 2 may be damaged due to traveling of the railway train 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, the sensor device 2 is provided at the main girder G of the superstructure 7.

[0039] The floor plate F, the main girders G, and the like of the superstructure 7 are deflected in a vertical direction due to a load of the railway train 6 traveling on the superstructure 7. Each sensor device 2 measures an acceleration of the deflection of the floor plate F or the main girder G caused by the load of the railway train 6 traveling on the superstructure 7.

(1-2) Deflection Model

[0040] Here, a model of deflection of a bridge when a railway train moves on one bridge will be described. Here, the model is information such as an equation indicating a correspondence relationship between predetermined information and a derivation result.

[0041] In the following, the number of railway vehicles (the number of railway vehicles) formed in the railway train moving on the bridge is defined as N. An entry time point, which is a time point at which the railway train enters the bridge, is defined as t.sub.i. Here, the entry of the railway train into the bridge means that the wheels of a first axle of a railway vehicle C.sub.1 (a first railway vehicle from the head of the railway train) have entered the bridge. In the following, an exit time point, which is a time point at which the railway train exits from the bridge, is defined as to. Here, the exit of the railway train from the bridge means that the wheels of a rearmost axle of a railway vehicle CN (the rearmost railway vehicle of the railway train) have exited from the bridge. In addition, in the following, a period during which the railway train passes through the bridge (a period from the time point t.sub.i to the time point t.sub.o) is defined as t.sub.s. Hereinafter, N, t.sub.i, t.sub.o, and t.sub.s are collectively referred to as observation information.

[0042] In the following description, a bridge length, which is a length of the bridge in the traveling direction of the railway train, is defined as LB. The bridge length is an example of a structure length. A distance from an end in a direction in which the railway train enters among ends of the bridge in the longitudinal direction to the observation point is defined as L.sub.x. FIG. 3 shows the lengths L.sub.B and L.sub.x. In the following description, the end in a direction in which the railway train enters among the ends of the bridge in the longitudinal direction is referred to as an entry end. In addition, in the following, an end in a direction in which the railway train exits among the ends of the bridge in the longitudinal direction is referred to as an exit end. A vehicle length, which is a length in the traveling direction of the mth railway vehicle from the head of the railway train, is defined as L.sub.c(m). The vehicle length is an example of a moving object length that is a length of the moving object in the traveling direction. In the following, the lengths L.sub.c(1) to L.sub.c(N) are collectively referred to as a length L.sub.c. The mth railway vehicle from the head of the railway train is defined as C.sub.m. The number of axles in the railway vehicle Cm is defined as a.sub.r(m). Hereinafter, a.sub.r(l) to a.sub.r(N) are collectively referred to as a.sub.r. Hereinafter, the a.sub.r(m) axles in the railway vehicle Cm are defined as a first axle, a second axle, a third axle, . . . , an a.sub.r(m)th axle in order from the head of the railway vehicle C.sub.m.

[0043] A distance from a front end of the railway vehicle C.sub.m in a traveling direction to the first axle is defined as L.sub.a(a.sub.w(m, 1)). Here, a.sub.w(α, β) indicates a βth axle from a head axle of the αth railway vehicle in the railway train. The number of axles in the railway vehicle Cm is defined as a.sub.r(m). A distance between a (n−1)th axle and a nth axle in the railway vehicle C.sub.m is defined as L.sub.a(a.sub.w(m, n)), n being an integer of 2 or more. That is, for two or more β, L.sub.a(a.sub.w(α, β) indicates a distance between the βth axle and the (β-1)th axle in a railway vehicle C.sub.α. In addition, L.sub.a(a.sub.w(α, 1)) indicates a distance between a first axle of the railway vehicle C.sub.α and the front end of the railway vehicle C.sub.α in the traveling direction. Hereinafter, L.sub.a(a.sub.w(1, 1)) to L.sub.a(a.sub.w(N, a.sub.r(N))) are collectively referred to as L.sub.a. Each L.sub.a indicates a position of the corresponding axle in the corresponding railway vehicle. For example, L.sub.a(a.sub.w(m, 1)) indicates that the first axle is present behind the front end of the railway vehicle C.sub.m by a distance of L.sub.a(a.sub.w(m, 1)). L.sub.a(a.sub.w(m, 2)) indicates that the second axle is present behind the first axle of the railway vehicle C.sub.m by a distance of L.sub.a(a.sub.w(m, 2)).

[0044] Here, a railway vehicle having a similar four-axle configuration is formed in the railway train. That is, a.sub.r(m) is 4, m being 1, 2, . . . , N. FIG. 4 shows the lengths L.sub.c(m), L.sub.a(a.sub.w (m, 1)), L.sub.a(a.sub.w (m, 2)) L.sub.a(a.sub.w (m, 3)), and L.sub.a(a.sub.w (m, 4)) in the railway vehicle C.sub.m.

[0045] Hereinafter, L.sub.B, L.sub.x, L.sub.c, a.sub.r, and L.sub.a are collectively referred to as environment information.

[0046] As shown in the following Equation (1), t.sub.s is obtained as a difference between to and t.sub.i.

t.sub.s=t.sub.o−t.sub.i  (1)

[0047] The total number T.sub.ar of wheels of the railway train is obtained by the following Equation (2).

[00001]$\begin{array}{cc}{T}_{a_r}=\underset{m=1}{\overset{N}{.Math.}}{a}_r\left(m\right)& \left(2\right)\end{array}$

[0048] A distance from the first axle of the railway vehicle C.sub.1 at the head of the railway vehicle to the nth axle of the mth railway vehicle C.sub.m of the railway vehicle is represented as D.sub.wa(a.sub.w(m, n)). D.sub.wa(a.sub.w(m, n)) is obtained from the following Equation (3).

[00002]$\begin{array}{cc}{D}_{\mathrm{wa}}\left(a_w\left(m,n\right)\right)=\underset{y=1}{\overset{m}{.Math.}}{L}_c\left(y\right)+\underset{x=1}{\overset{n}{.Math.}}{L}_a\left(a_w\left(m,x\right)\right)-L_c\left(m\right)-L_a\left(a_w\left(1,1\right)\right)& \left(3\right)\end{array}$

[0049] The distance from the first axle of the railway vehicle C.sub.1 at the head of the railway train to the last axle a.sub.r(N) of the rearmost railway vehicle C.sub.N is D.sub.wa(a.sub.w(N, a.sub.r(N))). By using D.sub.wa(a.sub.w(N, a.sub.r(N))), an average velocity v.sub.a of the railway train passing through the bridge is represented by the following Equation (4).

[00003]$\begin{array}{cc}{v}_a=\frac{L_B}{t_s}+\frac{D_{\mathrm{wa}}\left(a_w\left(N,a_r\left(N\right)\right)\right)}{t_s}& \left(4\right)\end{array}$

[0050] From Equation (3) and Equation (4), the following Equation (5) is established.

[00004]$\begin{array}{cc}{v}_a=\frac{L_B}{t_S}+\frac{1}{t_S}\left[\underset{m=1}{\overset{N}{.Math.}}{L}_c\left(m\right)+\underset{n=1}{\overset{a_r\left(N\right)}{.Math.}}{L}_a\left(a_w\left(N,n\right)\right)-L_c\left(N\right)-L_a\left(a_w\left(1,1\right)\right)\right]& \left(5\right)\end{array}$

[0051] Next, deflection generated in the bridge when a load is applied to the bridge will be described.

[0052] FIG. 5 is a schematic view of the bridge. FIG. 5 shows a situation in which a load P is applied to the bridge. Here, a distance between a position of the bridge to which the load P is applied and the entry end is represented by a. A distance between the position of the bridge to which the load P is applied and the exit end is represented by b. In this case, the bending moment at the position of the bridge to which the load P is applied is represented by the following Equation (6).

[00005]$\begin{array}{cc}M=\frac{\mathrm{abP}}{L_B}& \left(6\right)\end{array}$

[0053] FIG. 6 shows the bending moment at each position of the bridge due to the load P. As shown in FIG. 6, the bending moment generated in the bridge due to the load P is 0 at the entry end, increases proportionally as the position approaches the position to which the load P is applied from the entry end, and becomes a value represented by Equation (6) at the position to which the load P is applied. The bending moment generated in the bridge due to the load P decreases proportionally as the position approaches the exit end from the position to which the load P is applied, and becomes 0 at the exit end. Therefore, the bending moment at an optional position X in the bridge is represented by the following Equation (7).

[00006]$\begin{array}{cc}M=\frac{\mathrm{bP}}{L_B}x-H_{a}P\left(x-a\right)& \left(7\right)\end{array}$

[0054] In Equation (7), x represents the distance from the entry end to the position X in the traveling direction of the railway train. Ha in Equation (7) is a value represented by the following Equation (8).

[00007]$\begin{array}{cc}{H}_a=\{\begin{array}{cc}x\le a:& 0\\ x>a:& 1\end{array}& \left(8\right)\end{array}$

[0055] Between the bending moment and a deflection w of the bridge at the optional position X, a relationship represented by the following Equation (9) is established.

[00008]$\begin{array}{cc}-M=\mathrm{EI}\frac{d\theta }{\mathrm{dx}}=\mathrm{EI}\frac{d^{2}w}{{\mathrm{dx}}^2}& \left(9\right)\end{array}$

[0056] θ in Equation (9) is an angle formed by a horizontal line and the deflected bridge at the position X. From Equation (7) and Equation (9), the following Equation (10) is established.

[00009]$\begin{array}{cc}\frac{d^{2}w}{{\mathrm{dx}}^2}=-\frac{1}{\mathrm{EI}}\left(\frac{\mathrm{bP}}{L_B}x-H_{a}P\left(x-a\right)\right)& \left(10\right)\end{array}$

[0057] By integrating both sides of Equation (10) twice with x, the following Equation (11) representing the deflection w at the position X is obtained.

[00010]$\begin{array}{cc}w=\frac{P}{6{\mathrm{EIL}}_B}\left\{-{\mathrm{bx}}^3+H_a{L_B\left(x-a\right)}^3+ℊ1x+ℊ2\right\}& \left(11\right)\end{array}$

[0058] In Equation (11), g1 and g2 are constant terms. Here, since the bridge is supported at the entry end and the exit end, no deflection is generated at the positions of the entry end and the exit end. That is, in Equation (11), when x=0 and x=L.sub.B, both sides are 0. Therefore, g1 and g2 are represented by the following Equation (12) and Equation (13).

g1=ab(a+2b)  (12)

g2=0  (13)

[0059] From Equation (11), Equation (12), and Equation (13), the following Equation (14) representing the deflection w at the position X is obtained.

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

[0060] When the load P is applied to a center of the bridge in the longitudinal direction, the maximum deflection among the deflection generated in the bridge due to the application of the load P is generated at the center of the bridge in the longitudinal direction. When this maximum deflection is w.sub.0.51, an equation representing w.sub.0.51 is obtained. When the load P is applied to the center of the bridge in the longitudinal direction, a=b=0.5L.sub.B. Since the position X of a target for which the deflection is to be obtained is the center of the bridge in the longitudinal direction, x=0.5L.sub.B. In this case, since x≤a, H.sub.a=0 is obtained from Equation (8). By substituting x=0.5L.sub.B, a=b=0.5L.sub.B, and H.sub.a=0 into Equation (14), the following Equation (15) representing the deflection w.sub.0.51 is obtained.

[00012]$\begin{array}{cc}{w}_{0.5l}=\frac{P}{48\mathrm{EI}}{L_B}^3& \left(15\right)\end{array}$

[0061] The deflection at an optional position in the bridge represented by Equation (14) is normalized using w.sub.0.51.

[0062] When the position of the load P exists on the entry end side with respect to the position X, that is, when x>a, H.sub.a=1 is obtained from Equation (8), and Equation (14) is represented as Equation (16) below.

[00013]$\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(16\right)\end{array}$

[0063] a=L.sub.Br. Here, r is a real number from 0 to 1. Since b=L.sub.B−a, b is represented as b=L.sub.B(1−r). When a=L.sub.Br and b=L.sub.B(1−r) are substituted into Equation (16) and the deflection is normalized by dividing by w.sub.0.51, the following Equation (17) representing normalized deflection w.sub.std at the position X when x>a is obtained.

[00014]$\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(17\right)\end{array}$

[0064] Similarly, when the position of the load P exists on the exit end side with respect to the position X, that is, when x≤a, H.sub.a=0 is obtained from Equation (8), and Equation (14) is represented as the following Equation (18).

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

[0065] a=L.sub.Br. Here, r is a real number from 0 to 1. Since b=L.sub.B−a, b is represented as b=L.sub.B(1−r). When a=L.sub.Br and b=L.sub.B(1−r) are substituted into Equation (18) and the deflection is normalized by dividing by w.sub.0.51, the following Equation (19) representing the normalized deflection w.sub.std at the position X when x≤a is obtained.

[00016]$\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{\mathrm{xr}}^2+\frac{x^3}{{L_B}^2}\right\}& \left(19\right)\end{array}$

[0066] By substituting L.sub.x for x in Equation (17) and Equation (19), the normalized deflection w.sub.std at the deflection observation point is represented as the following Equation (20) as a function of r.

[00017]$\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(20\right)\end{array}$

[0067] The function R(r) in Equation (20) is a function represented by the following Equation (21).

[00018]$\begin{array}{cc}R\left(r\right)=\{\begin{array}{cc}{L}_x>L_{B}r:& L_B{r}^3+\frac{3{L_x}^2}{L_B}r\\ L_x\le L_{B}r:& 3L_x{r}^2+\frac{{L_x}^3}{{L_B}^2}\end{array}& \left(21\right)\end{array}$

[0068] Here, using Equation (20) and Equation (21), a function indicating a temporal change in deflection generated at an observation point due to a load applied to the bridge via a wheel of any one axle a.sub.w(m, n) is obtained. First, a period required for the wheel of one axle of the railway train to reach the observation point from the entry end is defined as t.sub.xn. t.sub.xn is obtained from L.sub.x and v.sub.a by the following Equation (22).

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

[0069] A period during which one wheel of the railway train crosses the bridge, that is, a period from the entry end to the exit end is defined by t.sub.ln. t.sub.ln is obtained from L.sub.B and v.sub.a by the following Equation (23).

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

[0070] A time point at which the wheel of the nth axle a.sub.w(m, n) of the mth railway vehicle of the railway train reaches the entry end is defined as t.sub.o(m, n). t.sub.o(m, n) is obtained from t.sub.i, v.sub.a, and D.sub.wa(a.sub.w(m, n)) by the following Equation (24).

[00021]$\begin{array}{cc}{t}_0\left(m,n\right)=t_i+\frac{D_{\mathrm{wa}}\left(a_w\left(m,n\right)\right)}{v_a}& \left(24\right)\end{array}$

[0071] From Equation (22), L.sub.x is represented as the following Equation (25).

L.sub.x=v.sub.at.sub.xn  (25)

[0072] From Equation (23), L.sub.B is represented as the following Equation (26).

L.sub.B=v.sub.at.sub.ln  (26)

[0073] The position of the axle a.sub.w(m, n) is a load position. Therefore, the position of the axle a.sub.w(m, n) is a position at a distance of a=L.sub.Br in the direction from the entry end to the exit end. When a variable indicating the time point is t, a distance from the entry end of a.sub.w(m, n) at the time point t is equal to a distance traveled by the railway vehicle from the time point t.sub.o(m, n) to the time point t. Therefore, the following Equation (27) is established.

L.sub.Br=v.sub.a(t−t.sub.0(m,n))  (27)

[0074] From Equation (27), r is represented as in the following Equation (28).

[00022]$\begin{array}{cc}r=\frac{v_a\left(t-t_0\left(m,n\right)\right)}{L_B}=\frac{v_a\left(t-t_0\left(m,n\right)\right)}{v_a{t}_{\ln }}=\frac{\left(t-t_0\left(m,n\right)\right)}{t_{\ln }}& \left(28\right)\end{array}$

[0075] By replacing L.sub.x, L.sub.B, and r in Equation (20) and Equation (21) using Equation (25), Equation (26), and Equation (28), a function w.sub.std(a.sub.w(m, n), t) in the following Equation (29) is obtained as a model indicating a temporal change in deflection generated at the observation point due to the load applied to the bridge via the wheel of the axle a.sub.w(m, n). A function R(t) in Equation (29) is a function represented by the following Equation (30).

[00023]$\begin{array}{cc}{w}_{\mathrm{std}}\left(a_w\left(m,n\right),t\right)=\{\begin{array}{c}\begin{array}{c}\begin{array}{c}t<t_0\left(m,n\right):0\\ t_0\left(m,n\right)\le t\le t_0\left(m,n\right)+t_{\ln }:\end{array}\\ \frac{8}{t_{\ln }}\left\{{t_{\mathrm{xn}}\left(\frac{t-t_0\left(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(m,n\right)}{t_{\ln }}\right)-R\left(t\right)\right\}\end{array}\\ t_0\left(m,n\right)+t_{\ln }<t:0\end{array}& \left(29\right)\end{array}$$\begin{array}{cc}R\left(t\right)=\{\begin{array}{c}\begin{array}{c}\begin{array}{c}t<t_0\left(m,n\right):0\\ \begin{array}{c}{t}_0\left(m,n\right)\le t\le t_0\left(m,n\right)+t_{\ln }.Math.t_{\mathrm{xn}}>t-t_0\left(m,n\right):\\ {t_{\ln }\left(\frac{t-t_0\left(m,n\right)}{t_{\ln }}\right)}^3+\frac{3t_{\mathrm{xn}}^2}{t_{\ln }}\left(\frac{t-t_0\left(m,n\right)}{t_{\ln }}\right)\end{array}\end{array}\\ t_0\left(m,n\right)\le t\le t_0\left(m,n\right)+t_{\ln }.Math.t_{\mathrm{xn}}\le t-t_0\left(m,n\right):\end{array}\\ 3{t_{\mathrm{xn}}\left(\frac{t-t_0\left(m,n\right)}{t_{\ln }}\right)}^2+\frac{t_{\mathrm{xn}}^3}{t_{\ln }^2}\\ t_0\left(m,n\right)+t_{\ln }<t:0\end{array}& \left(30\right)\end{array}$

[0076] When the observation information and the environment information (t.sub.i, t.sub.o, N, L.sub.B, L.sub.x, L.sub.c(1) to L.sub.c(N), a.sub.r(1) to a.sub.r(N), and L.sub.a(a.sub.w(1, 1)) to L.sub.a(a.sub.w(N.sub.r a.sub.r(N)))) are known, w.sub.std(a.sub.w(m, n), t) is obtained using the information. For example, t.sub.s is obtained from t.sub.i and t.sub.o using Equation (1). From t.sub.s, N, a.sub.r, L.sub.a, and L.sub.c, v.sub.a is obtained using Equation (5). From v.sub.a, L.sub.B, and L.sub.x, t.sub.xn and tin are obtained using Equation (22) and Equation (23). From L.sub.a, L.sub.c, and t.sub.i, t.sub.o(m, n) is obtained using Equation (3) and Equation (24). Then, by substituting the obtained t.sub.xn, t.sub.ln, and t.sub.o(m, n) into the Equation (29) and Equation (30), the function w.sub.std(a.sub.w(m, n), t) of t is obtained.

[0077] An example of a change in the deflection amount at the observation point indicated by w.sub.std(a.sub.w(m, n), t) is shown in FIG. 7. In a graph of FIG. 7, a horizontal axis represents time, and a vertical axis represents the deflection amount. In accordance with the movement of one railway vehicle C.sub.m, a set of wheels for each of the a.sub.r(m) axles moves on the bridge. Therefore, a function C.sub.std(m, t) serving as a model indicating a temporal change in the deflection amount generated at the observation point due to the movement of one railway vehicle C.sub.m is obtained as the sum of w.sub.std(a.sub.w(m, n), t) for the respective axles as in the following Equation (31).

[00024]$\begin{array}{cc}{C}_{\mathrm{std}}\left(m,t\right)=\underset{n=1}{\overset{a_r\left(m\right)}{.Math.}}{w}_{\mathrm{std}}\left(a_w\left(m,n\right),t\right)& \left(31\right)\end{array}$

[0078] FIG. 8 shows how the deflection amount changes at the observation point indicated by the function C.sub.std(m, t) when a.sub.r(m) is 4, that is, when the railway vehicle C.sub.m has a four-axle configuration. In a graph of FIG. 8, a horizontal axis represents time, and a vertical axis represents the deflection amount. A solid line in the graph of FIG. 8 indicates C.sub.std(m, t), and each dotted line in the graph indicates w.sub.std(a.sub.w(m, n), t) for each axle.

[0079] In accordance with the movement of the railway train, the N railway vehicles move on the bridge. Therefore, a function T.sub.std(t) serving as a model indicating a temporal change in the deflection amount generated at the observation point due to movement of one railway train is obtained as the sum of C.sub.std(m, t) for the respective railway vehicles as in the following Equation (32).

[00025]$\begin{array}{cc}{T}_{\mathrm{std}}\left(t\right)=\underset{m=1}{\overset{N}{.Math.}}{C}_{\mathrm{std}}\left(m,t\right)& \left(32\right)\end{array}$

[0080] FIG. 9 shows how the deflection amount changes at the observation point indicated by the function T.sub.std(t) when N is 16, that is, when 16 railway vehicles are formed in the railway train. In a graph of FIG. 9, a horizontal axis represents time, and a vertical axis represents the deflection amount. A solid line in the graph of FIG. 9 indicates T.sub.std(t), and each dotted line in the graph indicates C.sub.std(m, t) for each railway vehicle. As shown in the graph of FIG. 9, the waveform is obtained by adding together the deflection of each passing railway vehicle, and it can be seen that vibration occurs in a cycle in which continuous vehicles pass through the bridge.

[0081] The deflection model of the bridge is described as above.

(1-3) Verification Experiment

[0082] Deflection amounts C.sub.std(1, t) to C.sub.std(N, t) and T.sub.std(t) generated when a railway train formed with the same railway vehicles passes through the bridge are considered under conditions in which the observation information and the environment information have the following values. That is, N=4, t.sub.i=7.21 [sec], t.sub.o=8.777 [sec], t.sub.s=1.567 [sec], L.sub.B=25 [in] L.sub.x=12.5 [In] L.sub.c=25 [In] a.sub.r=4, L.sub.a(a.sub.w (m, 1))=2.5 [m] for each of m=1 t.sub.o N, L.sub.a(a.sub.w(m, 2))=2.5 [m] for each of m=1 t.sub.o N, L.sub.a(a.sub.w(m, 3))=15 [m] for each of m=1 t.sub.o N, and L.sub.a(a.sub.w(m, 4))=2.5 [m] for each of m=1 to N.

[0083] The deflection amounts C.sub.std(1, t) to C.sub.std(4, t) by four railway vehicles included in the railway train at this time are shown in FIG. 10. A cycle of vibration generated in the bridge when the railway vehicles continuously pass through the bridge is defined as T.sub.f. The vibration generated in the bridge when the railway vehicles continuously pass through the bridge is vibration generated when the continuous railway vehicles pass through the bridge. For this reason, the cycle T.sub.f is a time difference between the entry time points of the continuous railway vehicles passing through the bridge into the bridge. Since the bridge is deflected by the railway vehicle from the time point when the railway vehicle enters the bridge, the time difference between a start time point of the deflection indicated by C.sub.std(m, t) and a start time point of the deflection indicated by C.sub.std(m+1, t) is the cycle T.sub.f. FIG. 10 shows the deflection generated in the bridge due to the passage of each railway vehicle of the railway train when the railway train passes through the bridge. In a graph of FIG. 10, a horizontal axis represents time, and a vertical axis represents the deflection amount. As shown in FIG. 10, the deflection due to the passing of the railway vehicles in front and rear occurs at a time difference of T.sub.f.

[0084] Since the cycle T.sub.f is a time difference between the entry time points of the continuous railway vehicles passing through the bridge into the bridge, as shown in the following Equation (33), the cycle T.sub.f can be regarded as a period during which the railway vehicle having a vehicle length L.sub.c(m) passes through at a velocity v.sub.a.

[00026]$\begin{array}{cc}{T}_f=\frac{L_C\left(m\right)}{v_a}& \left(33\right)\end{array}$

[0085] A period during which the railway vehicle C.sub.m of the railway train passes through the bridge is defined as t.sub.c(m). t.sub.c(m) is an example of a moving object passing period, which is a period during which the railway vehicle C.sub.m which is a moving object passes through a bridge which is a structure. t.sub.c(m) is a period from the time point when the first axle of the railway vehicle C.sub.m reaches the entry end to the time point when the a.sub.r(m)th axle of the railway vehicle C.sub.m reaches the exit end. That is, t.sub.c(m) is a period during which the railway vehicle C.sub.m moves by a total distance of the bridge length L.sub.B and the distance from the first axle which is the foremost axle of the railway vehicle C.sub.m to the a.sub.r(m)th axle which is the rearmost axle of the railway vehicle C.sub.m. Therefore, t.sub.c(m) is represented by the following Equation (34).

[00027]$\begin{array}{cc}{t}_c\left(m\right)=\left\{L_B+\underset{n=1}{\overset{a_r\left(m\right)}{.Math.}}{L}_a\left(a_w\left(m,n\right)\right)-L_a\left(a_w\left(m,1\right)\right)\right\}/v_a& \left(34\right)\end{array}$

[0086] When the railway train passes through the bridge, the number of railway vehicles for which the subsequent railway vehicle is present among the railway vehicles formed in the railway train is defined as C.sub.Tn. Among the railway vehicles formed in the railway train, for the railway vehicles other than the rearmost railway vehicle, there is a subsequent railway vehicle. Therefore, C.sub.Tn is a number smaller than N by 1. That is, the following Equation (35) is established.

t.sub.s=C.sub.TnT.sub.f+t.sub.c(m)  (35)

[0087] FIG. 11 shows C.sub.std(1, t) to C.sub.std(N, t) and T.sub.std(t). In a graph of FIG. 11, a horizontal axis represents time, and a vertical axis represents the deflection amount. A solid line in the graph of FIG. 11 indicates T.sub.std(t), and dotted lines in the graph indicate C.sub.std(1, t) to C.sub.std(4, t). As shown in FIG. 11, the passing period t.sub.s is the sum of C.sub.Tn T.sub.fs and a period t.sub.c(m) during which one railway vehicle C.sub.m passes through the bridge. That is, the following Equation (36) is established.

N=C.sub.Tn+1  (36)

[0088] From Equation (35) and Equation (36), the number N of railway vehicles formed in the railway train is represented by the following Equation (37).

[00028]$\begin{array}{cc}N=\frac{t_s-t_c\left(m\right)}{T_f}+1& \left(37\right)\end{array}$

[0089] T.sub.f is also a period required for the railway train to move by the vehicle length of one railway vehicle. Therefore, the distance by which the railway train travels during the passing period t.sub.s is the sum of the length of the (N−1) railway vehicles and the distance by which the railway train travels during the period t.sub.c(m) at the velocity v.sub.a. Therefore, the following Equation (38) is established.

[00029]$\begin{array}{cc}{v}_a=\frac{\left(N-1\right)L_C\left(m\right)+t_c\left(m\right)v_a}{t_s}& \left(38\right)\end{array}$

[0090] From Equation (38), the following Equation (39) is established. From Equation (39), it can be confirmed that Equation (37) is established.

[00030]$\begin{array}{cc}\left(N-1\right)=\frac{v_a{t}_s-t_c\left(m\right)v_a}{L_C\left(m\right)}=\frac{t_s-t_c\left(m\right)}{\frac{L_C\left(m\right)}{v_a}}=\frac{t_s-t_c\left(m\right)}{T_f}& \left(39\right)\end{array}$

[0091] It is considered that the deflection amount T.sub.std(t) generated in the bridge when the railway train passes through the bridge includes, as a component of the fundamental frequency, a component of vibration generated in the bridge in accordance with the movement of the continuous railway vehicles. Here, the fundamental frequency is a frequency of a component having the lowest frequency included in the signal. This fundamental frequency is defined as F.sub.f. F.sub.f is also the frequency of the vibration generated in the bridge in accordance with the movement of the continuous railway vehicles, and thus can be represented as a reciprocal of t.sub.f as shown in the following Equation (40).

[00031]$\begin{array}{cc}{F}_f=\frac{1}{T_f}& \left(40\right)\end{array}$

[0092] From Equation (33) and Equation (40), the velocity v.sub.a is represented by the product of F.sub.f and L.sub.c(m) as in the following Equation (41).

v.sub.a=F.sub.fL.sub.c(m)  (41)

[0093] Therefore, t.sub.c(m) represented by Equation (34) is a value obtained by dividing the total distance of the bridge length L.sub.B and the distance from the first axle which is the foremost axle of the railway vehicle C.sub.m to the a.sub.r(m)th axle which is the rearmost axle by the product of F.sub.f and L.sub.c(m).

[0094] From Equation (37) and Equation (40), the number N of railway vehicles formed in the railway train is represented as a value obtained by adding one to the product of the fundamental frequency F.sub.f and a value obtained by subtracting the passing period t.sub.c(m) of one railway vehicle C.sub.m passing through the bridge from the passing period t.sub.s of the railway train passing through the bridge, and is represented by the following Equation (42).

N=(t.sub.s−t.sub.c(m))F.sub.f+1  (42)

[0095] The inventors found that the average velocity v.sub.a of the railway train is represented by the product of the fundamental frequency F.sub.f and the length of one railway vehicle C.sub.m included in the railway train, as represented by Equation (41). In addition, the inventors found that the period t.sub.c(m) during which one railway vehicle C.sub.m passes through the bridge is represented as a period during which the railway vehicle C.sub.m moves a length of the sum of the length L.sub.B of the bridge and the distance from the first axle of the railway vehicle C.sub.m to the a.sub.r(m)th axle at the velocity v.sub.a, as represented by Equation (34). Further, the inventors found that the number N of railway vehicles formed in the railway train is represented as a value obtained by adding one to the product of the fundamental frequency F.sub.f and a value obtained by subtracting t.sub.c(m) from t.sub.s, as represented by Equation (42).

[0096] Then, the inventors conceived a method of deriving the number of railway vehicles formed in the railway train by using time-series data of a displacement at an observation point set on a bridge on which the railway train moves.

[0097] The method conceived by the inventors is as follows.

[0098] The time-series data of the displacement at the observation point set on the bridge on which the railway train moves is acquired. Hereinafter, the time-series data of the displacement at the observation point set on the bridge is defined as u(t). u(t) is data of discrete values of the displacement measured in a predetermined cycle, and is data in which each discrete value is associated with a measurement time point. In addition, the lengths L.sub.B, L.sub.c, and L.sub.a are acquired as the environmental information. Then, based on the time-series data u(t), the fundamental frequency F.sub.f of u(t) is acquired as the frequency of the vibration generated in the bridge due to the passage of the continuous railway vehicles formed in the railway train. Based on u(t), the period t.sub.s during which the railway train passes through the bridge is derived. Then, based on the lengths L.sub.B, L.sub.c and L.sub.a, the fundamental frequency F.sub.f, and the passing period t.sub.s, the number of railway vehicles included in the railway train is derived using the relationships represented by Equation (34), Equation (41), and Equation (42).

[0099] The inventors conducted an experiment for confirming the effectiveness of the conceived method. This experiment will be described below.

[0100] The inventors acquired time-series data of the displacement at the observation point by measuring the displacement of the bridge at the observation point set on the bridge in a predetermined cycle ΔT when the railway train formed with 16 railway vehicles passes through a predetermined bridge. This predetermined bridge is referred to as a first bridge. In the first bridge, with the passage of the railway train, resonance with vibration generated in the first bridge due to the passage of the railway train formed with the continuous railway vehicles did not occur. Here, the resonance means that the vibration generated in the bridge due to the passage of the railway train and a natural frequency of the bridge approximate each other, and thus the natural frequency or harmonic vibration of the natural frequency is excited in the bridge due to the vibration generated in the bridge due to the passage of the railway train. Hereinafter, the vibration generated by resonance is referred to as a dynamic response.

[0101] The environmental information at this time is known and is as follows. L.sub.B=25 [m], L.sub.x=12.5 [m], L.sub.c=25 [m], a.sub.r=4, L.sub.a(a.sub.w(m, 1))=2.5 [m] for each of m=1 to N, L.sub.a(a.sub.w(m, 2))=2.5 [m] for each of m=1 to N, L.sub.a(a.sub.w(m, 3))=15 [m] for each of m=1 to N, and L.sub.a(a.sub.w(m, 4))=2.5 [m] for each of m=1 to N.

[0102] The observation information is not known except for N. That is, it is assumed that t.sub.i, t.sub.o, and t.sub.s are not known in advance. The number N of railway vehicles formed in the railway train is 16.

[0103] FIG. 12 shows time-series data u(t) of the displacement of the observation point in the first bridge acquired by the inventors. In a graph of FIG. 12, a horizontal axis represents time, and a vertical axis represents the deflection amount.

[0104] The inventors performed fast Fourier transform (FFT) on u(t). A result of the FFT of u(t) is shown in FIG. 13. In a graph of FIG. 13, a horizontal axis represents the frequency, and a vertical axis represents the intensity of the corresponding frequency component. Then, the inventors derived the fundamental frequency F.sub.f of u(t) from the result of the FFT of u(t) as the frequency of the vibration generated in the first bridge in accordance with the movement of the continuous railway vehicles. Specifically, the inventors specified a peak corresponding to the lowest frequency from the result of the FFT of u(t) except for a side lobe generated due to an influence of a window function used in the FFT, and obtained the specified peak as the fundamental frequency. In an example of FIG. 13, the inventors obtained a value of 3.01 Hz as the fundamental frequency F.sub.f from the graph of FIG. 13.

[0105] The inventors also obtained the cycle T.sub.f of the vibration generated in the first bridge due to the passage of the railway train by obtaining a reciprocal of the obtained fundamental frequency F.sub.f. The inventors performed low-pass filter processing for attenuating a component of a frequency equal to or higher than the fundamental frequency on u(t) by performing a moving average on u(t) in the cycle T.sub.f. u(t) subjected to the low-pass filter processing is defined as u.sub.lp(t)=u.sub.lp(kΔT). Here, k is a variable indicating what number of observations when the deflection amount is observed in a cycle at the observation point. That is, when a data cycle (time resolution) of the observation of the deflection amount is ΔT, t=kΔT. Hereinafter, the low-pass filter processing will be described in detail.

[0106] As shown in the following Equation (43), a moving average interval k.sub.mf adjusted to the time resolution of the data is obtained from the cycles T.sub.f and ΔT.

[00032]$\begin{array}{cc}{k}_{\mathrm{mf}}=2.Math.\frac{T_f}{2\Delta T}.Math.+1& \left(43\right)\end{array}$

[0107] By using k.sub.mf, u.sub.lp(t) is obtained by the following Equation (44).

[00033]$\begin{array}{cc}{u}_{\mathrm{lp}}\left(t\right)=u_{\mathrm{lp}}\left(k\Delta T\right)=\frac{1}{k_{\mathrm{mf}}}\underset{n=k-\frac{k_{\mathrm{mf}}-1}{2}}{\overset{k+\frac{k_{\mathrm{mf}}-1}{2}}{.Math.}}u\left(n\Delta T\right)& \left(44\right)\end{array}$

[0108] The inventors obtained k.sub.mf from the cycles T.sub.f and ΔT using Equation (43), and obtained u.sub.lp(t) from the obtained k.sub.mf and u(t) using Equation (44). The above is the details of the low-pass filter processing. The low-pass filter processing may be processing of applying an FIR filter under a condition of attenuating a component of a frequency equal to or higher than the fundamental frequency F.sub.f.

[0109] Then, the inventors specified, from ulp(t), two consecutive pieces of data between which a predetermined threshold C.sub.L related to the deflection amount is. Here, the fact that the threshold C.sub.L is between the two consecutive pieces of data of u.sub.lp(t) means that the threshold CL is included in a range between the values of the two pieces of consecutively measured data of the displacement included in u.sub.lp(t), that is, a range from the smaller value of these pieces of data of the displacement to the larger value of these pieces of data of the displacement. The threshold C.sub.L is a value of deflection generated in the bridge in accordance with the entry of the railway train into the bridge, and is, for example, a value of deflection at the observation point of the bridge when the railway vehicle is disposed such that the wheel of the first axle at the head railway vehicle is placed in the vicinity of the entry end. The threshold C.sub.L may be another value as long as the threshold C.sub.L can used for detecting the entry of the railway train into the bridge, and may be, for example, a deflection amount of the observation point of the bridge when a predetermined weight is applied to the vicinity of the entry end. The threshold C.sub.L may be a value of a predetermined ratio (for example, 10%, 1%, or the like) of the maximum value of the deflection amount at the observation point of the bridge when the railway train passes through the bridge. The threshold C.sub.L may be set to a value of any data included in u.sub.lp(t).

[0110] FIG. 14 shows u.sub.lp(t) and the threshold C.sub.L. In a graph of FIG. 14, a horizontal axis represents time (t=kΔT), and a vertical axis represents the deflection amount. A solid line in the graph of FIG. 14 indicates u.sub.lp(t), and a dotted line in the graph indicates u(t). In portions surrounded by dotted circles in FIG. 14, u.sub.lp(t) and the threshold C.sub.L intersect with each other. FIG. 15 shows an enlarged view of a portion where u.sub.lp(t) and C.sub.L intersect with each other (a portion of a dotted circle on the left side in the graph of FIG. 14). In a graph of FIG. 15, a horizontal axis represents time, and a vertical axis represents the deflection amount. Black dots in FIG. 15 indicate data of discrete values included in u.sub.lp(t). In an example of FIG. 15, the threshold C.sub.L is between data k−1 included in u.sub.lp(t) and data k.

[0111] The inventors specified a later one of two time points corresponding to the two consecutive pieces of data between which the specified C.sub.L is. In the example of FIG. 15, a time point kΔT corresponding to the data k is specified.

[0112] In an example of FIG. 14, the inventors specified two pieces of data in a portion of a dotted circle on the right side in FIG. 14 as two consecutive pieces of data between which C.sub.L is, and specified the later one of two time points corresponding to the specified two pieces of data.

[0113] Then, the inventors derived the earlier one of the specified time points as the entry time point t.sub.i of the railway train entering into the bridge. The inventors also derived the later one of the specified time points as the exit time point t.sub.o of the railway train exiting from the bridge. In the example of FIG. 14, the inventors derived that the entry time point t.sub.i=7.2 [s] and the exit time point t.sub.o=12.795 [s]. Then, the inventors derived t.sub.s as t.sub.o−t.sub.i=12.795−7.2=5.595 [sec].

[0114] The inventors derived, using Equation (41), the average velocity v.sub.a of the railway train as 25 [m]×3.01 [Hz]=75.25 [m/sec] based on L.sub.c(m)=25 [m] and the obtained F.sub.f=3.01 [Hz]. The inventors also derived, using Equation (34), the period t.sub.c(m) during which one railway vehicle C.sub.m passes through the first bridge as (25+2.5+2.5+15+2.5−2.5)/75.25=0.5980 [seconds] based on L.sub.B, L.sub.a, and the derived v.sub.a. In addition, the inventors derived, using Equation (42), the number N of railway vehicles formed in the railway train as (5.595−0.5980)×3.01+1=16.04097 based on the derived F.sub.f, t.sub.s, and t.sub.c(m). This estimated value is approximated to 16, which is the number N of railway vehicles. Therefore, the inventors confirmed that the number of railway vehicles included in the railway train can be derived with high accuracy by using the conceived method.

[0115] In addition, the inventors conducted an experiment of deriving the number of railway vehicles formed in the railway train by using the conceived method when a similar railway train passes through a second bridge different from the first bridge. In the second bridge, with the passage of the railway train, resonance with vibration generated in the bridge due to the passage of the railway train formed with the continuous railway vehicles occurred.

[0116] The environmental information at this time is known and is as follows. L.sub.B=25 [m], L.sub.x=12.5 [m], L.sub.c=25 [m], a.sub.r=4, L.sub.a(a.sub.w(m, 1))=2.5 [m] for each of m=1 to N, L.sub.a(a.sub.w(m, 2))=2.5 [m] for each of m=1 to N, L.sub.a(a.sub.w(m, 3))=15 [m] for each of m=1 to N, and L.sub.a(a.sub.w(m, 4))=2.5 [m] for each of m=1 to N.

[0117] The observation information is not known except for N. That is, it is assumed that t.sub.i, t.sub.o, and t.sub.s are not known in advance. The number N of railway vehicles formed in the railway vehicle is 16.

[0118] The inventors acquired time-series data of the displacement at the observation point by measuring the displacement of the bridge at the observation point set on the second bridge in a predetermined cycle ΔT when the railway train passes through the second bridge. The time-series data of the displacement of the observation point in the second bridge acquired by the inventors is defined as u(t). FIG. 16 shows this u(t). In a graph of FIG. 16, a horizontal axis represents time, and a vertical axis represents the deflection amount. The graph of FIG. 16 shows a dynamic response generated in the bridge as compared with the graph of FIG. 12.

[0119] The inventors performed the FFT on u(t). The result of the FFT of u(t) is shown in FIG. 17. In a graph of FIG. 17, a horizontal axis represents the frequency, and a vertical axis represents the intensity of the corresponding frequency component. The dynamic response is generated by resonance with vibration generated in the second bridge due to the passage of the railway train. Therefore, the component of the dynamic response is a component having substantially the same frequency as the frequency of the vibration generated in the second bridge due to the passage of the railway train, and a component having a frequency that is an integral multiple of the frequency. Therefore, the inventors found that even when u(t) includes an influence of the dynamic response, the fundamental frequency of u(t) indicates the frequency of the vibration generated in the second bridge due to the passage of the railway train.

[0120] Therefore, the inventors derived the fundamental frequency F.sub.f of u(t) from the result of the FFT of u(t) as the frequency of the vibration generated in the second bridge in accordance with the movement of the continuous railway vehicles. In an example of FIG. 17, the inventors obtained a value of 2.8 Hz as the fundamental frequency F.sub.f from the graph of FIG. 17.

[0121] The inventors also obtained the cycle T.sub.f of the vibration generated in the second bridge due to the passage of the railway train by obtaining the reciprocal of the obtained fundamental frequency F.sub.f. The inventors performed the low-pass filter processing for attenuating a component of a frequency equal to or higher than the fundamental frequency on u(t) by performing a moving average on u(t), and determined u.sub.lp(t). Specifically, the inventors obtained k.sub.mf from the cycles T.sub.f and ΔT using Equation (43), and obtained u.sub.lp(t) from the obtained k.sub.mf and u(t) using Equation (44). Therefore, the influence of the dynamic response on u(t) is reduced. FIG. 18 shows u.sub.lp(t) superimposed on u(t). In a graph of FIG. 18, a horizontal axis represents time, and a vertical axis represents the deflection amount. A solid line in the graph of FIG. 18 indicates u.sub.lp(t), and a dotted line in the graph indicates u(t). FIG. 18 shows a state in which the influence of the dynamic response is reduced.

[0122] Then, the inventors specified, from u.sub.lp(t), two sets of two consecutive pieces of data between which the predetermined threshold C.sub.L related to the deflection amount is. For each of the two specified sets, the inventors specified a later one of two time points corresponding to the two consecutive pieces of data between which the specified C.sub.L is. The inventors derived the earlier one of the specified time points as the entry time point t.sub.i, and derived the later one as the exit time point t.sub.o. The inventors derived t.sub.i as 5.984 [sec] and to as 12.284 [sec]. The inventors derived to as t.sub.o−t.sub.i=12.284-5.984=6.3 [seconds]. Since the influence of the dynamic response is reduced by using u.sub.lp(t) subjected to the low-pass filter processing, the entry time point t.sub.i and the exit time point t.sub.o can be derived more accurately.

[0123] The inventors derived, using Equation (41), the average velocity v.sub.a of the railway train as 25 [m]×2.8 [Hz]=70 [m/sec] based on L.sub.c(m)=25 [m] and the obtained F.sub.f=2.8 [Hz]. The inventors also derived, using Equation (34), the period t.sub.c(m) during which one railway vehicle C.sub.m passes through the bridge as (46.77+2.5+2.5+15+2.5-2.5)/70=0.954 [seconds] based on L.sub.B=46.77 [m], L.sub.a, and the derived v.sub.a. In addition, the inventors derived, using Equation (42), the number N of railway vehicles formed in the railway train as (6.3−0.954)×2.8+1=15.9692 based on the derived F.sub.f, t.sub.s, and t.sub.c(m). This estimated value is approximated to 16, which is the number N of railway vehicles. Therefore, the inventors confirmed that the number of railway vehicles included in the railway train can be derived with high accuracy by using the conceived method.

[0124] Based on the method conceived by the inventors, the derivation system 10 of the present embodiment derives the number of railway vehicles formed in the railway train 6 based on the time-series data of the displacement at the observation point in the unit bridge girder of the bridge 5.

(1-4) Details of Elements

[0125] Here, each of the measurement device 1, the sensor device 2, and the server device 3 of the derivation system 10 will be described in detail with reference to FIG. 19.

[0126] In the present embodiment, the derivation system 10 derives the observation information (the number N of railway vehicles formed in the railway train 6, the time point t.sub.i at which the railway train 6 enters the unit bridge girder, the time point t.sub.o at which the railway train 6 exits the unit bridge girder, and the period t.sub.s during which the railway train 6 passes through the unit bridge girder) based on data measured by the measurement device 1.

[0127] The measurement device 1 measures the deflection at the observation point via the sensor device 2. In the present embodiment, the measurement device 1 is installed on the bridge abutment 8b, but may be installed at another position. The measurement device 1 includes a control unit 100, a storage unit 110, and a communication unit 120. The control unit 100 includes a processor such as a CPU (Central Processing Unit), a ROM (Read Only Memory), a RAM (Random Access Memory), and the like. The control unit 100 implements each function of the measurement device 1 by loading various programs recorded in the ROM or the like in the RAM and executing the programs via the CPU. The storage unit 110 stores various programs, measured deflection data, and the like. The communication unit 120 includes a circuit used for wired or wireless communication with an external device.

[0128] The sensor device 2 detects the acceleration as the predetermined physical quantity at the observation point. The sensor device 2 includes a control unit 200, an acceleration sensor 210, a storage unit 220, and a communication unit 230. The control unit 200 includes a processor such as a CPU, a ROM, a RAM, and the like. The control unit 200 implements each function of the sensor device 2 by loading various programs recorded in the ROM or the like in the RAM and executing the programs via the CPU.

[0129] The acceleration sensor 210 is an acceleration sensor such as a quartz acceleration sensor or a MEMS acceleration sensor capable of detecting an acceleration generated in each axial direction of three axes orthogonal to one another. In the present embodiment, the acceleration sensor 210 is disposed such that one axis is parallel to the vertical direction in order to more accurately detect the acceleration in the vertical direction. However, an installation location of the sensor device 2 in the superstructure 7 may be inclined. Even when one of the three detection axes of the acceleration sensor 210 is not installed in alignment with the vertical direction, the measurement device 1 combines the accelerations of the three axes and detects the acceleration in the vertical direction.

[0130] The control unit 200 of the sensor device 2 detects an acceleration in a cycle in the vertical direction at the observation point on the bridge 5 via the acceleration sensor 210, and transmits the detected acceleration data to the measurement device 1. The control unit 100 of the measurement device 1 measures the deflection of the bridge 5 in the vertical direction at the observation point at an acceleration detection time point based on the acceleration data transmitted from the sensor device 2. In the present embodiment, the control unit 100 obtains the deflection of the bridge 5 in the vertical direction at the observation point by integrating the acceleration indicated by the data transmitted from the sensor device 2 twice with respect to time. Then, the control unit 100 transmits the measured deflection data to the server device 3. In the present embodiment, the sensor device 2 detects the acceleration in a predetermined cycle ΔT. Therefore, the measurement device 1 measures time-series data of the deflection in the ΔT cycle. That is, the measured time-series data is data of discrete values of displacement measured in the ΔT cycle, and is data in which each discrete value is associated with the measurement time point.

[0131] The server device 3 derives the number of railway vehicles included in the railway train 6 based on the deflection of the observation point measured by the measurement device 1. The server device 3 is an example of a derivation device. The server device 3 includes a control unit 300, a storage unit 310, and a communication unit 320. The control unit 300 includes a processor such as a CPU, a ROM, a RAM, and the like. The control unit 300 implements functions of an acquisition unit 301, an environment information acquisition unit 302, a fundamental frequency derivation unit 303, a passing period derivation unit 304, and a number derivation unit 305, by loading various programs recorded in the ROM or the like into the RAM and executing the programs via the CPU. The storage unit 310 stores various programs, the detected deflection data, and the like. The communication unit 320 includes a circuit used for wired or wireless communication with an external device.

[0132] The acquisition unit 301 has a function of acquiring the time-series data of the deflection generated at the observation point as a response caused by the movement of the railway train 6 on each bridge in the bridge 5. With the function of the acquisition unit 301, the control unit 300 acquires the time-series data u(t) of the deflection generated at the observation point from the measurement device 1.

[0133] The environment information acquisition unit 302 has a function of acquiring environment information including information on the bridge length of the unit bridge girder, the vehicle length of the railway vehicles formed in the railway train 6, and the position of the axle that is the installation position of the wheel in the railway vehicle. With the function of the environment information acquisition unit 302, the control unit 300 acquires information on the bridge length L.sub.B of the unit bridge girder, the vehicle length L.sub.c of each railway vehicle of the railway train 6, and the distance L.sub.a indicating the position of each railway vehicle of the railway train 6 as the environment information. In the present embodiment, the environment information is stored in advance in the storage unit 310, and the control unit 300 acquires the environment information from the storage unit 310. However, the control unit 300 may acquire the environment information by using another method such as receiving the environment information from an external device.

[0134] The fundamental frequency derivation unit 303 has a function of deriving the fundamental frequency F.sub.f of u(t), which indicates the frequency of the vibration generated in the unit bridge girder due to the passage of the railway train 6, based on the time-series data u(t). The control unit 300 executes the FFT on u(t) with the function of the fundamental frequency derivation unit 303. The control unit 300 detects peaks from the FFT result. The control unit 300 specifies, among the detected peaks, a peak corresponding to a minimum frequency obtained by excluding a peak of a side lobe generated due to an influence of a window function used in the FFT. The control unit 300 derives the frequency corresponding to the specified peak as the fundamental frequency F.sub.f of u(t).

[0135] The passing period derivation unit 304 has a function of deriving the passing period t.sub.s during which the railway train 6 passes through the unit bridge girder based on the time-series data u(t). With the function of the passing period derivation unit 304, the control unit 300 derives the reciprocal of the fundamental frequency F.sub.f acquired with the function of the fundamental frequency derivation unit 303, thereby deriving the cycle T.sub.f of the vibration generated in the unit bridge girder caused by the passage of the railway train 6. The control unit 300 derives the interval k.sub.mf using Equation (43) based on the derived T.sub.f and ΔT which is a predetermined cycle. The control unit 300 applies a low-pass filter to u(t) and obtains u.sub.lp(t) which is u(t) subjected to the low-pass filter processing by taking a moving average in the derived interval k.sub.mf for each value of u(t). Specifically, the control unit 300 derives u.sub.lp(t) using Equation (44) based on the derived interval k.sub.mf.

[0136] Then, the control unit 300 specifies, from u.sub.lp(t), two consecutive pieces of data between which a predetermined threshold C.sub.L related to the deflection amount is. The control unit 300 specifies a later one of two time points corresponding to the two consecutive pieces of data between which the specified C.sub.L is.

[0137] The control unit 300 derives the earlier one of the specified time points as the entry time point t.sub.i of the railway train 6 entering into the unit bridge girder. In addition, the control unit 300 derives the later one of the specified time points as the exit time point t.sub.o of the railway train 6 exiting from the unit bridge girder.

[0138] As described above, in the present embodiment, the control unit 300 derives the later one of the two time points corresponding to the two consecutive pieces of data between which C.sub.L is, which are included in u.sub.lp(t), as the entry time point t.sub.i and the exit time point t.sub.o. However, the control unit 300 may derive other time points as the entry time point t.sub.i and the exit time point t.sub.o. For example, the control unit 300 may specify, from u.sub.lp(t), two consecutive pieces of data between which the predetermined threshold C.sub.L related to the deflection amount is, and derive, as the entry time point t.sub.i and the exit time point t.sub.o, a time point that is during a period after one time point of the time points corresponding to the two specified pieces of data and before the other time point. In the example of FIG. 15, the control unit 300 may derive, as the entry time point t.sub.i, a time point after a time point (k−1)ΔT corresponding to the data k−1 and before the time point kΔT corresponding to the data k (for example, time point (k−1)ΔT, a time point corresponding to a point where u.sub.lp(t) and C.sub.L intersect with each other, or the like). In addition, the control unit 300 may obtain a curve obtained by interpolating data included in u.sub.lp(t), and obtain time points corresponding to intersection points of the obtained curve and C.sub.L as t.sub.i and t.sub.o.

[0139] It is conceivable that one of two consecutive pieces of data between which C.sub.L included in u.sub.lp(t) is present is equal to C.sub.L. For example, in the example of FIG. 15, the value of the data k may be equal to C.sub.L. In this case, the control unit 300 may select any one of two sets, that is, a set of data equal to C.sub.L and data preceding the data and a set of data equal to C.sub.L and data following the data, as two consecutive pieces of data between which C.sub.L is. In the example of FIG. 15, when the data k is equal to the C.sub.L, the control unit 300 selects any one of the two sets of the set of the data k−1 and the data k and the set of the data k and the data k+1 as two consecutive pieces of data between which C.sub.L is. The control unit 300 may derive a time point in a period between two time points corresponding to two pieces of data included in the selected set as t.sub.i or t.sub.o.

[0140] In the present embodiment, the control unit 300 derives the time point associated with any data included in u.sub.lp(t) as the entry time point t.sub.i and the exit time point t.sub.o. As a result, the control unit 300 can easily acquire and utilize the data of u.sub.lp(t) corresponding to each measurement time point of ΔT interval including the entry time point t.sub.i and the exit time point t.sub.o by referring to u.sub.lp(t). On the other hand, when deriving the time point not associated with any data included in u.sub.lp(t) as the entry time point t.sub.i and the exit time point t.sub.o, the control unit 300 obtains the data of u.sub.lp(t) corresponding to each measurement time point of the ΔT interval including t.sub.i and to by resampling from the original u.sub.lp(t) or the like, which increases time and effort of processing.

[0141] The control unit 300 derives the entry time point and the exit time point by using u.sub.lp(t) in which a vibration component of a frequency equal to or higher than the fundamental frequency is attenuated, thereby reducing an influence of the vibration component of a frequency equal to or higher than the fundamental frequency and more accurately deriving the entry time point and the exit time point.

[0142] However, the control unit 300 may not derive u.sub.lp(t). In this case, the control unit 300 may derive t.sub.i and to using u(t) instead of u.sub.lp(t), for example.

[0143] Then, the control unit 300 derives the passing period t.sub.s as t.sub.o−t.sub.i.

[0144] The number derivation unit 305 has a function of deriving the number of railway vehicles included in the railway train 6 based on the environmental information acquired with the function of the environmental information acquisition unit 302, the fundamental frequency F.sub.f derived with the function of the fundamental frequency derivation unit 303, and the passing period t.sub.s derived with the function of the passing period derivation unit 304.

[0145] With the function of the number derivation unit 305, the control unit 300 derives the average velocity v.sub.a of the railway train using Equation (41) based on the fundamental frequency F.sub.f, and the vehicle length L.sub.c(m) of the railway vehicles of the railway train 6 indicated by the environmental information. In the present embodiment, since the railway train 6 is formed with the same railway vehicle, L.sub.c(m) has the same value for each of m=1 to N.

[0146] Based on the derived v.sub.a and L.sub.B and L.sub.a that are indicated by the environment information, the control unit 300 derives the period t.sub.c(m) during which one railway vehicle passes through the bridge using Equation (34). Then, based on the derived F.sub.f, t.sub.s, and t.sub.c(m), the control unit 300 derives the number N of railway vehicles formed in the railway train using Equation (42).

[0147] As described above, according to the configuration of the present embodiment, the derivation system 10 can derive the number of railway vehicles based on the time-series data of the deflection at the observation point and the environmental information. In addition, the derivation system 10 derives the fundamental frequency F.sub.f from the result of the FFT of u(t), derives the average velocity v.sub.a of the railway train using Equation (41) based on the vehicle length L.sub.c(m) and F.sub.f, derives the period t.sub.c(m) during which one railway vehicle passes through the bridge using Equation (34) based on the derived v.sub.a, L.sub.B, and L.sub.a, and derives the number N of railway vehicles formed in the railway train using Equation (42) based on the derived F.sub.f, t.sub.s and t.sub.c (m). As described above, the derivation system 10 can obtain the number of railway vehicles formed in the railway train 6 with a smaller amount of calculation and a lower load than in the case of obtaining the number of railway vehicles formed in the railway train 6 with an inverse analysis method. Thus, the derivation system 10 can obtain the number of railway vehicles formed in the railway train 6 with a lower load.

[0148] In addition, since the number of railway vehicles is derived with higher accuracy by the derivation system 10, it is possible to diagnose the bridge with higher accuracy using the derived number. For example, based on the number N derived in the processing of the present embodiment, T.sub.std(t) represented by Equation (32) is derived using the deflection model. Since the number N is derived with higher accuracy, the deflection amount T.sub.std(t) is also derived with higher accuracy. In addition, by using the deflection amount T.sub.std(t) derived in this manner, it is possible to calculate an impact coefficient of impact on the bridge with high accuracy. By using such an impact coefficient, it is possible to diagnose the bridge with high accuracy. Therefore, in dynamic design including risk factors including the resonance phenomenon of the railroad bridge, accurate measurement and derivation of the number of vehicles are effective.

(2) Derivation Processing

[0149] Processing of deriving the number of railway vehicles of the railway train 6 executed by the server device 3 will be described with reference to FIG. 20. The server device 3 starts processing in FIG. 20 in response to the fact that the data of the displacement at the observation point is transmitted from the measurement device 1, but may start the processing in FIG. 20 at any timing such as a designated timing.

[0150] In S100, the control unit 300 acquires the time-series data u(t) of the deflection generated at the observation point from the measurement device 1 by the function of the acquisition unit 301. S100 is an example of an acquisition step.

[0151] In S105, with the function of the environment information acquisition unit 302, the control unit 300 acquires information on the bridge length L.sub.B of the unit bridge girder, the vehicle length L.sub.c of each railway vehicle of the railway train 6, and the distance L.sub.a indicating the position of each railway vehicle of the railway train 6 as the environment information. S105 is an example of an environment information acquisition step.

[0152] In S110, the control unit 300 executes the FFT on u(t) acquired in S100 with the function of the fundamental frequency derivation unit 303. The control unit 300 detects peaks from the FFT result. The control unit 300 specifies, among the detected peaks, a peak corresponding to a minimum frequency obtained by excluding a peak of a side lobe generated due to an influence of a window function used in the FFT. The control unit 300 derives the frequency corresponding to the specified peak as the fundamental frequency F.sub.f of u(t). S110 is an example of a fundamental frequency derivation step.

[0153] In S115, the control unit 300 derives the cycle T.sub.f by deriving the reciprocal of the fundamental frequency F.sub.f with the function of the passing period derivation unit 304. The control unit 300 derives k.sub.mf using Equation (43) based on the derived T.sub.f and ΔT which is a predetermined cycle. Then, the control unit 300 derives u.sub.lp(t) from the derived k.sub.mf and u(t) using Equation (44). The control unit 300 specifies an intersection point between the derived u.sub.lp(t) and the predetermined threshold C.sub.L related to the deflection amount, that is, two points at which u.sub.lp(t)=C.sub.L, derives the earlier one of the time points corresponding to the specified points as the entry time point t.sub.i, and derives the later one as the exit time point t.sub.o. Then, the control unit 300 derives the passing period t.sub.s as t.sub.o−t.sub.i. S115 is an example of a passing period derivation step.

[0154] In S120, the control unit 300 derives the average velocity v.sub.a of the railway train using Equation (41) based on the fundamental frequency F.sub.f and the vehicle length L.sub.c (m) of the railway vehicles of the railway train 6 indicated by the environmental information with the function of the number derivation unit 305. Based on the derived v.sub.a and L.sub.B and L.sub.a that are indicated by the environment information, the control unit 300 derives the period t.sub.c(m) during which one railway vehicle passes through the bridge using Equation (34). Then, based on the derived F.sub.f, t.sub.s, and t.sub.c(m), the control unit 300 derives the number N of railway vehicles formed in the railway train using Equation (42). S120 is an example of a number derivation step. Processing of deriving t.sub.c(m) in S120 is an example of a moving object passing period derivation step.

(3) Other Embodiments

[0155] The above embodiments are examples for carrying out the present disclosure, and various other embodiments can be adopted. The method of deriving the number of railway vehicles of the railway train from the displacement at the observation point as in the above embodiment can also be implemented as an invention of a program and an invention of a method.

[0156] Further, a configuration in which the function of the server device 3 is implemented by a plurality of devices may be adopted. The functions of the server device 3 may be distributed and implemented in a plurality of devices. In addition, the functions of the server device 3 may be implemented in another device. For example, the functions of the acquisition unit 301, the environment information acquisition unit 302, the fundamental frequency derivation unit 303, the passing period derivation unit 304, and the number derivation unit 305 may be implemented in the measurement device 1. The server device 3 may be distributed in a plurality of devices, or the like. Further, the above embodiments are examples, and an embodiment in which a part of the configuration is omitted or another configuration is added may be adopted.

[0157] In the above embodiment, the derivation system 10 derives the number of railway vehicles included in the railway train 6 formed with one or more railway vehicles that are moving objects. However, the derivation system 10 may derive the number of moving objects included in another formation moving object. The derivation system 10 may derive the number of moving objects included in another formation moving object. For example, the derivation system 10 may derive the number of trucks included in a formation truck in which one or more trucks are coupled.

[0158] Further, in the above embodiment, the derivation system 10 derives the number of moving objects included in the formation moving object that moves on the bridge 5. However, the derivation system 10 may derive the number of moving objects included in the formation moving object that moves on a structure different from the bridge, such as a base that supports a railroad track.

[0159] In addition, in the above embodiment, the number of sensor devices 2 included in the derivation system 10 is two, but may be one or three or more.

[0160] In the above embodiment, the control unit 300 acquires, as the time-series data u(t), the data of displacement (deflection) measured from the acceleration detected via the acceleration sensor 210. However, the control unit 300 may acquire, as u(t), displacement data of the bridge derived from a physical quantity detected via a sensor such as an impact sensor, a pressure-sensitive sensor, a strain gauge, an image measuring device, a load cell, or a displacement gauge. For example, the control unit 300 may detect the displacement of the observation point and acquire the detected displacement data by the image measuring device capturing an image of a predetermined object in a cycle, the predetermined object being disposed at the observation point of the bridge 5. The control unit 300 may acquire data of a physical quantity different from the displacement of the bridge as u(t). For example, the control unit 300 may acquire, as u(t), the number of pixels indicating the displacement amount of the predetermined object disposed at the observation point of the bridge 5 in the image captured via the image measuring device.

[0161] Further, in the above embodiment, the control unit 300 specifies the peak corresponding to the lowest frequency, except for the side lobe generated due to the influence of the window function used in the FFT, from the result of the FFT on the time-series data u(t) acquired by the function of the acquisition unit 301, and obtains the specified peak as the fundamental frequency F.sub.f. However, in consideration of an influence of a noise generated in the result of the FFT on u(t), the control unit 300 may obtain the fundamental frequency F.sub.f. For example, the control unit 300 may specify a peak equal to or greater than a predetermined threshold corresponding to the lowest frequency, except for the side lobe generated due to the influence of the window function used in the FFT, from the result of the FFT on u(t), and obtain the specified peak as the fundamental frequency F.sub.f.

[0162] The time-series data may be data acquired at a data rate of twice or more the frequency of vibration assumed to occur in the structure due to the movement of the formation moving object.

[0163] Further, the present disclosure can also be applied as a program executed by a computer or a method. In addition, the program and method as described above may be implemented as a single device or may be implemented by using components included in a plurality of devices, and includes various aspects. In addition, it is possible to appropriately change the configuration such that a part of the configuration is software and a part of the configuration is hardware. Further, the present disclosure is also applicable to a recording medium of a program. As a matter of course, the recording medium of the program may be a magnetic recording medium, a semiconductor memory, or the like, and any recording medium to be developed in the future can be considered in the same manner.