INSPECTION SYSTEM, INSPECTION METHOD, AND PROGRAM

Description

TECHNICAL FIELD

[0001] The present invention relates to an inspection system, an inspection method, and a program, and in particular, is ones to be suitable when used for inspecting a track of a railway vehicle.

BACKGROUND ART

[0002] When a railway vehicle travels on a track, the position of the track changes due to a load from the railway vehicle. When such a change in track occurs, the railway vehicle is liable to exhibit abnormal behavior. Thus, detection of abnormality of the track has been performed conventionally by causing a railway vehicle to travel on the track.

[0003] Patent Literature 1 describes that angular displacements of wheel sets in a yawing direction, state variables derived by a Kalman filter, and a forward-and-backward-direction force are substituted into motion equations that describe yawings of the wheel sets to estimate an alignment irregularity amount.

CITATION LIST

PATENT LITERATURE

[0004] 

Patent Literature 1: International Publication Pamphlet No. WO 2017/164133

Patent Literature 2: Japanese Laid-open Patent Publication No. 2017-53773

SUMMARY OF INVENTION

TECHNICAL PROBLEM

[0005] The present inventors have learned that in the technique described in Patent Literature 1, when a disturbance that is not considered in the motion equations occurs, the error in an estimated value of the alignment irregularity amount becomes large.

[0006] The present invention has been made in consideration of the above problem, and an object thereof is to enable accurate detection of track irregularity of a railway vehicle without using a special measuring apparatus.

SOLUTION TO PROBLEM

[0007] An inspection system of the present invention includes: a data acquisition means configured to acquire measured data being time-series data of measured values to be measured by causing a railway vehicle including a vehicle body, a bogie, and a wheel set to travel on a track; a first track state calculation means configured to calculate an estimated value of a first physical quantity; a correction amount calculation means configured to calculate a correction amount for the estimated value of the first physical quantity based on the estimated value of the first physical quantity calculated by the first track state calculation means and an actual value of the first physical quantity;
a second track state calculation means configured to calculate an estimated value of the first physical quantity after the correction amount is calculated; and a track state correction means configured to correct the estimated value of the first physical quantity calculated by the second track state calculation means by using the correction amount, in which the measured data contain a measured value of a forward-and-backward-direction force, the forward-and-backward-direction force is a force in a forward and backward direction to occur in a member disposed between the wheel set and the bogie on which the wheel set is provided, the member is a member for supporting an axle box, the forward and backward direction is a direction along a traveling direction of the railway vehicle, the first physical quantity is a physical quantity reflecting a state of the track, the first track state calculation means and the second track state calculation means are configured to use a relational expression representing the relation between the first physical quantity at a position of the wheel set and the forward-and-backward-direction force and a measured value of the forward-and-backward-direction force to calculate the estimated value of the first physical quantity, the measured value of the forward-and-backward-direction force used in the first track state calculation means is contained in the measured data acquired by the data acquisition means before the correction amount is calculated, and the measured value of the forward-and-backward-direction force used in the second track state calculation means is contained in the measured data acquired by the data acquisition means after the correction amount is calculated.

[0008] An inspection method of the present invention, carried out by an inspection system, includes: a data acquisition step of acquiring measured data being time-series data of measured values to be measured by causing a railway vehicle including a vehicle body, a bogie, and a wheel set to travel on a track; a first track state calculation step of calculating an estimated value of a first physical quantity; a correction amount calculation step of calculating a correction amount for the estimated value of the first physical quantity based on the estimated value of the first physical quantity calculated by the first track state calculation step and an actual value of the first physical quantity; a second track state calculation step of calculating an estimated value of the first physical quantity after the correction amount is calculated; and a track state correction step of correcting the estimated value of the first physical quantity calculated by the second track state calculation step by using the correction amount, in which the measured data contain a measured value of a forward-and-backward-direction force, the forward-and-backward-direction force is a force in a forward and backward direction to occur in a member disposed between the wheel set and the bogie on which the wheel set is provided, the member is a member for supporting an axle box, the forward and backward direction is a direction along a traveling direction of the railway vehicle, the first physical quantity is a physical quantity reflecting a state of the track, the first track state calculation step and the second track state calculation step use a relational expression representing the relation between the first physical quantity at a position of the wheel set and the forward-and-backward-direction force and a measured value of the forward-and-backward-direction force to calculate the estimated value of the first physical quantity, the measured value of the forward-and-backward-direction force used in the first track state calculation step is contained in the measured data acquired by the data acquisition step before the correction amount is calculated, and the measured value of the forward-and-backward-direction force used in the second track state calculation step is contained in the measured data acquired by the data acquisition step after the correction amount is calculated.

[0009] A computer program of the present invention comprising instructions which, when the program is executed by a computer, cause the computer to execute steps including: a data acquisition step of acquiring measured data being time-series data of measured values to be measured by causing a railway vehicle including a vehicle body, a bogie, and a wheel set to travel on a track; a first track state calculation step of calculating an estimated value of a first physical quantity; a correction amount calculation step of calculating a correction amount for the estimated value of the first physical quantity based on the estimated value of the first physical quantity calculated by the first track state calculation step and an actual value of the first physical quantity; a second track state calculation step of calculating an estimated value of the first physical quantity after the correction amount is calculated; and a track state correction step of correcting the estimated value of the first physical quantity calculated by the second track state calculation step by using the correction amount, in which the measured data contain a measured value of a forward-and-backward-direction force, the forward-and-backward-direction force is a force in a forward and backward direction to occur in a member disposed between the wheel set and the bogie on which the wheel set is provided, the member is a member for supporting an axle box, the forward and backward direction is a direction along a traveling direction of the railway vehicle, the first physical quantity is a physical quantity reflecting a state of the track, the first track state calculation step and the second track state calculation step use a relational expression representing the relation between the first physical quantity at a position of the wheel set and the forward-and-backward-direction force and a measured value of the forward-and-backward-direction force to calculate the estimated value of the first physical quantity, the measured value of the forward-and-backward-direction force used in the first track state calculation step is contained in the measured data acquired by the data acquisition step before the correction amount is calculated, and the measured value of the forward-and-backward-direction force used in the second track state calculation step is contained in the measured data acquired by the data acquisition step after the correction amount is calculated.

BRIEF DESCRIPTION OF DRAWINGS

[0010] 

[Fig. 1] Fig. 1 is a view illustrating one example of an outline of a railway vehicle.

[Fig. 2] Fig. 2 is a view conceptually illustrating directions of main motions of components of the railway vehicle.

[Fig. 3A] Fig. 3A is a view illustrating one example of an alignment irregularity amount in a linear track.

[Fig. 3B] Fig. 3B is a view illustrating one example of an alignment irregularity amount in a curved track.

[Fig. 4] Fig. 4 is a diagram illustrating one example of a functional configuration of an inspection apparatus.

[Fig. 5] Fig. 5 is a diagram illustrating one example of a hardware configuration of the inspection apparatus.

[Fig. 6] Fig. 6 is a flowchart illustrating one example of first preprocessing.

[Fig. 7] Fig. 7 is a flowchart illustrating one example of second preprocessing.

[Fig. 8] Fig. 8 is a flowchart illustrating one example of main processing.

[Fig. 9] Fig. 9 is a view illustrating one example of a distribution of eigenvalues of an autocorrelation matrix.

[Fig. 10] Fig. 10 is a view illustrating one example of time-series data of a measured value of a forward-and-backward-direction force (a measured value) and time-series data of a predicted value of the forward-and-backward-direction force (a calculated value).

[Fig. 11] Fig. 11 is a view illustrating one example of time-series data of a high-frequency component of the forward-and-backward-direction force.

[Fig. 12A] Fig. 12A is a view illustrating a first example of a relation between an estimated value of the alignment irregularity amount, an actual value of the alignment irregularity amount, a traveling velocity of the railway vehicle, and a curvature of a rail and a distance from a starting point of the railway vehicle.

[Fig. 12B] Fig. 12B is a view illustrating a second example of the relation between the estimated value of the alignment irregularity amount, the actual value of the alignment irregularity amount, the traveling velocity of the railway vehicle, and the curvature of the rail and the distance from the starting point of the railway vehicle.

[Fig. 13A] Fig. 13A is a view illustrating a third example of the relation between the estimated value of the alignment irregularity amount, the actual value of the alignment irregularity amount, the traveling velocity of the railway vehicle, and the curvature of the rail and the distance from the starting point of the railway vehicle.

[Fig. 13B] Fig. 13B is a view illustrating a fourth example of the relation between the estimated value of the alignment irregularity amount, the actual value of the alignment irregularity amount, the traveling velocity of the railway vehicle, and the curvature of the rail and the distance from the starting point of the railway vehicle.

[Fig. 14A] Fig. 14A is a view illustrating a fifth example of the relation between the estimated value of the alignment irregularity amount, the actual value of the alignment irregularity amount, the traveling velocity of the railway vehicle, and the curvature of the rail and the distance from the starting point of the railway vehicle.

[Fig. 14B] Fig. 14B is a view illustrating a sixth example of the relation between the estimated value of the alignment irregularity amount, the actual value of the alignment irregularity amount, the traveling velocity of the railway vehicle, and the curvature of the rail and the distance from the starting point of the railway vehicle.

[Fig. 15] Fig. 15 is a view explaining one example of flange contact.

[Fig. 16A] Fig. 16A is a view illustrating a first example of a relation between a second correction amount and the distance from the starting point of the railway vehicle.

[Fig. 16B] Fig. 16B is a view illustrating a second example of the relation between the second correction amount and the distance from the starting point of the railway vehicle.

[Fig. 16C] Fig. 16C is a view illustrating a third example of the relation between the second correction amount and the distance from the starting point of the railway vehicle.

[Fig. 17A] Fig. 17A is a view illustrating a first example of a relation between a corrected estimated value of the alignment irregularity amount and the distance from the starting point of the railway vehicle.

[Fig. 17B] Fig. 17B is a view illustrating a second example of the relation between the corrected estimated value of the alignment irregularity amount and the distance from the starting point of the railway vehicle.

[Fig. 18A] Fig. 18A is a view illustrating a third example of the relation between the corrected estimated value of the alignment irregularity amount and the distance from the starting point of the railway vehicle.

[Fig. 18B] Fig. 18B is a view illustrating a fourth example of the relation between the corrected estimated value of the alignment irregularity amount and the distance from the starting point of the railway vehicle.

[Fig. 19A] Fig. 19A is a view illustrating a fifth example of the relation between the corrected estimated value of the alignment irregularity amount and the distance from the starting point of the railway vehicle.

[Fig. 19B] Fig. 19B is a view illustrating a sixth example of the relation between the corrected estimated value of the alignment irregularity amount and the distance from the starting point of the railway vehicle.

[Fig. 20] Fig. 20 is a view illustrating one example of a constitution of an inspection system.

DESCRIPTION OF EMBODIMENTS

[0011] Hereinafter, there will be explained embodiments of the present invention with reference to the drawings.

(Outline)

[0012] First, there will be explained the outline of the embodiment of the present invention.

[0013] Fig. 1 is a view illustrating one example of an outline of a railway vehicle. Incidentally, in Fig. 1, the railway vehicle is set to proceed in the positive direction of the x axis (the x axis is an axis along a traveling direction of the railway vehicle). Further, the z axis is set to a direction perpendicular to a track 16 (the ground) (a height direction of the railway vehicle). The y axis is set to a horizontal direction perpendicular to the traveling direction of the railway vehicle (a direction perpendicular to both the traveling direction and the height direction of the railway vehicle). Further, the railway vehicle is set to a commercial vehicle. Incidentally, in each of the drawings, the mark of ● added inside O indicates the direction from the far side of the sheet toward the near side, and the mark of × added inside ○ indicates the direction from the near side of the sheet toward the far side.

[0014] As illustrated in Fig. 1, in this embodiment, the railway vehicle includes a vehicle body 11, bogies 12a, 12b, and wheel sets 13a to 13d. As above, in this embodiment, the railway vehicle including the single vehicle body 11 provided with the two bogies 12a, 12b and four sets of the wheel sets 13a to 13d will be explained as an example. The wheel sets 13a to 13d have axles 15a to 15d and wheels 14a to 14d provided on both ends of the axles 15a to 15d respectively. In this embodiment, the case of the bogies 12a, 12b each being a bolsterless bogie will be explained as an example. Incidentally, in Fig. 1, for convenience of illustration, only the wheels 14a to 14d on one side of the wheel sets 13a to 13d are illustrated, but wheels are also provided on the other side of the wheel sets 13a to 13d (in the example illustrated in Fig. 1, there are eight wheels in total). Further, the railway vehicle includes components other than the components illustrated in Fig. 1 (components and so on to be explained in later-described motion equations), but for convenience of illustration, illustrations of these components are omitted in Fig. 1. For example, the bogies 12a, 12b have bogie frames, bolster springs, and so on. Further, an axle box is disposed on the both sides of each of the wheel sets 13a to 13d in the direction along the y axis. Further, the bogie frame and the axle box are coupled to each other by an axle box suspension. The axle box suspension is a device (suspension) to be disposed between the axle box and the bogie frame. The axle box suspension absorbs vibration to be conveyed to the railway vehicle from the track 16. Further, the axle box suspension supports the axle box in a state where the position of the axle box relative to the bogie frame is restricted, so as to prevent the axle box from moving in a direction along the x axis and a direction along the y axis relative to the bogie frame (so as to prevent these movements from occurring preferably). The axle box suspension is disposed on the both sides of each of the wheel sets 13a to 13d in the direction along the y axis. Incidentally, the railway vehicle itself can be fabricated by a well-known technique, and thus its detailed explanation is omitted here.

[0015] When the railway vehicle travels on the track 16, acting force (creep force) between the wheels 14a to 14d and the track 16 becomes a vibration source and the vibration sequentially propagates to the wheel sets 13a to 13d, the bogies 12a, 12b, and the vehicle body 11. Fig. 2 is a view conceptually illustrating directions of the main motions of the components (the wheel sets 13a to 13d, the bogies 12a, 12b, and the vehicle body 11) of the railway vehicle. The x axis, the y axis, and the z axis illustrated in Fig. 2 correspond to the x axis, the y axis, and the z axis illustrated in Fig. 1 respectively.

[0016] In this embodiment as illustrated in Fig. 2, there will be explained, as an example, the case where the wheel sets 13a to 13d, the bogies 12a, 12b, and the vehicle body 11 perform pivoting motion about the x axis as a pivot axis, pivoting motion about the z axis as a pivot axis, and motion in the direction along the y axis. In the following explanation, the pivoting motion about the x axis as a pivot axis is referred to as rolling as necessary, the pivoting direction about the x axis as a pivot axis is referred to as a rolling direction as necessary, and the direction along the x axis is referred to as the forward and backward direction as necessary. Incidentally, the forward and backward direction is the traveling direction of the railway vehicle. In this embodiment, the direction along the x axis is the traveling direction of the railway vehicle. Further, the pivoting motion about the z axis as a pivot axis is referred to as yawing as necessary, the pivoting direction about the z axis as a pivot axis is referred to as a yawing direction as necessary, and the direction along the z axis is referred to as the up and down direction as necessary. Incidentally, the up and down direction is a direction perpendicular to the track 16. Further, the motion in the direction along the y axis is referred to as a transversal vibration as necessary, and the direction along the y axis is referred to as the right and left direction as necessary. Incidentally, the right and left direction is a direction perpendicular to both the forward and backward direction (the traveling direction of the railway vehicle) and the up and down direction (the direction perpendicular to the track 16).
Further, the railway vehicle performs motions other than these, but in each of the embodiments, these motions are not considered in order to simplify the explanation. However, these motions may be considered.

[0017] As described in Patent Literature 1, the present inventors devised a method of calculating, as one example of a first physical quantity reflecting track irregularity (appearance failure of the track 16), an alignment irregularity amount by using a measured value of force in the forward and backward direction to occur in a member disposed between the wheel sets 13a to 13b (13c to 13d) and the bogie 12a (12b) on which these wheel sets 13a to 13b (13c to 13d) are provided. In the following explanation, the force in the forward and backward direction to occur in the member is referred to as a forward-and-backward-direction force as necessary.

[0018] The alignment irregularity amount is calculated by using an equation representing the relation between the alignment irregularity amount and the forward-and-backward-direction force, which is an equation based on a motion equation describing motion when the railway vehicle travels on a linear track. The track 16 includes a linear portion and a curved portion. In the following explanation, the linear portion of the track 16 is referred to as a linear track as necessary and the curved portion of the track 16 is referred to as a curved track as necessary.

[0019] When a state equation is constituted by using a motion equation that describes motion of the railway vehicle traveling on the curved track in the case of performing filtering with a filter (Kalman filter) performing later-described data assimilation, state variables may diverge. Therefore, the state equation in the case of performing filtering with the filter (Kalman filter) to perform data assimilation is constituted by using a motion equation that describes motion of the railway vehicle traveling on the linear track.

[0020] It is necessary to consider centrifugal force or the like that the railway vehicle receives when traveling in the motion equation describing the motion of the railway vehicle traveling on the curved track. Accordingly, the motion equation describing the motion of the railway vehicle traveling on the curved track includes a term containing a curvature radius of the rail. Therefore, when the state variables are derived using the filter (Kalman filter) that performs data assimilation constituted by using the motion equation that describes the motion of the railway vehicle traveling on the linear track when the railway vehicle is traveling on the curved track, there is a risk that it becomes impossible to derive the state variables with high accuracy.

[0021] The present inventors focused attention on the fact that the measured value of the forward-and-backward-direction force when the railway vehicle travels on the curved track has a certain bias relative to that when traveling on the linear track. Thus, the present inventors thought that by reducing a low-frequency component (behavior of the aforementioned bias) from time-series data of the measured value of the forward-and-backward-direction force, the low-frequency component due to the railway vehicle traveling on the curved track can be reduced from an estimated value of the state variable even where the filter (Kalman filter) performing the later-described data assimilation is constituted by using an equation based on the motion equation that describes the motion of the railway vehicle when traveling on the linear track. From this, the present inventors devised calculating the alignment irregularity amount by giving the time-series data of the value of the forward-and-backward-direction force from which the low-frequency component has been reduced to an equation representing the relation between the alignment irregularity amount and the forward-and-backward-direction force, which is an equation based on the motion equation describing the motion of the railway vehicle when traveling on the linear track. The alignment irregularity amount is calculated as above, thereby making it possible to calculate the alignment irregularity amount in the curved track regardless of using the equation based on the motion equation describing the motion of the railway vehicle when traveling on the linear track. Further, the calculating equation of the alignment irregularity amount results in the same calculating equation even on the curved track or the linear track.

[0022] Further, the present inventors found out that depending on at least any of a traveling state of the railway vehicle and an installation state of the track 16, the calculation accuracy of the alignment irregularity amount may decrease because the disturbance not considered in the motion equations describing the motions of the railway vehicle affects the measured value of the forward-and-backward-direction force. Examples of such a traveling state of the railway vehicle where the disturbance is likely to occur include a state where the railway vehicle is traveling at low speed, a state where the railway vehicle is decelerating rapidly, a state where the railway vehicle is accelerating rapidly, a state where the railway vehicle is traveling in contact with flanges, and a state where the railway vehicle is traveling on a seam of the rail. Further, examples of such an installation state of the track 16 where the disturbance is likely to occur include a state where the rail has a sharp curve (state where the rail has a large curvature), a state where the track 16 is installed at a place of a specific structure, a state where the rail has a seam, and a state where the track 16 is a ballastless track. Examples of the specific structure include station platforms, bridges, tunnels, turnouts, railroad crossings, and guardrails.

[0023] Such a disturbance is represented by the difference between an estimated value and an actual measured value of the alignment irregularity amount. In the case of the same railway vehicle, measured data do not vary significantly due to characteristics inherent in the railway vehicle. Examples of the characteristics inherent in the railway vehicle include individual differences in the vehicle body 11, individual differences in the bogies 12a, 12b, individual differences in the wheel sets 13a to 13d, and individual differences in strain gauges that measure the forward-and-backward-direction force. Further, connection states of these are cited as the characteristics inherent in the railway vehicle. Further, in the case of the same railway vehicle, a traveling velocity at each position of the track 16 does not vary significantly. From this, the present inventors found out that the previously-described difference between the estimated value and the actual measured value of the alignment irregularity amount does not vary significantly depending on the traveling date and time of the railway vehicle in the case where the same railway vehicle is traveling on the same position. Thus, the previously-described difference between the estimated value and the actual measured value of the alignment irregularity amount is derived beforehand as a correction amount for the estimated value of the alignment irregularity amount at each position of the track 16 on which the railway vehicle travels. Thereafter, the railway vehicle is made to travel on the track 16, to thereby obtain the estimated value of the alignment irregularity amount again at each position of the track 16. The estimated value of the alignment irregularity amount obtained as above is corrected by the correction amount at the position of the track 16 where the estimated value has been obtained. In this manner, the alignment irregularity amount at each position of the track 16 is obtained. In this embodiment, the alignment irregularity amount after correction is set as a final alignment irregularity amount.

(Motion equation)

[0024] Next, there will be explained one example of the motion equation that describes the motion when the railway vehicle travels on the linear track. In this embodiment, there will be explained, as an example, the case where the railway vehicle has 21 degrees of freedom while taking the motion equations described in Patent Literature 1 as an example. That is, it is set that the wheel sets 13a to 13d perform the motion in the right and left direction (transversal vibration) and the motion in the yawing direction (yawing) (2 × 4 sets = eight degrees of freedom). Further, it is set that the bogies 12a, 12b perform the motion in the right and left direction (transversal vibration), the motion in the yawing direction (yawing), and the motion in the rolling direction (rolling) (3 × 2 sets = six degrees of freedom). Further, the vehicle body 11 performs the motion in the right and left direction (transversal vibration), the motion in the yawing direction (yawing), and the motion in the rolling direction (rolling) (3 × 1 sets = three degrees of freedom). Further, it is set that air springs (the bolster springs) each provided on the bogies 12a, 12b perform the motion in the rolling direction (rolling) (1 × 2 sets = two degrees of freedom). Further, it is set that yaw dampers each provided on the bogies 12a, 12b perform the motion in the yawing direction (yawing) (1 × 2 sets = two degrees of freedom).

[0025] Incidentally, the degree of freedom is not limited to 21 degrees of freedom. When the degree of freedom increases, calculation accuracy improves, but a calculation load becomes high. Further, there is a risk that a later-described Kalman filter no longer operates stably. It is possible to appropriately determine the degree of freedom considering these points. Further, the following motion equations can be achieved by representing actions in the respective directions (the right and left direction, the yawing direction, and the rolling direction) of the respective components (the vehicle body 11, the bogies 12a, 12b, and the wheel sets 13a to 13d) based on the descriptions of Patent Literature 1, for example. Thus, outlines of these motion equations will be explained here, and their detailed explanations are omitted. Incidentally, in each of the following equations, the term containing the curvature radius (curvature) of the track 16 (rail) does not exist. That is, each of the following equations is an equation expressing the railway vehicle traveling on the linear track. The equation expressing the railway vehicle traveling on the linear track can be obtained by setting the curvature radius of the track 16 (rail) to be infinite (the curvature to 0 (zero)) in the equation expressing the railway vehicle traveling on the curved track.

[0026] In each of the following equations, each subscript w indicates the wheel sets 13a to 13d. Variables to which (only) the subscript w is added indicate that they are common to the wheel sets 13a to 13d. Subscripts w1, w2, w3, and w4 indicate the wheel sets 13a, 13b, 13c, and 13d respectively.

[0027] Subscripts t, T indicate the bogies 12a, 12b. Variables to which (only) the subscripts t, T are added indicate that they are common to the bogies 12a, 12b. Subscripts t1, t2 indicate the bogies 12a, 12b respectively.

[0028] Subscripts b, B indicate the vehicle body 11.

[0029] A subscript x indicates the forward and backward direction or the rolling direction, and a subscript y indicates the right and left direction, and a subscript z indicates the up and down direction or the yawing direction.

[0030] Further, " · · " and " · " each added above a variable indicate a second-order time differential and a first-order time differential respectively.

[0031] Incidentally, when the following motion equations are explained, explanations of the already-explained variables are omitted as necessary. Further, the motion equations themselves are the same as those described in Patent Literature 1.

[0032] [Transversal vibration of the wheel set]

[0033] The motion equations that describe the transversal vibrations of the wheel sets 13a to 13d (motion in the right and left direction) are expressed by (1) Equation to (4) Equation below.
[Mathematical equation 1]

[0034] m_w is the mass of the wheel sets 13a to 13d. y_w1 · · is acceleration of the wheel set 13a in the right and left direction (in the equation, · · is added above y_w1 (the same is true of the other variables below)). f₂ is a lateral creep coefficient (incidentally, the lateral creep coefficient f₂ may be given for each of the wheel sets 13a to 13d). v is a traveling velocity of the railway vehicle. y_w1 · is a velocity of the wheel set 13a in the right and left direction (in the equation, · is added above y_w1 (the same is true of the other variables below)). C_wy is a damping constant of the axle box suspension coupling the axle box and the wheel set in the right and left direction. y_t1 · is a velocity of the bogie 12a in the right and left direction. a represents 1/2 of each distance between the wheel sets 13a and 13b and between the wheel sets 13c and 13d in the forward and backward direction, which are provided on the bogies 12a, 12b (the distance between the wheel sets 13a and 13b and the distance between the wheel sets 13c and 13d, which are provided on the bogies 12a, 12b, each become 2a). φ_t1 · is an angular velocity of the bogie 12a in the yawing direction. h₁ is a distance between the middle of the axle and the center of gravity of the bogie 12a in the up and down direction. φ_t1 · is an angular velocity of the bogie 12a in the rolling direction. φ _w1 is a pivot amount (angular displacement) of the wheel set 13a in the yawing direction. K_wy is a spring constant of the axle box suspension in the right and left direction. y_w1 is a displacement of the wheel set 13a in the right and left direction. y_t1 is a displacement of the bogie 12a in the right and left direction. φ_t1 is a pivot amount (angular displacement) of the bogie 12a in the yawing direction. φ_t1 is a pivot amount (angular displacement) of the bogie 12a in the rolling direction. Incidentally, respective variables in (2) Equation to (4) Equation are represented by being replaced with the variables in (1) Equation according to the meanings of the aforementioned subscripts.

[Yawing of the wheel set]

[0035] The motion equations that describe the yawings of the wheel sets 13a to 13d are expressed by (5) Equation to (8) Equation below.
[Mathematical equation 2]

[0036] I_wz is a moment of inertia of the wheel sets 13a to 13d in the yawing direction. φ_w1 · · is angular acceleration of the wheel set 13a in the yawing direction. f₁ is a longitudinal creep coefficient. b is a distance in the right and left direction between contacts between the two wheels, which are attached to each of the wheel sets 13a to 13d, and the track 16 (rail). φ_w1 · is an angular velocity of the wheel set 13a in the yawing direction. C_wx is a damping constant of the axle box suspension in the forward and backward direction. b₁ represents the length of 1/2 of the interval between the axle box suspensions in the right and left direction (the interval of the two axle box suspensions, which are provided on the right and left sides of the single wheel set, in the right and left direction becomes 2b₁). γ is a tread slope. r is a radius of the wheels 14a to 14d. y_R1 is an alignment irregularity amount at the position of the wheel set 13a. s_a is an offset from the middle of the axles 15a to 15d to an axle box suspension spring in the forward and backward direction. y_t1 is a displacement of the bogie 12a in the right and left direction. K_wx is a spring constant of the axle box suspension in the forward and backward direction. Incidentally, respective variables in (6) Equation to (8) Equation are represented by being replaced with the variables in (5) Equation according to the meanings of the aforementioned subscripts. However, y_R2, y_R3, and y_R4 are alignment irregularity amounts at the positions of the wheel sets 13b, 13c, and 13d respectively.

[0037] Here, the alignment irregularity is a lateral displacement of a rail in a longitudinal direction as described in Japan Industrial Standard (JIS E 1001: 2001). The alignment irregularity amount is an amount of the displacement. Fig. 3A and Fig. 3B each illustrate one example of the alignment irregularity amount y_R1 at the position of the wheel set 13a. In Fig. 3A, the case of the track 16 being the linear track will be explained as an example. In Fig. 3B, the case of the track 16 being the curved track will be explained as an example. In Fig. 3A and Fig. 3B, 16a denotes a rail and 16b denotes a crosstie. In Fig. 3A, it is set that the wheel 14a of the wheel set 13a is in contact with the rail 16a at a position 301. In Fig. 3B, it is set that the wheel 14a of the wheel set 13a is in contact with the rail 16a at a position 302. The alignment irregularity amount y_R1 at the position of the wheel set 13a is a distance in the right and left direction between the contact position between the wheel 14a of the wheel set 13a and the rail 16a and the position of the rail 16a in the case where this position is assumed as a regular state. The position of the wheel set 13a is the contact position between the wheel 14a of the wheel set 13a and the rail 16a. The alignment irregularity amounts y_R2, y_R3, and y_R4 at the positions of the wheel sets 13b, 13c, and 13d are also defined in the same manner as the alignment irregularity amount y_R1 at the position of the wheel set 13a.

[Transversal vibration of the bogie]

[0038] The motion equations that describe the transversal vibrations of the bogies 12a, 12b (motion in the right and left direction) are expressed by (9) Equation and (10) Equation below.
[Mathematical equation 3]

[0039] m_T is the mass of the bogies 12a, 12b. y_t1 · · is acceleration of the bogie 12a in the right and left direction. c'₂ is a damping constant of a lateral movement damper. h₄ is a distance between the center of gravity of the bogie 12a and the lateral movement damper in the up and down direction. y_b · is a velocity of the vehicle body 11 in the right and left direction. L represents 1/2 of the interval between the center of the bogie 12a and the center of the bogie 12b in the forward and backward direction (the interval between the center of the bogie 12a and the center of the bogie 12b in the forward and backward direction becomes 2L). φ_b · is an angular velocity of the vehicle body 11 in the yawing direction. h₅ is a distance between the lateral movement damper and the center of gravity of the vehicle body 11 in the up and down direction. φ_b · is an angular velocity of the vehicle body 11 in the rolling direction. y_w2 · is a velocity of the wheel set 13b in the right and left direction. k'₂ is a spring constant of the air spring (bolster spring) in the right and left direction. h₂ is a distance between the center of gravity of each of the bogies 12a, 12b and the center of the air spring (bolster spring) in the up and down direction. y_b is a displacement of the vehicle body 11 in the right and left direction. φ_b is a pivot amount (angular displacement) of the vehicle body 11 in the yawing direction. h₃ is a distance between the center of the air spring (bolster spring) and the center of gravity of the vehicle body 11 in the up and down direction. φ_b is a pivot amount (angular displacement) of the vehicle body 11 in the rolling direction. Incidentally, respective variables in (10) Equation are represented by being replaced with the variables in (9) Equation according to the meanings of the aforementioned subscripts.

[Yawing of the bogie]

[0040] The motion equations that describe the yawings of the bogies 12a, 12b are expressed by (11) Equation and (12) Equation below.
[Mathematical equation 4]

[0041] I_Tz is a moment of inertia of the bogies 12a, 12b in the yawing direction. φ_t1 · · is angular acceleration of the bogie 12a in the yawing direction. φ_w2 · is an angular velocity of the wheel set 13b in the yawing direction. φ_w2 is a pivot amount (angular displacement) of the wheel set 13b in the yawing direction. y_w2 is a displacement of the wheel set 13b in the right and left direction. k'₀ is stiffness of a rubber bush of the yaw damper. b'₀ represents 1/2 of the interval between the two yaw dampers, which are disposed on the right and left sides of each of the bogies 12a, 12b, in the right and left direction (the interval between the two yaw dampers, which are disposed on the right and left sides of each of the bogies 12a, 12b, in the right and left direction becomes 2b'₀). φ_y1 is a pivot amount (angular displacement) of the yaw damper disposed on the bogie 12a in the yawing direction. k"₂ is a spring constant of the air spring (bolster spring) in the right and left direction. b₂ represents 1/2 of the interval between the two air springs (bolster springs), which are disposed on the right and left sides of each of the bogies 12a, 12b, in the right and left direction (the interval between the two air springs (bolster springs), which are disposed on the right and left sides of each of the bogies 12a, 12b, in the right and left direction becomes 2b₂). Incidentally, respective variables in (12) Equation are represented by being replaced with the variables in (11) Equation according to the meanings of the aforementioned subscripts.

[Rolling of the bogie]

[0042] The motion equations that describe the rollings of the bogies 12a, 12b are expressed by (13) Equation and (14) Equation below.
[Mathematical equation 5]

[0043] I_Tx is a moment of inertia of the bogies 12a, 12b in the rolling direction. φ_t1 · · is angular acceleration of the bogie 12a in the rolling direction. c₁ is a damping constant of an axle damper in the up and down direction. b'₁ represents 1/2 of the interval between the two axle dampers, which are disposed on the right and left sides of each of the bogies 12a, 12b, in the right and left direction (the interval between the two axle dampers, which are disposed on the right and left sides of each of the bogies 12a, 12b, in the right and left direction becomes 2b'₁). c₂ is a damping constant of the air spring (bolster spring) in the up and down direction. φ_a1 · is an angular velocity of the air spring (bolster spring) disposed on the bogie 12a in the rolling direction. k₁ is a spring constant of an axle spring in the up and down direction. λ is a value obtained by dividing the volume of the air spring (bolster spring) main body by the volume of an auxiliary air chamber. k₂ is a spring constant of the air spring (bolster spring) in the up and down direction. φ_a1 is a pivot amount (angular displacement) of the air spring (bolster spring) disposed on the bogie 12a in the rolling direction. k₃ is equivalent stiffness by a change in effective pressure receiving area of the air spring (bolster spring). Incidentally, respective variables in (14) Equation are represented by being replaced with the variables in (13) Equation according to the meanings of the aforementioned subscripts. However, φ_a2 is a pivot amount (angular displacement) of the air spring (bolster spring) disposed on the bogie 12b in the rolling direction.

[Transversal vibration of the vehicle body]

[0044] The motion equation that describes the transversal vibration of the vehicle body 11 (motion in the right and left direction) is expressed by (15) Equation below.
[Mathematical equation 6]

[0045] m_B is the mass of the bogies 12a, 12b. y_b · · is acceleration of the vehicle body 11 in the right and left direction. y_t2 · is a velocity of the bogie 12b in the right and left direction. φ_t2 · is an angular velocity of the bogie 12b in the rolling direction. y_t2 is a displacement of the bogie 12b in the right and left direction. φ_t2 is a pivot amount (angular displacement) of the bogie 12b in the rolling direction.

[Yawing of the vehicle body]

[0046] The motion equation that describes the yawing of the vehicle body 11 is expressed by (16) Equation below.
[Mathematical equation 7]

[0047] I_Bz is a moment of inertia of the vehicle body 11 in the yawing direction. φ_b · · is angular acceleration of the vehicle body 11 in the yawing direction. c₀ is a damping constant of the yaw damper in the forward and backward direction. φ_y1 · is an angular velocity of the yaw damper disposed on the bogie 12a in the yawing direction. φ_y2 · is an angular velocity of the yaw damper disposed on the bogie 12b in the yawing direction. φ_t2 is a pivot amount (angular displacement) of the bogie 12b in the yawing direction.

[Rolling of the vehicle body]

[0048] The motion equation that describes the rolling of the vehicle body 11 is expressed by (17) Equation below.
[Mathematical equation 8]

[0049] I_Bx is a moment of inertia of the vehicle body 11 in the rolling direction. φ_b · · is angular acceleration of the vehicle body 11 in the rolling direction.

[Yawing of the damper]

[0050] The motion equations that describe the yawing of the yaw damper disposed on the bogie 12a and the yawing of the yaw damper disposed on the bogie 12b are expressed by (18) Equation and (19) Equation below respectively.
[Mathematical equation 9]

[0051] φ_y2 is a pivot amount (angular displacement) of the yaw damper disposed on the bogie 12b in the yawing direction.

[Rolling of the air spring (bolster spring)]

[0052] The motion equations that describe the rolling of the air spring (bolster spring) disposed on the bogie 12a and the rolling of the air spring (bolster spring) disposed on the bogie 12b are expressed by (20) Equation and (21) Equation below respectively.
[Mathematical equation 10]

[0053] φ_a2 · is an angular velocity of the air spring (bolster spring) disposed on the bogie 12b in the rolling direction.

(Forward-and-backward-direction force)

[0054] Next, the forward-and-backward-direction force will be explained. Incidentally, the forward-and-backward-direction force itself is the same as that described in Patent Literature 1.

[0055] In-phase components of the longitudinal creep force in one wheel of right and left wheels in one wheel set and the longitudinal creep force in the other wheel are components corresponding to a braking force and a driving force. Accordingly, the forward-and-backward-direction force is preferably determined so as to correspond to an opposite-phase component of the longitudinal creep force. The opposite-phase component of the longitudinal creep force is a component to be opposite in phase to each other between the longitudinal creep force in one wheel of the right and left wheels in one wheel set and the longitudinal creep force in the other wheel. That is, the opposite-phase component of the longitudinal creep force is a component, of the longitudinal creep force, in the direction in which the axle is twisted. In this case, the forward-and-backward-direction force becomes a component opposite in phase to each other out of forward-and-backward-direction components of forces to occur in the aforementioned two members attached to both the right and left sides of one wheel set.

[0056] Hereinafter, there will be explained concrete examples of the forward-and-backward-direction force in the case where the forward-and-backward-direction force is determined so as to correspond to the opposite-phase component of the longitudinal creep force.

[0057] In the case of the axle box suspension being a mono-link type axle box suspension, the axle box suspension includes a link, and the axle box and the bogie frame are coupled by the link. A rubber bush is attached to both ends of the link. In this case, the forward-and-backward-direction force becomes, out of forward-and-backward-direction components of loads that two links, which are attached to right and left ends of one wheel set one by one, receive, the component to be opposite in phase to each other. Further, due to arrangement and constitution of the links, the link mainly receives, out of loads in the forward and backward direction, the right and left direction, and the up and down direction, the load in the forward and backward direction. Accordingly, one strain gauge only needs to be attached to each link, for example. By using a measured value of the strain gauge, the forward-and-backward-direction component of the load that this link receives is derived, to thereby obtain a measured value of the forward-and-backward-direction force. Further, in place of applying such a design, a forward-and-backward-direction displacement of the rubber bush attached to the link may be measured by a displacement meter. In this case, the product of a measured displacement and a spring constant of this rubber bush is set as the measured value of the forward-and-backward-direction force. In the case of the axle box suspension being the mono-link type axle box suspension, the previously-described member for supporting the axle box becomes the link or the rubber bush.

[0058] Incidentally, in the load measured by the strain gauge attached to the link, not only the component in the forward and backward direction, but also at least one component of a component in the right and left direction and a component in the up and down direction is sometimes contained. However, even in such a case, due to the structure of the axle box suspension, the load of the component in the right and left direction and the load of the component in the up and down direction that the link receives are sufficiently smaller than the load of the component in the forward and backward direction. Accordingly, only attaching one strain gauge to each link makes it possible to obtain a measured value of the forward-and-backward-direction force, which has accuracy to be required practically. In this manner, the components other than the component in the forward and backward direction are sometimes included in the measured value of the forward-and-backward-direction force. Thus, three or more strain gauges may be attached to each link so as to cancel the strains in the up and down direction and the right and left direction. This makes it possible to improve the accuracy of the measured value of the forward-and-backward-direction force.

[0059] In the case of the axle box suspension being an axle beam type axle box suspension, the axle box suspension includes an axle beam, and the axle box and the bogie frame are coupled by the axle beam. The axle beam may be formed integrally with the axle box. A rubber bush is attached to a bogie frame-side end of the axle beam. In this case, the forward-and-backward-direction force becomes, out of forward-and-backward-direction components of loads that two axle beams, which are attached to right and left ends of one wheel set one by one, receive, the component to be opposite in phase to each other. Further, due to arrangement and constitution of the axle beams, the axle beam is likely to receive, out of loads in the forward and backward direction, the right and left direction, and the up and down direction, the load in the right and left direction, in addition to the load in the forward and backward direction. Accordingly, two or more strain gauges are attached to each axle beam so as to cancel the strain in the right and left direction, for example. By using measured values of these strain gauges, the forward-and-backward-direction component of the load that the axle beam receives is derived, to thereby obtain a measured value of the forward-and-backward-direction force. Further, in place of applying such a design, a forward-and-backward-direction displacement of the rubber bush attached to the axle beam may be measured by a displacement meter. In this case, the product of a measured displacement and a spring constant of this rubber bush is set as the measured value of the forward-and-backward-direction force. In the case of the axle box suspension being the axle beam type axle box suspension, the previously-described member for supporting the axle box becomes the axle beam or the rubber bush.

[0060] Incidentally, in the load measured by the strain gauge attached to the axle beam, not only the components in the forward and backward direction and the right and left direction, but also the component in the up and down direction is sometimes included. However, even in such a case, due to the structure of the axle box suspension, the load of the component in the up and down direction that the axle beam receives is sufficiently smaller than the load of the component in the forward and backward direction and the load of the component in the right and left direction. Accordingly, unless the strain gauge is attached so as to cancel the load of the component in the up and down direction that the axle beam receives, a measured value of the forward-and-backward-direction force, which has accuracy to be required practically, can be obtained. In this manner, the components other than the component in the forward and backward direction are sometimes included in the measured forward-and-backward-direction force, and three or more strain gauges may be attached to each axle beam so as to cancel the strain in the up and down direction as well as the strain in the right and left direction. This makes it possible to improve the accuracy of the measured value of the forward-and-backward-direction force.

[0061] In the case of the axle box suspension being a leaf spring type axle box suspension, the axle box suspension includes a leaf spring, and the axle box and the bogie frame are coupled by the leaf spring. A rubber bush is attached to ends of the leaf spring. In this case, the forward-and-backward-direction force becomes, out of forward-and-backward-direction components of loads that two leaf springs, which are attached to right and left ends of one wheel set one by one, receive, the component to be opposite in phase to each other. Further, due to arrangement and constitution of the leaf springs, the leaf spring is likely to receive, out of loads in the forward and backward direction, the right and left direction, and the up and down direction, the load in the right and left direction and the load in the up and down direction, in addition to the load in the forward and backward direction. Accordingly, three or more strain gauges are attached to each leaf spring so as to cancel the strains in the right and left direction and the up and down direction, for example. By using measured values of these strain gauges, the forward-and-backward-direction component of the load that the leaf spring receives is derived, to thereby obtain a measured value of the forward-and-backward-direction force. Further, in place of applying such a design, a forward-and-backward-direction displacement of the rubber bush attached to the leaf spring may be measured by a displacement meter. In this case, the product of a measured displacement and a spring constant of this rubber bush is set as the measured value of the forward-and-backward-direction force. In the case of the axle box suspension being the leaf spring type axle box suspension, the previously-described member for supporting the axle box becomes the leaf spring or the rubber bush.

[0062] Incidentally, as the previously-described displacement meter, a well-known laser displacement meter or eddy current displacement meter can be used.

[0063] Further, the forward-and-backward-direction force has been explained here by taking the case of the system of the axle box suspension being a mono-link type, an axle beam type, and a leaf spring type as an example. However, the system of the axle box suspension is not limited to the mono-link type, the axle beam type, and the leaf spring type. In conformity with the system of the axle box suspension, the forward-and-backward-direction force can be determined in the same manner as in the mono-link type, the axle beam type, and the leaf spring type.

[0064] Further, the case where a measured value of a single forward-and-backward-direction force can be obtained in one wheel set will be explained as an example, in order to simplify the explanation below. That is, the railway vehicle illustrated in Fig. 1 has the four wheel sets 13a to 13d. Accordingly, it is possible to obtain measured values of four forward-and-backward-direction forces T₁ to T₄.

(First embodiment)

[0065] Next, there will be explained a first embodiment of the present invention.

inspection apparatus 400>

[0066] Fig. 4 is a diagram illustrating one example of a functional configuration of an inspection apparatus 400. Fig. 5 is a diagram illustrating one example of a hardware configuration of the inspection apparatus 400. Fig. 6 is a flowchart illustrating one example of first preprocessing in the inspection apparatus 400. The first preprocessing is processing for setting state equations and observation equations used in second preprocessing and main processing. Fig. 7 is a flowchart illustrating one example of the second preprocessing in the inspection apparatus 400. The second preprocessing is processing to derive a correction amount for an estimated value of the previously-described alignment irregularity amount after the first preprocessing is finished. Fig. 8 is a flowchart illustrating one example of the main processing in the inspection apparatus 400. The main processing is processing to derive an estimated value of a final alignment irregularity amount after the first preprocessing and the second preprocessing are finished. In this embodiment, as illustrated in Fig. 1, the case where the inspection apparatus 400 is mounted on the railway vehicle will be explained as an example. In the following explanation, the railway vehicle is set the same as the railway vehicle with the inspection apparatus 400 mounted thereon.

[0067] In Fig. 4, the inspection apparatus 400 includes, as its functions, a state equation storage unit 401, an observation equation storage unit 402, a data acquisition unit 403, a first frequency adjustment unit 404, a filter operation unit 405, a second frequency adjustment unit 406, a first track state calculation unit 407, an actual value acquisition unit 408, a correction amount calculation unit 409, a correction amount storage unit 410, a second track state calculation unit 411, a track state correction unit 412, and an output unit 413.

[0068] In Fig. 5, the inspection apparatus 400 includes a CPU 501, a main memory 502, an auxiliary memory 503, a communication circuit 504, a signal processing circuit 505, an image processing circuit 506, an I/F circuit 507, a user interface 508, a display 509, and a bus 510.

[0069] The CPU 501 overall controls the entire inspection apparatus 400. The CPU 501 uses the main memory 502 as a work area to execute a program stored in the auxiliary memory 503. The main memory 502 stores data temporarily. The auxiliary memory 503 stores various data, in addition to programs to be executed by the CPU 501. The auxiliary memory 503 stores state equations, observation equations, and correction amounts (a first correction amount, a second correction amount) to be described later. The state equation storage unit 401, the observation equation storage unit 402, and the correction amount storage unit 410 are fabricated by using the CPU 501 and the auxiliary memory 503, for example.

[0070] The communication circuit 504 is a circuit intended for performing communication with the outside of the inspection apparatus 400. The communication circuit 504 receives information of the measured value of the forward-and-backward-direction force and pieces of information of measured values of accelerations of the vehicle body 11, the bogies 12a, 12b, and the wheel sets 13a to 13d in the right and left direction, for example. The communication circuit 504 may perform radio communication or wire communication with the outside of the inspection apparatus 400. The communication circuit 504 is connected to an antenna provided on the railway vehicle in the case of performing radio communication.

[0071] The signal processing circuit 505 performs various pieces of signal processing on signals received in the communication circuit 504 and signals input according to the control by the CPU 501. The data acquisition unit 403 and the actual value acquisition unit 408 are fabricated by using the CPU 501, the communication circuit 504, and the signal processing circuit 505, for example. Further, the first frequency adjustment unit 404, the filter operation unit 405, the second frequency adjustment unit 406, the first track state calculation unit 407, the correction amount calculation unit 409, the second track state calculation unit 411, and the track state correction unit 412 are fabricated by using the CPU 501 and the signal processing circuit 505, for example.

[0072] The image processing circuit 506 performs various pieces of image processing on signals input according to the control by the CPU 501. The signal that has been subjected to the image processing is output on the display 509.

[0073] The user interface 508 is a part through which an operator gives an instruction to the inspection apparatus 400. The user interface 508 includes buttons, switches, dials, and so on, for example. Further, the user interface 508 may include a graphical user interface using the display 509.

[0074] The display 509 displays an image based on a signal output from the image processing circuit 506. The I/F circuit 507 exchanges data with a device connected to the I/F circuit 507. In Fig. 5, as the device to be connected to the I/F circuit 507, the user interface 508 and the display 509 are illustrated. However, the device to be connected to the I/F circuit 507 is not limited to these. For example, a portable storage medium may be connected to the I/F circuit 507. Further, at least a part of the user interface 508 and the display 509 may be provided outside the inspection apparatus 400.

[0075] The output unit 413 is fabricated by using the communication circuit 504, the signal processing circuit 505, and at least any one of the image processing circuit 506, the I/F circuit 507, and the display 509, for example.

[0076] Incidentally, the CPU 501, the main memory 502, the auxiliary memory 503, the signal processing circuit 505, the image processing circuit 506, and the I/F circuit 507 are connected to the bus 510. Communication among these components is performed via the bus 510. Further, the hardware of the inspection apparatus 400 is not limited to the one illustrated in Fig. 5 as long as it can perform later-described functions of the inspection apparatus 400.

[State equation storage unit 401, S601]

[0077] The state equation storage unit 401 stores state equations. In this embodiment, the case of using the state equations described in Patent Literature 1 will be explained as an example. As described previously, in this embodiment, the motion equations that describe the yawings of the wheel sets 13a to 13d of (5) Equation to (8) Equation are not included in the state equation, and the state equation is constituted as follows.

[0078] First, the motion equations that describe the transversal vibrations of the bogies 12a, 12b (motion in the right and left direction) of (9) Equation and (10) Equation, the motion equations that describe the rollings of the bogies 12a, 12b of (13) Equation and (14) Equation, the motion equation that describes the transversal vibration of the vehicle body 11 (motion in the right and left direction) of (15) Equation, the motion equation that describes the yawing of the vehicle body 11 of (16) Equation, the motion equation that describes the rolling of the vehicle body 11 of (17) Equation, the motion equations that describe the yawings of the yaw damper disposed on the bogie 12a and the yaw damper disposed on the bogie 12b of (18) Equation and (19) Equation, and the motion equations that describe the rollings of the air spring (bolster spring) disposed on the bogie 12a and the air spring (bolster spring) disposed on the bogie 12b of (20) Equation and (21) Equation are used as they are to constitute the state equation.

[0079] In the meantime, in the motion equations that describe the transversal vibrations of the wheel sets 13a to 13d (motion in the right and left direction) of (1) Equation to (4) Equation and the motion equations that describe the yawings of the bogies 12a, 12b of (11) Equation and (12) Equation, the pivot amounts (angular displacements) φ_w1 to φ_w4 and the angular velocities φ_w1 · to φ_w4 · of the wheel sets 13a to 13d in the yawing direction are included. Results obtained after eliminating these variables from (1) Equation to (4) Equation, (11) Equation, and (12) Equation are used to constitute the state equation.

[0080] First, the forward-and-backward-direction forces T₁ to T₄ of the wheel sets 13a to 13d are expressed by (22) Equation to (25) Equation below. In this manner, the forward-and-backward-direction forces T₁ to T₄ are determined according to the differences between the angular displacements φ_w1 to φ_w4 of the wheel sets in the yawing direction and the angular displacements φ_t1 to φ_t2 of the bogies on which these wheel sets are provided in the yawing direction.
[Mathematical equation 11]

[0081] Transformation variables e₁ to e₄ are defined as in (26) Equation to (29) Equation below. As above, the transformation variables e₁ to e₄ are defined by the differences between the angular displacements φ_t1 to φ_t2 of the bogies in the yawing direction and the angular displacements φ_w1 to φ_w4 of the wheel sets in the yawing direction. The transformation variables e₁ to e₄ are variables for performing mutual transformation between the angular displacements φ_t1 to φ_t2 of the bogies in the yawing direction and the angular displacements φ_w1 to φ_w4 of the wheel sets in the yawing direction.
[Mathematical equation 12]

[0082] When (26) Equation to (29) Equation are modified, (30) Equation to (33) Equation below are obtained.
[Mathematical equation 13]

[0083] When (30) Equation to (33) Equation are substituted into the motion equations that describe the transversal vibrations of the wheel sets 13a to 13d (motion in the right and left direction) of (1) Equation to (4) Equation, (34) Equation to (37) Equation below are obtained.
[Mathematical equation 14]

[0084] As above, the motion equations that describe the transversal vibrations of the wheel sets 13a to 13d (motion in the right and left direction) of (1) Equation to (4) Equation are expressed by using the transformation variables e₁ to e₄, thereby making it possible to eliminate the pivot amounts (angular displacements) φ_w1 to φ_w4 of the wheel sets 13a to 13d in the yawing direction that are included in these motion equations.

[0085] When (22) Equation to (25) Equation are substituted into the motion equations that describe the yawings of the bogies 12a, 12b of (11) Equation and (12) Equation, (38) Equation and (39) Equation below are obtained.
[Mathematical equation 15]

[0086] As above, the motion equations that describe the yawings of the bogies 12a, 12b of (11) Equation and (12) Equation are expressed by using the forward-and-backward-direction forces T₁ to T₄, thereby making it possible to eliminate the angular displacements φ_w1 to φ_w4 and the angular velocities φ_w1 · to φ_w4 · of the wheel sets 13a to 13d in the yawing direction that are included in these motion equations.

[0087] Further, when (26) Equation to (29) Equation are substituted into (22) Equation to (25) Equation, (40) Equation to (43) Equation below are obtained.
[Mathematical equation 16]

[0088] As above, in this embodiment, as in (34) Equation to (37) Equation, the motion equations that describe the transversal vibrations of the wheel sets 13a to 13d (motion in the right and left direction) are expressed, and at the same time, as in (38) Equation and (39) Equation, the motion equations that describe the yawings of the bogies 12a, 12b are expressed, and by using these, the state equation is constituted. Further, (40) Equation to (43) Equation are ordinary differential equations, and actual values of the transformation variables e₁ to e₄, which are solutions of the equations, can be derived by using the values of the forward-and-backward-direction forces T₁ to T₄ in the wheel sets 13a to 13d. Here, the values of the forward-and-backward-direction forces T₁ to T₄ are that a signal strength of a low-frequency component to be generated due to the railway vehicle traveling on the curved portion of the track is reduced from the time-series data of the measured value of the forward-and-backward-direction force by the later-described first frequency adjustment unit 404.

[0089] The actual values of the transformation variables e₄ to e₄ derived as above are given to (34) Equation to (37) Equation. Further, the values of the forward-and-backward-direction forces T₁ to T₄ in the wheel sets 13a to 13d are given to (38) Equation and (39) Equation. Here, the values of the forward-and-backward-direction forces T₁ to T₄ are that a signal strength of a low-frequency component to be generated due to the railway vehicle traveling on the curved portion of the track is reduced from the time-series data of the measured value of the forward-and-backward-direction force by the later-described first frequency adjustment unit 404.

[0090] In this embodiment, variables illustrated in (44) Equation below are set as the state variables, and by using the motion equations of (9) Equation, (10) Equation, (13) Equation to (21) Equation, and (34) Equation to (39) Equation, the state equation is constituted.
[Mathematical equation 17]

[0091] The state equation storage unit 401 receives the state equation constituted as above, for example, based on the operation of the user interface 508 by an operator and stores it.

[Observation equation storage unit 402, S602]

[0092] The observation equation storage unit 402 stores observation equations. In this embodiment, the acceleration of the vehicle body 11 in the right and left direction, the accelerations of the bogies 12a, 12b in the right and left direction, and the accelerations of the wheel sets 13a to 13d in the right and left direction are set to observation variables. These observation variables are observation variables of filtering by a later-described Kalman filter. In this embodiment, the motion equations that describe the transversal vibrations of (34) Equation to (37) Equation, (9) Equation, (10) Equation, and (15) Equation are used to constitute an observation equation. The observation equation storage unit 402 receives the observation equation constituted in this manner, for example, based on the operation of the user interface 508 by an operator and stores it.

[0093] After the state equation and the observation equation are stored in the inspection apparatus 400 as above, the data acquisition unit 403, the first frequency adjustment unit 404, the filter operation unit 405, the second frequency adjustment unit 406, the first track state calculation unit 407, the actual value acquisition unit 408, the correction amount calculation unit 409, and the correction amount storage unit 410 start. That is, after the first preprocessing by the flowchart in Fig. 6 is finished, the second preprocessing by the flowchart in Fig. 7 starts.

[Data acquisition unit 403, S701]

[0094] The data acquisition unit 403 acquires measured data with a predetermined sampling period.

[0095] In this embodiment, the data acquisition unit 403 acquires, as the measured data, time-series data of a measured value of the acceleration of the vehicle body 11 in the right and left direction, time-series data of measured values of the accelerations of the bogies 12a, 12b in the right and left direction, and time-series data of measured values of the accelerations of the wheel sets 13a to 13d in the right and left direction. The respective accelerations are measured by using strain gauges attached to, for example, the vehicle body 11, the bogies 12a, 12b, and the wheel sets 13a to 13d respectively and an arithmetic device that calculates the accelerations by using measured values of these strain gauges. Incidentally, the measurement of the accelerations can be performed by a well-known technique, and thus its detailed explanation is omitted.

[0096] Further, the data acquisition unit 403 acquires, as the measured data, time-series data of the measured value of the forward-and-backward-direction force. The method of measuring the forward-and-backward-direction force is as described previously.

[0097] The data acquisition unit 403 can acquire the measured data by performing communication with the previously-described arithmetic device, for example. At Step S701, the data acquisition unit 403 acquires measured data in the entire traveling section of the railway vehicle.

[First frequency adjustment unit 404, S702]

[0098] The first frequency adjustment unit 404 reduces (preferably removes) the signal strength of the low-frequency component contained in the time-series data of the measured value of the forward-and-backward-direction force (a second physical quantity) out of the measured data acquired by the data acquisition unit 403. A signal of this low-frequency component is a signal that is not measured when the railway vehicle is traveling on the linear track, but is measured when the railway vehicle is traveling on the curved track. That is, the signal measured when the railway vehicle is traveling on the curved track can be regarded as a signal obtained by superimposing the signal of this low-frequency component on the signal measured when the railway vehicle is traveling on the linear track.

[0099] The present inventors devised a model in which an AR (Auto-regressive) model is corrected. Then, the present inventors devised reducing the signal strength of the low-frequency component contained in the time-series data of the measured value of the forward-and-backward-direction force by using this model. In the following explanation, the model devised by the present inventors is referred to as a corrected AR model. In contrast to this, the well-known AR model is referred to as an AR model simply. Hereinafter, there will be explained one example of the corrected AR model.

[0100] A value of time-series data y of a physical quantity at a time k (1 is k is M) is set to y_k. M is a number indicating, as the time-series data y of the physical quantity, data until when are contained, and is preset. In the following explanation, the time-series data of the physical quantity will be abbreviated to data y as necessary. The AR model approximating the value y_k of the data y is as in (45) Equation below, for example. The AR model is, as illustrated in (45) Equation, an equation expressing a predicted value y^_k of the physical quantity at the time k (m + 1 ≦ k ≦ M) in the data y by using an actual value y _{k - l} of the physical quantity at a time k - l (1 ≦ l ≦ m) prior to the time k in the data y. Incidentally, y^_k is expressed by adding ^ above y_k in (45) Equation.
[Mathematical equation 18]

[0101] In (45) Equation, α is a coefficient of the AR model. m is a number of the value of the data y to be used for approximating the value y_k of the data y at the time k in the AR model, and is a number among values y_{k - l} to y_{k - m} of the data y at continuous times k - l to k - m prior to the time k. m is an integer less than M. As m, for example, 1500 can be used.

[0102] Then, there is derived a conditional expression for approximating the predicted value y^_k of the physical quantity at the time k by the AR model to the value y_k by using a least square method. As the condition for approximating the predicted value y^_k of the physical quantity at the time k by the AR model to the value y_k, it is possible to employ a condition that minimizes a square error between the predicted value y^_k of the physical quantity at the time k by the AR model and the value y_k, for example. That is, the least square method is used in order to approximate the predicted value y^_k of the physical quantity at the time k by the AR model to the value y_k. (46) Equation below is a conditional expression for minimizing the square error between the predicted value y^_k of the physical quantity at the time k by the AR model and the value y_k.
[Mathematical equation 19]

[0103] The relation of (47) Equation below is established by (46) Equation.
[Mathematical equation 20]

[0104] Further, (47) Equation is modified (expressed in matrix notation form), and thereby (48) Equation below is obtained.
[Mathematical equation 21]

[0105] R_j1 in (48) Equation is called autocorrelation of the data y, and is a value defined by (49) Equation below. |j - l| at this time is referred to as a time lag.
[Mathematical equation 22]

[0106] Based on (48) Equation, (50) Equation below is considered. (50) Equation is an equation derived from a condition that minimizes the error between the predicted value y^_k of the physical quantity at the time k by the AR model and the value y_k of the physical quantity at the time k corresponding to the predicted value y^_k. (50) Equation is called a Yule-Walker equation. Further, (50) Equation is a linear equation in which a vector composed of coefficients of the AR model is set to a variable vector. A constant vector on the left side in (50) Equation is a vector whose component is the autocorrelation of the data y with a time lag of 1 to m. In the following explanation, the constant vector on the left side in (50) Equation is referred to as an autocorrelation vector as necessary. Further, a coefficient matrix on the right side in (50) Equation is a matrix whose component is the autocorrelation of the data y with a time lag of 0 to m - 1. In the following explanation, the coefficient matrix on the right side in (50) Equation is referred to as an autocorrelation matrix as necessary.
[Mathematical equation 23]

[0107] Further, the autocorrelation matrix on the right side in (50) Equation (a matrix of m × m composed of R_j1) is described as an autocorrelation matrix R as in (51) Equation below.
[Mathematical equation 24]

[0108] In general, when deriving the coefficient of the AR model, a method of solving a coefficient α of (50) Equation is used. In (50) Equation, the coefficient α is derived so as to make the predicted value y^_k of the physical quantity at the time k derived by the AR model come close to the value y_k of the physical quantity at the time k as much as possible. Therefore, frequency characteristics of the AR model include a large number of frequency components contained in the value y_k of the data y at each time.

[0109] Thus, the present inventors focused on the autocorrelation matrix R to be multiplied by the coefficient α of the AR model and earnestly examined it. As a result, the present inventors found out that it is possible to reduce the effect of a high-frequency component contained in the data y by using a part of eigenvalues of the autocorrelation matrix R. That is, the present inventors found out that it is possible to rewrite the autocorrelation matrix R so that the low-frequency component is emphasized.

[0110] There will be explained a concrete example of the above below.

[0111] The autocorrelation matrix R is subjected to singular value decomposition. Elements of the autocorrelation matrix R are symmetric. Thus, when the autocorrelation matrix R is subjected to singular value decomposition, as in (52) Equation below, the result becomes the product of an orthogonal matrix U, a diagonal matrix Σ, and a transposed matrix of the orthogonal matrix U.
[Mathematical equation 25]

[0112] The diagonal matrix Σ in (52) Equation is a matrix whose diagonal component is the eigenvalues of the autocorrelation matrix R as illustrated in (53) Equation below. The diagonal component of the diagonal matrix Σ is set to σ₁₁, σ₂₂, ···, σ_mm. Further, the orthogonal matrix U is a matrix in which each column component vector is an eigenvector of the autocorrelation matrix R. The column component vector of the diagonal matrix U is set to u₁, u₂, ···, u_m. There is a correspondence relation in which the eigenvalue of the autocorrelation matrix R responsive to an eigenvector u_j is σ_jj. The eigenvalue of the autocorrelation matrix R is a variable reflecting the strength of each frequency component included in a time waveform of the predicted value y^_k of the physical quantity at the time k by the AR model.
[Mathematical equation 26]

[0113] The values of σ₁₁, σ₂₂, ···, σ_mm being the diagonal components of the diagonal matrix Σ obtained by the result of the singular value decomposition of the autocorrelation matrix R are set in descending order in order to simplify the illustration of the mathematical equation. A matrix R' is defined as in (54) Equation below by using, out of the eigenvalues of the autocorrelation matrix R illustrated in (53) Equation, s pieces of the eigenvalues, which are chosen from the largest. s is a number that is 1 or more and less than m. In this embodiment, s is preset. The matrix R' is a matrix resulting from approximating the autocorrelation matrix R by using s pieces of the eigenvalues out of the eigenvalues of the autocorrelation matrix R.
[Mathematical equation 27]

[0114] A matrix U_s in (54) Equation is a matrix of m × s composed of s pieces of the column component vectors (eigenvectors corresponding to the eigenvalues to be used), which are chosen from the left of the orthogonal matrix U of (52) Equation. That is, the matrix U_s is a submatrix composed of the left elements of m × s cut out from the orthogonal matrix U. Further, U_s^T in (54) Equation is a transposed matrix of U_s. U_s^T is a matrix of s × m composed of s pieces of row component vectors, which are chosen from the top of the matrix U^T in (52) Equation. The matrix Σ_s in (54) Equation is a matrix of s × s composed of s pieces of columns, which are chosen from the left, and s pieces of rows, which are chosen from the top, of the diagonal matrix Σ in (52) Equation. That is, the matrix Σ_s is a submatrix composed of the top and left elements of s × s cut out from the diagonal matrix Σ.

[0115] When the matrix Σ_s and the matrix U_s are expressed by the matrix elements, (55) Equation below is obtained.
[Mathematical equation 28]

[0116] By using the matrix R' in place of the autocorrelation matrix R, the relational expression of (50) Equation is rewritten into (56) Equation below.
[Mathematical equation 29]

[0117] (56) Equation is modified, and thereby (57) Equation below is obtained as the equation deriving the coefficient α . The model that calculates the predicted value y^_k of the physical quantity at the time k from (45) Equation while using the coefficient α derived by (57) Equation is the "corrected AR model."
[Mathematical equation 30]

[0118] The case where the values of σ₁₁, σ₂₂, ···, σ_mm being the diagonal components of the diagonal matrix Σ are set in descending order has been explained here as an example. However, it is not necessary to set the diagonal components of the diagonal matrix Σ in descending order during a process of calculating the coefficient α. In this case, the matrix U_s is not the submatrix composed of the left elements of m × s cut out from the orthogonal matrix U, but becomes a submatrix composed of the cut out column component vectors corresponding to the eigenvalues to be used (the eigenvectors). Further, the matrix Σ_s is not the submatrix composed of the top and left elements of s × s cut out from the diagonal matrix Σ, but becomes a submatrix to be cut out so as to make the eigenvalues used for determining the coefficient of the corrected AR model become the diagonal components.

[0119] (57) Equation is an equation to be used for determining the coefficient of the corrected AR model. The matrix U_s in (57) Equation is a matrix (a third matrix) in which the eigenvectors corresponding to the eigenvalues used for determining the coefficient of the corrected AR model are set to the column component vectors, which is the submatrix of the orthogonal matrix U obtained by the singular value decomposition of the autocorrelation matrix R. Further, the matrix Σ_s in (57) Equation is a matrix (a second matrix) in which the eigenvalues used for determining the coefficient of the corrected AR model are set to the diagonal components, which is the submatrix of the diagonal matrix obtained by the singular value decomposition of the autocorrelation matrix R. The matrix U_sΣ_sU_s^T in (57) Equation is a matrix (a first matrix) derived from the matrix Σ_s and the matrix U_s.

[0120] The right side of (57) Equation is calculated, and thereby the coefficient α of the corrected AR model is derived. One example of the method of deriving the coefficient α of the corrected AR model has been explained above. Here, as the method of deriving the coefficient of the AR model to be the base of the corrected AR model, the method of using the least square method for the predicted value y^_k of the physical quantity at the time k has been set in order to make the method understandable intuitively. However, there has been known a method of defining the AR model by using the concept of a stochastic process and deriving its coefficient generally. In this case, the autocorrelation is expressed by autocorrelation of the stochastic process (a population). This autocorrelation of the stochastic process is expressed as a function of a time lag. Thus, the autocorrelation of the data y in this embodiment may be replaced with a value calculated by another calculating formula as long as it approximates the autocorrelation of the stochastic process. For example, R₂₂ to R_mm are autocorrelation with a time lag of 0 (zero), but they may be replaced with R₁₁.

[0121] The number s of the eigenvalues extracted from the autocorrelation matrix R illustrated in (53) Equation can be determined from a distribution of the eigenvalues of the autocorrelation matrix R, for example.

[0122] As the physical quantity in the explanation of the previously-described corrected AR model, the forward-and-backward-direction force is applied here. The value of the forward-and-backward-direction force varies according to the state of the railway vehicle. Thus, the railway vehicle is first made to travel on the track 16 to obtain the data y of the measured value of the forward-and-backward-direction force. The autocorrelation matrix R is derived by using (49) Equation and (51) Equation for each of the obtained data y. The autocorrelation matrix R is subjected to singular value decomposition expressed by (52) Equation, to thereby derive the eigenvalues of the autocorrelation matrix R. Fig. 9 is a view illustrating one example of the distribution of the eigenvalues of the autocorrelation matrix R. In Fig. 9, eigenvalues σ₁₁ to σ_mm, which are obtained by the autocorrelation matrix R in each of the data y of the measured value of the forward-and-backward-direction force T₁ in the wheel set 13a being subjected to singular value decomposition, are aligned in ascending order and are plotted. In Fig. 9, the horizontal axis is an index of the eigenvalue and the vertical axis is the value of the eigenvalue.

[0123] The example illustrated in Fig. 9 includes one eigenvalue having a value significantly higher than the others. Further, the example includes two eigenvalues that have a value relatively larger as compared to the others and are not regarded as 0 (zero), which are not as high as the aforementioned eigenvalue having a significantly higher value. This makes it possible to employ, for example, two or three as the number s of the eigenvalues to be extracted from the autocorrelation matrix R illustrated in (53) Equation. There is no significant difference in the results whichever number is employed.

[0124] The first frequency adjustment unit 404 performs the following processing by using the value y_k of the data y of the measured value of the forward-and-backward-direction force at the time k acquired in the data acquisition unit 403.

[0125] First, the first frequency adjustment unit 404 generates the autocorrelation matrix R using (49) Equation and (51) Equation based on the data y of the measured value of the forward-and-backward-direction force and preset numbers M, m.

[0126] Then, the first frequency adjustment unit 404 performs singular value decomposition on the autocorrelation matrix R, to thereby derive the orthogonal matrix U and the diagonal matrix Σ of (52) Equation, and derives the eigenvalues σ₁₁ to σ_mm of the autocorrelation matrix R from the diagonal matrix Σ.

[0127] Then, the first frequency adjustment unit 404 chooses s pieces of the eigenvalues σ₁₁ to σ_ss from the largest from among the plural eigenvalues σ₁₁ to σ_mm of the autocorrelation matrix R as the eigenvalues of the autocorrelation matrix R to be used for deriving the coefficient α of the corrected AR model.

[0128] Then, the first frequency adjustment unit 404 determines the coefficient α of the corrected AR model using (57) Equation based on the data y of the measured value of the forward-and-backward-direction force, the eigenvalues σ₁₁ to σ_ss, and the orthogonal matrix U obtained by the singular value decomposition of the autocorrelation matrix R.

[0129] Then, the first frequency adjustment unit 404 derives the predicted value y^_k of the data y of the measured value of the forward-and-backward-direction force at the time k from (45) Equation based on the coefficient α of the corrected AR model and the data y of the measured value of the forward-and-backward-direction force. Time-series data of the predicted value y^_k of the forward-and-backward-direction force result in the time-series data from which the low-frequency component contained in the data y of the measured value of the forward-and-backward-direction force has been extracted.

[0130] Fig. 10 is a view illustrating one example of the time-series data of the measured value of the forward-and-backward-direction force (the measured value) and the time-series data of the predicted value of the forward-and-backward-direction force (the calculated value). Incidentally, in this embodiment, the measured values of the four forward-and-backward-direction forces T₁ to T₄ are obtained. That is, four pieces of the data y of the forward-and-backward-direction force are obtained. In Fig. 10, the measured value and the calculated value of each of four pieces of the data y are illustrated. The horizontal axis in Fig. 10 indicates a measuring time and a calculating time of the forward-and-backward-direction forces T₁ to T₄, each of which is an elapsed time (second) from a reference time when the reference time is set to 0 (zero). The vertical axis indicates the forward-and-backward-direction forces T₁ to T₄ (Nm).

[0131] In Fig. 10, the calculated value of the forward-and-backward-direction force T₁ in the wheel set 13a is biased at about 15 seconds to 35 seconds (namely, a value larger than that at another time is exhibited). This period corresponds to the period when the wheel set 13a passes through the curved track. The calculated value of the forward-and-backward-direction force T₂ in the wheel set 13b, the calculated value of the forward-and-backward-direction force T₃ in the wheel set 13c, and the calculated value of the forward-and-backward-direction force T₄ in the wheel set 13d are also biased during the period when the wheel sets 13b, 13c, and 13d pass through the curved track similarly to the calculated value of the forward-and-backward-direction force T₁ in the wheel set 13a.

[0132] Accordingly, in Fig. 10, removal of the calculated values from the measured values of the forward-and-backward-direction forces T₁ to T₄ in the wheel sets 13a to 13d makes it possible to remove the low-frequency components, which are due to the wheel sets 13a to 13d passing through the curved track, from the signals of the forward-and-backward-direction forces T₁ to T₄. That is, in Fig. 10, when the calculated values are removed from the measured values of the forward-and-backward-direction forces T₁ to T₄ in the wheel sets 13a to 13d, as the forward-and-backward-direction forces T₁ to T₄ when the wheel sets 13a to 13d have passed through the curved track, the forward-and-backward-direction forces equivalent to those when the wheel sets 13a to 13d have passed through the linear track can be obtained.

[0133] Thus, the first frequency adjustment unit 404 subtracts the time-series data of the predicted value y^_k of the forward-and-backward-direction force from the time-series data (the data y) of the measured value y_k of the forward-and-backward-direction force. In the following explanation, the time-series data resulting from the subtraction of the time-series data of the predicted value y^_k of the forward-and-backward-direction force from the time-series data (the data y) of the measured value y_k of the forward-and-backward-direction force are referred to as time-series data of a high-frequency component of the forward-and-backward-direction force as necessary. Further, a value of the time-series data of the high-frequency component of the forward-and-backward-direction force at each sampling time is referred to as a value of the high-frequency component of the forward-and-backward-direction force as necessary.

[0134] Fig. 11 is a view illustrating one example of the time-series data of the high-frequency component of the forward-and-backward-direction force. The vertical axis in Fig. 11 indicates the time-series data of the high-frequency components of the forward-and-backward-direction forces T₁, T₂, T₃, and T₄. That is, the high-frequency components of the forward-and-backward-direction forces T₁, T₂, T₃, and T₄ illustrated on the vertical axis in Fig. 11 are ones obtained by subtracting the calculated values from the measured values of the forward-and-backward-direction forces T₁, T₂, T₃, and T₄ in the wheel sets 13a, 13b, 13c, and 13d that are illustrated in Fig. 10 respectively. Further, the horizontal axis in Fig. 11 indicates a measuring time and a calculating time of the forward-and-backward-direction forces T₁ to T₄, each of which is an elapsed time (second) from a reference time when the reference time is set to 0 (zero), similarly to the horizontal axis in Fig. 10.

[0135] The first frequency adjustment unit 404 derives the time-series data of the high-frequency components of the forward-and-backward-direction forces T₁ to T₄ as above.

[Filter operation unit 405, S703]

[0136] The filter operation unit 405 sets the observation equation as the observation equation stored by the observation equation storage unit 402, sets the state equation as the state equation stored by the state equation storage unit 401, and determines estimated values of the state variables illustrated in (44) Equation by the Kalman filter. At this time, the filter operation unit 405 uses, out of pieces of the measured data acquired in the data acquisition unit 403, the measured data excluding the forward-and-backward-direction forces T₁ to T₄ and the time-series data of the high-frequency components of the forward-and-backward-direction forces T₁ to T₄ generated in the first frequency adjustment unit 404. As described previously, in this embodiment, in the measured data, the measured value of the acceleration of the vehicle body 11 in the right and left direction, the measured values of the accelerations of the bogies 12a, 12b in the right and left direction, and the measured values of the accelerations of the wheel sets 13a to 13d in the right and left direction are contained. As for the forward-and-backward-direction forces T₁ to T₄ in the wheel sets 13a to 13d, the time-series data of the high-frequency components of the forward-and-backward-direction forces T₁ to T₄ generated in the first frequency adjustment unit 404 are used without using the measured data (the measured values) acquired in the data acquisition unit 403.

[0137] The Kalman filter is one of the methods of performing data assimilation. That is, the Kalman filter is one example of the method to determine an estimated value of an unobserved variable (state variable) so as to make the difference between, of an observable variable (observation variable), a measured value and an estimated value small (minimum). The filter operation unit 405 derives a Kalman gain at which the difference between, of the observation variable, the measured value and the estimated value becomes small (minimum) and derives the estimated value of the unobserved variable (state variable) at that time. In the Kalman filter, the following observation equation of (58) Equation and the following state equation of (59) Equation are used.

[0138] In (58) Equation, Y is a vector in which the measured value of the observation variable is stored. H is an observation model. X is a vector in which the state variable is stored. V is observation noise. In (59) Equation, X · indicates a time differentiation of X. Φ is a linear model. W is system noise. Incidentally, the Kalman filter itself can be fabricated by a well-known technique, and thus its detailed explanation is omitted.

[0139] The filter operation unit 405 determines the estimated values of the state variables illustrated in (44) Equation with a predetermined sampling period, to thereby generate time-series data of the estimated values of the state variables illustrated in (44) Equation.

[Second frequency adjustment unit 406, S704]

[0140] Unless the signal strength of the low-frequency component contained in the time-series data of the measured value of the forward-and-backward-direction force is removed sufficiently by the first frequency adjustment unit 404, the signal of the low-frequency component due to the railway vehicle traveling on the curved track may be left in the time-series data of the estimated values of the state variables generated by the filter operation unit 405. Thus, the second frequency adjustment unit 406 reduces (preferably removes) the signal strength of the low-frequency component contained in the time-series data of the estimated values of the state variables (the second physical quantity) generated by the filter operation unit 405. Incidentally, in the case where it is possible to determine the number s of the eigenvalues to be extracted from the autocorrelation matrix R illustrated in (53) Equation so as to sufficiently remove the signal strength of the low-frequency component contained in the time-series data of the measured value of the forward-and-backward-direction force by the first frequency adjustment unit 404, the processing of the second frequency adjustment unit 406 is no longer required.

[0141] In this embodiment, the second frequency adjustment unit 406 uses the corrected AR model to reduce the signal strength of the low-frequency component contained in the time-series data of the estimated values of the state variables similarly to the first frequency adjustment unit 404.

[0142] The second frequency adjustment unit 406 performs the following processing for each state variable with a predetermined sampling period.

[0143] As the physical quantity in the explanation of the previously-described corrected AR model, the state variable is applied here. That is, the data y of the state variable result in the time-series data of the estimated values of the state variables generated by the filter operation unit 405. The estimated values of the state variables each vary according to the state of the railway vehicle.

[0144] First, the second frequency adjustment unit 406 generates the autocorrelation matrix R using (49) Equation and (51) Equation based on the data y of the estimated values of the state variables and the preset numbers M and m.

[0145] Then, the second frequency adjustment unit 406 performs singular value decomposition on the autocorrelation matrix R, to thereby derive the orthogonal matrix U and the diagonal matrix Σ of Equation (52), and derives the eigenvalues σ₁₁ to σ_mm of the autocorrelation matrix R from the diagonal matrix Σ.

[0146] Then, the second frequency adjustment unit 406 chooses s pieces of the eigenvalues σ₁₁ to σ_ss from the largest from among the plural eigenvalues σ₁₁ to σ_mm of the autocorrelation matrix R as the eigenvalues of the autocorrelation matrix R to be used for deriving the coefficient α of the corrected AR model. s is preset for each state variable. For example, the railway vehicle is made to travel on the track 16, to then obtain the data y of the estimated value of each of the state variables in such a manner as explained so far. Then, a distribution of the eigenvalues of the autocorrelation matrix R is made individually for each state variable. From the distributions of the eigenvalues of the autocorrelation matrix R, the number s of the eigenvalues to be extracted from the autocorrelation matrix R illustrated in (53) Equation is determined for each of the state variables.

[0147] Then, the second frequency adjustment unit 406 determines the coefficient α of the corrected AR model using (57) Equation based on the data y of the estimated value of the state variable, the eigenvalues σ₁₁ to σ_ss, and the orthogonal matrix U obtained by the singular value decomposition of the autocorrelation matrix R.

[0148] Then, the second frequency adjustment unit 406 derives the predicted value y^_k of the data y of the estimated value of the state variable at the time k from (45) Equation based on the coefficient α of the corrected AR model and the data y of the estimated value of the state variable. Time-series data of the predicted value y^_k of the state variable result in the time-series data from which the low-frequency component contained in the data y of the estimated value of the state variable has been extracted.

[0149] Then, the second frequency adjustment unit 406 subtracts the time-series data of the predicted value y^_k of the state variable from the data y of the estimated value of the state variable. In the following explanation, time-series data resulting from the subtraction of the time-series data of the predicted value y^_k of the state variable from the data y of the estimated value of the state variable are referred to as time-series data of a high-frequency component of the state variable as necessary.

[First track state calculation unit 407, S705]

[0150] When (22) Equation to (25) Equation are substituted into the motion equations that describe the yawings of the wheel sets 13a to 13d of (5) Equation to (8) Equation, (60) Equation to (63) Equation below are obtained.
[Mathematical equation 31]

[0151] In this embodiment, as illustrated in (60) Equation to (63) Equation, relational expressions representing the relations between the forward-and-backward-direction forces T₁ to T₄ and the alignment irregularity amounts y_R1 to y_R4 at the positions of the wheel sets 13a to 13d are set.

[0152] The first track state calculation unit 407 calculates estimated values of the pivot amounts (angular displacements) φ_w1 to φ_w4 of the wheel sets 13a to 13d in the yawing direction by (30) Equation to (33) Equation. Then, the first track state calculation unit 407 gives the estimated values of the pivot amounts (angular displacements) φ_w1 to φ_w4 of the wheel sets 13a to 13d in the yawing direction, the values of the high-frequency components of the state variables generated in the second frequency adjustment unit 406, and the values of the high-frequency components of the forward-and-backward-direction forces T₁ to T₄ generated in the first frequency adjustment unit 404 to (60) Equation to (63) Equation, to thereby calculate the alignment irregularity amounts y_R1 to y_R4 at the positions of the wheel sets 13a to 13d. The state variables to be used here are the displacements y_t1 to y_t2 of the bogies 12a, 12b in the right and left direction, the velocities y_t1 · to y_t2 · of the bogies 12a, 12b in the right and left direction, the displacements y_w1 to y_w4 of the wheel sets 13a to 13d in the right and left direction, and the velocities y_w1 · to y_w4 · of the wheel sets 13a to 13d in the right and left direction. The first track state calculation unit 407 performs such a calculation of the alignment irregularity amounts y_R1 to y_R4 as above with a predetermined sampling period, to thereby obtain the time-series data of the alignment irregularity amounts y_R1 to Y_R4.

[0153] Then, the first track state calculation unit 407 calculates an alignment irregularity amount y_R from the alignment irregularity amounts y_R1 to y_R4. For example, the first track state calculation unit 407 matches phases of the time-series data of the alignment irregularity amounts y_R2 to y_R4 to a phase of the time-series data of the alignment irregularity amount y_R1. That is, the first track state calculation unit 407 calculates, from the distance between the wheel set 13a and the wheel sets 13b to 13d in the forward and backward direction and the velocity of the railway vehicle, a delay time between the time when the wheel set 13a passes through a certain position and the time when the wheel sets 13b to 13d pass through the certain position. The first track state calculation unit 407 displaces the phases of the time-series data of the alignment irregularity amounts y_R2 to y_R4 by this delay time.

[0154] The first track state calculation unit 407 calculates an arithmetic mean value of the sum of the values of the alignment irregularity amounts y_R1 to y_R4 whose phases are matched at the same sampling time as the alignment irregularity amount y_R at this sampling time. The first track state calculation unit 407 performs such a calculation at each sampling time, to thereby obtain time-series data of the alignment irregularity amount y_R. The phases of the alignment irregularity amounts y_R2 to y_R4 are matched to the phase of the alignment irregularity amount y_R1, thereby making it possible to cancel disturbance factors existing in common in the time-series data of the alignment irregularity amounts y_R1 to y_R4.

[0155] Incidentally, the first track state calculation unit 407 may find a moving average of each of the alignment irregularity amounts y_R1 to y_R4 whose phases are matched (namely, pass each of the alignment irregularity amounts y_R1 to y_R4 through a low-pass filter) and calculate the alignment irregularity amount y_R from the alignment irregularity amounts y_R1 to y_R4 whose moving averages have been found.

[0156] Further, the first track state calculation unit 407 may calculate, as the alignment irregularity amount y_R, an arithmetic mean value of two of the values of the alignment irregularity amounts y_R1 to y_R4 whose phases are matched at the same sampling time from which the maximum value and the minimum value are removed.

[0157] The inspection apparatus 400 uses the measured data at each sampling time acquired by the data acquisition unit 403 while the railway vehicle is traveling in the entire traveling section of the railway vehicle, to execute pieces of the processing of the first frequency adjustment unit 404, the filter operation unit 405, the second frequency adjustment unit 406, and the first track state calculation unit 407.

[0158] In this manner, the first track state calculation unit 407 can obtain the alignment irregularity amount y_R at each sampling time while the railway vehicle is traveling in the entire traveling section. The first track state calculation unit 407 calculates a traveling position of the railway vehicle at each sampling time based on, for example, a traveling velocity of the railway vehicle and an elapsed time from the time when the railway vehicle starts to travel. In this embodiment, the case where the traveling position of the railway vehicle is the position of the wheel set 13a is explained as an example. The first track state calculation unit 407 calculates the alignment irregularity amount y_R at each traveling position of the railway vehicle based on the alignment irregularity amount y_R at each sampling time and the traveling position of the railway vehicle at each sampling time. In the following explanation, the value calculated in this manner is referred to as an estimated value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle or an estimated value of the alignment irregularity amount as necessary.

[0159] Incidentally, the first track state calculation unit 407 does not always need to calculate the traveling position of the railway vehicle at each sampling time as described previously. For example, the first track state calculation unit 407 may derive the traveling position of the railway vehicle at each sampling time by using a GPS (Global Positioning System) .

[Actual value acquisition unit 408, S706]

[0160] The actual value acquisition unit 408 acquires an actual measured value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle. The actual measured value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle is set to be measured before the second preprocessing is started. The timing for acquiring the actual measured value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle is not limited to a period between Step S705 and Step S707. The timing for acquiring the actual measured value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle may be any timing as long as it is the timing prior to Step S707. For example, the actual value acquisition unit 408 may acquire the actual measured value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle before the flowchart in Fig. 7 is started. In the following explanation, the actual measured value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle is referred to as an actual measured value of the alignment irregularity amount or an actual measured value as necessary.

[0161] The actual measured value of the alignment irregularity amount is a value to be obtained by directly measuring the alignment irregularity amount. The actual measured value of the alignment irregularity amount can be obtained as follows, for example. A test vehicle provided with a sensor that directly measures the alignment irregularity amount is made to travel. The direct measurement of the alignment irregularity amount by the sensor during traveling of the test vehicle is performed repeatedly in a predetermined cycle, to thereby obtain the alignment irregularity amount in the entire traveling section of the railway vehicle. Besides, the actual measured value of the alignment irregularity amount can also be obtained by using the measuring device described in Patent Literature 2, for example. As above, the actual measured value of the alignment irregularity amount can be obtained by a well-known technique. Thus, its detailed explanation is omitted here.

[0162] Fig. 12A and Fig. 12B to Fig. 14A and Fig. 14B are views illustrating first to sixth examples of the relation between the estimated value of the alignment irregularity amount (y_R), the actual value of the alignment irregularity amount (y_R), the traveling velocity of the railway vehicle (v), and the curvature of the track 16 (rail) (1/R) and the distance from the starting point of the railway vehicle respectively. Incidentally, the estimated value of the alignment irregularity amount is the one calculated by the first track state calculation unit 407. The actual value of the alignment irregularity amount is the one acquired by the actual value acquisition unit 408. Further, in Fig. 12A and Fig. 12B to Fig. 14A and Fig. 14B, for convenience of notation, illustration of data of a portion where the distance from the starting point of the railway vehicle is small is omitted.

[0163] In Fig. 12A and Fig. 12B to Fig. 14A and Fig. 14B, graphs 1211, 1221, 1311, 1321, 1411, and 1421 each indicate the estimated value of the alignment irregularity amount calculated by the first track state calculation unit 407. Graphs 1212, 1222, 1312, 1322, 1412, and 1422 each indicate the actual value of the alignment irregularity amount acquired by the actual value acquisition unit 408. Graphs 1213, 1223, 1313, 1323, 1413, and 1423 each indicate the traveling velocity of the railway vehicle. Graphs 1214, 1224, 1314, 1324, 1414, and 1424 each indicate the curvature 1/R of the track 16 (rail).

[0164] In Fig. 12A and Fig. 12B to Fig. 14A and Fig. 14B, the fact that the curvature 1/R is 0 (zero) indicates the linear track, and the fact that the curvature 1/R is a value other than 0 (zero) indicates the curved track.

[0165] Fig. 12A and Fig. 12B illustrate that the graphs 1214 and 1224 are the same and the traveling sections are the same. Fig. 12A and Fig. 12B illustrate that the traveling velocities of the railway vehicle are different as illustrated in the graphs 1213, 1223. This reveals that due to the different traveling velocities of the railway vehicle in the same traveling section, as illustrated in the graphs 1211, 1221, the estimated values of the alignment irregularity amount are different but the difference is not so large.

[0166] Further, the graphs 1212 and 1222 (the actual values of the alignment irregularity amount) are the same. As illustrated in the graphs 1211, 1212, it becomes clear that there is a difference between the estimated value of the alignment irregularity amount calculated by the first track state calculation unit 407 and the actual value of the alignment irregularity amount acquired by the actual value acquisition unit 408. The same is true of the graphs 1221, 1222.

[0167] As above, it becomes clear that the estimation accuracy of the alignment irregularity amount decreases depending on the traveling state of the railway vehicle and the installation state of the track 16.

[0168] As illustrated in the graphs 1214, 1224, the traveling sections illustrated in Fig. 12A and Fig. 12B are a sharp curve with a curvature radius R of 171 m. Therefore, the railway vehicle is in flange contact with the track 16.

[0169] Here, the flange contact is explained. Fig. 15 is a view explaining one example of the flange contact. Fig. 15 illustrates a cross section of the case where left and right rails of a track 16 and a single wheel set 13 are cut perpendicularly to the traveling direction of a railway vehicle (x axis direction). Further, Fig. 15 illustrates a state of the wheel set 13 in the case where the track 16 (rail) is curved to the right (in the negative direction of the y axis) and the railway vehicle is traveling while curving to the right. Incidentally, Fig. 15 also illustrates a lateral creep force F_y^L_i and a normal load N^L_i in a left wheel 14L and a lateral creep force F_y^R_i and a normal load N^R_i in a right wheel 14R.

[0170] As illustrated in Fig. 15, when traveling on the rail curved to the right, the railway vehicle receives an acting force in the left direction (the positive direction of the y axis) and the wheel set 13 moves to the left, and thereby reaction forces in the left and right direction from the contact position between the wheel 14L and the rail and the contact position between the wheel 14R and the rail increase to reach a force balance point. When this acting force further increases, the wheel set 13 further moves to the left, and when a contact angle α^L becomes the same as a flange angle a^L of the left wheel 14L, the left wheel 14L is to come into contact with the rail at a flange, as illustrated in Fig. 15. Such a contact is referred to as the flange contact. In the meantime, in this state, the right wheel 14R comes into contact with the rail at a tread.

[0171] Fig. 13A and Fig. 13B illustrate that the graphs 1314 and 1324 are the same and the traveling sections are the same. As illustrated in the graphs 1314, 1324, the curvature 1/R is 0 (zero), and thus it becomes clear that the traveling sections illustrated in Fig. 13A and Fig. 13B are a linear track. Fig. 13A and Fig. 13B illustrate that the traveling velocities of the railway vehicle are different as illustrated in the graphs 1313, 1323. This reveals that due to the different traveling velocities of the railway vehicle in the same traveling section, as illustrated in the graphs 1311, 1321, the estimated values of the alignment irregularity amount are different but the difference is not so large. Further, in Fig. 13A and Fig. 13B, an S/N ratio of the measured value of the forward-and-backward-direction force decreases because the traveling velocity of the railway vehicle decreases to 30 km/h or less. Thus, as illustrated in the graphs 1311, 1321, a high-frequency noise is mixed into the estimated value of the alignment irregularity amount. However, it becomes clear that a feature amount of the alignment irregularity amount (such as a way how the graph changes) is captured in the graphs 1311, 1321.

[0172] Further, the graphs 1312 and 1322 (the actual values of the alignment irregularity amount) are the same. As illustrated in the graphs 1311, 1312, it becomes clear that there is a difference between the estimated value of the alignment irregularity amount calculated by the first track state calculation unit 407 and the actual value of the alignment irregularity amount acquired by the actual value acquisition unit 408. The same is true of the graphs 1321, 1322.

[0173] As above, it becomes clear that the estimation accuracy of the alignment irregularity amount decreases depending on the traveling state of the railway vehicle.

[0174] Fig. 14A and Fig. 14B illustrate that the graphs 1414 and 1424 are the same and the traveling sections are the same. Fig. 14A and Fig. 14B illustrate that the traveling velocities of the railway vehicle are different as illustrated in the graphs 1413, 1423. This reveals that due to the different traveling velocities of the railway vehicle in the same traveling section, as illustrated in the graphs 1411, 1421, the estimated values of the alignment irregularity amount are different but the difference is not so large. Further, the graphs 1412 and 1422 (the actual values of the alignment irregularity amount) are the same. As illustrated in the graphs 1411, 1412, it becomes clear that there is a difference between the estimated value of the alignment irregularity amount calculated by the first track state calculation unit 407 and the actual value of the alignment irregularity amount acquired by the actual value acquisition unit 408. The same is true of the graphs 1421, 1422.

[0175] As illustrated in the graphs 1414, 1424, the traveling sections illustrated in Fig. 14A and Fig. 14B are a gentle curve with the curvature radius R of 993 m, and the railway vehicle does not come into flange contact with the track 16.

[0176] As above, it becomes clear that the estimation accuracy of the alignment irregularity amount decreases depending on the installation state of the track 16.

[Correction amount calculation unit 409, correction amount storage unit 410, S707 to S711]

[0177] The correction amount calculation unit 409 calculates, when the estimated value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle is calculated by the first track state calculation unit 407, a correction amount at each position in the entire traveling section of the railway vehicle. The correction amount at each position in the entire traveling section of the railway vehicle is the correction amount for the estimated value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle, which is calculated by the later-described second track state calculation unit 411.

[0178] The correction amount calculation unit 409 calculates the correction amount at each position in the entire traveling section of the railway vehicle based on the estimated value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle that is calculated by the first track state calculation unit 407 and the actual measured value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle that is acquired by the actual value acquisition unit 408.

[0179] In this embodiment, the correction amount calculation unit 409 calculates the correction amount at each position in the entire traveling section of the railway vehicle as follows.

[0180] The correction amount calculation unit 409 extracts a pair of values at the same position, which is a pair of the estimated value of the alignment irregularity amount that is calculated by the first track state calculation unit 407 and the actual measured value of the alignment irregularity amount that is acquired by the actual value acquisition unit 408. The correction amount calculation unit 409 calculates, as the correction amount at the position, a value obtained by subtracting the extracted actual measured value of the alignment irregularity amount from the extracted estimated value of the alignment irregularity amount. The correction amount calculation unit 409 performs such a calculation of the correction amount by using the estimated value of the alignment irregularity amount and the actual measured value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle. In this manner, the correction amount at each position in the entire traveling section of the railway vehicle is calculated.

[0181] In this embodiment, when the railway vehicle travels through the entire traveling section once, the correction amount calculation unit 409 calculates one set of correction amounts at all the positions in the entire traveling section of the railway vehicle (Step S707).

[0182] Incidentally, the correction amount calculation unit 409 performs interpolation processing on the correction amount at each position in the entire traveling section of the railway vehicle, thereby being able to calculate the correction amounts at all the positions of the entire traveling section of the railway vehicle.

[0183] In the following explanation, the correction amount at each position in the entire traveling section of the railway vehicle, which is obtained in this manner by the railway vehicle traveling through the entire traveling section once, is referred to as a first correction amount at each position in the entire traveling section of the railway vehicle or a first correction amount as necessary.

[0184] The first correction amount at each position in the entire traveling section of the railway vehicle may be set to the correction amount for the estimated value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle. However, in this embodiment, the correction amount calculation unit 409 calculates, as the first correction amount at a certain position in the entire traveling section of the railway vehicle, the correction amount for the estimated value of the alignment irregularity amount at the certain position by using a plurality of first correction amounts. This is because it is possible to improve the accuracy of the correction amount for the estimated value of the alignment irregularity amount.

[0185] In this embodiment, as one example, an addition mean value of a plurality of the first correction amounts is set to the correction amount for the estimated value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle. In the following explanation, the correction amount for the estimated value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle, which is obtained in this manner by using a plurality of the first correction amounts, is referred to as a second correction amount at each position in the entire traveling section of the railway vehicle, or a second correction amount as necessary.

[0186] The correction amount calculation unit 409 stores, when the first correction amount at each position in the entire traveling section of the railway vehicle is obtained, this first correction amount at each position in the entire traveling section of the railway vehicle temporarily (Step S708).

[0187] Then, the correction amount calculation unit 409 determines whether or not a predetermined number of first correction amounts necessary for calculating the addition mean value have been obtained (Step S709). The predetermined number may be any number as long as it is two or more. As a result of this determination, in the case where a predetermined number of first correction amounts necessary for calculating the addition mean value have not been obtained (NO at Step S709), the inspection apparatus 400 performs Steps S701 to S708 described previously when the railway vehicle is traveling in the entire traveling section again, and stores a new first correction amount.

[0188] When a predetermined number of the first correction amounts necessary for calculating the addition mean value are obtained as above (YES at Step S709), the correction amount calculation unit 409 calculates the addition mean value of a predetermined number of the first correction amounts as the second correction amount (Step S710). The correction amount storage unit 410 stores the second correction amount (Step S711). As described previously, the second correction amount is the correction amount for the estimated value of the alignment irregularity amount, and is used in the later-described track state correction unit 412.

[0189] After the second correction amount is stored in the correction amount storage unit 410 as above, the data acquisition unit 403, the first frequency adjustment unit 404, the filter operation unit 405, the second frequency adjustment unit 406, the second track state calculation unit 411, the track state correction unit 412, and the output unit 413 start. That is, after the second preprocessing by the flowchart in Fig. 7 is finished, the main processing by the flowchart in Fig. 8 is started. At the time of the main processing, the first track state calculation unit 407, the actual value acquisition unit 408, and the correction amount calculation unit 409 do not start. Further, the flowchart in Fig. 8 is executed repeatedly every time the sampling time arrives.

[0190] Fig. 16A to Fig. 16C are views illustrating first to third examples of the relation between a second correction amount M and the distance from the starting point of the railway vehicle respectively. Fig. 16A illustrates the second correction amount M obtained from the results illustrated in Fig. 12A and Fig. 12B. Fig. 16B illustrates the second correction amount M obtained from the results illustrated in Fig. 13A and Fig. 13B. Fig. 16C illustrates the second correction amount M obtained from the results illustrated in Fig. 14A and Fig. 14B.

[Data acquisition unit 403, S801]

[0191] The data acquisition unit 403 acquires measured data with a predetermined sampling period. At Step S801, the data acquisition unit 403 acquires one set of measured data at a sampling time. Incidentally, the measured data to be acquired by the data acquisition unit 403 are the same in objects to be measured as the measured data to be acquired at Step S701, and thus, their detailed explanation is omitted here.

[First frequency adjustment unit 404, S802]

[0192] The first frequency adjustment unit 404 reduces (preferably removes) the signal strength of the low-frequency component contained in the time-series data of the measured value of the forward-and-backward-direction force out of the measured data acquired by the data acquisition unit 403. Incidentally, the processing at Step S802 is the same as that at Step S702, and thus its detailed explanation is omitted here .

[0193] However, at Step S702, the first frequency adjustment unit 404 derives the time-series data from which the low-frequency component contained in the data y of the measured value of the forward-and-backward-direction force has been extracted after the measured data in the entire traveling section of the railway vehicle are obtained. In contrast to this, at Step S802, the first frequency adjustment unit 404 derives the time-series data from which the low-frequency component contained in the data y of the measured value of the forward-and-backward-direction force has been extracted every time the data acquisition unit 403 acquires the value y_k of the data y of the measured value of the forward-and-backward-direction force at the time k with a predetermined sampling period.

[Filter operation unit 405, S803]

[0194] The filter operation unit 405 sets the observation equation as the observation equation stored by the observation equation storage unit 402, sets the state equation as the state equation stored by the state equation storage unit 401, and determines estimated values of the state variables illustrated in (44) Equation by the Kalman filter. Incidentally, the processing at Step S803 is the same as that at Step S703, and thus its detailed explanation is omitted here .

[Second frequency adjustment unit 406, S804]

[0195] The second frequency adjustment unit 406 reduces (preferably removes) the signal strength of the low-frequency component contained in the time-series data of the estimated values of the state variables generated by the filter operation unit 405. Incidentally, the processing at Step S804 is the same as that at Step S704, and thus its detailed explanation is omitted here.

[Second track state calculation unit 411, S805]

[0196] The second track state calculation unit 411 calculates the alignment irregularity amounts y_R1 to y_R4 and calculates, as the estimated value of the alignment irregularity amount, the alignment irregularity amount y_R from the alignment irregularity amounts y_R1 to y_R4. The processing at Step S805 is the same as that at Step S705, and thus its detailed explanation is omitted here. However, at Step S705, the second track state calculation unit 411 calculates the estimated value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle. In contrast to this, at Step S805, the second track state calculation unit 411 calculates the estimated value of the alignment irregularity amount at a traveling position of the railway vehicle corresponding to the current sampling time.

[Track state correction unit 412, S806]

[0197] The track state correction unit 412 reads the second correction amount at the traveling position of the railway vehicle corresponding to the current sampling time from the correction amount storage unit 410. The track state correction unit 412 uses the second correction amount at the traveling position of the railway vehicle corresponding to the current sampling time, which is read from the correction amount storage unit 410, to correct the estimated value of the alignment irregularity amount at the traveling position of the railway vehicle corresponding to the current sampling time, which is calculated by the second track state calculation unit 411.

[0198] In this embodiment, the track state correction unit 412 subtracts the second correction amount at the traveling position of the railway vehicle corresponding to the current sampling time that is read from the correction amount storage unit 410 from the estimated value of the alignment irregularity amount at the traveling position of the railway vehicle corresponding to the current sampling time that is calculated by the second track state calculation unit 411, to thereby correct the estimated value of the alignment irregularity amount at the traveling position of the railway vehicle corresponding to the current sampling time that is calculated by the second track state calculation unit 411. In the following explanation, the estimated value of the alignment irregularity amount at the traveling position of the railway vehicle corresponding to the current sampling time, which is corrected in this manner, is referred to as a corrected estimated value of the alignment irregularity amount as necessary. The corrected estimated value of the alignment irregularity amount becomes the estimated value of the final alignment irregularity amount.

[0199] Incidentally, in the case where the correction amount calculation unit 409 sets the second correction amount to the value obtained by subtracting the estimated value of the alignment irregularity amount from the actual measured value of the alignment irregularity amount, the track state correction unit 412 corrects the estimated value of the alignment irregularity amount at the traveling position of the railway vehicle corresponding to the current sampling time that is calculated by the second track state calculation unit 411 as follows. That is, the track state correction unit 412 adds the estimated value of the alignment irregularity amount at the traveling position of the railway vehicle corresponding to the current sampling time that is calculated by the second track state calculation unit 411 and the second correction amount at the traveling position of the railway vehicle corresponding to the current sampling time that is read from the correction amount storage unit 410 together, to thereby correct the estimated value of the alignment irregularity amount at the traveling position of the railway vehicle corresponding to the current sampling time that is calculated by the second track state calculation unit 411.

[0200] Fig. 17A and Fig. 17B to Fig. 19A and Fig. 19B are views illustrating first to sixth examples of the relation between the corrected estimated value of the alignment irregularity amount and the distance from the starting point of the railway vehicle respectively. Incidentally, here, for simplicity, it is assumed that the estimated value of the alignment irregularity amount calculated by the second track state calculation unit 411 is the same as the estimated value of the alignment irregularity amount calculated by the first track state calculation unit 407.

[0201] That is, graphs 1711 and 1721 in Fig. 17A and Fig. 17B illustrate corrected estimated values of the alignment irregularity amount obtained by correcting the estimated values of the alignment irregularity amount (the graphs 1211, 1221) illustrated in Fig. 12A and Fig. 12B by the correction amount M illustrated in Fig. 16A respectively. Further, the graphs 1212, 1222, 1712, and 1722 (the actual values of the alignment irregularity amount) are the same.

[0202] Graphs 1811 and 1821 in Fig. 18A and Fig. 18B illustrate corrected estimated values of the alignment irregularity amount obtained by correcting the estimated values of the alignment irregularity amount (the graphs 1311, 1321) illustrated in Fig. 13A and Fig. 13B by the correction amount M illustrated in Fig. 16B respectively. Further, the graphs 1312, 1322, 1812, and 1822 (the actual values of the alignment irregularity amount) are the same.

[0203] Graphs 1911 and 1921 in Fig. 19A and Fig. 19B illustrate corrected estimated values of the alignment irregularity amount obtained by correcting the estimated values of the alignment irregularity amount (the graphs 1411, 1421) illustrated in Fig. 14A and Fig. 14B by the correction amount M illustrated in Fig. 16C respectively. Further, the graphs 1412, 1422, 1912, and 1922 (the actual values of the alignment irregularity amount) are the same.

[0204] As illustrated in Fig. 17A and Fig. 17B to Fig. 19A and Fig. 19B, it becomes clear that in all the cases, the corrected estimated value of the alignment irregularity amount agrees with the actual measured value with high accuracy.

[Output unit 413, S807]

[0205] The output unit 413 outputs information of the corrected estimated value of the alignment irregularity amount that is calculated by the track state correction unit 412. At this time, the output unit 413 may output information indicating that the track 16 is abnormal in the where the corrected estimated value of the alignment irregularity amount is larger than a preset value. As a form of output, it is possible to employ at least any one of displaying the information on a computer display, transmitting the information to an external device, and storing the information in an internal or external storage medium, for example.

<Summary>

[0206] In this embodiment as above, the inspection apparatus 400 acquires the measured values of the forward-and-backward-direction forces T₁ to T₄ by making the railway vehicle travel. The inspection apparatus 400 uses the measured values of the forward-and-backward-direction forces T₁ to T₄ and the relational expression between the forward-and-backward-direction forces T₁ to T₄ and the alignment irregularity amounts y_R1 to y_R4 at the positions of the wheel sets 13a to 13d to derive the estimated value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle. The inspection apparatus 400 uses the estimated value and the actual measured value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle to calculate the second correction amount as the correction amount for the estimated value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle. Thereafter, the inspection apparatus 400 derives the estimated value of the alignment irregularity amount at the traveling position of the railway vehicle as described previously by making the railway vehicle travel. The inspection apparatus 400 corrects the estimated value of the alignment irregularity amount at the traveling position of the railway vehicle, which is derived in this manner, by the second correction amount at this traveling position. Accordingly, it is possible to detect the irregularity in the track 16 of the railway vehicle with high accuracy without using a special measuring apparatus.

[0207] Further, in this embodiment, the inspection apparatus 400 reduces the signal strength of the low-frequency components contained in the time-series data of the measured values of the forward-and-backward-direction forces T₁ to T₄ and generates the time-series data of the high-frequency components of the forward-and-backward-direction forces T₁ to T₄. The inspection apparatus 400 gives the time-series data of the high-frequency components of the forward-and-backward-direction forces T₁ to T₄ to the relational expression between the forward-and-backward-direction forces T₁ to T₄ and the alignment irregularity amounts y_R1 to y_R4 at the positions of the wheel sets 13a to 13d, to thereby calculate the alignment irregularity amounts y_R1 to y_R4 at the positions of the wheel sets 13a to 13d. This relational expression is an expression based on the motion equations that describe the motions of the railway vehicle when traveling on the linear track (namely, the equations not including the curvature radius R of the track 16 (the rail)). Accordingly, it is possible to detect the irregularity in the curved track with high accuracy without using a special measuring apparatus.

[0208] Further, in this embodiment, the inspection apparatus 400 generates the autocorrelation matrix R from the data y of the measured value of the forward-and-backward-direction force, and by using s pieces of the eigenvalues from the largest chosen from the eigenvalues obtained by the singular value decomposition of the autocorrelation matrix R, determines the coefficient α of the corrected AR model approximating the data y of the measured value of the forward-and-backward-direction force. Accordingly, it is possible to determine the coefficient α so as to make the signal of the low-frequency component contained in the data y of the measured value of the forward-and-backward-direction force remain and prevent the high-frequency component from remaining. The inspection apparatus 400 calculates the predicted value y^_k of the forward-and-backward-direction force at the time k by giving the data y of the measured value of the forward-and-backward-direction force at the time k - 1 (1 ≦ 1 ≦ m), which is prior to the time k, to the corrected AR model whose coefficient α is determined in this manner. Accordingly, it is possible to reduce the signal of the low-frequency component, which is due to the railway vehicle traveling on the curved track, from the data y of the measured value of the forward-and-backward-direction force without estimating a cutoff frequency beforehand.

[0209] Further, in this embodiment, the inspection apparatus 400 gives, out of the measured data acquired in the data acquisition unit 403, the measured data excluding the forward-and-backward-direction forces T₁ to T₄ and the time-series data of the high-frequency components of the forward-and-backward-direction forces T₁ to T₄ generated in the first frequency adjustment unit 404 to the Kalman filter, to derive the state variables (y_w1 · to y_w4 ·, y_w1 to y_w4, y_t1 · to y_t2 ·, y_t1 to y_t2, φ_t1 · to φ_t2 ·, φ_t1 to φ_t2, φ_t1 · to φ_t2 ·, φ_t1 to φ_t2, y_b ·, y_b, φ_b ·, φ_b, φ_b · , φ_b, φ_y1, φ_y2, φ_a1, and φ_a2). Next, the inspection apparatus 400 reduces (preferably removes) the signal strength of the low-frequency components contained in the time-series data of the estimated values of the state variables, to thereby calculate the values of the high-frequency components of the state variables. Next, the inspection apparatus 400 uses the values of the high-frequency components of the pivot amounts (angular displacements) φ_t1 to φ_t2 of the bogies 12a, 12b in the yawing direction and the actual values of the transformation variables e₁ to e₄, to derive the pivot amounts (angular displacements) φ_w1 to φ_w4 of the wheel sets 13a to 13d in the yawing direction. Next, the inspection apparatus 400 substitutes the pivot amounts (angular displacements) φ_w1 to φ_w4 of the wheel sets 13a to 13d in the yawing direction, the values of the high-frequency components of the state variables, and the values of the high-frequency components of the forward-and-backward-direction forces T₁ to T₄ into the motion equations that describe the yawings of the wheel sets 13a to 13d, to calculate the alignment irregularity amounts y_R1 to y_R4 at the positions of the wheel sets 13a to 13d. Then, the inspection apparatus 400 calculates the alignment irregularity amount y_R from the alignment irregularity amounts y_R1 to y_R4 at the positions of the wheel sets 13a to 13d. Accordingly, it becomes no longer necessary to constitute the state equation by using the motion equations including the alignment irregularity amounts y_R1 to y_R4 at the positions of the wheel sets 13a to 13d as variables as the motion equations that describe the yawings of the wheel sets 13a to 13d. Thereby, it becomes no longer necessary to create a model of the track 16, and at the same time, it is possible to reduce the number of state variables. In this embodiment, it is possible to reduce the degrees of freedom of the model to 17 degrees of freedom from 21 degrees of freedom, and at the same time, it is possible to reduce the number of state variables to 30 from 38. Further, the number of measured values to be used in the Kalman filter increases by the forward-and-backward-direction forces T₁ to T₄.

[0210] On the other hand, when the motion equations that describe the yawings of the wheel sets 13a to 13d of (5) Equation to (8) Equation are included in the state equation without using the forward-and-backward-direction forces T₁ to T₄, the calculation sometimes becomes unstable, failing to obtain the estimated results. That is, unless the state variables are selected, the calculation sometimes becomes unstable, failing to obtain the estimated results. Further, even if the estimated results are obtained, the accuracy to detect the irregularity of the track 16 becomes higher in the method of this embodiment as compared to the method of not selecting the state variables. This is because in this embodiment, not including the motion equations that describe the yawings of the wheel sets 13a to 13d in the state equation and using the measured value of the forward-and-backward-direction force are achieved.

[0211] Further, the strain gauges can be used as sensors in this embodiment, thus not requiring special sensors. Accordingly, it is possible to accurately detect the irregularity of the track 16 (track irregularity) inexpensively. Further, since it is not necessary to use special sensors, the strain gauges are attached to a commercial vehicle and the inspection apparatus 400 is mounted on the commercial vehicle, thereby making it possible to detect the irregularity of the track 16 in real time during traveling of the commercial vehicle. Accordingly, it is possible to detect the irregularity of the track 16 without traveling of an inspection car. However, the strain gauges may be attached to the inspection car and the inspection apparatus 400 may be mounted on the inspection car.

<Modified example>

[0212] The case where the addition mean value of a plurality of the first correction amounts is set to the correction amount for the estimated value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle has been explained as an example in this embodiment. However, it is not always necessary to derive the correction amount for the estimated value of the alignment irregularity amount at each position in the entire traveling section of the railway vehicle in this manner.

[0213] For example, the inspection apparatus 400 calculates a plurality of first correction amounts as the first correction amount at the same position in a state where the traveling velocities of the railway vehicle are mutually different. The inspection apparatus 400 performs a regression analysis using a plurality of these first correction amounts and calculates coefficients of a regression formula. An objective variable of the regression formula is the second correction amount. An explanatory variable of the regression formula includes the traveling velocity of the railway vehicle. The inspection apparatus 400 derives such a regression formula at each position in the entire traveling section of the railway vehicle. Thereafter, the inspection apparatus 400 (track state correction unit 412) reads the regression formula corresponding to the traveling position of the railway vehicle corresponding to the current sampling time from the correction amount storage unit 410. Then, the inspection apparatus 400 (track state correction unit 412) substitutes the traveling velocity of the railway vehicle corresponding to the current sampling time into the regression formula to calculate the second correction amount.

[0214] Further, it is not necessary to calculate the correction amount for the estimated value of the alignment irregularity amount at a certain position in the entire traveling section of the railway vehicle by using a plurality of the first correction amounts. In this case, the correction amount for the estimated value of the alignment irregularity amount at a certain position in the entire traveling section of the railway vehicle is determined by the single first correction amount at the certain position. If the above is performed, the accuracy of the correction amount for the estimated value of the alignment irregularity amount may decrease. However, it becomes no longer necessary to make the railway vehicle travel for a plurality of times in the second preprocessing. For example, it is possible to determine which method to employ based on the combination of the accuracy of the correction amount for the estimated value of the alignment irregularity amount and the time and effort of the second preprocessing.

[0215] Further, the case where the data acquisition unit 403 acquires the measured data in the entire traveling section of the railway vehicle at Step S701 in the flowchart in Fig. 7 has been explained as an example. However, it is not always necessary to constitute this embodiment as above. For example, similarly to the flowchart in Fig. 8, at Step S701, the data acquisition unit 403 may acquire one set of measured data at the sampling time. In this case, the processing at Step S701 to the processing at Step S708 are performed repeatedly for each traveling position of the railway vehicle corresponding to the sampling time. This processing is performed repeatedly until the correction amount at each position in the entire traveling section of the railway vehicle (the first correction amount) is obtained.

[0216] Further, in this embodiment, the case where the measured values of the forward-and-backward-direction force used in the first track state calculation unit 407 and the second track state calculation unit 411 are the measured values in the same railway vehicle has been explained as an example. In this case, the inspection apparatus 400 that calculates the estimated value of the alignment irregularity amount and the inspection apparatus 400 that calculates the second correction amount result in the inspection apparatus 400 mounted on the same railway vehicle. This is preferable because it is possible to prevent errors caused by the characteristics inherent in the railway vehicle from being contained in the second correction amount. However, it is not always necessary to constitute this embodiment as above. The same second correction amount may be used for a plurality of railway vehicles of the same type that travel for the same traveling section, for example. Further, the same second correction amount may be used for a plurality of railway vehicles with the same route name.

[0217] Further, the case of using the corrected AR model has been explained as an example in this embodiment. However, it is not always necessary to use the corrected AR model and reduce the signal of the low-frequency component, which is due to the railway vehicle traveling on the curved track, from the data y of the measured value of the forward-and-backward-direction force. For example, in the case where it is possible to specify a frequency band due to the railway vehicle traveling on the curved track, a high-pass filter may be used to reduce the signal of the low-frequency component, which is due to the railway vehicle traveling on the curved track, from the data y of the measured value of the forward-and-backward-direction force.

[0218] Further, it is not always necessary to reduce the signal of the low-frequency component, which is due to the railway vehicle traveling on the curved track, from the data y of the measured value of the forward-and-backward-direction force. In the case of calculating the alignment irregularity amount of the linear track, for example, it is not necessary to do the above. In this case, the first frequency adjustment unit 404 and the second frequency adjustment unit 406 are no longer required.

[0219] Further, the case where the wheel set to be a standard when phase matching is the wheel set 13a has been explained as an example in this embodiment. However, the wheel set to be a standard may be the wheel set 13b, 13c, or 13d other than the wheel set 13a.

[0220] In this embodiment, the case of using the Kalman filter has been explained as an example. However, it is not always necessary to use the Kalman filter as long as a filter that derives the estimated values of the state variables so that the error between, of the observation variable, the measured value and the estimated value becomes minimum or the expected value of this error becomes minimum (that is, a filter that performs data assimilation) is used. For example, a particle filter may be used. Incidentally, as the error between, of the observation variable, the measured value and the estimated value, for example, a square error between, of the observation variable, the measured value and the estimated value is cited.

[0221] Further, in this embodiment, the case of deriving the alignment irregularity amount has been explained as an example. However, it is not always necessary to derive the alignment irregularity amount as long as a physical quantity that reflects the track irregularity (appearance failure of the track 16) is derived as the physical quantity (first physical quantity) reflecting the state of the track 16. For example, in addition to or in place of the alignment irregularity amount, calculations of (64) Equation to (67) Equation may be performed, to thereby derive a lateral force to occur when the railway vehicle travels on the linear track (stress in the right and left direction between the wheel and the rail). However, Q₁, Q₂, Q₃, and Q₄ are lateral forces in the wheels 14a, 14b, 14c, and 14d respectively. f₃ represents a spin creep coefficient.
[Mathematical equation 32]

[0222] Further, in this embodiment, the case of including the state variables that represent the state of the vehicle body 11 has been explained as an example. However, the vehicle body 11 is a part into which vibrations by acting forces between the wheels 14a to 14d and the track 16 (creep force) propagate finally. Accordingly, it is not necessary to include the state variables representing the state of the vehicle body 11 in the case where the effect by the propagation in the vehicle body 11 is judged to be small, for example. In such a case, out of the motion equations of (1) Equation to (21) Equation, the motion equations that describe the transversal vibration, the yawing, and the rolling of the vehicle body 11 of (15) Equation to (17) Equation and the motion equations that describe the yawing of the yaw damper disposed on the bogie 12a and the yawing of the yaw damper disposed on the bogie 12b of (18) Equation and (19) Equation are no longer required. Further, in the motion equations of (1) Equation to (21) Equation, the value inside {} that includes the state amount relating to the vehicle body (state amount including the subscript of b) and the state amount relating to the vehicle body (state amount including the subscript of b) (for example, the third term {φ_a2 - φ_b} on the left side of (21) Equation)) is set to 0 (zero).

[0223] Further, in this embodiment, the case of the bogies 12a, 12b each being a bolsterless bogie has been explained as an example. However, the bogies 12a, 12b are not limited to the bolsterless bogie. Besides, according to the components of the railway vehicle, forces that the railway vehicle receives, the directions of the motions of the railway vehicle, or the like, the motion equations are rewritten appropriately. That is, the motion equations are not limited to the ones explained in this embodiment as an example.

(Second embodiment)

[0224] Next, there will be explained a second embodiment.

[0225] In the first embodiment, the case where the inspection apparatus 400 mounted on the railway vehicle calculates and corrects the estimated value of the alignment irregularity amount has been explained as an example. In contrast to this, in this embodiment, a data processing device in which some functions of the inspection apparatus 400 are mounted is disposed in an operation center. The data processing device receives measured data transmitted from the railway vehicle and calculates and corrects the estimated value of the alignment irregularity amount by using the received measured data. In this manner, in this embodiment, the functions that the inspection apparatus 400 in the first embodiment has are shared and executed by the railway vehicle and the operation center. Constitutions and processing due to this are mainly different between this embodiment and the first embodiment. Accordingly, in the explanation of this embodiment, the same reference numerals and symbols as those added to Fig. 1 to Fig. 19A and Fig. 19B are added to the same parts as those in the first embodiment, or the like, and their detailed explanations are omitted.

[0226] Fig. 20 is a view illustrating one example of a configuration of an inspection system. In Fig. 20, the inspection system includes data collecting devices 2010a, 2010b, and a data processing device 2020. In Fig. 20, one example of functional configurations of the data collecting devices 2010a, 2010b and the data processing device 2020 is also illustrated. Incidentally, each hardware of the data collecting devices 2010a, 2010b and the data processing device 2020 can be fabricated by the one illustrated in Fig. 5, for example. Accordingly, detailed explanations of the hardware configurations of the data collecting devices 2010a, 2010b and the data processing device 2020 are omitted.

[0227] The data collecting devices 2010a, 2010b are mounted on each railway vehicle one by one. The data processing device 2020 is disposed at the operation center. The operation center centrally manages operations of a plurality of railway vehicles, for example.

<Data collecting devices 2010a, 2010b>

[0228] The data collecting devices 2010a, 2010b can be fabricated by the same components. The data collecting devices 2010a, 2010b include data acquisition units 2011a, 2011b and data transmission units 2012a, 2012b.

[Data acquisition units 2011a, 2011b]

[0229] The data acquisition units 2011a, 2011b have the same function as that of the data acquisition unit 403. That is, the data acquisition units 2011a, 2011b acquire the same measured data as those acquired in the data acquisition unit 403. Concretely, the data acquisition units 2011a, 2011b acquire, as the measured data, a measured value of acceleration of the vehicle body 11 in the right and left direction, measured values of accelerations of the bogies 12a, 12b in the right and left direction, measured values of accelerations of the wheel sets 13a to 13d in the right and left direction, and a measured value of the forward-and-backward-direction force. Strain gauges and an arithmetic device for obtaining these measured values are the same as those explained in the first embodiment.

[Data transmission units 2012a, 2012b]

[0230] The data transmission units 2012a, 2012b transmit the measured data acquired in the data acquisition units 2011a, 2011b to the data processing device 2020. In this embodiment, the data transmission units 2012a, 2012b transmit the measured data acquired in the data acquisition units 2011a, 2011b to the data processing device 2020 by radio. At this time, the data transmission units 2012a, 2012b add identification numbers of the railway vehicles in which the data collecting devices 2010a, 2010b are mounted to the measured data acquired in the data acquisition units 2011a, 2011b. In this manner, the data transmission units 2012a, 2012b transmit the measured data with the identification numbers of the railway vehicles added thereto.

<Data processing device 2020>

[Data reception unit 2021]

[0231] A data reception unit 2021 receives the measured data transmitted by the data transmission units 2012a, 2012b. To the measured data, the identification numbers of the railway vehicles, which are transmission sources of the measured data, have been added.

[Data storage unit 2022]

[0232] A data storage unit 2022 stores the measured data received in the data reception unit 2021. The data storage unit 2022 stores the measured data every identification number of the railway vehicle. The data storage unit 2022 specifies the traveling position of the railway vehicle at the time of receipt of the measured data based on the current operation situation of the railway vehicle and the time of receipt of the measured data, and stores information of the specified traveling position and the measured data in association with each other. Incidentally, the data collecting devices 2010a, 2010b may collect the information of the current traveling position of the railway vehicle and contain the collected information in the measured data.

[Data reading unit 2023]

[0233] A data reading unit 2023 reads the measured data stored in the data storage unit 2022. The data reading unit 2023 can read, out of the measured data stored in the data storage unit 2022, the measured data designated by an operator. Further, the data reading unit 2023 can also read the measured data matching a preset condition at a preset timing. In this embodiment, the measured data read by the data reading unit 2023 are determined based on at least any one of the identification number and the traveling position of the railway vehicle, for example.

[0234] A state equation storage unit 401, an observation equation storage unit 402, a first frequency adjustment unit 404, a filter operation unit 405, a second frequency adjustment unit 406, a first track state calculation unit 407, an actual value acquisition unit 408, a correction amount calculation unit 409, a correction amount storage unit 410, a second track state calculation unit 411, a track state correction unit 412, and an output unit 413 are the same as those explained in the first embodiment. Accordingly, their detailed explanations are omitted here. Incidentally, the filter operation unit 405 uses the measured data read by the data reading unit 2023 in place of using the measured data acquired in the data acquisition unit 403, and determines the estimated values of the state variables illustrated in (44) Equation.

<Summary>

[0235] As above, in this embodiment, the data collecting devices 2010a, 2010b mounted on the railway vehicles collect the measured data to transmit them to the data processing device 2020. The data processing device 2020 disposed at the operation center stores the measured data received from the data collecting devices 2010a, 2010b and uses the stored measured data to calculate the estimated value of the alignment irregularity amount. Accordingly, in addition to the effects explained in the first embodiment, for example, the following effects are exhibited. That is, the data processing device 2020 can calculate the final alignment irregularity amount y_R at an arbitrary timing by reading the measured data at an arbitrary timing. Further, the data processing device 2020 can output time-series variation of the estimated value of the final alignment irregularity amount at the same position. Further, the data processing device 2020 can output the estimated values of the alignment irregularity amounts in a plurality of routes for each route.

<Modified example>

[0236] In this embodiment, the case where the data collecting devices 2010a, 2010b directly transmit the measured data to the data processing device 2020 has been explained as an example. However, it is not always necessary to constitute this embodiment as above. An inspection system may be built by using cloud computing, for example.

[0237] Besides, it is possible to employ the various modified examples explained in the first embodiment also for this embodiment.

[0238] Further, in the first embodiment, the case where the state equation storage unit 401, the observation equation storage unit 402, the data acquisition unit 403, the first frequency adjustment unit 404, the filter operation unit 405, the second frequency adjustment unit 406, the first track state calculation unit 407, the actual value acquisition unit 408, the correction amount calculation unit 409, the correction amount storage unit 410, the second track state calculation unit 411, the track state correction unit 412, and the output unit 413 are included in one apparatus has been explained as an example. However, it is not always necessary to constitute the first embodiment as above. Functions of the state equation storage unit 401, the observation equation storage unit 402, the data acquisition unit 403, the first frequency adjustment unit 404, the filter operation unit 405, the second frequency adjustment unit 406, the first track state calculation unit 407, the actual value acquisition unit 408, the correction amount calculation unit 409, the correction amount storage unit 410, the second track state calculation unit 411, the track state correction unit 412, and the output unit 413 may be fabricated by a plurality of apparatuses. In this case, the inspection system is constituted by using these plural apparatuses.

(Another embodiment)

[0239] The embodiments of the present invention explained above can be fabricated by causing a computer to execute a program. Further, a computer-readable recording medium in which the aforementioned program is recorded and a computer program product such as the aforementioned program can also be applied as the embodiment of the present invention. As the recording medium, it is possible to use a flexible disk, a hard disk, an optical disk, a magneto-optic disk, a CD-ROM, a magnetic tape, a nonvolatile memory card, a ROM, or the like, for example.

[0240] Further, the embodiments of the present invention explained above merely illustrate concrete examples of implementing the present invention, and the technical scope of the present invention is not to be construed in a restrictive manner by these embodiments. That is, the present invention may be implemented in various forms without departing from the technical features as defined in the appended claims.

INDUSTRIAL APPLICABILITY

[0241] The present invention can be utilized for inspecting tracks of railway vehicles.

Claims

1. An inspection system (400), comprising:

a data acquisition means (403, 2011a, 2011b) configured to acquire measured data being time-series data of measured values to be measured by causing a railway vehicle including a vehicle body (11), a bogie (12a, 12b), and a wheel set (13, 13a-13d) to travel on a track (16);

a first track state calculation means (407) configured to calculate an estimated value of a first physical quantity;

a correction amount calculation means (409) configured to calculate a correction amount for the estimated value of the first physical quantity based on the estimated value of the first physical quantity calculated by the first track state calculation means (407) and an actual value of the first physical quantity;

a second track state calculation means (411) configured to calculate an estimated value of the first physical quantity after the correction amount is calculated; and

a track state correction means (412) configured to correct the estimated value of the first physical quantity calculated by the second track state calculation means (411) by using the correction amount, wherein

the measured data contain a measured value of a forward-and-backward-direction force (T₁-T₄),

the forward-and-backward-direction force (T₁-T₄) is a force in a forward and backward direction (x) to occur in a member disposed between the wheel set (13, 13a-13d) and the bogie (12a, 12b) on which the wheel set (13, 13a-13d) is provided,

the member is a member for supporting an axle box,

the forward and backward direction (x) is a direction along a traveling direction of the railway vehicle,

the first physical quantity is a physical quantity reflecting a state of the track (16),

the first track state calculation means (407) and the second track state calculation means (411) are configured to use a relational expression representing the relation between the first physical quantity at a position of the wheel set (13, 13a-13d) and the forward-and-backward-direction force (T₁-T₄) and a measured value of the forward-and-backward-direction force (T₁-T₄) to calculate the estimated value of the first physical quantity,

the measured value of the forward-and-backward-direction force (T₁-T₄) used in the first track state calculation means (407) is contained in the measured data acquired by the data acquisition means (403, 2011a, 2011b) before the correction amount is calculated, and

the measured value of the forward-and-backward-direction force (T₁-T₄) used in the second track state calculation means (411) is contained in the measured data acquired by the data acquisition means (403, 2011a, 2011b) after the correction amount is calculated.

2. The inspection system (400) according to claim 1, wherein
the correction amount calculation means (409) is configured to calculate, based on estimated values of a plurality of the first physical quantities calculated by the first track state calculation means (407) using measured values of the forward-and-backward-direction force (T₁-T₄) when the railway vehicle is traveling at the same position and actual values of the first physical quantities, a correction amount at the position as the correction amount for the estimated value of the first physical quantity.

3. The inspection system (400) according to claim 1 or 2, wherein

the correction amount calculation means (409) is configured to calculate, based on estimated values of a plurality of the first physical quantities calculated by the first track state calculation means (407) using measured values of the forward-and-backward-direction force (T₁-T₄) when the railway vehicle is traveling at the same position, a traveling velocity of the railway vehicle traveling at the position, and actual values of the first physical quantities, a correction amount for the estimated value of the first physical quantity, and

the correction amount for the estimated value of the first physical quantity is a correction amount according to the position and the traveling velocity of the railway vehicle.

4. The inspection system (400) according to any one of claims 1 to 3, wherein
the measured values of the forward-and-backward-direction force (T₁-T₄) used by the first track state calculation means (407) and the second track state calculation means (411) are measured values in the same railway vehicle.

5. The inspection system (400) according to any one of claims 1 to 4, further comprising:

a frequency adjustment means (404, 406) configured to reduce a signal strength of a low-frequency component to be generated due to the railway vehicle traveling on a curved portion of the track (16) from time-series data of a second physical quantity, wherein

the second physical quantity is a physical quantity to vary in a value according to a state of the railway vehicle,

the frequency adjustment means (404, 406) includes a first frequency adjustment means (404) configured to reduce a signal strength of a low-frequency component to be generated due to the railway vehicle traveling on a curved portion of the track (16) from time-series data of a measured value of the forward-and-backward-direction force (T₁-T₄) being one of the second physical quantity,

the first track state calculation means (407) and the second track state calculation means (411) are configured to use the relational expression and the value of the forward-and-backward-direction force (T₁-T₄) from which the signal strength of the low-frequency component has been reduced by the first frequency adjustment means (404) to calculate the estimated value of the first physical quantity, and

the relational expression is an expression not including a curvature radius (R) of a rail (16a).

6. The inspection system (400) according to claim 5, wherein

the frequency adjustment means (404, 406) is configured to use the time-series data of the second physical quantity and determines a coefficient in a corrected AR model, and uses the corrected AR model whose coefficient is determined and the time-series data of the second physical quantity and reduces a signal strength of a low-frequency component to be generated due to the railway vehicle traveling on a curved portion of the track (16) from the time-series data of the second physical quantity,

the corrected AR model is an expression representing a predicted value of the second physical quantity by using a value of the second physical quantity and the coefficient responsive to the value,

the frequency adjustment means (404, 406) is configured to use an equation in which a first matrix is set to a coefficient matrix and an autocorrelation vector is set to a constant vector and determines the coefficient,

the autocorrelation vector is a vector whose component is autocorrelation of the time-series data of the second physical quantity with a time lag of 1 to m, m being a number of the measured value used in the corrected AR model,

the first matrix is a matrix U_sΣ_sU_s^T derived from a second matrix Σ_s that is derived from s pieces of eigenvalues of an autocorrelation matrix, s being a number set to be 1 or more and less than m, and a diagonal matrix Σ and a third matrix U_s that is derived from s pieces of the eigenvalues and an orthogonal matrix U,

the autocorrelation matrix is a matrix whose component is autocorrelation of the time-series data of the second physical quantity with a time lag of 0 to m - 1,

the diagonal matrix is a matrix whose diagonal component is eigenvalues of the autocorrelation matrix that are derived by singular value decomposition of the autocorrelation matrix,

the orthogonal matrix is a matrix in which an eigenvector of the autocorrelation matrix is set to a column component vector,

the second matrix is a submatrix of the diagonal matrix and is a matrix whose diagonal component is s pieces of the eigenvalues, and

the third matrix is a submatrix of the orthogonal matrix and is a matrix in which eigenvectors corresponding to s pieces of the eigenvalues are set to column component vectors.

7. The inspection system (400) according to claim 6, wherein
out of eigenvalues of the autocorrelation matrix, the largest eigenvalue is included in s pieces of the eigenvalues.

8. The inspection system (400) according to any one of claims 1 to 7, further comprising:

a filter operation means (405) configured to perform an operation using a filter performing data assimilation by using the measured data, a state equation, and an observation equation, and thereby determines estimated values of state variables being variables to determine estimated values in the state equation, wherein

the measured data further contain measured values of accelerations of the bogie (12a, 12b) and the wheel set (13, 13a-13d) in a right and left direction (y),

the right and left direction (y) is a direction perpendicular to both the forward and backward direction (x) and an up and down direction (z) being a direction perpendicular to the track (16),

the forward-and-backward-direction force (T₁-T₄) is a force to be determined according to a difference between an angular displacement of the wheel set (13, 13a-13d) in a yawing direction and an angular displacement of the bogie (12a, 12b) on which the wheel set (13, 13a-13d) is provided in the yawing direction,

the yawing direction is a pivoting direction with the up and down direction (z) set as a pivot axis,

the state equation is an equation described by using the state variables, the forward-and-backward-direction force (T₁-T₄), and a transformation variable,

the state variables include a displacement and a velocity of the bogie (12a, 12b) in the right and left direction (y), an angular displacement and an angular velocity of the bogie (12a, 12b) in the yawing direction, an angular displacement and an angular velocity of the bogie (12a, 12b) in a rolling direction, a displacement and a velocity of the wheel set (13, 13a-13d) in the right and left direction (y), and an angular displacement of an air spring attached to the railway vehicle in the rolling direction, and do not include an angular displacement or an angular velocity of the wheel set (13, 13a-13d) in the yawing direction,

the rolling direction is a pivoting direction with the forward and backward direction (x) set as a pivot axis,

the transformation variable is a variable that performs mutual transformation between the angular displacement of the wheel set (13, 13a-13d) in the yawing direction and the angular displacement of the bogie (12a, 12b) in the yawing direction,

the observation equation is an equation described by using an observation variable and the transformation variable,

the observation variable includes the accelerations of the bogie (12a, 12b) and the wheel set (13, 13a-13d) in the right and left direction (y),

the filter operation means (405) is configured to use a measured value of the observation variable, the state equation into which the measured value of the forward-and-backward-direction force (T₁-T₄) and an actual value of the transformation variable are substituted, and the observation equation into which the actual value of the transformation variable is substituted, and determines estimated values of the state variables when an error between, of the observation variable, the measured value and an estimated value or an expected value of the error becomes minimum,

the first track state calculation means (407) and the second track state calculation means (411) are configured to use an estimated value of the angular displacement of the bogie (12a, 12b) in the yawing direction, which is one of the state variables determined by the filter operation means (405), and the actual value of the transformation variable to calculate an estimated value of the angular displacement of the wheel set (13, 13a-13d) in the yawing direction, and use the estimated value of the angular displacement of the wheel set (13, 13a-13d) in the yawing direction, the value of the forward-and-backward-direction force (T₁-T₄), and the relational expression to calculate the estimated value of the first physical quantity,

the relational expression is an expression in which a motion equation describing motion of the wheel set (13, 13a-13d) in the yawing direction is expressed by using the forward-and-backward-direction force (T₁-T₄), and

the actual value of the transformation variable is derived by using the measured value of the forward-and-backward-direction force (T₁-T₄).

9. The inspection system (400) according to claim 8, wherein

the state equation is constituted by using a motion equation that describes motion of the wheel set (13, 13a-13d) in the right and left direction (y), a motion equation that describes motion of the bogie (12a, 12b) in the right and left direction (y), a motion equation that describes motion of the bogie (12a, 12b) in the yawing direction, a motion equation that describes motion of the bogie (12a, 12b) in the rolling direction, and a motion equation that describes motion of the air spring in the rolling direction,

the motion equation that describes the motion of the wheel set (13, 13a-13d) in the right and left direction (y) is a motion equation described by using the transformation variable in place of using the angular displacement of the wheel set (13, 13a-13d) in the yawing direction,

the motion equation that describes the motion of the bogie (12a, 12b) in the yawing direction is a motion equation described by using the forward-and-backward-direction force (T₁-T₄) in place of using the angular displacement and the angular velocity of the wheel set (13, 13a-13d) in the yawing direction, and

the transformation variable is represented by a difference between the angular displacement of the bogie (12a, 12b) in the yawing direction and the angular displacement of the wheel set (13, 13a-13d) in the yawing direction.

10. The inspection system (400) according to claim 8 or 9, wherein

the data acquisition means (403, 2011a, 2011b) is configured to further acquire a measured value of acceleration of the vehicle body (11) in the right and left direction (y),

the observation variable further includes the acceleration of the vehicle body (11) in the right and left direction (y),

the state variables further include a displacement and a velocity of the vehicle body (11) in the right and left direction (y), an angular displacement and an angular velocity of the vehicle body (11) in the yawing direction, an angular displacement and an angular velocity of the vehicle body (11) in the rolling direction, and an angular displacement of a yaw damper attached to the railway vehicle in the yawing direction, and

the filter operation means (405) is configured to determine the state variables when differences between, of the accelerations of the vehicle body (11), the bogie (12a, 12b), and the wheel set (13, 13a-13d) in the right and left direction (y), measured values and calculated values become minimum.

11. The inspection system (400) according to claim 10, wherein
the state equation is constituted by further using a motion equation that describes motion of the vehicle body (11) in the right and left direction (y), a motion equation that describes motion of the vehicle body (11) in the yawing direction, a motion equation that describes motion of the vehicle body (11) in the rolling direction, and a motion equation that describes motion of the yaw damper in the yawing direction.

12. The inspection system (400) according to any one of claims 8 to 11, wherein

the observation equation is constituted by further using the motion equation that describes the motion of the wheel set (13, 13a-13d) in the right and left direction (y) and the motion equation that describes the motion of the bogie (12a, 12b) in the right and left direction (y), and

the motion equation that describes the motion of the wheel set (13, 13a-13d) in the right and left direction (y) is a motion equation described by using the transformation variable in place of using the angular displacement of the wheel set (13, 13a-13d) in the yawing direction.

13. The inspection system (400) according to claim 12, wherein
the observation equation is constituted by further using the motion equation that describes the motion of the vehicle body (11) in the right and left direction (y).

14. The inspection system (400) according to any one of claims 8 to 13, wherein

the first track state calculation means (407) and the second track state calculation means (411) are configured to derive an alignment irregularity amount (y_R1-y_R4) of the track (16) as the estimated value of the first physical quantity based on the displacement and the velocity of the bogie (12a, 12b) in the right and left direction (y) being the state variables determined by the filter operation means (405), the displacement and the velocity of the wheel set (13, 13a-13d) in the right and left direction (y) being the state variables determined by the filter operation means (405), the estimated value of the angular displacement of the wheel set (13, 13a-13d) in the yawing direction, the measured value of the forward-and-backward-direction force (T₁-T₄), and the motion equation that describes the motion of the wheel set (13, 13a-13d) in the yawing direction, and

the motion equation that describes the motion of the wheel set (13, 13a-13d) in the yawing direction includes the forward-and-backward-direction force (T₁-T₄) and the alignment irregularity amount (y_R1-y_R4) of the track (16) as variables.

15. The inspection system (400) according to any one of claims 8 to 14, further comprising:

a frequency adjustment means (404, 406) configured to reduce a signal strength of a low-frequency component to be generated due to the railway vehicle traveling on a curved portion of the track (16) from time-series data of a second physical quantity, wherein

the second physical quantity is a physical quantity to vary in a value according to a state of the railway vehicle, and

the frequency adjustment means includes a second frequency adjustment means (406) configured to reduce a signal strength of a low-frequency component to be generated due to the railway vehicle traveling on a curved portion of the track (16) from time-series data of the estimated values of the state variables being one of the second physical quantity.

16. The inspection system (400) according to any one of claims 8 to 15, wherein
the filter is a Kalman filter.

17. The inspection system (400) according to any one of claims 1 to 16, wherein

the first physical quantity is an alignment irregularity amount (y_R1-y_R4) of the track (16) or a lateral force (Q₁-Q₄) that is a stress in the right and left direction (y) between a wheel (14L, 14R, 14a-14d) provided on the wheel set (13, 13a-13d) and the track (16), and

the right and left direction (y) is a direction perpendicular to both the forward and backward direction (x) and the up and down direction (z) being a direction perpendicular to the track (16).

18. The inspection system (400) according to any one of claims 1 to 17, wherein

the forward-and-backward-direction force (T₁-T₄) is, out of components in the forward and backward direction (x) of forces to occur in the two members attached to both sides of the single wheel set (13, 13a-13d) in the right and left direction (y), the component opposite in phase to each other, and

the right and left direction (y) is a direction perpendicular to both the forward and backward direction (x) and the up and down direction (z) being a direction perpendicular to the track (16).

19. An inspection method carried out by an inspection system, the inspection method comprising:

a data acquisition step (S701) of acquiring measured data being time-series data of measured values to be measured by causing a railway vehicle including a vehicle body (11), a bogie (12a, 12b), and a wheel set (13, 13a-13d) to travel on a track (16);

a first track state calculation step (S705) of calculating an estimated value of a first physical quantity;

a correction amount calculation step (S707) of calculating a correction amount for the estimated value of the first physical quantity based on the estimated value of the first physical quantity calculated by the first track state calculation step (S705) and an actual value of the first physical quantity;

a second track state calculation step of calculating an estimated value of the first physical quantity after the correction amount is calculated (S707); and

a track state correction step of correcting the estimated value of the first physical quantity calculated by the second track state calculation step by using the correction amount, wherein

the measured data contain a measured value of a forward-and-backward-direction force (T₁-T₄),

the forward-and-backward-direction force (T₁-T₄) is a force in a forward and backward direction (x) to occur in a member disposed between the wheel set (13, 13a-13d) and the bogie (12a, 12b) on which the wheel set (13, 13a-13d) is provided,

the member is a member for supporting an axle box,

the forward and backward direction (x) is a direction along a traveling direction of the railway vehicle,

the first physical quantity is a physical quantity reflecting a state of the track (16),

the first track state calculation step (S705) and the second track state calculation step use a relational expression representing the relation between the first physical quantity at a position of the wheel set (13, 13a-13d) and the forward-and-backward-direction force (T₁-T₄) and a measured value of the forward-and-backward-direction force (T₁-T₄) to calculate the estimated value of the first physical quantity,

the measured value of the forward-and-backward-direction force (T₁-T₄) used in the first track state calculation step (S705) is contained in the measured data acquired by the data acquisition step (S701) before the correction amount is calculated (S707), and

the measured value of the forward-and-backward-direction force (T₁-T₄) used in the second track state calculation step is contained in the measured data acquired by the data acquisition step after the correction amount is calculated (S707).

20. A computer program comprising instructions which, when the program is executed by a computer, cause the computer to execute steps comprising:

a data acquisition step (S701) of acquiring measured data being time-series data of measured values to be measured by causing a railway vehicle including a vehicle body (11), a bogie (12a, 12b), and a wheel set (13, 13a-13d) to travel on a track (16);

a first track state calculation step (S705) of calculating an estimated value of a first physical quantity;

a correction amount calculation step (S707) of calculating a correction amount for the estimated value of the first physical quantity based on the estimated value of the first physical quantity calculated by the first track state calculation step (S705) and an actual value of the first physical quantity;

a second track state calculation step of calculating an estimated value of the first physical quantity after the correction amount is calculated (S707); and

a track state correction step of correcting the estimated value of the first physical quantity calculated by the second track state calculation step by using the correction amount, wherein

the measured data contain a measured value of a forward-and-backward-direction force (T₁-T₄),

the forward-and-backward-direction force (T₁-T₄) is a force in a forward and backward direction (x) to occur in a member disposed between the wheel set (13, 13a-13d) and the bogie (12a, 12b) on which the wheel set (13, 13a-13d) is provided,

the member is a member for supporting an axle box,

the forward and backward direction (x) is a direction along a traveling direction of the railway vehicle,

the first physical quantity is a physical quantity reflecting a state of the track (16),

the first track state calculation step (S705) and the second track state calculation step use a relational expression representing the relation between the first physical quantity at a position of the wheel set (13, 13a-13d) and the forward-and-backward-direction force (T₁-T₄) and a measured value of the forward-and-backward-direction force (T₁-T₄) to calculate the estimated value of the first physical quantity,

the measured value of the forward-and-backward-direction force (T₁-T₄) used in the first track state calculation step (S705) is contained in the measured data acquired by the data acquisition step (S701) before the correction amount is calculated (S707), and

the measured value of the forward-and-backward-direction force (T₁-T₄) used in the second track state calculation step is contained in the measured data acquired by the data acquisition step after the correction amount is calculated (S707).

Ansprüche

1. Inspektionssystem (400), umfassend:

ein Datenerfassungsmittel (403, 2011a, 2011b), das so konfiguriert ist, dass es Messdaten erfasst, bei denen es sich um Zeitreihendaten von Messwerten handelt, die zu messen sind, indem ein Schienenfahrzeug, das einen Fahrzeugkörper (11), ein Drehgestell (12a, 12b) und einen Radsatz (13, 13a-13d) aufweist, veranlasst wird, auf einem Gleis (16) zu fahren;

ein erstes Gleiszustandsberechnungsmittel (407), das konfiguriert ist, um einen geschätzten Wert einer ersten physikalischen Größe zu berechnen;

ein Korrekturbetragsberechnungsmittel (409), das so konfiguriert ist, dass es einen Korrekturbetrag für den geschätzten Wert der ersten physikalischen Größe basierend auf dem geschätzten Wert der ersten physikalischen Größe, der von dem ersten Gleiszustandsberechnungsmittel (407) berechnet wurde, und einem tatsächlichen Wert der ersten physikalischen Größe berechnet;

eine zweites Gleiszustandsberechnungsmittel (411), das konfiguriert ist, um einen geschätzten Wert der ersten physikalischen Größe zu berechnen, nachdem der Korrekturbetrag berechnet wurde; und

ein Gleiszustandskorrekturmittel (412), das so konfiguriert ist, dass es den geschätzten Wert der ersten physikalischen Größe, der von dem zweiten Gleiszustandsberechnungsmittel (411) berechnet wurde, unter Verwendung des Korrekturbetrags korrigiert, wobei

die gemessenen Daten einen gemessenen Wert einer vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) enthalten,

die vorwärts- und rückwärtsgerichtete Kraft (T₁-T₄) eine Kraft in einer Vorwärts- und Rückwärtsrichtung (x) ist, die in einem Element auftritt, das zwischen dem Radsatz (13, 13a-13d) und dem Drehgestell (12a, 12b) angeordnet ist, auf dem der Radsatz (13, 13a-13d) vorgesehen ist,

das Element ein Element zum Tragen eines Achslagers ist,

die Vorwärts- und Rückwärtsrichtung (x) eine Richtung entlang einer Fahrtrichtung des Schienenfahrzeugs ist,

die erste physikalische Größe eine physikalische Größe ist, die einen Zustand des Gleises (16) wiedergibt,

das erste Gleiszustandsberechnungsmittel (407) und das zweite Gleiszustandsberechnungsmittel (411) so konfiguriert sind, dass sie einen relationalen Ausdruck, der die Beziehung zwischen der ersten physikalischen Größe an einer Position des Radsatzes (13, 13a-13d) und der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) darstellt, und einen gemessenen Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) verwenden, um den geschätzten Wert der ersten physikalischen Größe zu berechnen,

der gemessene Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄), der in dem ersten Gleiszustandsberechnungsmittel (407) verwendet wird, in den gemessenen Daten enthalten ist, die durch das Datenerfassungsmittel (403, 2011a, 2011b) erfasst werden, bevor der Korrekturbetrag berechnet wird, und

der gemessene Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄), der in dem zweiten Gleiszustandsberechnungsmittel (411) verwendet wird, in den gemessenen Daten enthalten ist, die durch das Datenerfassungsmittel (403, 2011a, 2011b) erfasst werden, nachdem der Korrekturbetrag berechnet wurde.

2. Inspektionssystem (400) nach Anspruch 1, wobei
das Korrekturbetragsberechnungsmittel (409) konfiguriert ist, um, basierend auf geschätzten Werten einer Vielzahl der ersten physikalischen Größen, die durch das erste Gleiszustandsberechnungsmittel (407) unter Verwendung gemessener Werte der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) berechnet werden, wenn das Schienenfahrzeug an der gleichen Position fährt, und tatsächlichen Werten der ersten physikalischen Größen, einen Korrekturbetrag an der Position als den Korrekturbetrag für den geschätzten Wert der ersten physikalischen Größe zu berechnen.

3. Inspektionssystem (400) nach Anspruch 1 oder 2, wobei

das Korrekturbetragsberechnungsmittel (409) konfiguriert ist, um, basierend auf geschätzten Werten einer Vielzahl der ersten physikalischen Größen, die durch das erste Gleiszustandsberechnungsmittel (407) unter Verwendung gemessener Werte der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄), wenn das Schienenfahrzeug an derselben Position fährt, einer Fahrgeschwindigkeit des an der Position fahrenden Schienenfahrzeugs und tatsächlichen Werten der ersten physikalischen Größen, einen Korrekturbetrag für den geschätzten Wert der ersten physikalischen Größe zu berechnen, und

der Korrekturbetrag für den geschätzten Wert der ersten physikalischen Größe ein Korrekturbetrag gemäß der Position und der Fahrgeschwindigkeit des Schienenfahrzeugs ist.

4. Inspektionssystem (400) nach einem der Ansprüche 1 bis 3, wobei
die von dem ersten Gleiszustandsberechnungsmittel (407) und dem zweiten Gleiszustandsberechnungsmittel (411) verwendeten Messwerte der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) Messwerte in demselben Schienenfahrzeug sind.

5. Inspektionssystem (400) nach einem der Ansprüche 1 bis 4, ferner umfassend:

ein Frequenzanpassungsmittel (404, 406), das so konfiguriert ist, dass es eine Signalstärke einer niederfrequenten Komponente, die aufgrund des auf einem gekrümmten Abschnitt des Gleises (16) fahrenden Schienenfahrzeugs zu erzeugen ist, aus Zeitreihendaten einer zweiten physikalischen Größe reduziert, wobei

die zweite physikalische Größe eine physikalische Größe ist, deren Wert gemäß einem Zustand des Schienenfahrzeugs variiert,

das Frequenzeinstellmittel (404, 406) ein erstes Frequenzeinstellmittel (404) aufweist, das so konfiguriert ist, dass es eine Signalstärke einer Niederfrequenzkomponente, die aufgrund des auf einem gekrümmten Abschnitt des Gleises (16) fahrenden Schienenfahrzeugs zu erzeugen ist, aus Zeitreihendaten eines gemessenen Werts der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄), die eine der zweiten physikalischen Größe ist, reduziert,

das erste Gleiszustandsberechnungsmittel (407) und das zweite Gleiszustandsberechnungsmittel (411) so konfiguriert sind, dass sie den relationalen Ausdruck und den Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) verwenden, von dem die Signalstärke der Niederfrequenzkomponente durch das erste Frequenzeinstellmittel (404) reduziert wurde, um den geschätzten Wert der ersten physikalischen Größe zu berechnen, und

der relationale Ausdruck ein Ausdruck ist, der einen Krümmungsradius (R) einer Schiene (16a) nicht aufweist.

6. Inspektionssystem (400) nach Anspruch 5, wobei

das Frequenzanpassungsmittel (404, 406) konfiguriert ist, um die Zeitreihendaten der zweiten physikalischen Größe zu verwenden und einen Koeffizienten in einem korrigierten AR-Modell zu bestimmen, und das korrigierte AR-Modell, dessen Koeffizient bestimmt ist, und die Zeitreihendaten der zweiten physikalischen Größe verwendet und eine Signalstärke einer Niederfrequenzkomponente, die aufgrund des auf einem gekrümmten Abschnitt des Gleises (16) fahrenden Schienenfahrzeugs zu erzeugen ist, aus den Zeitreihendaten der zweiten physikalischen Größe reduziert,

das korrigierte AR-Modell ein Ausdruck ist, der einen vorhergesagten Wert der zweiten physikalischen Größe darstellt, indem ein Wert der zweiten physikalischen Größe und der Koeffizient, der auf den Wert reagiert, verwendet werden,

das Frequenzanpassungsmittel (404, 406) so konfiguriert ist, dass sie eine Gleichung verwendet, in der eine erste Matrix auf eine Koeffizientenmatrix und ein Autokorrelationsvektor auf einen konstanten Vektor festgelegt ist und den Koeffizienten bestimmt,

der Autokorrelationsvektor ein Vektor ist, dessen Komponente die Autokorrelation der Zeitreihendaten der zweiten physikalischen Größe mit einer Zeitverzögerung von 1 bis m ist, wobei m eine Zahl des in dem korrigierten AR-Modell verwendeten Messwertes ist,

die erste Matrix eine Matrix U_sΣ_sU_s^T ist, die aus einer zweiten Matrix Σ_s, die aus s Teilen von Eigenwerten einer Autokorrelationsmatrix abgeleitet ist, wobei s eine Zahl ist, die auf 1 oder mehr und weniger als m festgelegt ist, und einer Diagonalmatrix Σ und einer dritten Matrix U_s, die aus s Teilen der Eigenwerte und einer orthogonalen Matrix U abgeleitet ist,

die Autokorrelationsmatrix eine Matrix ist, deren Komponente die Autokorrelation der Zeitreihendaten der zweiten physikalischen Größe mit einer Zeitverzögerung von 0 bis m - 1 ist,

die Diagonalmatrix eine Matrix ist, deren Diagonalkomponente Eigenwerte der Autokorrelationsmatrix sind, die durch Singulärwertzerlegung der Autokorrelationsmatrix abgeleitet werden,

die orthogonale Matrix ist eine Matrix, in der ein Eigenvektor der Autokorrelationsmatrix auf einen Spaltenkomponentenvektor gesetzt ist,

die zweite Matrix ist eine Untermatrix der Diagonalmatrix und ist eine Matrix, deren Diagonalkomponente s Stücke der Eigenwerte sind, und

die dritte Matrix ist eine Untermatrix der orthogonalen Matrix und ist eine Matrix, in der Eigenvektoren, die s Teilen der Eigenwerte entsprechen, auf Spaltenkomponentenvektoren gesetzt sind.

7. Inspektionssystem (400) nach Anspruch 6, wobei
von den Eigenwerten der Autokorrelationsmatrix der größte Eigenwert in s Stücken der Eigenwerte enthalten ist.

8. Inspektionssystem (400) nach einem der Ansprüche 1 bis 7, ferner umfassend:

ein Filteroperationsmittel (405), das so konfiguriert ist, dass es eine Operation unter Verwendung eines Filters durchführt, das eine Datenassimilation unter Verwendung der gemessenen Daten, einer Zustandsgleichung und einer Beobachtungsgleichung durchführt, und dadurch geschätzte Werte von Zustandsvariablen bestimmt, die Variablen zum Bestimmen geschätzter Werte in der Zustandsgleichung sind, wobei

die gemessenen Daten ferner gemessene Werte von Beschleunigungen des Drehgestells (12a, 12b) und des Radsatzes (13, 13a-13d) in einer rechten und linken Richtung (y) enthalten,

die rechte und linke Richtung (y) eine Richtung senkrecht zur Vorwärts- und Rückwärtsrichtung (x) ist und eine Aufwärts- und Abwärtsrichtung (z), die eine Richtung senkrecht zum Gleis (16) ist,

die vorwärts- und rückwärtsgerichtete Kraft (T₁-T₄) eine Kraft ist, die gemäß einer Differenz zwischen einer Winkelverschiebung des Radsatzes (13, 13a-13d) in einer Gierrichtung und einer Winkelverschiebung des Drehgestells (12a, 12b), auf dem der Radsatz (13, 13a-13d) in der Gierrichtung vorgesehen ist, zu bestimmen ist,

die Gierrichtung eine Schwenkrichtung ist, wobei die Auf- und Abwärtsrichtung (z) als Schwenkachse festgelegt ist,

die Zustandsgleichung eine Gleichung ist, die durch Verwendung der Zustandsvariablen, der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) und einer Transformationsvariablen beschrieben wird,

die Zustandsvariablen eine Verschiebung und eine Geschwindigkeit des Drehgestells (12a, 12b) in der rechten und linken Richtung (y), eine Winkelverschiebung und eine Winkelgeschwindigkeit des Drehgestells (12a, 12b) in der Gierrichtung, eine Winkelverschiebung und eine Winkelgeschwindigkeit des Drehgestells (12a, 12b) in einer Rollrichtung, eine Verschiebung und eine Geschwindigkeit des Radsatzes (13, 13a-13d) in der rechten und linken Richtung (y) und eine Winkelverschiebung einer an dem Schienenfahrzeug angebrachten Luftfeder in der Rollrichtung umfassen, und nicht eine Winkelverschiebung oder eine Winkelgeschwindigkeit des Radsatzes (13, 13a-13d) in der Gierrichtung umfassen,

die Rollrichtung eine Schwenkrichtung ist, wobei die Vorwärts- und Rückwärtsrichtung (x) als Schwenkachse festgelegt ist,

die Transformationsvariable eine Variable ist, die eine gegenseitige Transformation zwischen der Winkelverschiebung des Radsatzes (13, 13a-13d) in der Gierrichtung und der Winkelverschiebung des Drehgestells (12a, 12b) in der Gierrichtung durchführt,

die Beobachtungsgleichung eine Gleichung ist, die durch Verwendung einer Beobachtungsvariablen und der Transformationsvariablen beschrieben wird,

die Beobachtungsvariable die Beschleunigungen des Drehgestells (12a, 12b) und des Radsatzes (13, 13a-13d) in der rechten und linken Richtung (y) umfasst,

das Filteroperationsmittel (405) so konfiguriert ist, dass es einen gemessenen Wert der Beobachtungsvariablen, die Zustandsgleichung, in die der gemessene Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) und ein tatsächlicher Wert der Transformationsvariablen eingesetzt werden, und die Beobachtungsgleichung, in die der tatsächliche Wert der Transformationsvariablen eingesetzt wird, verwendet und geschätzte Werte der Zustandsvariablen bestimmt, wenn ein Fehler zwischen der Beobachtungsvariablen, dem gemessenen Wert und einem geschätzten Wert oder einem erwarteten Wert des Fehlers minimal wird,

das erste Gleiszustandsberechnungsmittel (407) und das zweite Gleiszustandsberechnungsmittel (411) so konfiguriert sind, dass sie einen geschätzten Wert der Winkelverschiebung des Drehgestells (12a, 12b) in der Gierrichtung, die eine der von dem Filteroperationsmittel (405) bestimmten Zustandsvariablen ist, und den aktuellen Wert der Transformationsvariablen verwenden, um einen geschätzten Wert der Winkelverschiebung des Radsatzes (13, 13a-13d) in der Gierrichtung zu berechnen, und den geschätzten Wert der Winkelverschiebung des Radsatzes (13, 13a-13d) in der Gierrichtung, den Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) und den relationalen Ausdruck zu verwenden, um den geschätzten Wert der ersten physikalischen Größe zu berechnen,

der relationale Ausdruck ein Ausdruck ist, in dem eine Bewegungsgleichung, die die Bewegung des Radsatzes (13, 13a-13d) in der Gierrichtung beschreibt, unter Verwendung der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) ausgedrückt wird, und

der tatsächliche Wert der Transformationsvariablen unter Verwendung des gemessenen Wertes der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) abgeleitet wird.

9. Inspektionssystem (400) nach Anspruch 8, wobei

die Zustandsgleichung durch Verwendung einer Bewegungsgleichung, die die Bewegung des Radsatzes (13, 13a-13d) in der rechten und linken Richtung (y) beschreibt, einer Bewegungsgleichung, die die Bewegung des Drehgestells (12a, 12b) in der rechten und linken Richtung (y) beschreibt, einer Bewegungsgleichung, die die Bewegung des Drehgestells (12a, 12b) in der Gierrichtung beschreibt, einer Bewegungsgleichung, die die Bewegung des Drehgestells (12a, 12b) in der Rollrichtung beschreibt, und einer Bewegungsgleichung, die die Bewegung der Luftfeder in der Rollrichtung beschreibt, gebildet wird,

die Bewegungsgleichung, die die Bewegung des Radsatzes (13, 13a-13d) in der rechten und linken Richtung (y) beschreibt, eine Bewegungsgleichung ist, die durch Verwenden der Transformationsvariablen anstelle der Verwendung der Winkelverschiebung des Radsatzes (13, 13a-13d) in der Gierrichtung beschrieben wird,

die Bewegungsgleichung, die die Bewegung des Drehgestells (12a, 12b) in der Gierrichtung beschreibt, eine Bewegungsgleichung ist, die beschrieben wird, indem die vorwärts- und rückwärtsgerichtete Kraft (T₁-T₄) anstelle der Winkelverschiebung und der Winkelgeschwindigkeit des Radsatzes (13, 13a-13d) in der Gierrichtung verwendet wird, und

die Transformationsvariable durch eine Differenz zwischen der Winkelverschiebung des Drehgestells (12a, 12b) in der Gierrichtung und der Winkelverschiebung des Radsatzes (13, 13a-13d) in der Gierrichtung dargestellt wird.

10. Inspektionssystem (400) nach Anspruch 8 oder 9, wobei

das Datenerfassungsmittel (403, 2011a, 2011b) so konfiguriert ist, dass es ferner einen Messwert der Beschleunigung des Fahrzeugkörper (11) in der rechten und linken Richtung (y) erfasst,

die Beobachtungsvariable ferner die Beschleunigung des Fahrzeugkörpers (11) in der rechten und linken Richtung (y) umfasst,

die Zustandsvariablen ferner eine Verschiebung und eine Geschwindigkeit des Fahrzeugkörpers (11) in der rechten und linken Richtung (y), eine Winkelverschiebung und eine Winkelgeschwindigkeit des Fahrzeugkörpers (11) in der Gierrichtung, eine Winkelverschiebung und eine Winkelgeschwindigkeit des Fahrzeugkörpers (11) in der Rollrichtung und eine Winkelverschiebung eines an dem Schienenfahrzeug angebrachten Gierdämpfers in der Gierrichtung umfassen, und

das Filteroperationsmittel (405) konfiguriert ist, um die Zustandsvariablen zu bestimmen, wenn Differenzen zwischen den Beschleunigungen des Fahrzeugkörpers (11), des Drehgestells (12a, 12b) und des Radsatzes (13, 13a-13d) in der rechten und linken Richtung (y), den gemessenen Werten und den berechneten Werten minimal werden.

11. Inspektionssystem (400) nach Anspruch 10, wobei
die Zustandsgleichung durch weitere Verwendung einer Bewegungsgleichung, die die Bewegung des Fahrzeugkörpers (11) in der rechten und linken Richtung (y) beschreibt, einer Bewegungsgleichung, die die Bewegung des Fahrzeugkörpers (11) in der Gierrichtung beschreibt, einer Bewegungsgleichung, die die Bewegung des Fahrzeugkörpers (11) in der Rollrichtung beschreibt, und einer Bewegungsgleichung, die die Bewegung des Gierdämpfers in der Gierrichtung beschreibt, gebildet wird.

12. Inspektionssystem (400) nach einem der Ansprüche 8 bis 11, wobei

die Beobachtungsgleichung durch weitere Verwendung der Bewegungsgleichung, die die Bewegung des Radsatzes (13, 13a-13d) in der rechten und linken Richtung (y) beschreibt, und der Bewegungsgleichung, die die Bewegung des Drehgestells (12a, 12b) in der rechten und linken Richtung (y) beschreibt, gebildet wird, und

die Bewegungsgleichung, die die Bewegung des Radsatzes (13, 13a-13d) in der rechten und linken Richtung (y) beschreibt, eine Bewegungsgleichung ist, die durch Verwendung der Transformationsvariablen anstelle der Verwendung der Winkelverschiebung des Radsatzes (13, 13a-13d) in der Gierrichtung beschrieben wird.

13. Inspektionssystem (400) nach Anspruch 12, wobei
die Beobachtungsgleichung durch weitere Verwendung der Bewegungsgleichung gebildet wird, die die Bewegung des Fahrzeugkörpers (11) in der rechten und linken Richtung (y) beschreibt.

14. Inspektionssystem (400) nach einem der Ansprüche 8 bis 13, wobei

das erste Gleiszustandsberechnungsmittel (407) und das zweite Gleiszustandsberechnungsmittel (411) konfiguriert sind, um einen Ausrichtungsunregelmäßigkeitsbetrag (y_R1-y_R4) des Gleises (16) als den geschätzten Wert der ersten physikalischen Größe basierend auf der Verschiebung und der Geschwindigkeit des Drehgestells (12a, 12b) in der rechten und linken Richtung (y) abzuleiten, die die Zustandsvariablen sind, die durch das Filteroperationsmittel (405) bestimmt werden, der Verschiebung und der Geschwindigkeit des Radsatzes (13, 13a-13d) in der rechten und linken Richtung (y), die die von dem Filteroperationsmittel (405) bestimmten Zustandsvariablen sind, der geschätzte Wert der Winkelverschiebung des Radsatzes (13, 13a-13d) in der Gierrichtung, der gemessene Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄), und die Bewegungsgleichung, die die Bewegung des Radsatzes (13, 13a-13d) in der Gierrichtung beschreibt, und

die Bewegungsgleichung, die die Bewegung des Radsatzes (13, 13a-13d) in der Gierrichtung beschreibt, die vorwärts- und rückwärtsgerichtete Kraft (T₁-T₄) und den Ausrichtungsunregelmäßigkeitsbetrag (y_R1-y_R4) des Gleises (16) als Variablen aufweist.

15. Inspektionssystem (400) nach einem der Ansprüche 8 bis 14, ferner umfassend:

ein Frequenzanpassungsmittel (404, 406), das so konfiguriert ist, dass es eine Signalstärke einer niederfrequenten Komponente, die aufgrund des auf einem gekrümmten Abschnitt des Gleises (16) fahrenden Schienenfahrzeugs zu erzeugen ist, aus Zeitreihendaten einer zweiten physikalischen Größe reduziert, wobei

die zweite physikalische Größe eine physikalische Größe ist, deren Wert gemäß einem Zustand des Schienenfahrzeugs variiert, und

das Frequenzeinstellmittel ein zweites Frequenzeinstellmittel (406) umfasst, das konfiguriert ist, um eine Signalstärke einer Niederfrequenzkomponente, die aufgrund des auf einem gekrümmten Abschnitt des Gleises (16) fahrenden Schienenfahrzeugs zu erzeugen ist, aus Zeitreihendaten der geschätzten Werte der Zustandsvariablen, die eine der zweiten physikalischen Größe sind, zu reduzieren.

16. Inspektionssystem (400) nach einem der Ansprüche 8 bis 15, wobei
der Filter ein Kalman-Filter ist.

17. Inspektionssystem (400) nach einem der Ansprüche 1 bis 16, wobei

die erste physikalische Größe ein Ausrichtungsunregelmäßigkeitsbetrag (y_R1-y_R4) des Gleises (16) oder eine Seitenkraft (Q₁-Q₄) ist, die eine Spannung in der rechten und linken Richtung (y) zwischen einem Rad (14L, 14R, 14a-14d), das an dem Radsatz (13, 13a-13d) vorgesehen ist, und des Gleises (16) ist, und

die rechte und linke Richtung (y) eine Richtung senkrecht sowohl zur Vorwärts- und Rückwärtsrichtung (x) als auch zur Auf- und Abwärtsrichtung (z) ist, die eine Richtung senkrecht zum Gleis (16) ist.

18. Inspektionssystem (400) nach einem der Ansprüche 1 bis 17, wobei

die vorwärts- und rückwärtsgerichtete Kraft (T₁-T₄) aus den Komponenten der in Vorwärts- und Rückwärtsrichtung (x) wirkenden Kräfte, die in den beiden an beiden Seiten des einzelnen Radsatzes (13, 13a-13d) befestigten Elementen in der rechten und linken Richtung (y) auftreten, die Komponente ist, die in der Phase zueinander entgegengesetzt ist, und

die rechte und linke Richtung (y) eine Richtung senkrecht sowohl zur Vorwärts- als auch zur Rückwärtsrichtung (x) ist und die Aufwärts- und Abwärtsrichtung (z) eine Richtung senkrecht zum Gleis (16) ist.

19. Inspektionsverfahren, das von einem Inspektionssystem ausgeführt wird, wobei das Inspektionsverfahren umfasst:

einen Datenerfassungsschritt (S701) zum Erfassen von Messdaten, bei denen es sich um Zeitreihendaten von Messwerten handelt, die zu messen sind, indem ein Schienenfahrzeug, das einen Fahrzeugkörper (11), ein Drehgestell (12a, 12b) und einen Radsatz (13, 13a-13d) aufweist, veranlasst wird, auf einem Gleis (16) zu fahren;

einen ersten Gleiszustandsberechnungsschritt (S705) des Berechnens eines geschätzten Wertes einer ersten physikalischen Größe;

einen Korrekturbetragsberechnungsschritt (S707) des Berechnens eines Korrekturbetrags für den geschätzten Wert der ersten physikalischen Größe basierend auf dem geschätzten Wert der ersten physikalischen Größe, der durch den ersten Gleiszustandsberechnungsschritt (S705) berechnet wurde, und einem tatsächlichen Wert der ersten physikalischen Größe;

einen zweiten Gleiszustandsberechnungsschritt des Berechnens eines geschätzten Wertes der ersten physikalischen Größe, nachdem der Korrekturbetrag berechnet ist (S707); und

einen Gleiszustandskorrekturschritt zum Korrigieren des geschätzten Wertes der ersten physikalischen Größe, der durch den zweiten Gleiszustandsberechnungsschritt berechnet wurde, unter Verwendung des Korrekturbetrags, wobei

die gemessenen Daten einen gemessenen Wert einer vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) enthalten,

die vorwärts- und rückwärtsgerichtete Kraft (T₁-T₄) eine Kraft in einer Vorwärts- und Rückwärtsrichtung (x) ist, die in einem Element auftritt, das zwischen dem Radsatz (13, 13a-13d) und dem Drehgestell (12a, 12b) angeordnet ist, auf dem der Radsatz (13, 13a-13d) vorgesehen ist,

das Element ein Element zum Tragen eines Achslagers ist,

die Vorwärts- und Rückwärtsrichtung (x) eine Richtung entlang einer Fahrtrichtung des Schienenfahrzeugs ist,

die erste physikalische Größe eine physikalische Größe ist, die einen Zustand des Gleises (16) wiedergibt,

der erste Gleiszustandsberechnungsschritt (S705) und der zweite Gleiszustandsberechnungsschritt einen relationalen Ausdruck, der die Beziehung zwischen der ersten physikalischen Größe an einer Position des Radsatzes (13, 13a-13d) und der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) darstellt, und einen gemessenen Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) verwenden, um den geschätzten Wert der ersten physikalischen Größe zu berechnen,

der gemessene Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄), der in dem ersten Gleiszustandsberechnungsschritt (S705) verwendet wird, in den gemessenen Daten enthalten ist, die durch den Datenerfassungsschritt (S701) erfasst werden, bevor der Korrekturbetrag berechnet wird (S707), und

der gemessene Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄), der in dem zweiten Gleiszustandsberechnungsschritt verwendet wird, in den gemessenen Daten enthalten ist, die durch den Datenerfassungsschritt erfasst wurden, nachdem der Korrekturbetrag berechnet wurde (S707).

20. Computerprogramm das Anweisungen umfasst, die, wenn das Programm von einem Computer ausgeführt wird, den Computer veranlassen, Schritte auszuführen, die umfassen:

einen Datenerfassungsschritt (S701) zum Erfassen von Messdaten, bei denen es sich um Zeitreihendaten von Messwerten handelt, die zu messen sind, indem ein Schienenfahrzeug, das einen Fahrzeugkörper (11), ein Drehgestell (12a, 12b) und einen Radsatz (13, 13a-13d) aufweist, veranlasst wird, auf einem Gleis (16) zu fahren;

einen ersten Gleiszustandsberechnungsschritt (S705) des Berechnens eines geschätzten Wertes einer ersten physikalischen Größe;

einen Korrekturbetragsberechnungsschritt (S707) des Berechnens eines Korrekturbetrags für den geschätzten Wert der ersten physikalischen Größe basierend auf dem geschätzten Werts der ersten physikalischen Größe, der durch den ersten Gleiszustandsberechnungsschritt (S705) berechnet wurde, und einem tatsächlichen Wert der ersten physikalischen Größe;

einen zweiten Gleiszustandsberechnungsschritt des Berechnens eines geschätzten Wertes der ersten physikalischen Größe, nachdem der Korrekturbetrag berechnet ist (S707); und

einen Gleiszustandskorrekturschritt zum Korrigieren des geschätzten Wertes der ersten physikalischen Größe, der durch den zweiten Gleiszustandsberechnungsschritt berechnet wurde, unter Verwendung des Korrekturbetrags, wobei

die gemessenen Daten einen gemessenen Wert einer vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) enthalten,

die vorwärts- und rückwärtsgerichtete Kraft (T₁-T₄) eine Kraft in einer Vorwärts- und Rückwärtsrichtung (x) ist, die in einem Element auftritt, das zwischen dem Radsatz (13, 13a-13d) und dem Drehgestell (12a, 12b) angeordnet ist, auf dem der Radsatz (13, 13a-13d) vorgesehen ist,

das Element ein Element zum Tragen eines Achslagers ist,

die Vorwärts- und Rückwärtsrichtung (x) eine Richtung entlang einer Fahrtrichtung des Schienenfahrzeugs ist,

die erste physikalische Größe eine physikalische Größe ist, die einen Zustand des Gleises (16) wiedergibt,

der erste Gleiszustandsberechnungsschritt (S705) und der zweite Gleiszustandsberechnungsschritt einen relationalen Ausdruck, der die Beziehung zwischen der ersten physikalischen Größe an einer Position des Radsatzes (13, 13a-13d) und der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) darstellt, und einen gemessenen Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄) verwenden, um den geschätzten Wert der ersten physikalischen Größe zu berechnen,

der gemessene Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄), der in dem ersten Gleiszustandsberechnungsschritt (S705) verwendet wird, in den gemessenen Daten enthalten ist, die durch den Datenerfassungsschritt (S701) erfasst werden, bevor der Korrekturbetrag berechnet wird (S707), und

der gemessene Wert der vorwärts- und rückwärtsgerichteten Kraft (T₁-T₄), der in dem zweiten Gleiszustandsberechnungsschritt verwendet wird, in den gemessenen Daten enthalten ist, die durch den Datenerfassungsschritt erfasst wurden, nachdem der Korrekturbetrag berechnet wurde (S707).

Revendications

1. Système d'inspection (400), comprenant :

un moyen d'acquisition de données (403, 2011a, 2011b) configuré pour acquérir des données mesurées qui sont des données chronologiques de valeurs mesurées à mesurer en provoquant le déplacement d'un véhicule ferroviaire comprenant un châssis de véhicule (11), un bogie (12a, 12b) et un train de roues (13, 13a-13d) sur une voie ferrée (16) ;

un premier moyen de calcul d'état de voie ferrée (407) configuré pour calculer une valeur estimée d'une première quantité physique ;

un moyen de calcul de quantité de correction (409) configuré pour calculer une quantité de correction pour la valeur estimée de la première quantité physique sur la base de la valeur estimée de la première quantité physique calculée par le premier moyen de calcul d'état de voie ferrée (407) et d'une valeur réelle de la première quantité physique ;

un second moyen de calcul d'état de voie ferrée (411) configuré pour calculer une valeur estimée de la première quantité physique après que la quantité de correction a été calculée ; et

un moyen de correction d'état de voie ferrée (412) configuré pour corriger la valeur estimée de la première quantité physique calculée par le second moyen de calcul d'état de voie ferrée (411) en utilisant la quantité de correction,

dans lequel

les données mesurées contiennent une valeur mesurée d'une force en marche avant et en marche arrière (T₁-T₄),

la force en marche avant et en marche arrière (T₁-T₄) est une force dans le sens de la marche avant et de la marche arrière (x) destinée à se produire dans un élément disposé entre le train de roues (13, 13a-13d) et le bogie (12a, 12b) sur lequel le train de roues (13, 13a-13d) est prévu,

l'élément est un élément destiné à supporter une boîte d'essieu,

le sens de la marche avant et de la marche arrière (x) est une direction le long d'une direction de déplacement du véhicule ferroviaire,

la première quantité physique est une quantité physique qui reflète un état de la voie ferrée (16),

le premier moyen de calcul d'état de voie ferrée (407) et le second moyen de calcul de voie ferrée (411) sont configurés pour utiliser une expression relationnelle qui représente la relation entre la première quantité physique à un emplacement du train de roues (13, 13a-13d) et la force en marche avant et en marche arrière (T₁-T₄) et une valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄) afin de calculer la valeur estimée de la première quantité physique,

la valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄) utilisée dans le premier moyen de calcul d'état de voie ferrée (407) est contenue dans les données mesurées acquises par le moyen d'acquisition de données (403, 2011a, 2011b) avant que la quantité de correction ne soit calculée, et

la valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄) utilisée dans le second moyen de calcul d'état de voie ferrée (411) est contenue dans les données mesurées acquises par le moyen d'acquisition de données (403, 2011a, 2011b) après que la quantité de correction a été calculée.

2. Système d'inspection (400) selon la revendication 1, dans lequel
le moyen de calcul de quantité de correction (409) est configuré pour calculer, sur la base des valeurs estimées d'une pluralité des premières quantités physiques calculées par le premier moyen de calcul d'état de voie ferrée (407) en utilisant les valeurs mesurées de la force en marche avant et en marche arrière (T₁-T₄) lorsque le véhicule ferroviaire se déplace au même emplacement et des valeurs réelles des premières quantités physiques, une quantité de correction à l'emplacement en guise de quantité de correction pour la valeur estimée de la première quantité physique.

3. Système d'inspection (400) selon la revendication 1 ou 2, dans lequel

le moyen de calcul de quantité de correction (409) est configuré pour calculer, sur la base des valeurs estimées d'une pluralité des premières quantités physiques calculées par le premier moyen de calcul d'état de voie ferrée (407) en utilisant les valeurs mesurées de la force en marche avant et en marche arrière (T₁-T₄) lorsque le véhicule ferroviaire se déplace au même emplacement, d'une vitesse de déplacement du véhicule ferroviaire qui se déplace à l'emplacement, et des valeurs réelles des premières quantités physiques, une quantité de correction pour la valeur estimée de la première quantité physique, et

la quantité de correction pour la valeur estimée de la première quantité physique est une quantité de correction selon l'emplacement et la vitesse de déplacement du véhicule ferroviaire.

4. Système d'inspection (400) selon l'une quelconque des revendications 1 à 3, dans lequel
les valeurs mesurées de la force en marche avant et en marche arrière (T₁-T₄) utilisées par le premier moyen de calcul d'état de voie ferrée (407) et le second moyen de calcul d'état de voie ferrée (411) sont des valeurs mesurées sur le même véhicule ferroviaire.

5. Système d'inspection (400) selon l'une quelconque des revendications 1 à 4, comprenant en outre :

un moyen de réglage de fréquence (404, 406) configuré pour réduire une intensité de signal d'une composante à basse fréquence générée en raison du véhicule ferroviaire qui se déplace sur une partie incurvée de la voie ferrée (16) à partir des données chronologiques d'une seconde quantité physique, dans lequel

la seconde quantité physique est une quantité physique dont la valeur varie selon un état du véhicule ferroviaire,

le moyen de réglage de fréquence (404, 406) comprend un premier moyen de réglage de fréquence (404) configuré pour réduire une intensité de signal d'une composante à basse fréquence générée en raison du véhicule ferroviaire qui se déplace sur une partie incurvée de la voie ferrée (16) à partir des données chronologiques d'une valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄) qui est l'une de la seconde quantité physique,

le premier moyen de calcul d'état de voie ferrée (407) et le second moyen de calcul d'état de voie ferrée (411) sont configurés pour utiliser l'expression relationnelle et la valeur de la force en marche avant et en marche arrière (T₁-T₄) à partir de laquelle l'intensité de signal de la composante à basse fréquence a été réduite par le premier moyen de réglage de fréquence (404) afin de calculer la valeur estimée de la première quantité physique, et

l'expression relationnelle est une expression qui n'inclut pas un rayon de courbure (R) d'un rail (16a).

6. Système d'inspection (400) selon la revendication 5, dans lequel

le moyen de réglage de fréquence (404, 406) est configuré pour utiliser les données chronologiques de la seconde quantité physique et détermine un coefficient dans un modèle autorégressif corrigé, et utilise le modèle autorégressif corrigé dont le coefficient est déterminé et les données chronologiques de la seconde quantité physique et réduit une intensité de signal d'une composante à basse fréquence générée en raison du véhicule ferroviaire qui se déplace sur une partie incurvée de la voie ferrée (16) à partir des données chronologiques de la seconde quantité physique,

le modèle autorégressif corrigé est une expression qui représente une valeur prédite de la seconde quantité physique en utilisant une valeur de la seconde quantité physique et le coefficient en réponse à la valeur,

le moyen de réglage de fréquence (404, 406) est configuré pour utiliser une équation dans laquelle une première matrice est définie comme une matrice à coefficient et un vecteur d'autocorrélation est défini comme un vecteur constant et détermine le coefficient,

le vecteur d'autocorrélation est un vecteur dont la composante est une autocorrélation des données chronologiques de la seconde quantité physique avec un retard de 1 à m, m étant un chiffre de la valeur mesurée utilisée dans le modèle autorégressif corrigé,

la première matrice est une matrice U_sΣ_sU_s^T dérivée d'une seconde matrice Σ_s qui est dérivée de s valeurs propres d'une matrice d'autocorrélation, s étant un chiffre défini pour être égal ou supérieur à 1 et inférieur à m, et d'une matrice diagonale Σ et d'une troisième matrice U_s qui est dérivée de s des valeurs propres et d'une matrice orthogonale U,

la matrice d'autocorrélation est une matrice dont la composante est une autocorrélation des données chronologiques de la seconde quantité physique avec un retard de 0 à m -1,

la matrice diagonale est une matrice dont la composante diagonale correspond aux valeurs propres de la matrice d'autocorrélation qui sont dérivées par la décomposition en valeurs singulières de la matrice d'autocorrélation,

la matrice orthogonale est une matrice dans laquelle un vecteur propre de la matrice d'autocorrélation est défini comme un vecteur de composante de colonne,

la seconde matrice est une sous-matrice de la matrice diagonale et est une matrice dont la composante diagonale correspond à s des valeurs propres, et

la troisième matrice est une sous-matrice de la matrice orthogonale et est une matrice dans laquelle les vecteurs propres qui correspondent à s des valeurs propres sont définis comme des vecteurs de composantes de colonnes.

7. Système d'inspection (400) selon la revendication 6, dans lequel
parmi les valeurs propres de la matrice d'autocorrélation, la valeur propre la plus élevée est incluse à s des valeurs propres.

8. Système d'inspection (400) selon l'une quelconque des revendications 1 à 7, comprenant en outre :

un moyen de fonctionnement de filtre (405) configuré pour effectuer une opération à l'aide d'un filtre qui effectue une assimilation de données en utilisant les données mesurées, une équation d'état et une équation d'observation, et détermine ainsi les valeurs estimées de variables d'état qui sont des variables destinées à déterminer les valeurs estimées dans l'équation d'état, dans lequel

les données mesurées contiennent en outre des valeurs mesurées d'accélérations du bogie (12a, 12b) et du train de roues (13, 13a-13d) dans la direction de gauche et de droite (y),

la direction de gauche et de droite (y) est une direction perpendiculaire au sens de la marche avant et de la marche arrière (x) et à une direction vers le haut et vers le bas (z) qui est une direction perpendiculaire à la voie ferrée (16),

la force en marche avant et en marche arrière (T₁-T₄) est une force à déterminer selon une différence entre un déplacement angulaire du train de roues (13, 13a-13d) dans une direction de lacet et un déplacement angulaire du bogie (12a, 12b) sur lequel le train de roues (13, 13a-13d) est prévu dans la direction de lacet,

la direction de lacet est une direction de pivotement avec la direction vers le haut et vers le bas (z) définie comme un axe de pivotement,

l'équation d'état est une équation décrite en utilisant les variables d'état, la force en marche avant et en marche arrière (T₁-T₄) et une variable de transformation,

les variables d'état comprennent un déplacement et une vitesse du bogie (12a, 12b) dans la direction de gauche et de droite (y), un déplacement angulaire et une vitesse angulaire du bogie (12a, 12b) dans la direction de lacet, un déplacement angulaire et une vitesse angulaire du bogie (12a, 12b) dans une direction de roulis, un déplacement et une vitesse du train de roues (13, 13a-13d) dans la direction de gauche et de droite (y), et un déplacement angulaire d'un ressort pneumatique fixé sur le véhicule ferroviaire dans la direction de roulis, et n'incluent pas un déplacement angulaire ou une vitesse angulaire du train de roues (13, 13a-13d) dans la direction de lacet,

la direction de roulis est une direction de pivotement avec le sens de la marche avant et de la marche arrière (x) défini comme un axe de pivotement,

la variable de transformation est une variable qui effectue une transformation mutuelle entre le déplacement angulaire du train de roues (13, 13a-13d) dans la direction de lacet et le déplacement angulaire du bogie (12a, 12b) dans la direction de lacet,

l'équation d'observation est une équation décrite en utilisant une variable d'observation et la variable de transformation,

la variable d'observation inclut les accélérations du bogie (12a, 12b) et du train de roues (13, 13a-13d) dans la direction de la droite et de la gauche (y),

le moyen de fonctionnement de filtre (405) est configuré pour utiliser une valeur mesurée de la variable d'observation, l'équation d'état dans laquelle la valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄) et une valeur réelle de la variable de transformation sont remplacées, et l'équation d'observation dans laquelle la valeur réelle de la variable de transformation est remplacée, et détermine les valeurs estimées des variables d'état lorsqu'une erreur entre la variable d'observation, la valeur mesurée et une valeur estimée ou une valeur prévue de l'erreur devient minimum,

le premier moyen de calcul d'état de voie ferrée (407) et le second moyen de calcul d'état de voie ferrée (411) sont configurés pour utiliser une valeur estimée du déplacement angulaire du bogie (12a, 12b) dans la direction de lacet, qui est l'une des variables d'état déterminées par le moyen de fonctionnement de filtre (405), et la valeur réelle de la variable de transformation pour calculer une valeur estimée du déplacement angulaire du train de roues (13, 13a-13d) dans la direction de lacet, et utilisent la valeur estimée du déplacement angulaire du train de roues (13, 13a-13d) dans la direction de lacet, la valeur de la force en marche avant et en marche arrière (T₁-T₄), et l'expression relationnelle pour calculer la valeur estimée de la première quantité physique,

l'expression relationnelle est une expression dans laquelle une équation de déplacement qui décrit le déplacement du train de roues (13, 13a-13d) dans la direction de lacet est exprimée en utilisant la force en marche avant et en marche arrière (T₁-T₄), et

la valeur réelle de la variable de transformation est dérivée en utilisant la valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄).

9. Système d'inspection (400) selon la revendication 8, dans lequel

l'équation d'état est constituée en utilisant une équation de déplacement qui décrit le déplacement du train de roues (13, 13a-13d) dans la direction de gauche et de droite (y), une équation de déplacement qui décrit le déplacement du bogie (12a, 12b) dans la direction de gauche et de droite (y), une équation de déplacement qui décrit le déplacement du bogie (12a, 12b) dans la direction de lacet, une équation de déplacement qui décrit le déplacement du bogie (12a, 12b) dans la direction de roulis, et une équation de déplacement qui décrit le déplacement du ressort pneumatique dans la direction de roulis,

l'équation de déplacement qui décrit le déplacement du train de roues (13, 13a-13d) dans la direction de gauche et de droite (y) est une équation de déplacement décrite en utilisant la variable de transformation au lieu d'utiliser la déplacement angulaire du train de roues (13, 13a-13d) dans la direction de lacet,

l'équation de déplacement qui décrit le déplacement du bogie (12a, 12b) dans la direction de lacet est une équation de déplacement décrite en utilisant la force en marche avant et en marche arrière (T₁-T₄) au lieu d'utiliser le déplacement angulaire et la vitesse angulaire du train de roues (13, 13a-13d) dans la direction de lacet, et

la variable de transformation est représentée par une différence entre le déplacement angulaire du bogie (12a, 12b) dans la direction de lacet et le déplacement angulaire du train de roues (13, 13a-13d) dans la direction de lacet.

10. Système d'inspection (400) selon la revendication 8 ou 9, dans lequel

le moyen d'acquisition de données (403, 2011a, 2011b) est configuré pour acquérir également une valeur mesurée d'accélération du châssis du véhicule (11) dans la direction de droite et de gauche (y),

la variable d'observation inclut en outre l'accélération du châssis du véhicule (11) dans la direction de droite et de gauche (y),

les variables d'état incluent en outre un déplacement et une vitesse du châssis du véhicule (11) dans la direction de droite et de gauche (y), un déplacement angulaire et une vitesse angulaire du châssis du véhicule (11) dans la direction de lacet, un déplacement angulaire et une vitesse angulaire du châssis du véhicule (11) dans la direction de roulis, et un déplacement angulaire d'un amortisseur de lacet fixé sur le véhicule ferroviaire dans la direction de lacet, et

le moyen de fonctionnement de filtre (405) est configuré pour déterminer les variables d'état lorsque les différences entre les accélérations du châssis du véhicule (11), du bogie (12a, 12b) et du train de roues (13, 13a-13d) dans la direction de droite et de gauche (y), les valeurs mesurées et les valeurs calculées deviennent minimum.

11. Système d'inspection (400) selon la revendication 10, dans lequel
l'équation d'état est constituée en utilisant en outre une équation de déplacement qui décrit le déplacement du châssis du véhicule (11) dans la direction de droite et de gauche (y), une équation de déplacement qui décrit le déplacement du châssis du véhicule (11) dans la direction de lacet, une équation de déplacement qui décrit le déplacement du châssis du véhicule (11) dans la direction de roulis, et une équation de déplacement qui décrit le déplacement de l'amortisseur de lacet dans la direction de lacet.

12. Système d'inspection (400) selon l'une quelconque des revendications 8 à 11, dans lequel

l'équation d'observation est constituée en utilisant en outre l'équation de déplacement qui décrit le déplacement du train de roues (13, 13a-13d) dans la direction de droite et de gauche (y) et l'équation de déplacement qui décrit le déplacement du bogie (12a, 12b) dans la direction de droite et de gauche (y), et

l'équation de déplacement qui décrit le déplacement du train de roues (13, 13a-13d) dans la direction de droite et de gauche (y) est une équation de déplacement décrite en utilisant la variable de transformation au lieu d'utiliser le déplacement angulaire du train de roues (13, 13a-13d) dans la direction de lacet.

13. Système d'inspection (400) selon la revendication 12, dans lequel
l'équation d'observation est constituée en utilisant en outre l'équation de déplacement qui décrit le déplacement du châssis du véhicule (11) dans la direction de droite et de gauche (y).

14. Système d'inspection (400) selon l'une quelconque des revendications 8 à 13, dans lequel

le premier moyen de calcul d'état de voie ferrée (407) et le second moyen de calcul d'état de voie ferrée (411) sont configurés pour dériver une quantité d'irrégularité d'alignement (y_R1-y_R4) de la voie ferrée (16) en guise de valeur estimée de la première quantité physique sur la base du déplacement et de la vitesse du bogie (12a, 12b) dans la direction de droite et de gauche (y) qui sont les variables d'état déterminées par le moyen de fonctionnement de filtre (405), du déplacement et de la vitesse du train de roues (13, 13a-13d) dans la direction de droite et de gauche (y) qui sont les variables d'état déterminées par le moyen de fonctionnement de filtre (405), de la valeur estimée du déplacement angulaire du train de roues (13, 13a-13d) dans la direction de lacet, de la valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄), et de l'équation de déplacement qui décrit le déplacement du train de roues (13, 13a-13d) dans la direction de lacet, et

l'équation de déplacement qui décrit le déplacement du train de roues (13, 13a-13d) dans la direction de lacet inclut la force en marche avant et en marche arrière (T₁-T₄) et la quantité d'irrégularité d'alignement (y_R1-y_R4) de la voie ferrée (16) comme des variables.

15. Système d'inspection (400) selon l'une quelconque des revendications 8 à 14, comprenant en outre :

un moyen de réglage de fréquence (404, 406) configuré pour réduire une intensité de signal d'une composante à basse fréquence générée en raison du véhicule ferroviaire qui se déplace sur une partie incurvée de la voie ferrée (16) à partir des données chronologiques d'une seconde quantité physique, dans lequel

la seconde quantité physique est une quantité physique qui varie selon un état du véhicule ferroviaire, et

le moyen de réglage de fréquence comprend un second moyen de réglage de fréquence (406) configuré pour réduire une intensité de signal d'une composante à basse fréquence générée en raison du véhicule ferroviaire qui se déplace sur une partie incurvée de la voie ferrée (16) à partir des données chronologiques des valeurs estimées des variables d'état qui sont l'une de la seconde quantité physique.

16. Système d'inspection (400) selon l'une quelconque des revendications 8 à 15, dans lequel
le filtre est un filtre de Kalman.

17. Système d'inspection (400) selon l'une quelconque des revendications 1 à 16, dans lequel

la première quantité physique est une quantité d'irrégularité d'alignement (y_R1-y_R4) de la voie ferrée (16) ou une force latérale (Q₁-Q₄) qui est une contrainte dans la direction de droite et de gauche (y) entre une roue (14L, 14R, 14a-14d) prévue sur le train de roues (13, 13a-13d) et la voie ferrée (16), et

la direction de droite et de gauche (y) est une direction perpendiculaire au sens de la marche avant et de la marche arrière (x) et à la direction vers le haut et vers le bas (z) qui est une direction perpendiculaire à la voie ferrée (16).

18. Système d'inspection (400) selon l'une quelconque des revendications 1 à 17, dans lequel

la force en marche avant et en marche arrière (T₁-T₄) est, parmi les composantes dans le sens de la marche avant et de la marche arrière (x) des forces qui se produisent dans les deux éléments fixés sur les deux côtés du train de roues (13, 13a-13d) dans la direction de droite et de gauche (y), la composante opposée en phase aux autres, et

la direction de droite et de gauche (y) est une direction perpendiculaire au sens de la marche avant et de la marche arrière (x) et à la direction vers le haut et vers le bas (z) qui est une direction perpendiculaire à la voie ferrée (16).

19. Procédé d'inspection exécuté par un système d'inspection, le procédé d'inspection comprenant :

une étape d'acquisition de données (S701) comprenant l'acquisition de données mesurées qui sont les données chronologiques de valeurs mesurées à mesurer en provoquant le déplacement d'un véhicule ferroviaire comprenant un châssis de véhicule (11), un bogie (12a, 12b) et un train de roues (13, 13a-13d) sur une voie ferrée (16) ;

une première étape de calcul d'état de voie ferrée (S705) comprenant le calcul d'une valeur estimée d'une première quantité physique ;

une étape de calcul de quantité de correction (S707) comprenant le calcul d'une quantité de correction pour la valeur estimée de la première quantité physique sur la base de la valeur estimée de la première quantité physique calculée par la première étape de calcul d'état de voie ferrée (S705) et d'une valeur réelle de la première quantité physique ;

une seconde étape de calcul d'état de voie ferrée comprenant le calcul d'une valeur estimée de la première quantité physique une fois que la quantité de correction a été calculée (S707) ; et

une étape de correction d'état de voie ferrée comprenant la correction de la valeur estimée de la première quantité physique calculée par la seconde étape de calcul d'état de voie ferrée à l'aide de la quantité de correction, dans lequel

les données mesurées contiennent une valeur mesurée d'une force en marche avant et en marche arrière (T₁-T₄),

la force en marche avant et en marche arrière (T₁-T₄) est une force dans le sens de la marche avant et de la marche arrière (x) destinée à se produire dans un élément disposé entre le train de roues (13, 13a-13d) et le bogie (12a, 12b) sur lequel le train de roues (13, 13a-13d) est prévu,

l'élément est un élément destiné à supporter une boîte d'essieu,

le sens de la marche avant et de la marche arrière (x) est une direction le long d'une direction de déplacement du véhicule ferroviaire,

la première quantité physique est une quantité physique qui reflète un état de la voie ferrée (16),

la première étape de calcul d'état de voie ferrée (S705) et la seconde étape de calcul de voie ferrée utilisent une expression relationnelle qui représente la relation entre la première quantité physique à un emplacement du train de roues (13, 13a-13d) et la force en marche avant et en marche arrière (T₁-T₄) et une valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄) afin de calculer la valeur estimée de la première quantité physique,

la valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄) utilisée à la première étape de calcul d'état de voie ferrée (S705) est contenue dans les données mesurées acquises par l'étape d'acquisition de données (S701) avant que la quantité de correction ne soit calculée (S707), et

la valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄) utilisée à la seconde étape de calcul d'état de voie ferrée est contenue dans les données mesurées acquises par l'étape d'acquisition de données après que la quantité de correction a été calculée (S707).

20. Programme informatique comprenant des instructions qui, lorsque le programme est exécuté par un ordinateur, permettent à l'ordinateur d'exécuter des étapes comprenant :

une étape d'acquisition de données (S701) comprenant l'acquisition de données mesurées qui sont les données chronologiques de valeurs mesurées à mesurer en provoquant le déplacement d'un véhicule ferroviaire comprenant un châssis de véhicule (11), un bogie (12a, 12b) et un train de roues (13, 13a-13d) sur une voie ferrée (16) ;

une première étape de calcul d'état de voie ferrée (S705) comprenant le calcul d'une valeur estimée d'une première quantité physique ;

une étape de calcul de quantité de correction (S707) comprenant le calcul d'une quantité de correction pour la valeur estimée de la première quantité physique sur la base de la valeur estimée de la première quantité physique calculée par la première étape de calcul d'état de voie ferrée (S705) et d'une valeur réelle de la première quantité physique ;

une seconde étape de calcul d'état de voie ferrée comprenant le calcul d'une valeur estimée de la première quantité physique une fois que la quantité de correction a été calculée (S707) ; et

une étape de correction d'état de voie ferrée comprenant la correction de la valeur estimée de la première quantité physique calculée par la seconde étape de calcul d'état de voie ferrée à l'aide de la quantité de correction, dans lequel

les données mesurées contiennent une valeur mesurée d'une force en marche avant et en marche arrière (T₁-T₄),

la force en marche avant et en marche arrière (T₁-T₄) est une force dans le sens de la marche avant et de la marche arrière (x) destinée à se produire dans un élément disposé entre le train de roues (13, 13a-13d) et le bogie (12a, 12b) sur lequel le train de roues (13, 13a-13d) est prévu,

l'élément est un élément destiné à supporter une boîte d'essieu,

le sens de la marche avant et de la marche arrière (x) est une direction le long d'une direction de déplacement du véhicule ferroviaire,

la première quantité physique est une quantité physique qui reflète un état de la voie ferrée (16),

la première étape de calcul d'état de voie ferrée (S705) et la seconde étape de calcul de voie ferrée utilisent une expression relationnelle qui représente la relation entre la première quantité physique à un emplacement du train de roues (13, 13a-13d) et la force en marche avant et en marche arrière (T₁-T₄) et une valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄) afin de calculer la valeur estimée de la première quantité physique,

la valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄) utilisée à la première étape de calcul d'état de voie ferrée (S705) est contenue dans les données mesurées acquises par l'étape d'acquisition de données (S701) avant que la quantité de correction ne soit calculée (S707), et

la valeur mesurée de la force en marche avant et en marche arrière (T₁-T₄) utilisée à la seconde étape de calcul d'état de voie ferrée est contenue dans les données mesurées acquises par l'étape d'acquisition de données après que la quantité de correction a été calculée (S707).

Drawing