EP3895955A1 - Inspection system, inspecting method, and program - Google Patents
Inspection system, inspecting method, and program Download PDFInfo
- Publication number
- EP3895955A1 EP3895955A1 EP19894582.6A EP19894582A EP3895955A1 EP 3895955 A1 EP3895955 A1 EP 3895955A1 EP 19894582 A EP19894582 A EP 19894582A EP 3895955 A1 EP3895955 A1 EP 3895955A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- equation
- backward
- state
- motion
- wheel set
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000007689 inspection Methods 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims description 52
- 230000001133 acceleration Effects 0.000 claims abstract description 64
- 230000009466 transformation Effects 0.000 claims abstract description 54
- 230000033001 locomotion Effects 0.000 claims description 190
- 238000006073 displacement reaction Methods 0.000 claims description 147
- 238000009795 derivation Methods 0.000 claims description 88
- 238000005096 rolling process Methods 0.000 claims description 61
- 230000008859 change Effects 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 2
- 230000002123 temporal effect Effects 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 description 117
- 238000012545 processing Methods 0.000 description 40
- 239000000725 suspension Substances 0.000 description 33
- 238000004364 calculation method Methods 0.000 description 21
- 238000005070 sampling Methods 0.000 description 17
- 239000013598 vector Substances 0.000 description 16
- 238000004891 communication Methods 0.000 description 14
- 241001247986 Calotropis procera Species 0.000 description 13
- 238000000354 decomposition reaction Methods 0.000 description 12
- 238000009826 distribution Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 230000005540 biological transmission Effects 0.000 description 8
- 230000014509 gene expression Effects 0.000 description 7
- 238000013016 damping Methods 0.000 description 6
- 238000013500 data storage Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 230000005484 gravity Effects 0.000 description 5
- 238000001914 filtration Methods 0.000 description 4
- 238000005309 stochastic process Methods 0.000 description 4
- 230000003247 decreasing effect Effects 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 206010000117 Abnormal behaviour Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 206010027175 memory impairment Diseases 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L23/00—Control, warning or like safety means along the route or between vehicles or trains
- B61L23/04—Control, warning or like safety means along the route or between vehicles or trains for monitoring the mechanical state of the route
- B61L23/042—Track changes detection
- B61L23/047—Track or rail movements
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61K—AUXILIARY EQUIPMENT SPECIALLY ADAPTED FOR RAILWAYS, NOT OTHERWISE PROVIDED FOR
- B61K9/00—Railway vehicle profile gauges; Detecting or indicating overheating of components; Apparatus on locomotives or cars to indicate bad track sections; General design of track recording vehicles
- B61K9/08—Measuring installations for surveying permanent way
Definitions
- the present invention relates to an inspection system, an inspection method, and a program, and is suitably used for inspecting a track of a railway vehicle, in particular.
- the present application is based upon and claims the benefit of priority of the prior Japanese Patent Application No. 2018-230834 , filed in Japan on December 10, 2018, the entire contents of which are incorporated herein by reference.
- Patent Literature 1 discloses a technique in which angular displacements of wheel sets in a yawing direction, state variables derived by a filter that performs data assimilation, and a measured value of a forward-and-backward-direction force being a force in a forward and backward direction that occurs in a member for supporting an axle box are substituted into motion equations that describe yawings of the wheel sets to derive track irregularity (alignment irregularity amount and the like) of a railway vehicle.
- Patent Literature 1 International Publication Pamphlet No. WO 2017/164133
- the measured value of the forward-and-backward-direction force and measured values of accelerations of the wheel set and a bogie (and a vehicle body in addition to the above, according to need), respectively, in a right and left direction are used when performing the data assimilation.
- These measured values can be obtained without using special sensors, but the number of sensors disposed in the railway vehicle is preferably as small as possible.
- the present invention has been made in view of the problems as described above, and an object thereof is to reduce the number of sensors which are used for detecting track irregularity of a railway vehicle.
- An inspection system of the present invention includes: a data acquisition means that acquires data of a measured value of a forward-and-backward-direction force as data of a measured value to be measured by causing a railway vehicle including a vehicle body, a bogie, and a wheel set to travel on a track; a state variable derivation means that derives state variables being variables to be determined in a state equation constituted by using motion equations that describe motions of the railway vehicle, by using the measured value of the forward-and-backward-direction force; and a track state derivation means that derives information reflecting a state of the track, in which the forward-and-backward-direction force is a force in a forward and backward direction that occurs in a member disposed between the wheel set and the bogie on which the wheel set is provided, and is a force to be determined according to a difference between an angular displacement of the wheel set in a yawing direction and an angular displacement of the bogie on which the wheel set is provided in the
- An inspection method of the present invention includes: a data acquisition step of acquiring data of a measured value of a forward-and-backward-direction force as data of a measured value to be measured by causing a railway vehicle including a vehicle body, a bogie, and a wheel set to travel on a track; a state variable derivation step of deriving state variables being variables to be determined in a state equation constituted by using motion equations that describe motions of the railway vehicle, by using the measured value of the forward-and-backward-direction force; and a track state derivation step of deriving information reflecting a state of the track, in which the forward-and-backward-direction force is a force in a forward and backward direction that occurs in a member disposed between the wheel set and the bogie on which the wheel set is provided, and is a force to be determined according to a difference between an angular displacement of the wheel set in a yawing direction and an angular displacement of the bogie on which the wheel set is provided in
- a program of the present invention causes a computer to execute: a data acquisition step of acquiring data of a measured value of a forward-and-backward-direction force as data of a measured value to be measured by causing a railway vehicle including a vehicle body, a bogie, and a wheel set to travel on a track; a state variable derivation step of deriving state variables being variables to be determined in a state equation constituted by using motion equations that describe motions of the railway vehicle, by using the measured value of the forward-and-backward-direction force; and a track state derivation step of deriving information reflecting a state of the track, in which the forward-and-backward-direction force is a force in a forward and backward direction that occurs in a member disposed between the wheel set and the bogie on which the wheel set is provided, and is a force to be determined according to a difference between an angular displacement of the wheel set in a yawing direction and an angular displacement of the bogie on which the wheel
- Fig. 1 is a view illustrating one example of an outline of a railway vehicle.
- 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).
- the z axis is set to a direction vertical to a track 16 (the ground) (a height direction of the railway vehicle).
- the y axis is set to a horizontal direction vertical to the traveling direction of the railway vehicle (a direction vertical to both the traveling direction and the height direction of the railway vehicle).
- the railway vehicle is set to a commercial vehicle.
- the mark of ⁇ added inside ⁇ indicates the direction from the far side of the sheet toward the near side
- the mark of ⁇ added inside ⁇ indicates the direction from the near side of the sheet toward the far side.
- the railway vehicle includes a vehicle body 11, bogies 12a, 12b, and wheel sets 13a to 13d.
- 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.
- a case where each of the bogies 12a, 12b is a bolsterless bogie will be explained as an example. Note that in Fig.
- 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 .
- the bogies 12a, 12b have bogie frames, bolster springs, and so on.
- an axle box is disposed on both sides of each of the wheel sets 13a to 13d in the direction along the y axis.
- 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 transmitted to the railway vehicle from the track 16.
- 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. Note that the railway vehicle itself can be fabricated by a well-known technique, and thus its detailed explanation is omitted here.
- Fig. 2 is a view conceptually illustrating directions of 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.
- 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, will be explained as an example.
- 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
- the direction along the x axis is referred to as the forward and backward direction as necessary.
- the forward and backward direction is the traveling direction of the railway vehicle.
- the direction along the x axis is set to the traveling direction of the railway vehicle.
- 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
- the direction along the z axis is referred to as the up and down direction as necessary.
- the up and down direction is a direction vertical to the track 16.
- 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.
- the right and left direction is a direction vertical to both the forward and backward direction (the traveling direction of the railway vehicle) and the up and down direction (the direction vertical to the track 16).
- 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.
- the state variables are derived by performing filtering with a filter (Kalman filter) that performs data assimilation by using, as observation variables, accelerations y w1 • •, y w2 • •, y w3 • •, y w4 • • of the wheel sets 13a, 13b, 13c, 13d in the right and left direction, accelerations y t1 • •, y t2 • • of the bogies 12a, 12b in the right and left direction, and in addition to the above, acceleration y b • • of the vehicle body 11 in the right and left direction according to need.
- a filter Kalman filter
- Fig. 3 illustrates a measured value and a calculated value of each of the acceleration y t1 • • of the bogie 12a in the right and left direction and the accelerations y w1 • •, y w2 • • of the wheel sets 13a, 13b in the right and left direction.
- the calculated value is an estimated value of the observation variable calculated by data assimilation.
- the horizontal axis in Fig. 3 indicates an elapsed time (second) from a reference time when the reference time is set to 0 (zero). Concretely, the horizontal axis in Fig.
- the estimated value of the state variable is derived so that an error between a value which is normally given as the measured value of the observation variable and the estimated value of the observation variable becomes minimum or an expected value of this error becomes minimum.
- the measured values of the accelerations y w1 • •, y w2 • • of the wheel sets 13a, 13b in the right and left direction, and the measured values of the accelerations y t1 • •, y t2 • • of the bogie 12a in the right and left direction include a lot of noises.
- the present inventors obtained a finding that, depending on a state of the track 16 (rail), the estimated value does not approximate the measured value and becomes substantially a fixed value even if the data assimilation is performed.
- the present inventors considered that there may be a chance that, when deriving the state variables, the alignment irregularity amounts y R1 , y R2 , y R3 , y R4 can be derived without greatly decreasing accuracy, even if the measured values of these accelerations are not used.
- the present inventors devised a method of calculating an alignment irregularity amount by using a measured value of force in the forward and backward direction that occurs in a member disposed between the wheel sets 13a, 13b (13c, 13d) and the bogie 12a (12b) on which these wheel sets 13a, 13b (13c, 13d) are provided.
- the force in the forward and backward direction that occurs in the member is referred to as a forward-and-backward-direction force as necessary.
- the alignment irregularity amount is calculated by using an equation representing a 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.
- 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.
- 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 data assimilation
- state variables may diverge. Therefore, the state equation in the case of performing filtering with the filter (Kalman filter) that performs data assimilation is constituted by using a motion equation that describes motion of the railway vehicle traveling on the linear track.
- the motion equation describing the motion of the railway vehicle traveling on the curved track includes a term including a curvature radius of the rail. Therefore, when the state variables are derived by 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.
- the filter Karl filter
- 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.
- the component itself of the forward-and-backward-direction force due to the alignment irregularity occurs in the same manner even on the curved track or the linear track.
- the present inventors thought that the alignment irregularity amount itself has nothing to do with the amount of the aforementioned bias, and 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 if the filter (Kalman filter) performing the 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.
- the present inventors devised calculating the alignment irregularity amount by using the time-series data of the value of the forward-and-backward-direction force from which the low-frequency component has been reduced.
- 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.
- the calculating equation of the alignment irregularity amount results in the same calculating equation even on the curved track or the linear track. Note that there is a chance that even a track, which is designed as a linear track, actually has a curvature which may exert an influence on estimation accuracy of the alignment irregularity amount.
- the reduction in the low-frequency component (behavior of the aforementioned bias) from the time-series data of the measured value of the forward-and-backward-direction force contributes to improvement of the estimation accuracy of the alignment irregularity amount on not only the curved track but also the linear track.
- explanation will be made by setting that even a track, which is designed as a linear track but actually has a curvature which may exert an influence on estimation accuracy of the alignment irregularity amount, is also regarded as a curved track.
- the degree of freedom is not limited to 21 degrees of freedom.
- 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 by taking these points into consideration.
- 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 Non-Patent Literatures 1, 2, for example. Therefore, outlines of the respective motion equations will be explained here, and their detailed explanations are omitted.
- 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.
- 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, w4 indicate the wheel sets 13a, 13b, 13c, 13d respectively.
- 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.
- Subscripts b, B indicate the vehicle body 11.
- a subscript x indicates the forward and backward direction or the rolling direction
- a subscript y indicates the right and left direction
- a subscript z indicates the up and down direction or the yawing direction.
- 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 applies to the other variables below)).
- f 2 is a lateral creep coefficient (note that the lateral creep coefficient f 2 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 applies to 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 1 is a distance between the center 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.
- 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 1 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 1 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 1 ).
- ⁇ 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 center 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.
- Kwx is a spring constant of the axle box suspension in the forward and backward direction.
- 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. 4A and Fig. 4B each illustrate one example of the alignment irregularity amount y R1 at the position of the wheel set 13a.
- Fig. 4A the case of the track 16 being the linear track will be explained as an example.
- Fig. 4B the case of the track 16 being the curved track will be explained as an example.
- 16a denotes a rail
- 16b denotes a crosstie.
- 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 , y R4 at the positions of the wheel sets 13b, 13c, 13d are also defined in the same manner as the alignment irregularity amount y R1 at the position of the wheel set 13a.
- 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' 2 is a damping constant of a lateral movement damper.
- h 4 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 5 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' 2 is a spring constant of the air spring (bolster spring) in the right and left direction.
- h 2 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 3 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.
- 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' 0 is stiffness of a rubber bush of the yaw damper.
- b' 0 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' 0 ).
- ⁇ y1 is a pivot amount (angular displacement) of the yaw damper disposed on the bogie 12a in the yawing direction.
- k" 2 is a spring constant of the air spring (bolster spring) in the forward and backward direction.
- b 2 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 2 ).
- respective variables in (12) Equation are represented by being replaced with the variables in (11) Equation according to the meanings of the subscripts described above.
- 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 1 is a damping constant of an axle damper in the up and down direction.
- b' 1 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' 1 ).
- c 2 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 1 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 2 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 3 is equivalent stiffness by a change in effective pressure receiving area of the air spring (bolster spring).
- ⁇ a2 is a pivot amount (angular displacement) of the air spring (bolster spring) disposed on the bogie 12b in the rolling direction.
- 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.
- 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 0 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.
- 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.
- ⁇ y2 is a pivot amount (angular displacement) of the yaw damper disposed on the bogie 12b in the yawing direction.
- ⁇ a2 • is an angular velocity of the air spring (bolster spring) disposed on the bogie 12b in the rolling direction.
- 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.
- the forward-and-backward-direction force becomes a component opposite in phase to each other out of forward-and-backward-direction components of forces that occur in the aforementioned two members attached to both the right and left sides of one wheel set.
- the axle box suspension being a mono-link type axle box suspension
- the axle box suspension includes a link
- the axle box and the bogie frame are coupled by the link.
- a rubber bush is attached to both ends of the link.
- 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.
- 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.
- 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.
- a forward-and-backward-direction displacement of the rubber bush attached to the link may be measured by a displacement meter.
- 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.
- 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.
- the load measured by the strain gauge attached to the link sometimes includes 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.
- 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.
- 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.
- the axle box suspension being an axle beam type axle box suspension
- the axle box suspension includes an axle beam
- 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.
- 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.
- 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.
- 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.
- 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.
- a forward-and-backward-direction displacement of the rubber bush attached to the axle beam may be measured by a displacement meter.
- 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.
- 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.
- the load measured by the strain gauge attached to the axle beam sometimes includes 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.
- 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, even if the strain gauge is not 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.
- the components other than the component in the forward and backward direction are sometimes included in the measured forward-and-backward-direction force, and thus 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.
- the axle box suspension being a leaf spring type axle box suspension
- the axle box suspension includes a leaf spring
- the axle box and the bogie frame are coupled by the leaf spring.
- a rubber bush is attached to ends of the leaf spring.
- 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.
- 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.
- 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.
- 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.
- a forward-and-backward-direction displacement of the rubber bush attached to the leaf spring may be measured by a displacement meter.
- 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.
- 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.
- 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.
- the system of the axle box suspension is not limited to the mono-link type, the axle beam type, and the leaf spring type.
- 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.
- 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 1 to T 4 .
- Fig. 5 is a view illustrating one example of a functional configuration of an inspection apparatus 500.
- Fig. 6 is a view illustrating one example of a hardware configuration of the inspection apparatus 500.
- Fig. 7 is a flowchart illustrating one example of processing in the inspection apparatus 500.
- the inspection apparatus 500 is mounted on the railway vehicle will be explained as an example.
- the inspection apparatus 500 includes, as its functions, a storage unit 501, a data acquisition unit 502, a first frequency adjustment unit 503, a state variable derivation unit 504, a second frequency adjustment unit 505, a track state derivation unit 506, and an output unit 507.
- the inspection apparatus 500 includes a CPU 601, a main memory 602, an auxiliary memory 603, a communication circuit 604, a signal processing circuit 605, an image processing circuit 606, an I/F circuit 607, a user interface 608, a display 609, and a bus 610.
- the CPU 601 overall controls the entire inspection apparatus 500.
- the CPU 601 uses the main memory 602 as a work area to execute a program stored in the auxiliary memory 603.
- the main memory 602 stores data temporarily.
- the auxiliary memory 603 stores various kinds of data, in addition to programs to be executed by the CPU 601.
- the auxiliary memory 603 stores later-described state equations and observation equations.
- the storage unit 501 is fabricated by using the CPU 601 and the auxiliary memory 603, for example.
- the communication circuit 604 is a circuit intended for performing communication with the outside of the inspection apparatus 500.
- the communication circuit 604 receives information of the measured value of the forward-and-backward-direction force, for example.
- the communication circuit 604 may perform radio communication or wire communication with the outside of the inspection apparatus 500.
- the communication circuit 604 is connected to an antenna provided on the railway vehicle in the case of performing radio communication.
- the signal processing circuit 605 performs various kinds of signal processing on signals received by the communication circuit 604 and signals input according to the control made by the CPU 601.
- the data acquisition unit 502 is fabricated by using the CPU 601, the communication circuit 604, and the signal processing circuit 605, for example.
- the first frequency adjustment unit 503, the state variable derivation unit 504, the second frequency adjustment unit 505, and the track state derivation unit 506 are fabricated by using the CPU 601 and the signal processing circuit 605, for example.
- the image processing circuit 606 performs various kinds of image processing on signals input according to the control made by the CPU 601.
- the signal after being subjected to the image processing is output to the display 609.
- the user interface 608 is a part through which an operator gives an instruction to the inspection apparatus 500.
- the user interface 608 includes buttons, switches, dials, and so on, for example. Further, the user interface 608 may include a graphical user interface using the display 609.
- the display 609 displays an image based on a signal output from the image processing circuit 606.
- the I/F circuit 607 exchanges data with a device connected to the I/F circuit 607.
- the user interface 608 and the display 609 are illustrated.
- the device to be connected to the I/F circuit 607 is not limited to these.
- a portable storage medium may be connected to the I/F circuit 607.
- at least a part of the user interface 608 and the display 609 may be provided outside the inspection apparatus 500.
- the output unit 507 is fabricated by using at least any one of a set including the communication circuit 604, and the signal processing circuit 605, and a set including the image processing circuit 606, the I/F circuit 607, and the display 609, for example.
- the CPU 601, the main memory 602, the auxiliary memory 603, the signal processing circuit 605, the image processing circuit 606, and the I/F circuit 607 are connected to the bus 610. Communication among these components is performed via the bus 610. Further, the hardware of the inspection apparatus 500 is not limited to the one illustrated in Fig. 6 as long as it can realize later-described functions of the inspection apparatus 500.
- the storage unit 501 stores equations which are used when the later-described state variable derivation unit 504 derives the state variables.
- the storage unit 501 stores state equations and observation equations.
- Equation to (8) Equation (the motion equations that describe the yawings of the wheel sets 13a to 13d) are not included in the state equation, and the state equation is constituted as follows.
- Equation and (10) Equation the motion equations that describe the transversal vibrations of the bogies 12a, 12b (motion in the right and left direction)), (13) Equation and (14) Equation (the motion equations that describe the rollings of the bogies 12a, 12b), (15) Equation (the motion equation that describes the transversal vibration of the vehicle body 11 (motion in the right and left direction)), (16) Equation (the motion equation that describes the yawing of the vehicle body 11), (17) Equation (the motion equation that describes the rolling of the vehicle body 11), (18) Equation and (19) 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), and (20) Equation and (21) Equation (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) are used as they
- Equation to (4) Equation the motion equations that describe the transversal vibrations of the wheel sets 13a to 13d (motion in the right and left direction)
- (11) Equation and (12) Equation the motion equations that describe the yawings of the bogies 12a, 12b
- 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.
- the forward-and-backward-direction forces T 1 to T 4 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 and ⁇ t2 of the bogies on which these wheel sets are provided in the yawing direction.
- T 1 C wx b 1 2 ⁇ ⁇ t1 ⁇ ⁇ ⁇ w1 + K wx b 1 2 ⁇ t1 ⁇ ⁇ w1
- T 2 C wx b 1 2 ⁇ ⁇ t1 ⁇ ⁇ ⁇ w2 + K wx b 1 2 ⁇ t1 ⁇ ⁇ w2
- T 3 C wx b 1 2 ⁇ ⁇ t2 ⁇ ⁇ w3 + K wx b 1 2 ⁇ t2 ⁇ w3
- T 4 C wx b 1 2 ⁇ ⁇ t2 ⁇ ⁇ ⁇ w4 + K wx b 1 2 ⁇ t2 ⁇ ⁇ w4
- Transformation variables e 1 to e 4 are defined as in (26) Equation to (29) Equation below.
- the transformation variables e 1 to e 4 are defined by the differences between the angular displacements ⁇ t1 and ⁇ 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 1 to e 4 are variables for performing mutual transformation between the angular displacements ⁇ t1 and ⁇ t2 of the bogies in the yawing direction and the angular displacements ⁇ w1 to ⁇ w4 of the wheel sets in the yawing direction.
- Equation to (33) Equation below are obtained.
- 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.
- Equation to (4) Equation the motion equations that describe the transversal vibrations of the wheel sets 13a to 13d (motion in the right and left direction)
- e 1 to e 4 the transformation variables that describe the transversal vibrations of the wheel sets 13a to 13d (motion in the right and left direction)
- e 1 to e 4 the transformation variables that describe the transversal vibrations of the wheel sets 13a to 13d (motion in the right and left direction)
- Equation and (12) Equation are expressed by using the forward-and-backward-direction forces T 1 to T 4 , 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.
- Equation to (43) Equation below are obtained.
- 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.
- the state equation is constituted by using (34) Equation to (39) Equation.
- (40) Equation to (43) Equation are ordinary differential equations. Actual values of the transformation variables e 1 to e 4 , which are solutions of the ordinary differential equations, can be derived by using the values of the forward-and-backward-direction forces T 1 to T 4 in the wheel sets 13a to 13d.
- the values of the forward-and-backward-direction forces T 1 to T 4 are obtained by reducing a signal strength of a low-frequency component to be generated due to the railway vehicle traveling on the curved portion of the track from the time-series data of the measured value of the forward-and-backward-direction force, with the use of the later-described first frequency adjustment unit 503.
- the values of the transformation variables e 1 to e 4 derived as above are given to (34) Equation to (37) Equation. Further, the values of the forward-and-backward-direction forces T 1 to T 4 in the wheel sets 13a to 13d are given to (38) Equation and (39) Equation.
- the values of the forward-and-backward-direction forces T 1 to T 4 are obtained by reducing a signal strength of a low-frequency component to be generated due to the railway vehicle traveling on the curved portion of the track from the time-series data of the measured value of the forward-and-backward-direction force, with the use of the later-described first frequency adjustment unit 503.
- 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.
- the storage unit 501 receives the state equation constituted as above, for example, based on the operation of the user interface 608 made by an operator and stores it.
- 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.
- (34) Equation to (37) Equation, (9) Equation, (10) Equation, and (15) Equation are used to constitute an observation equation.
- the storage unit 501 receives the observation equation constituted in this manner, for example, based on the operation of the user interface 608 made by an operator and stores it.
- the data acquisition unit 502, the first frequency adjustment unit 503, the state variable derivation unit 504, the second frequency adjustment unit 505, the track state derivation unit 506, and the output unit 507 start.
- the processing according to the flowchart in Fig. 7 starts after the state equation and the observation equation are stored in the inspection apparatus 500.
- the data acquisition unit 502 acquires 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.
- the data acquisition unit 502 can acquire the time-series data of the measured value of the forward-and-backward-direction force by performing communication with an arithmetic device that calculates the forward-and-backward-direction force by using a measured value of a strain gauge for measuring the forward-and-backward-direction force, for example.
- the data acquisition unit 502 does not acquire 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 first frequency adjustment unit 503 reduces (preferably removes) the signal strength of the low-frequency component included in the time-series data of the measured value of the forward-and-backward-direction force acquired by the data acquisition unit 502.
- a signal of this low-frequency component is a signal that is not measured when the railway vehicle is traveling on the linear track whose curvature is 0 (zero), 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 whose curvature is 0 (zero).
- the present inventors devised a model in which an AR (Auto-regressive) model is corrected. Further, the present inventors came up with an idea of reducing the signal strength of the low-frequency component included in the time-series data of the measured value of the forward-and-backward-direction force by using this model.
- the model devised by the present inventors is referred to as a corrected AR model.
- the well-known AR model is referred to as an AR model simply.
- a value of time-series data of a physical quantity y at a time k (1 ⁇ k ⁇ M) is set to y k .
- M is a number indicating, as the time-series data of the physical quantity y, data until when is contained, and is preset.
- 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.
- y ⁇ k is expressed by adding ⁇ above y k 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-1 to y k-m of the data y at continuous times k-1 to k-m prior to the time k.
- m is an integer of less than M. As m, for example, 1500 can be used.
- 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 .
- 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 .
- Equation 21 Equation 21
- Equation 21 Equation 21
- R j0 R j1 R j2 ⁇ R jm ⁇ 1 ⁇ 2 ⁇ ⁇ m , 1 ⁇ j ⁇ m
- R jl in (48) Equation is called autocorrelation of the data y, and is a value defined by (49) Equation below.
- at this time is referred to as a time lag.
- 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 .
- Equation is called a Yule-Walker equation.
- (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.
- 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.
- the coefficient matrix on the right side in (50) Equation is referred to as an autocorrelation matrix as necessary.
- 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 of this, the present inventors found out that it is possible to reduce the influence of a high-frequency component included 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.
- 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 ⁇ 11 , ⁇ 22 , ⁇ ⁇ ⁇ , ⁇ mm.
- 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 orthogonal matrix U is set to u 1 , u 2 , ⁇ , u m .
- 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.
- ⁇ ⁇ 11 0 0 0 ⁇ 0 0 ⁇ mm
- U u 1 u 2 ⁇ u m
- ⁇ 11 , ⁇ 22 , ⁇ , ⁇ 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 the present embodiment, s is preset.
- a matrix U s in (54) Equation is a matrix of m ⁇ s composed of s pieces of the column component vectors, which are chosen from the left of the orthogonal matrix U of (52) Equation (eigenvectors corresponding to the eigenvalues to be used). 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 ⁇ .
- U S u 1 u 2 ⁇ u S
- Equation is modified, and thereby (57) Equation below is obtained as the equation for 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] ⁇ 1 ⁇ 2 ⁇ ⁇ m U S ⁇ S ⁇ 1 U S T R 10 R 20 ⁇ R m0
- 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).
- 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.
- the matrix U s in (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.
- 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 ⁇ s U s T in (57) Equation is a matrix (a first matrix) derived from the matrix ⁇ s and the matrix U s .
- 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.
- 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.
- This autocorrelation of the stochastic process is expressed as a function of a time lag. Therefore, the autocorrelation of the data y in the present embodiment may be replaced with a value calculated by another calculating formula as long as it approximates the autocorrelation of the stochastic process.
- R 22 to R mm are autocorrelation with a time lag of 0 (zero), but they may be replaced with R 11 .
- 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.
- the forward-and-backward-direction force is applied here.
- the value of the forward-and-backward-direction force varies according to the state and the like of the railway vehicle.
- 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. 8 is a view illustrating one example of the distribution of the eigenvalues of the autocorrelation matrix R. In Fig.
- eigenvalues ⁇ 11 to ⁇ mm which are obtained by the autocorrelation matrix R in each of the time-series data of the measured value y of the forward-and-backward-direction force T 1 in the wheel set 13a being subjected to singular value decomposition, are aligned in ascending order and are plotted.
- the horizontal axis is an index of the eigenvalue and the vertical axis is the value of the eigenvalue.
- the number of eigenvalue having a value which is significantly higher than the other values may change depending on the configuration of the railway vehicle, the configuration of the track, and so on. Therefore, the number s of eigenvalues to be extracted from the autocorrelation matrix R is not limited to these values as long as it is one or more.
- the first frequency adjustment unit 503 performs the following processing every time the data acquisition unit 502 acquires the value y k of the time-series data of the measured value y of the forward-and-backward-direction force at the time k at a predetermined sampling period.
- the first frequency adjustment unit 503 generates the autocorrelation matrix R by using (49) Equation and (51) Equation based on the time-series data of the measured value y of the forward-and-backward-direction force and preset numbers M, m.
- the first frequency adjustment unit 503 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 ⁇ 11 to ⁇ mm of the autocorrelation matrix R from the diagonal matrix ⁇ .
- the first frequency adjustment unit 503 chooses s pieces of the eigenvalues ⁇ 11 to ⁇ ss from the largest from among the plural eigenvalues ⁇ 11 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.
- the first frequency adjustment unit 503 determines the coefficient ⁇ of the corrected AR model by using (57) Equation based on the time-series data of the measured value y of the forward-and-backward-direction force, the eigenvalues ⁇ 11 to ⁇ ss , and the orthogonal matrix U obtained by the singular value decomposition of the autocorrelation matrix R.
- the first frequency adjustment unit 503 derives the predicted value y ⁇ k of the time-series data of the measured value y 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 time-series data of the measured value y of the forward-and-backward-direction force.
- Time-series data of the predicted value y ⁇ k of the forward-and-backward-direction force results in the time-series data from which the low-frequency component included in the time-series data of the measured value y of the forward-and-backward-direction force has been extracted.
- Fig. 9 is a view illustrating one example of time-series data of measured values of forward-and-backward-direction forces (measured values) and time-series data of predicted values of the forward-and-backward-direction forces (calculated values). Note that in the present embodiment, the measured values of the four forward-and-backward-direction forces T 1 to T 4 are obtained. That is, four pieces of the data y of the forward-and-backward-direction force are obtained. In Fig. 9 , the measured value and the calculated value of each of the four pieces of data y are illustrated. The horizontal axis in Fig.
- the vertical axis indicates the forward-and-backward-direction forces T 1 to T 4 (Nm).
- the calculated value of the forward-and-backward-direction force T 1 in the wheel set 13a is biased at about 15 seconds to 35 seconds. Specifically, the calculated value of the forward-and-backward-direction force T 1 in the wheel set 13a exhibits, at about 15 seconds to 35 seconds, a value larger than that at another time. 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 2 in the wheel set 13b, the calculated value of the forward-and-backward-direction force T 3 in the wheel set 13c, and the calculated value of the forward-and-backward-direction force T 4 in the wheel set 13d are also biased during the period when the wheel sets 13b, 13c, 13d pass through the curved track, similarly to the calculated value of the forward-and-backward-direction force T 1 in the wheel set 13a.
- the first frequency adjustment unit 503 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.
- 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 is referred to as time-series data of a high-frequency component of the forward-and-backward-direction force as necessary.
- 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.
- Fig. 10 is a view illustrating one example of the time-series data of the high-frequency components of the forward-and-backward-direction forces.
- the vertical axis in Fig. 10 indicates the time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 , T 2 , T 3 , T 4 . That is, the high-frequency components of the forward-and-backward-direction forces T 1 , T 2 , T 3 , T 4 illustrated on the vertical axis in Fig.
- Fig. 10 are ones obtained by subtracting the calculated values from the measured values of the forward-and-backward-direction forces T 1 , T 2 , T 3 , T 4 in the wheel sets 13a, 13b, 13c, 13d that are illustrated in Fig. 9 respectively.
- the horizontal axis in Fig. 10 indicates a measuring time and a calculating time of the forward-and-backward-direction forces T 1 to T 4 , 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. 9 .
- the first frequency adjustment unit 503 derives the time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 to T 4 as above.
- the state variable derivation unit 504 sets the observation equation as the observation equation stored by the storage unit 501, sets the state equation as the state equation stored by the storage unit 501, and determines estimated values of the state variables illustrated in (44) Equation by using the Kalman filter. At this time, the state variable derivation unit 504 uses time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 to T 4 generated by the first frequency adjustment unit 503.
- the time-series data of the measured value of the acceleration of the vehicle body 11 in the right and left direction is not used when determining the estimated values of the state variables.
- 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 reduce (minimize) the difference between, of an observable variable (observation variable), a measured value and an estimated value.
- the state variable derivation unit 504 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.
- the following observation equation of (58) Equation and the following state equation of (59) Equation are used.
- 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. Note that the Kalman filter itself can be fabricated by a well-known technique, and thus its detailed explanation is omitted.
- the measured values (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 used as they are.
- the measured value not the measured value but a preset fixed value is given as the value to be normally given as the measured value of the observation variable when performing the data assimilation.
- an average value of the pieces of time-series data of the accelerations is assumed to be 0 (zero), and all of fixed values to be given as the observation variables are set to 0 (zero). Therefore, in the present embodiment, when performing the data assimilation, the state variable derivation unit 504 derives the estimated value of the state variable so as to minimize an error of the estimated value of the observation variable with respect to the fixed value (0 (zero) in this case) or minimize an expected value of the error.
- the state variable derivation unit 504 determines the estimated values of the state variables illustrated in (44) Equation at a predetermined sampling period, to thereby generate time-series data of the estimated values of the state variables illustrated in (44) Equation.
- the signal strength of the low-frequency component included in the time-series data of the measured value of the forward-and-backward-direction force is not removed sufficiently by the first frequency adjustment unit 503
- the signal of the low-frequency component due to the railway vehicle traveling on the curved track may remain in the time-series data of the estimated values of the state variables generated by the state variable derivation unit 504.
- the second frequency adjustment unit 505 reduces (preferably removes) the signal strength of the low-frequency component included in the time-series data of the estimated values of the state variables generated by the state variable derivation unit 504.
- the second frequency adjustment unit 505 uses the corrected AR model to reduce the signal strength of the low-frequency component included in the time-series data of the estimated values of the state variables similarly to the first frequency adjustment unit 503.
- the second frequency adjustment unit 505 performs the following processing for each state variable at a predetermined sampling period.
- the state variable is applied here. That is, the data y of the state variable results in the time-series data of the estimated values of the state variables generated by the state variable derivation unit 504. Each of the estimated values of the state variables varies according to the state of the railway vehicle.
- the second frequency adjustment unit 505 generates the autocorrelation matrix R by 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.
- the second frequency adjustment unit 505 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 ⁇ 11 to ⁇ mm of the autocorrelation matrix R from the diagonal matrix ⁇ .
- the second frequency adjustment unit 505 chooses s pieces of the eigenvalues ⁇ 11 to ⁇ ss from the largest from among the plural eigenvalues ⁇ 11 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.
- the data y of the estimated value of each state variable can be obtained in a manner as explained so far in a state of causing the railway vehicle to travel on the track 16, for example.
- a distribution of the eigenvalues of the autocorrelation matrix R is made by the second frequency adjustment unit 505 individually for each state variable.
- the second frequency adjustment unit 505 determines the number s of the eigenvalues to be extracted from the autocorrelation matrix R illustrated in (53) Equation for each of the state variables.
- the second frequency adjustment unit 505 determines the coefficient ⁇ of the corrected AR model by using (57) Equation based on the data y of the estimated value of the state variable, the eigenvalues ⁇ 11 to ⁇ ss , and the orthogonal matrix U obtained by the singular value decomposition of the autocorrelation matrix R.
- the second frequency adjustment unit 505 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 results in the time-series data from which the low-frequency component included in the data y of the estimated value of the state variable has been extracted.
- the second frequency adjustment unit 505 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.
- a value of 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 at each sampling time is referred to as a value of a high-frequency component of the state variable as necessary.
- the track state derivation unit 506 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. Subsequently, the track state derivation unit 506 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 by the second frequency adjustment unit 505, and the values of the high-frequency components of the forward-and-backward-direction forces T 1 to T 4 generated by the first frequency adjustment unit 503 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 and y t2 of the bogies 12a, 12b in the right and left direction, the velocities y t1 ⁇ and 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 track state derivation unit 506 performs such a calculation of the alignment irregularity amounts y R1 to y R4 as above at a predetermined sampling period, to thereby obtain the time-series data of the alignment irregularity amounts y R1 to y R4 .
- the track state derivation unit 506 calculates a final alignment irregularity amount y R from the alignment irregularity amounts y R1 to y R4 .
- the track state derivation unit 506 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 track state derivation unit 506 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 track state derivation unit 506 displaces the phases of the time-series data of the alignment irregularity amounts y R2 to y R4 by this delay time.
- the track state derivation unit 506 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 final alignment irregularity amount y R at this sampling time.
- the track state derivation unit 506 performs such a calculation at each sampling time, to thereby obtain time-series data of the final 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 .
- the track state derivation unit 506 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 final alignment irregularity amount y R from the alignment irregularity amounts y R1 to y R4 whose moving averages have been found.
- the track state derivation unit 506 may calculate, as the final 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.
- the inspection apparatus 500 uses the time-series data of the measured value of the forward-and-backward-direction force at each sampling time acquired by the data acquisition unit 502 while the railway vehicle is traveling in the traveling section being the target of deriving the alignment irregularity amount, to execute pieces of the processing of the first frequency adjustment unit 503, the state variable derivation unit 504, the second frequency adjustment unit 505, and the track state derivation unit 506.
- the track state derivation unit 506 can obtain the alignment irregularity amount y R at each sampling time while the railway vehicle is traveling in the traveling section being the target of deriving the alignment irregularity amount.
- the track state derivation unit 506 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.
- the traveling position of the railway vehicle can be set to the position of the wheel set 13a, for example.
- the track state derivation unit 506 derives the final 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.
- the track state derivation unit 506 does not always need to calculate the traveling position of the railway vehicle at each sampling time as described previously.
- the track state derivation unit 506 may derive the traveling position of the railway vehicle at each sampling time by using a GPS (Global Positioning System).
- the output unit 507 outputs information of the final alignment irregularity amount y R that is calculated by the track state derivation unit 506. At this time, the output unit 507 may output information indicating that the track 16 is abnormal in the case where the final alignment irregularity amount y R is larger than a preset value.
- 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.
- the inspection apparatus 500 gives the measured values of the forward-and-backward-direction forces T 1 to T 4 and the actual values of the transformation variables e 1 to e 4 to the Kalman filter, to derive the state variables (y w1 ⁇ to y w4 ⁇ , y w1 to y w4 , y t1 ⁇ and y t2 ⁇ , y t1 and y t2 , ⁇ t1 ⁇ and ⁇ t2 ⁇ , ⁇ t1 and ⁇ t2 , ⁇ t1 ⁇ and ⁇ t2 ⁇ , ⁇ t1 and ⁇ t2 , y t1 and ⁇ t2 , y t1 ⁇ ⁇ t2 , y t1 ⁇ ⁇ t2 , ⁇ t1 and ⁇ t2 , y b ⁇ , y b , ⁇ b ⁇ , ⁇ b ⁇ , ⁇ b
- the inspection apparatus 500 uses the pivot amounts (angular displacements) ⁇ t1 and ⁇ t2 of the bogies 12a, 12b in the yawing direction included in the state variables and the actual values of the transformation variables e 1 to e 4 , to derive the pivot amounts (angular displacements) ⁇ w1 to ⁇ w4 of the wheel sets 13a to 13d in the yawing direction.
- the inspection apparatus 500 substitutes the pivot amounts (angular displacements) ⁇ w1 to ⁇ w4 of the wheel sets 13a to 13d in the yawing direction, the state variables, and the measured values of the forward-and-backward-direction forces T 1 to T 4 into the motion equations that describe the yawings of the wheel sets 13a to 13d, to derive the alignment irregularity amounts y R1 to y R4 at the positions of the wheel sets 13a to 13d. Subsequently, the inspection apparatus 500 derives the final alignment irregularity amount y R from the alignment irregularity amounts y R1 to y R4 .
- the alignment irregularity amounts y R1 to y R4 (the final alignment irregularity amount y R ) without greatly decreasing accuracy, with no use of the measured values of the accelerations of the vehicle body 11, the bogies 12a, 12b, and the wheel sets 13a to 13d in the right and left direction. Therefore, it is possible to reduce the number of sensors which are used for deriving the alignment irregularity amounts y R1 to y R4 (the final alignment irregularity amount y R ).
- the inspection apparatus 500 generates the autocorrelation matrix R from the time-series data of the measured value y of the forward-and-backward-direction force, and determines, 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, the coefficient ⁇ of the corrected AR model approximating the time-series data of the measured value y of the forward-and-backward-direction force. Therefore, it is possible to determine the coefficient ⁇ so as to make the signal of the low-frequency component included in the time-series data of the measured value y of the forward-and-backward-direction force remain and prevent the high-frequency component from remaining.
- the inspection apparatus 500 calculates the predicted value y ⁇ k of the forward-and-backward-direction force at the time k by giving the time-series data of the measured value y of the forward-and-backward-direction force at the time k-l (1 ⁇ 1 ⁇ m), which is prior to the time k, to the corrected AR model whose coefficient ⁇ is determined in this manner. Therefore, 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 time-series data of the measured value y of the forward-and-backward-direction force without estimating a cutoff frequency beforehand.
- the inspection apparatus 500 reduces the signal strength of the low-frequency components included in the time-series data of the measured values of the forward-and-backward-direction forces T 1 to T 4 in this manner and generates the time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 to T 4 .
- the inspection apparatus 500 gives the time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 to T 4 to the relational expression between the forward-and-backward-direction forces T 1 to T 4 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)).
- the preset fixed value is given as the value to be normally given as the measured value of the observation variable when performing the data assimilation.
- This fixed value is not limited to 0 (zero).
- the time-series data of the measured value of the acceleration of the vehicle body 11 in the right and left direction, the time-series data of the measured values of the accelerations of the bogies 12a, 12b in the right and left direction, and the time-series data of the measured values of the accelerations of the wheel sets 13a to 13d in the right and left direction when the railway vehicle having the inspection apparatus 500 mounted thereon or a railway vehicle equivalent to this railway vehicle (a railway vehicle whose structure is the same as that of this railway vehicle) travels on the track 16 being the target of deriving the alignment irregularity amounts y R1 to y R4 (the final alignment irregularity amount y R ) are obtained, and an average value of the respective pieces of time-series data is used as the fixed value.
- time-series data of a predicted value of acceleration y ⁇ k of the vehicle body 11 in the right and left direction, time-series data of predicted values of accelerations y ⁇ k of the bogies 12a, 12b in the right and left direction, and time-series data of predicted values y ⁇ k of accelerations of the wheel sets 13a to 13d in the right and left direction are derived by the above-described corrected AR model. Subsequently, an average value thereof may also be used as the fixed value.
- the case of using the corrected AR model has been explained as an example.
- 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 time-series data of the measured value y of the forward-and-backward-direction force.
- the processing of the first frequency adjustment unit 503 and the processing of the second frequency adjustment unit 505 are no longer required.
- the wheel set to be a standard when performing the phase matching is the wheel set 13a
- the wheel set to be a standard may be the wheel set 13b, 13c, or 13d other than the wheel set 13a.
- the case of using the Kalman filter has been explained as an example.
- a filter that derives the estimated values of the state variables so that the error of the estimated value of the observation variable with respect to the fixed value becomes minimum or the expected value of this error becomes minimum that is, a filter that performs data assimilation
- a particle filter may be used. Note that as the error of the estimated value of the observation variable with respect to the fixed value, for example, a square error between the estimated value of the observation variable and the fixed value is cited.
- 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 information that reflects the track irregularity (appearance failure of the track 16) is derived as information reflecting the state of the track 16.
- the following calculations of (64) Equation to (67) Equation may be performed, to thereby derive a lateral force that occurs when the railway vehicle travels on the linear track (stress in the right and left direction between the wheel and the rail).
- Q 1 , Q 2 , Q 3 , Q 4 are lateral forces in the wheels 14a, 14b, 14c, 14d respectively.
- f 3 represents a spin creep coefficient.
- the case of including the state variables that represent the state of the vehicle body 11 has been explained as an example.
- 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 influence by the propagation in the vehicle body 11 is judged to be small, for example.
- Equation to (21) Equation (15) Equation to (17) Equation (the motion equations that describe the transversal vibration, the yawing, and the rolling of the vehicle body 11) and (18) Equation and (19) Equation (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 no longer required.
- the state amount relating to the vehicle body (state amount including the subscript of b) and the value inside ⁇ that includes the state amount relating to the vehicle body (state amount including the subscript of b) are set to 0 (zero) .
- the bogies 12a, 12b each being a bolsterless bogie has been explained as an example.
- the bogies 12a, 12b are not limited to the bolsterless bogies.
- the motion equations are rewritten appropriately. That is, the motion equations are not limited to the ones exemplified in the present embodiment.
- the state equation includes a term that represents this external force.
- the state variables are derived by using the filter (Kalman filter) that performs the data assimilation by setting the values to be normally given as the measured values of the observation variables (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) to the fixed value (0 (zero)), when performing the data assimilation, has been described as an example.
- the filter Karl filter
- the present embodiment and the first embodiment are different mainly in the method of deriving the state variables (the function of the state variable derivation unit 504). Therefore, in the explanation of the present embodiment, the same reference numerals and symbols as those added to Fig. 1 to Fig. 10 are added to the same parts as those in the first embodiment, or the like, and their detailed explanations are omitted.
- the storage unit 501 does not store the state equation ((58) Equation) and the observation equation ((59) Equation), but stores the following motion equation of (68) Equation.
- X ⁇ c ⁇ X
- (68) Equation is one example of equation obtained by changing an equation representing the motion equations of (9) Equation, (10) Equation, (13) Equation to (21) Equation, and (34) Equation to (39) Equation by using the state variables illustrated in (44) Equation (an equation corresponding to (68) Equation in which c is set to 1), so that a temporal change in the state variables becomes smaller than that of this equation.
- (68) Equation is one obtained in a manner that, in an equation representing the motion equations of (9) Equation, (10) Equation, (13) Equation to (21) Equation, and (34) Equation to (39) Equation by using the state variables illustrated in (44) Equation, a term which is connected to a term of first-order time differential (X ⁇ ) of the state variable by an equal sign is multiplied by a forgetting factor c.
- (68) Equation corresponds to the state equation of (59) Equation in which the forgetting factor c is introduced, to thereby set the system noise W to 0 (zero).
- the forgetting factor c is a preset value, and is (theoretically) a value of greater than 0 and 1 or less (0 ⁇ c ⁇ 1).
- the forgetting factor c functions in a manner that the smaller the value thereof, the greater the degree of forgetfulness of a past observation value.
- the smaller the value of the forgetting factor c the smaller the influence of the measured value of the forward-and-backward-direction force on the estimated value (solution) of the state variable. For this reason, from a viewpoint of correctly obtaining the estimated value (solution) of the state variable, the value of the forgetting factor c is desirably close to 1.
- (68) Equation is directly solved without using the filter that performs the data assimilation. Accordingly, there is a need to suppress the divergence of the estimated value (solution) of the state variable.
- the value of the forgetting factor c is determined based on the viewpoint as described above.
- the forgetting factor c for example, a value of greater than 0.0 and 1.0 or less (0.0 ⁇ c ⁇ 1.0), preferably a value of 0.90 or more and 1.0 or less (0.90 ⁇ c ⁇ 1.0), more preferably a value of 0.95 or more and 1.0 or less (0.95 ⁇ c ⁇ 1.0), still more preferably a value of 0.99 or more and 1.0 or less (0.99 ⁇ c ⁇ 1.0), and most preferably 1.0 is selected among them.
- the forgetting factor c is selected so as to prevent the divergence of the estimated value (solution) of the state variable obtained by solving (68) Equation. If the divergence of the estimated value (solution) of the state variable obtained by solving (68) Equation does not occur, the estimated value (solution) of this state variable when the value of the forgetting factor is 1.0 becomes the most accurate solution. However, when the value of the forgetting factor c is 1.0, it is highly likely that the estimated value (solution) of the state variable obtained by solving (68) Equation diverges (the solution is not determined).
- the forgetting factor c may be selected by setting an upper limit value of the forgetting factor c to less than 1.0.
- the value of the forgetting factor c can be selected from, for example, a value of greater than 0.0 and less than 1.0 (0.0 ⁇ c ⁇ 1.0), preferably a value of 0.90 or more and less than 1.0 (0.90 ⁇ c ⁇ 1.0), more preferably a value of 0.95 or more and less than 1.0 (0.95 ⁇ c ⁇ 1.0), and still more preferably a value of 0.99 or more and less than 1.0 (0.99 ⁇ c ⁇ 1.0).
- Equation corresponds to the motion equations themselves of (9) Equation, (10) Equation, (13) Equation to (21) Equation, and (34) Equation to (39) Equation (an equation only representing the motion equations by using the state variables).
- the state variable derivation unit 504 uses the time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 to T 4 generated by the first frequency adjustment unit 503 to derive the actual values of the transformation variables e 1 to e 4 and substitutes the actual values into (34) Equation to (37) Equation, and substitutes the time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 to T 4 generated by the first frequency adjustment unit 503 into (38) Equation and (39) Equation as the measured values of the forward-and-backward-direction forces T 1 to T 4 to solve the equation of (68) Equation, thereby determining the estimated values of the state variables illustrated in (44) Equation.
- the method of solving the equation of (68) Equation can be realized by a well-known numerical solution (Euler method or the like), for example. Therefore, when deriving the estimated values of the state variables, the state variable derivation unit 504 does not use the time-series data of the measured value of the acceleration of the vehicle body 11 in the right and left direction, the time-series data of the measured values of the accelerations of the bogies 12a, 12b in the right and left direction, and the time-series data of the measured values of the accelerations of the wheel sets 13a to 13d in the right and left direction. Besides, the observation equation is also not used.
- the inspection apparatus 500 gives the measured values of the forward-and-backward-direction forces T 1 to T 4 and the actual values of the transformation variables e 1 to e 4 to the equation being the state equation with the system noise W set to 0 (zero) in which the term other than the time differential term X ⁇ of the state variable is multiplied by the forgetting factor c, to thereby derive the state variables (y w1 ⁇ to y w4 ⁇ , y w1 to y w4 , y t1 ⁇ and y t2 ⁇ , y t1 and y t2 , ⁇ t1 ⁇ and ⁇ t2 ⁇ , ⁇ t1 and ⁇ t2 , ⁇ t1 ⁇ and ⁇ t2 ⁇ , ⁇ t1 and ⁇ t2 , ⁇ t1 and ⁇ t2 ⁇ , ⁇ t1 and ⁇ t2 , y t1 and ⁇ t2 ⁇ ,
- G is a vector in which the term that does not depend on the state equation is stored.
- f is a matrix corresponding to the vector G.
- the inspection apparatus 500 mounted on the railway vehicle calculates the final alignment irregularity amount y R
- a data processing device in which some functions of the inspection apparatus 500 are mounted is disposed in an operation center.
- the data processing device receives measured data transmitted from the railway vehicle and calculates the final alignment irregularity amount y R by using the received measured data.
- the functions possessed by the inspection apparatus 500 in the first and second embodiments are shared and executed by the railway vehicle and the operation center. Constitutions and processing due to this are mainly different between the present embodiment and the first and second embodiments. Accordingly, in the explanation of the present embodiment, the same reference numerals and symbols as those added to Fig. 1 to Fig. 10 are added to the same parts as those in the first and second embodiments, or the like, and their detailed explanations are omitted. Note that the present embodiment can be applied to any of the first and second embodiments.
- Fig. 11 is a view illustrating one example of a configuration of an inspection system.
- the inspection system includes data collecting devices 1110a, 1110b, and a data processing device 1120.
- Fig. 11 one example of functional configurations of the data collecting devices 1110a, 1110b and the data processing device 1120 is also illustrated. Note that each hardware of the data collecting devices 1110a, 1110b and the data processing device 1120 can be fabricated by the one illustrated in Fig. 6 , for example. Accordingly, detailed explanations of the hardware configurations of the data collecting devices 1110a, 1110b and the data processing device 1120 are omitted.
- the data collecting devices 1110a, 1110b are mounted on each railway vehicle one by one.
- the data processing device 1120 is disposed at the operation center.
- the operation center centrally manages operations of a plurality of railway vehicles, for example.
- the data collecting devices 1110a, 1110b can be fabricated by the same components.
- the data collecting devices 1110a, 1110b include data acquisition units 1111a, 1111b and data transmission units 1112a, 1112b.
- the data acquisition units 1111a, 1111b have the same function as that of the data acquisition unit 502. That is, the data acquisition units 1111a, 1111b acquire the time-series data of the measured values of the forward-and-backward-direction forces, similarly to the data acquisition unit 502.
- the configuration for obtaining the measured value of the forward-and-backward-direction force is the same as that explained in the first embodiment.
- the data transmission units 1112a, 1112b transmit the time-series data of the measured values of the forward-and-backward-direction forces acquired by the data acquisition units 1111a, 1111b to the data processing device 1120.
- the data transmission units 1112a, 1112b transmit the time-series data of the measured values of the forward-and-backward-direction forces acquired by the data acquisition units 1111a, 1111b to the data processing device 1120 by radio.
- the data transmission units 1112a, 1112b add identification numbers of the railway vehicles on which the data collecting devices 1110a, 1110b are mounted to the time-series data of the measured values of the forward-and-backward-direction forces acquired by the data acquisition units 1111a, 1111b. In this manner, the data transmission units 1112a, 1112b transmit the time-series data of the measured values of the forward-and-backward-direction forces with the identification numbers of the railway vehicles added thereto.
- a data reception unit 1121 receives the time-series data of the measured values of the forward-and-backward-direction forces transmitted by the data transmission units 1112a, 1112b. To the time-series data of the measured values of the forward-and-backward-direction forces, the identification numbers of the railway vehicles, which are transmission sources of the time-series data of the measured values of the forward-and-backward-direction forces, have been added.
- a data storage unit 1122 stores the time-series data of the measured values of the forward-and-backward-direction forces received by the data reception unit 1121.
- the data storage unit 1122 stores the time-series data of the measured value of the forward-and-backward-direction force for each identification number of the railway vehicle.
- the data storage unit 1122 specifies the position of the railway vehicle at the time of receipt of the time-series data of the measured value of the forward-and-backward-direction force based on the current operation situation of the railway vehicle and the time of receipt of the time-series data of the measured value of the forward-and-backward-direction force, and stores information of the specified position and the time-series data of the measured value of the forward-and-backward-direction force in association with each other.
- the data collecting devices 1110a, 1110b may collect the information of the current positions of the railway vehicles and include the collected information in the time-series data of the measured values of the forward-and-backward-direction forces.
- a data reading unit 1123 reads the time-series data of the measured value of the forward-and-backward-direction force stored in the data storage unit 1122.
- the data reading unit 1123 can read, out of the time-series data of the measured values of the forward-and-backward-direction forces stored in the data storage unit 1122, data designated by an operator. Further, the data reading unit 1123 can also read the time-series data of the measured value of the forward-and-backward-direction force matching a preset condition at a preset timing.
- the time-series data of the measured value of the forward-and-backward-direction force read by the data reading unit 1123 is determined based on at least any one of the identification number and the position of the railway vehicle, for example.
- a storage unit 501, a first frequency adjustment unit 503, a state variable derivation unit 504, a second frequency adjustment unit 505, a track state derivation unit 506, and an output unit 507 are the same as those explained in the first embodiment. Accordingly, their detailed explanations are omitted here.
- the first frequency adjustment unit 503 uses the time-series data of the measured value of the forward-and-backward-direction force read by the data reading unit 1123 in place of using the time-series data of the measured value of the forward-and-backward-direction force acquired by the data acquisition unit 502, and generates the time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 to T 4 .
- the data collecting devices 1110a, 1110b mounted on the railway vehicles collect the time-series data of the measured values of the forward-and-backward-direction forces to transmit the data to the data processing device 1120.
- the data processing device 1120 disposed at the operation center stores the time-series data of the measured values of the forward-and-backward-direction forces received from the data collecting devices 1110a, 1110b and uses the stored time-series data of the measured values of the forward-and-backward-direction forces to calculate the final alignment irregularity amount y R . Accordingly, in addition to the effects explained in the first and second embodiments, for example, the following effects are exhibited.
- the data processing device 1120 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 1120 can output time-series variation of the final alignment irregularity amount y R at the same position. Further, the data processing device 1120 can output the final alignment irregularity amounts y R in a plurality of routes for each route.
- the storage unit 501, the data acquisition unit 502, the first frequency adjustment unit 503, the state variable derivation unit 504, the second frequency adjustment unit 505, the track state derivation unit 506, and the output unit 507 are included in one apparatus.
- Functions of the storage unit 501, the data acquisition unit 502, the first frequency adjustment unit 503, the state variable derivation unit 504, the second frequency adjustment unit 505, the track state derivation unit 506, and the output unit 507 may be fabricated by a plurality of apparatuses. In this case, the inspection system is constituted by using these plural apparatuses.
- derivation of a final alignment irregularity amount y R was performed based on a method in which as the values to be given as the measured values of the observation variables (the measured values of the accelerations of the vehicle body 11, the bogies 12a, 12b, and the wheel sets 13a to 13d in the right and left direction), not the preset fixed value but the measured values are given as they are (namely, the method described in Patent Literature 1).
- Fig. 12 illustrates the present calculation example, and is a view illustrating the curvature 1/R of the track 16 being the target of deriving the alignment irregularity amount and the traveling velocity v of the railway vehicle.
- a graph 1201 indicates the traveling velocity of the railway vehicle
- a graph 1202 indicates the curvature 1/R of the track 16.
- the horizontal axis in Fig. 12 indicates an elapsed time (second) from a reference time when the reference time is set to 0 (zero).
- FIG. 13A and Fig. 13B illustrates the present calculation example, and is a view illustrating a distribution of eigenvalues of an autocorrelation matrix R.
- Fig. 13A illustrates a distribution of eigenvalues of an autocorrelation matrix R with respect to the forward-and-backward-direction force T 1 in the wheel set 13a
- Fig. 13B illustrates a distribution of eigenvalues of an autocorrelation matrix R with respect to the forward-and-backward-direction force T 2 in the wheel set 13b.
- Fig. 14 illustrates the present calculation example, and is a view illustrating time-series data of the measured values y of the forward-and-backward-direction forces T 1 , T 2 and time-series data of predicted values y ⁇ k of the forward-and-backward-direction forces T 1 , T 2 (time-series data obtained by extracting low-frequency components included in the time-series data of the measured values y of the forward-and-backward-direction forces).
- the measured value indicates the time-series data of the measured value y of the forward-and-backward-direction force
- bias indicates the time-series data of the predicted value y ⁇ k of the forward-and-backward-direction force.
- the horizontal axis in Fig. 14 indicates a measuring time and a calculating time of the forward-and-backward-direction forces T 1 to T 4 , each of which is an elapsed time (second) from a reference time when the reference time is set to 0 (zero).
- Fig. 15 illustrates the present calculation example, and is a view illustrating the time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 , T 2 .
- the time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 , T 2 can be obtained by subtracting the time-series data of the predicted values y ⁇ k of the forward-and-backward-direction forces T 1 , T 2 from the time-series data of the measured values y of the forward-and-backward-direction forces T 1 , T 2 illustrated in Fig. 14 . Note that the horizontal axis in Fig.
- Fig. 16A and Fig. 16B are views illustrating the alignment irregularity amounts y R derived by the method of the first embodiment and the method described in Patent Literature 1 by using the time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 , T 2 illustrated in Fig. 15 .
- the calculated value indicates the alignment irregularity amount y R derived by the method described in Patent Literature 1
- the measured value indicates the measured value of the alignment irregularity amount y R .
- the calculated value indicates the alignment irregularity amount y R derived by the method of the first embodiment
- the measured value indicates the measured value of the alignment irregularity amount yR.
- the calculated value of the alignment irregularity amount y R an average value of the alignment irregularity amount y R at the position of the wheel set 13a and the alignment irregularity amount y R2 at the position of the wheel set 13b was used.
- the measured value illustrated in Fig. 16A and the measured value illustrated in Fig. 16B are the same.
- the horizontal axis in each of Fig. 16A and Fig. 16B indicates an elapsed time (second) from a reference time when the reference time is set to 0 (zero) and indicates a time corresponding to the position where the alignment irregularity amount y R exists.
- the illustration of data of a part at which a distance from a starting point of the railway vehicle is small is omitted for convenience of illustration.
- Fig. 17A is a view illustrating the alignment irregularity amount y R derived by the method of the second embodiment by using the time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 , T 2 illustrated in Fig. 15 .
- Fig. 17B is a view illustrating the alignment irregularity amount y R derived by the method described in Patent Literature 1 by using the time-series data of the high-frequency components of the forward-and-backward-direction forces T 1 , T 2 illustrated in Fig. 15 .
- the calculated value indicates the alignment irregularity amount y R derived by the method described in Patent Literature 1
- the measured value indicates the measured value of the alignment irregularity amount y R .
- Fig. 17A the calculated value indicates the alignment irregularity amount y R derived by the method described in Patent Literature 1
- the measured value indicates the measured value of the alignment irregularity amount y R .
- the calculated value indicates the alignment irregularity amount y R derived by the method of the second embodiment
- the measured value indicates the measured value of the alignment irregularity amount y R .
- the calculated value of the alignment irregularity amount y R an average value of the alignment irregularity amount y R1 at the position of the wheel set 13a and the alignment irregularity amount y R2 at the position of the wheel set 13b was used.
- the measured value illustrated in Fig. 17A and the measured value illustrated in Fig. 17B are the same (these measured values are also the same as the measured values illustrated in Fig. 16A and Fig. 16B ). Note that the horizontal axis in each of Fig. 17A and Fig.
- FIG. 17B indicates an elapsed time (second) from a reference time when the reference time is set to 0 (zero) and indicates a time corresponding to the position where the alignment irregularity amount y R exists. Further, in Fig. 17A and Fig. 17B , the illustration of data of a part at which a distance from a starting point of the railway vehicle is small is omitted for convenience of illustration.
- the alignment irregularity amount y R derived by the method of the first embodiment matches the alignment irregularity amount y R derived by the method described in Patent Literature 1 with good accuracy. Further, it can be understood that the calculated value and the measured value also match with good accuracy.
- the alignment irregularity amount y R derived by the method of the second embodiment matches the alignment irregularity amount y R derived by the method described in Patent Literature 1 with good accuracy. Further, it can be understood that the calculated value and the measured value also match with good accuracy.
- the embodiments of the present invention explained above can be fabricated by causing a computer to execute a program.
- 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.
- 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.
- Patent Literature 1 Note that the entire contents of the description and drawings of Patent Literature 1 can be incorporated herein by reference.
- the present invention can be utilized for inspecting railway vehicles, for example.
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Vehicle Body Suspensions (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Machines For Laying And Maintaining Railways (AREA)
Abstract
Description
- The present invention relates to an inspection system, an inspection method, and a program, and is suitably used for inspecting a track of a railway vehicle, in particular. The present application is based upon and claims the benefit of priority of the prior
Japanese Patent Application No. 2018-230834 - When a railway vehicle travels on a track, a position of the track changes due to a load from the railway vehicle. When such a change in track occurs, there is a possibility that the railway vehicle exhibits abnormal behavior. Accordingly,
Patent Literature 1 discloses a technique in which angular displacements of wheel sets in a yawing direction, state variables derived by a filter that performs data assimilation, and a measured value of a forward-and-backward-direction force being a force in a forward and backward direction that occurs in a member for supporting an axle box are substituted into motion equations that describe yawings of the wheel sets to derive track irregularity (alignment irregularity amount and the like) of a railway vehicle. - Patent Literature 1:
International Publication Pamphlet No. WO 2017/164133 - However, in the technique described in
Patent Literature 1, the measured value of the forward-and-backward-direction force and measured values of accelerations of the wheel set and a bogie (and a vehicle body in addition to the above, according to need), respectively, in a right and left direction are used when performing the data assimilation.
These measured values can be obtained without using special sensors, but the number of sensors disposed in the railway vehicle is preferably as small as possible. - The present invention has been made in view of the problems as described above, and an object thereof is to reduce the number of sensors which are used for detecting track irregularity of a railway vehicle.
- An inspection system of the present invention includes: a data acquisition means that acquires data of a measured value of a forward-and-backward-direction force as data of a measured value to be measured by causing a railway vehicle including a vehicle body, a bogie, and a wheel set to travel on a track; a state variable derivation means that derives state variables being variables to be determined in a state equation constituted by using motion equations that describe motions of the railway vehicle, by using the measured value of the forward-and-backward-direction force; and a track state derivation means that derives information reflecting a state of the track, in which the forward-and-backward-direction force is a force in a forward and backward direction that occurs in a member disposed between the wheel set and the bogie on which the wheel set is provided, and is a force to be determined according to a difference between an angular displacement of the wheel set in a yawing direction and an angular displacement of the bogie on which the wheel set is provided in the yawing direction, 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 yawing direction is a pivoting direction with an up and down direction being a direction vertical to the track set as a pivot axis, the state equation is an equation described by using the state variables, the forward-and-backward-direction force, and a transformation variable, the state variables include a displacement and a velocity of the bogie in a right and left direction, an angular displacement and an angular velocity of the bogie in the yawing direction, an angular displacement and an angular velocity of the bogie in a rolling direction, a displacement and a velocity of the wheel set in the right and left direction, and an angular displacement of an air spring attached to the railway vehicle in the rolling direction, and do not include an angular displacement and an angular velocity of the wheel set in the yawing direction, the rolling direction is a pivoting direction with the forward and backward direction set as a pivot axis, the transformation variable is a variable that performs mutual transformation between the angular displacement of the wheel set in the yawing direction and the angular displacement of the bogie in the yawing direction, the track state derivation means uses the angular displacement of the bogie in the yawing direction, which is one of the state variables derived by the state variable derivation means, and an actual value of the transformation variable to derive an estimated value of the angular displacement of the wheel set in the yawing direction, and uses the derived estimated value of the angular displacement of the wheel set in the yawing direction to derive the information reflecting the state of the track, the actual value of the transformation variable is derived by using the measured value of the forward-and-backward-direction force, and the state variable derivation means derives the state variables without using measured values of accelerations of the bogie, the wheel set, and the vehicle body in the right and left direction during a period in which the measured value of the forward-and-backward-direction force has been obtained.
- An inspection method of the present invention includes: a data acquisition step of acquiring data of a measured value of a forward-and-backward-direction force as data of a measured value to be measured by causing a railway vehicle including a vehicle body, a bogie, and a wheel set to travel on a track; a state variable derivation step of deriving state variables being variables to be determined in a state equation constituted by using motion equations that describe motions of the railway vehicle, by using the measured value of the forward-and-backward-direction force; and a track state derivation step of deriving information reflecting a state of the track, in which the forward-and-backward-direction force is a force in a forward and backward direction that occurs in a member disposed between the wheel set and the bogie on which the wheel set is provided, and is a force to be determined according to a difference between an angular displacement of the wheel set in a yawing direction and an angular displacement of the bogie on which the wheel set is provided in the yawing direction, 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 yawing direction is a pivoting direction with an up and down direction being a direction vertical to the track set as a pivot axis, the state equation is an equation described by using the state variables, the forward-and-backward-direction force, and a transformation variable, the state variables include a displacement and a velocity of the bogie in a right and left direction, an angular displacement and an angular velocity of the bogie in the yawing direction, an angular displacement and an angular velocity of the bogie in a rolling direction, a displacement and a velocity of the wheel set in the right and left direction, and an angular displacement of an air spring attached to the railway vehicle in the rolling direction, and do not include an angular displacement and an angular velocity of the wheel set in the yawing direction, the rolling direction is a pivoting direction with the forward and backward direction set as a pivot axis, the transformation variable is a variable that performs mutual transformation between the angular displacement of the wheel set in the yawing direction and the angular displacement of the bogie in the yawing direction, the track state derivation step uses the angular displacement of the bogie in the yawing direction, which is one of the state variables derived by the state variable derivation step, and an actual value of the transformation variable to derive an estimated value of the angular displacement of the wheel set in the yawing direction, and uses the derived estimated value of the angular displacement of the wheel set in the yawing direction to derive the information reflecting the state of the track, the actual value of the transformation variable is derived by using the measured value of the forward-and-backward-direction force, and the state variable derivation step derives the state variables without using measured values of accelerations of the bogie, the wheel set, and the vehicle body in the right and left direction during a period in which the measured value of the forward-and-backward-direction force has been obtained.
- A program of the present invention causes a computer to execute: a data acquisition step of acquiring data of a measured value of a forward-and-backward-direction force as data of a measured value to be measured by causing a railway vehicle including a vehicle body, a bogie, and a wheel set to travel on a track; a state variable derivation step of deriving state variables being variables to be determined in a state equation constituted by using motion equations that describe motions of the railway vehicle, by using the measured value of the forward-and-backward-direction force; and a track state derivation step of deriving information reflecting a state of the track, in which the forward-and-backward-direction force is a force in a forward and backward direction that occurs in a member disposed between the wheel set and the bogie on which the wheel set is provided, and is a force to be determined according to a difference between an angular displacement of the wheel set in a yawing direction and an angular displacement of the bogie on which the wheel set is provided in the yawing direction, 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 yawing direction is a pivoting direction with an up and down direction being a direction vertical to the track set as a pivot axis, the state equation is an equation described by using the state variables, the forward-and-backward-direction force, and a transformation variable, the state variables include a displacement and a velocity of the bogie in a right and left direction, an angular displacement and an angular velocity of the bogie in the yawing direction, an angular displacement and an angular velocity of the bogie in a rolling direction, a displacement and a velocity of the wheel set in the right and left direction, and an angular displacement of an air spring attached to the railway vehicle in the rolling direction, and do not include an angular displacement and an angular velocity of the wheel set in the yawing direction, the rolling direction is a pivoting direction with the forward and backward direction set as a pivot axis, the transformation variable is a variable that performs mutual transformation between the angular displacement of the wheel set in the yawing direction and the angular displacement of the bogie in the yawing direction, the track state derivation step uses the angular displacement of the bogie in the yawing direction, which is one of the state variables derived by the state variable derivation step, and an actual value of the transformation variable to derive an estimated value of the angular displacement of the wheel set in the yawing direction, and uses the derived estimated value of the angular displacement of the wheel set in the yawing direction to derive the information reflecting the state of the track, the actual value of the transformation variable is derived by using the measured value of the forward-and-backward-direction force, and the state variable derivation step derives the state variables without using measured values of accelerations of the bogie, the wheel set, and the vehicle body in the right and left direction during a period in which the measured value of the forward-and-backward-direction force has been obtained.
-
- [
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. 3] Fig. 3 is a view illustrating a measured value and a calculated value of each of acceleration of a bogie in a right and left direction and accelerations of wheel sets in the right and left direction. - [
Fig. 4A] Fig. 4A is a view illustrating one example of an alignment irregularity amount in a linear track. - [
Fig. 4B] Fig. 4B is a view illustrating one example of an alignment irregularity amount in a curved track. - [
Fig. 5] Fig. 5 is a view illustrating one example of a functional configuration of an inspection apparatus. - [
Fig. 6] Fig. 6 is a view illustrating one example of a hardware configuration of the inspection apparatus. - [
Fig. 7] Fig. 7 is a flowchart illustrating one example of processing in the inspection apparatus. - [
Fig. 8] Fig. 8 is a view illustrating one example of a distribution of eigenvalues of an autocorrelation matrix. - [
Fig. 9] Fig. 9 is a view illustrating one example of time-series data of measured values of forward-and-backward-direction forces (measured values) and time-series data of predicted values of the forward-and-backward-direction forces (calculated values). - [
Fig. 10] Fig. 10 is a view illustrating one example of time-series data of high-frequency components of the forward-and-backward-direction forces. - [
Fig. 11] Fig. 11 is a view illustrating one example of a configuration of an inspection system. - [
Fig. 12] Fig. 12 is a view illustrating a calculation example, and is a view illustrating acurvature 1/R of a track being a target of deriving an alignment irregularity amount and a traveling velocity of the railway vehicle. - [
Fig. 13A] Fig. 13A is a view illustrating a calculation example, and is a view illustrating a first example of a distribution of eigenvalues of an autocorrelation matrix R. - [
Fig. 13B] Fig. 13B is a view illustrating a calculation example, and is a view illustrating a second example of a distribution of eigenvalues of the autocorrelation matrix R. - [
Fig. 14] Fig. 14 is a view illustrating a calculation example, and is a view illustrating time-series data of measured values of forward-and-backward-direction forces and time-series data of predicted values of the forward-and-backward-direction forces. - [
Fig. 15] Fig. 15 is a view illustrating a calculation example, and is a view illustrating time-series data of high-frequency components of the forward-and-backward-direction forces. - [
Fig. 16A] Fig. 16A is a view illustrating a calculation example based on a method of a first embodiment, and is a view illustrating a first example of an alignment irregularity amount yR. - [
Fig. 16B] Fig. 16B is a view illustrating a calculation example based on the method of the first embodiment, and is a view illustrating a second example of the alignment irregularity amount yR. - [
Fig. 17A] Fig. 17A is a view illustrating a calculation example based on a method of a second embodiment, and is a view illustrating a first example of the alignment irregularity amount yR. - [
Fig. 17B] Fig. 17B is a view illustrating a calculation example based on the method of the second embodiment, and is a view illustrating the first example of the alignment irregularity amount yR. - Hereinafter, embodiments of the present invention will be explained while referring to the drawings.
- First, there will be explained a conception obtained by the present inventors when realizing embodiments of the present invention.
-
Fig. 1 is a view illustrating one example of an outline of a railway vehicle. Note that inFig. 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 vertical to a track 16 (the ground) (a height direction of the railway vehicle). The y axis is set to a horizontal direction vertical to the traveling direction of the railway vehicle (a direction vertical to both the traveling direction and the height direction of the railway vehicle). Further, the railway vehicle is set to a commercial vehicle. Note that in the respective drawings, the mark of ● added inside ○ 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. - As illustrated in
Fig. 1 , in the present embodiment, the railway vehicle includes avehicle body 11,bogies wheel sets 13a to 13d. As described above, in the present embodiment, the railway vehicle including thesingle vehicle body 11 provided with the twobogies axles 15a to 15d andwheels 14a to 14d provided on both ends of theaxles 15a to 15d respectively. In the present embodiment, a case where each of thebogies Fig. 1 , for convenience of illustration, only thewheels 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 inFig. 1 , there are eight wheels in total). Further, the railway vehicle includes components other than the components illustrated inFig. 1 (components and so on to be explained in later-described motion equations), but for convenience of illustration, illustrations of these components are omitted inFig. 1 . For example, thebogies 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. Note that the railway vehicle itself can be fabricated by a well-known technique, and thus its detailed explanation is omitted here. - When the railway vehicle travels on the
track 16, acting force (creep force) between thewheels 14a to 14d and thetrack 16 becomes a vibration source and the vibration sequentially propagates to the wheel sets 13a to 13d, thebogies vehicle body 11.Fig. 2 is a view conceptually illustrating directions of main motions of the components (the wheel sets 13a to 13d, thebogies Fig. 2 correspond to the x axis, the y axis, and the z axis illustrated inFig. 1 respectively. - As illustrated in
Fig. 2 , in the present embodiment, a case where the wheel sets 13a to 13d, thebogies 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, will be explained as an example. 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. Note that the forward and backward direction is the traveling direction of the railway vehicle. In the present embodiment, the direction along the x axis is set to 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. Note that the up and down direction is a direction vertical to thetrack 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. Note that the right and left direction is a direction vertical to both the forward and backward direction (the traveling direction of the railway vehicle) and the up and down direction (the direction vertical 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. - In the technique described in
Patent Literature 1, the state variables are derived by performing filtering with a filter (Kalman filter) that performs data assimilation by using, as observation variables, accelerations yw1 • •, yw2 • •, yw3 • •, yw4 • • of thewheel sets bogies vehicle body 11 in the right and left direction according to need. -
Fig. 3 illustrates a measured value and a calculated value of each of the acceleration yt1 • • of thebogie 12a in the right and left direction and the accelerations yw1 • •, yw2 • • of thewheel sets Fig. 3 indicates an elapsed time (second) from a reference time when the reference time is set to 0 (zero). Concretely, the horizontal axis inFig. 3 indicates a measuring time and a calculating time of the acceleration yt1 • • of thebogie 12a in the right and left direction and the accelerations yw1 • •, yw2 • • of thewheel sets - As illustrated in
Fig. 3 , the measured values of the accelerations yw1 • •, yw2• • of thewheel sets bogie 12a in the right and left direction include a lot of noises. Based on this, the present inventors obtained a finding that, depending on a state of the track 16 (rail), the estimated value does not approximate the measured value and becomes substantially a fixed value even if the data assimilation is performed. The same applies to the accelerations yw3 • •, yw4 • • of the wheel sets 13c, 13d in the right and left direction, the acceleration yt2 • • of thebogie 12b in the right and left direction, and the acceleration yb • • of thevehicle body 11 in the right and left direction. From the above, the present inventors considered that there may be a chance that, when deriving the state variables, the alignment irregularity amounts yR1, yR2, yR3, yR4 can be derived without greatly decreasing accuracy, even if the measured values of these accelerations are not used. - As described in
Patent Literature 1, the present inventors devised a method of calculating an alignment irregularity amount by using a measured value of force in the forward and backward direction that occurs in a member disposed between thewheel sets bogie 12a (12b) on which thesewheel sets - The alignment irregularity amount is calculated by using an equation representing a 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 thetrack 16 is referred to as a linear track as necessary and the curved portion of thetrack 16 is referred to as a curved track as necessary. - 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 data assimilation, state variables may diverge. Therefore, the state equation in the case of performing filtering with the filter (Kalman filter) that performs data assimilation is constituted by using a motion equation that describes motion of the railway vehicle traveling on the linear track.
- 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 including a curvature radius of the rail. Therefore, when the state variables are derived by 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.
- 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. The component itself of the forward-and-backward-direction force due to the alignment irregularity occurs in the same manner even on the curved track or the linear track. Accordingly, the present inventors thought that the alignment irregularity amount itself has nothing to do with the amount of the aforementioned bias, and 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 if the filter (Kalman filter) performing the 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 using the time-series data of the value of the forward-and-backward-direction force from which the low-frequency component has been reduced. 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. Note that there is a chance that even a track, which is designed as a linear track, actually has a curvature which may exert an influence on estimation accuracy of the alignment irregularity amount. Therefore, the reduction in the low-frequency component (behavior of the aforementioned bias) from the time-series data of the measured value of the forward-and-backward-direction force contributes to improvement of the estimation accuracy of the alignment irregularity amount on not only the curved track but also the linear track. Hereinafter, explanation will be made by setting that even a track, which is designed as a linear track but actually has a curvature which may exert an influence on estimation accuracy of the alignment irregularity amount, is also regarded as a curved track.
- Next, there will be explained one example of a motion equation that describes the motion of the railway vehicle. In the present embodiment, there will be explained, as an example, a 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 thebogies 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 (bolster springs) each provided on thebogies bogies - Note that 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 by taking these points into consideration. 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, thebogies Non-Patent Literatures - 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, w4 indicate the
wheel sets - Subscripts t, T indicate the
bogies bogies bogies - Subscripts b, B indicate the
vehicle body 11. - A subscript x indicates the forward and backward direction or the rolling direction, a subscript y indicates the right and left direction, and a subscript z indicates the up and down direction or the yawing direction.
- Further, " • • " and "•" each added above a variable indicate a second-order time differential and a first-order time differential respectively.
- Note that 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. -
- mw is the mass of the wheel sets 13a to 13d. yw1 • • is acceleration of the
wheel set 13a in the right and left direction (in the equation, • • is added above yw1 (the same applies to the other variables below)). f2 is a lateral creep coefficient (note that the lateral creep coefficient f2 may be given for each of the wheel sets 13a to 13d). v is a traveling velocity of the railway vehicle. yw1 • is a velocity of thewheel set 13a in the right and left direction (in the equation, • is added above yw1 (the same applies to the other variables below)). Cwy is a damping constant of the axle box suspension coupling the axle box and the wheel set in the right and left direction. yt1 • is a velocity of thebogie 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 thebogies bogies bogie 12a in the yawing direction. h1 is a distance between the center of the axle and the center of gravity of thebogie 12a in the up and down direction. φ t1 • is an angular velocity of thebogie 12a in the rolling direction. φ w1 is a pivot amount (angular displacement) of thewheel set 13a in the yawing direction. Kwy is a spring constant of the axle box suspension in the right and left direction. yw1 is a displacement of thewheel set 13a in the right and left direction. yt1 is a displacement of thebogie 12a in the right and left direction. φ t1 is a pivot amount (angular displacement) of thebogie 12a in the yawing direction. φ t1 is a pivot amount (angular displacement) of thebogie 12a in the rolling direction. Note that 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 subscripts described above. -
- Iwz 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. f1 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 thewheel set 13a in the yawing direction. Cwx is a damping constant of the axle box suspension in the forward and backward direction. b1 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 2b1). γ is a tread slope. r is a radius of thewheels 14a to 14d. yR1 is an alignment irregularity amount at the position of thewheel set 13a. sa is an offset from the center of theaxles 15a to 15d to an axle box suspension spring in the forward and backward direction. yt1 is a displacement of thebogie 12a in the right and left direction. Kwx is a spring constant of the axle box suspension in the forward and backward direction. Note that 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 subscripts described above. However, yR2, yR3, yR4 are alignment irregularity amounts at the positions of the wheel sets 13b, 13c, 13d respectively. - 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. 4A and Fig. 4B each illustrate one example of the alignment irregularity amount yR1 at the position of thewheel set 13a. InFig. 4A , the case of thetrack 16 being the linear track will be explained as an example. InFig. 4B , the case of thetrack 16 being the curved track will be explained as an example. InFig. 4A and Fig. 4B ,16a denotes a rail and 16b denotes a crosstie. InFig. 4A , it is set that thewheel 14a of thewheel set 13a is in contact with therail 16a at aposition 401. InFig. 4B , it is set that thewheel 14a of thewheel set 13a is in contact with therail 16a at aposition 402. The alignment irregularity amount yR1 at the position of thewheel set 13a is a distance in the right and left direction between the contact position between thewheel 14a of thewheel set 13a and therail 16a and the position of therail 16a in the case where this position is assumed as a regular state. The position of thewheel set 13a is the contact position between thewheel 14a of thewheel set 13a and therail 16a. The alignment irregularity amounts yR2, yR3, yR4 at the positions of the wheel sets 13b, 13c, 13d are also defined in the same manner as the alignment irregularity amount yR1 at the position of thewheel set 13a. -
- mT is the mass of the
bogies bogie 12a in the right and left direction. c'2 is a damping constant of a lateral movement damper. h4 is a distance between the center of gravity of thebogie 12a and the lateral movement damper in the up and down direction. yb • is a velocity of thevehicle body 11 in the right and left direction. L represents 1/2 of the interval between the center of thebogie 12a and the center of thebogie 12b in the forward and backward direction (the interval between the center of thebogie 12a and the center of thebogie 12b in the forward and backward direction becomes 2L). φ b • is an angular velocity of thevehicle body 11 in the yawing direction. h5 is a distance between the lateral movement damper and the center of gravity of thevehicle body 11 in the up and down direction. φ b • is an angular velocity of thevehicle body 11 in the rolling direction. yw2 • is a velocity of thewheel set 13b in the right and left direction. k'2 is a spring constant of the air spring (bolster spring) in the right and left direction. h2 is a distance between the center of gravity of each of thebogies vehicle body 11 in the right and left direction. φ b is a pivot amount (angular displacement) of thevehicle body 11 in the yawing direction. h3 is a distance between the center of the air spring (bolster spring) and the center of gravity of thevehicle body 11 in the up and down direction. φ b is a pivot amount (angular displacement) of thevehicle body 11 in the rolling direction. Note that respective variables in (10) Equation are represented by being replaced with the variables in (9) Equation according to the meanings of the subscripts described above. -
- ITz is a moment of inertia of the
bogies bogie 12a in the yawing direction. φ w2 • is an angular velocity of thewheel set 13b in the yawing direction. φ w2 is a pivot amount (angular displacement) of thewheel set 13b in the yawing direction. yw2 is a displacement of thewheel set 13b in the right and left direction. k'0 is stiffness of a rubber bush of the yaw damper. b'0 represents 1/2 of the interval between the two yaw dampers, which are disposed on the right and left sides of each of thebogies bogies bogie 12a in the yawing direction. k"2 is a spring constant of the air spring (bolster spring) in the forward and backward direction. b2 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 thebogies bogies -
- ITx is a moment of inertia of the
bogies bogie 12a in the rolling direction. c1 is a damping constant of an axle damper in the up and down direction. b'1 represents 1/2 of the interval between the two axle dampers, which are disposed on the right and left sides of each of thebogies bogies bogie 12a in the rolling direction. k1 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. k2 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 thebogie 12a in the rolling direction. k3 is equivalent stiffness by a change in effective pressure receiving area of the air spring (bolster spring). Note that respective variables in (14) Equation are represented by being replaced with the variables in (13) Equation according to the meanings of the subscripts described above. Note that φ a2 is a pivot amount (angular displacement) of the air spring (bolster spring) disposed on thebogie 12b in the rolling direction. -
- mB is the mass of the
bogies vehicle body 11 in the right and left direction. yt2 • is a velocity of thebogie 12b in the right and left direction. φ t2 • is an angular velocity of thebogie 12b in the rolling direction. yt2 is a displacement of thebogie 12b in the right and left direction. φ t2 is a pivot amount (angular displacement) of thebogie 12b in the rolling direction. -
- IBz is a moment of inertia of the
vehicle body 11 in the yawing direction. φ b • • is angular acceleration of thevehicle body 11 in the yawing direction. c0 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 thebogie 12a in the yawing direction. φ y2 • is an angular velocity of the yaw damper disposed on thebogie 12b in the yawing direction. φ t2 is a pivot amount (angular displacement) of thebogie 12b in the yawing direction. -
- IBx is a moment of inertia of the
vehicle body 11 in the rolling direction. φ b • • is angular acceleration of thevehicle body 11 in the rolling direction. -
- φ y2 is a pivot amount (angular displacement) of the yaw damper disposed on the
bogie 12b in the yawing direction. -
- φ a2 • is an angular velocity of the air spring (bolster spring) disposed on the
bogie 12b in the rolling direction. - Next, the forward-and-backward-direction force will be explained. Note that the forward-and-backward-direction force itself is the same as that described in
Patent Literature 1. - 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 that occur in the aforementioned two members attached to both the right and left sides of one wheel set.
- 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.
- 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.
- Note that the load measured by the strain gauge attached to the link sometimes includes 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. 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. Therefore, 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.
- 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.
- Note that the load measured by the strain gauge attached to the axle beam sometimes includes 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. 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, even if the strain gauge is not 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 thus 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.
- 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.
- Note that as the previously-described displacement meter, a well-known laser displacement meter or eddy current displacement meter can be used.
- 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.
- 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 fourwheel sets 13a to 13d. Accordingly, it is possible to obtain measured values of four forward-and-backward-direction forces T1 to T4. - Next, a first embodiment of the present invention will be described.
-
Fig. 5 is a view illustrating one example of a functional configuration of aninspection apparatus 500.Fig. 6 is a view illustrating one example of a hardware configuration of theinspection apparatus 500.Fig. 7 is a flowchart illustrating one example of processing in theinspection apparatus 500. In the present embodiment, as illustrated inFig. 1 , the case where theinspection apparatus 500 is mounted on the railway vehicle will be explained as an example. - In
Fig. 5 , theinspection apparatus 500 includes, as its functions, astorage unit 501, adata acquisition unit 502, a firstfrequency adjustment unit 503, a statevariable derivation unit 504, a secondfrequency adjustment unit 505, a trackstate derivation unit 506, and anoutput unit 507. - In
Fig. 6 , theinspection apparatus 500 includes aCPU 601, amain memory 602, anauxiliary memory 603, acommunication circuit 604, asignal processing circuit 605, animage processing circuit 606, an I/F circuit 607, auser interface 608, adisplay 609, and abus 610. - The
CPU 601 overall controls theentire inspection apparatus 500. TheCPU 601 uses themain memory 602 as a work area to execute a program stored in theauxiliary memory 603. Themain memory 602 stores data temporarily. Theauxiliary memory 603 stores various kinds of data, in addition to programs to be executed by theCPU 601. Theauxiliary memory 603 stores later-described state equations and observation equations. Thestorage unit 501 is fabricated by using theCPU 601 and theauxiliary memory 603, for example. - The
communication circuit 604 is a circuit intended for performing communication with the outside of theinspection apparatus 500. Thecommunication circuit 604 receives information of the measured value of the forward-and-backward-direction force, for example. Thecommunication circuit 604 may perform radio communication or wire communication with the outside of theinspection apparatus 500.
Thecommunication circuit 604 is connected to an antenna provided on the railway vehicle in the case of performing radio communication. - The
signal processing circuit 605 performs various kinds of signal processing on signals received by thecommunication circuit 604 and signals input according to the control made by theCPU 601. Thedata acquisition unit 502 is fabricated by using theCPU 601, thecommunication circuit 604, and thesignal processing circuit 605, for example. Further, the firstfrequency adjustment unit 503, the statevariable derivation unit 504, the secondfrequency adjustment unit 505, and the trackstate derivation unit 506 are fabricated by using theCPU 601 and thesignal processing circuit 605, for example. - The
image processing circuit 606 performs various kinds of image processing on signals input according to the control made by theCPU 601. The signal after being subjected to the image processing is output to thedisplay 609. - The
user interface 608 is a part through which an operator gives an instruction to theinspection apparatus 500. Theuser interface 608 includes buttons, switches, dials, and so on, for example. Further, theuser interface 608 may include a graphical user interface using thedisplay 609. - The
display 609 displays an image based on a signal output from theimage processing circuit 606. The I/F circuit 607 exchanges data with a device connected to the I/F circuit 607. InFig. 6 , as the device to be connected to the I/F circuit 607, theuser interface 608 and thedisplay 609 are illustrated. However, the device to be connected to the I/F circuit 607 is not limited to these. For example, a portable storage medium may be connected to the I/F circuit 607. Further, at least a part of theuser interface 608 and thedisplay 609 may be provided outside theinspection apparatus 500. - The
output unit 507 is fabricated by using at least any one of a set including thecommunication circuit 604, and thesignal processing circuit 605, and a set including theimage processing circuit 606, the I/F circuit 607, and thedisplay 609, for example. - Note that the
CPU 601, themain memory 602, theauxiliary memory 603, thesignal processing circuit 605, theimage processing circuit 606, and the I/F circuit 607 are connected to thebus 610. Communication among these components is performed via thebus 610. Further, the hardware of theinspection apparatus 500 is not limited to the one illustrated inFig. 6 as long as it can realize later-described functions of theinspection apparatus 500. - The
storage unit 501 stores equations which are used when the later-described statevariable derivation unit 504 derives the state variables. - In the present embodiment, the
storage unit 501 stores state equations and observation equations. - In the present embodiment, the case of using the state equations and the observation equations described in
Patent Literature 1 will be explained as an example. - First, the state equation will be described.
- In the present embodiment, (5) Equation to (8) Equation (the motion equations that describe the yawings of the wheel sets 13a to 13d) are not included in the state equation, and the state equation is constituted as follows.
- First, (9) Equation and (10) Equation (the motion equations that describe the transversal vibrations of the
bogies bogies bogie 12a and the yaw damper disposed on thebogie 12b), and (20) Equation and (21) Equation (the motion equations that describe the rollings of the air spring (bolster spring) disposed on thebogie 12a and the air spring (bolster spring) disposed on thebogie 12b) are used as they are to constitute the state equation. - Meanwhile, in (1) Equation to (4) Equation (the motion equations that describe the transversal vibrations of the wheel sets 13a to 13d (motion in the right and left direction)) and (11) Equation and (12) Equation (the motion equations that describe the yawings of the
bogies - First, the forward-and-backward-direction forces T1 to T4 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 T1 to T4 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 and φ t2 of the bogies on which these wheel sets are provided in the yawing direction.
[Mathematical equation 11] - Transformation variables e1 to e4 are defined as in (26) Equation to (29) Equation below. As above, the transformation variables e1 to e4 are defined by the differences between the angular displacements φ t1 and φ 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 e1 to e4 are variables for performing mutual transformation between the angular displacements φ t1 and φ 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] -
-
- As above, (1) Equation to (4) 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 by using the transformation variables e1 to e4, 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.
-
- As above, (11) Equation and (12) Equation (the motion equations that describe the yawings of the
bogies -
- As above, in the present 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 frequency adjustment unit 503. - The actual values of the transformation variables e1 to e4 derived as above are given to (34) Equation to (37) Equation. Further, the values of the forward-and-backward-direction forces T1 to T4 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 T1 to T4 are obtained by reducing a signal strength of a low-frequency component to be generated due to the railway vehicle traveling on the curved portion of the track from the time-series data of the measured value of the forward-and-backward-direction force, with the use of the later-described first
frequency adjustment unit 503. -
- The
storage unit 501 receives the state equation constituted as above, for example, based on the operation of theuser interface 608 made by an operator and stores it. - Next, the observation equation will be described.
- In the present embodiment, the acceleration of the
vehicle body 11 in the right and left direction, the accelerations of thebogies - The
storage unit 501 receives the observation equation constituted in this manner, for example, based on the operation of theuser interface 608 made by an operator and stores it. - After the state equation and the observation equation are stored in the
inspection apparatus 500 as above, thedata acquisition unit 502, the firstfrequency adjustment unit 503, the statevariable derivation unit 504, the secondfrequency adjustment unit 505, the trackstate derivation unit 506, and theoutput unit 507 start. Specifically, the processing according to the flowchart inFig. 7 starts after the state equation and the observation equation are stored in theinspection apparatus 500. - The
data acquisition unit 502 acquires 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. Thedata acquisition unit 502 can acquire the time-series data of the measured value of the forward-and-backward-direction force by performing communication with an arithmetic device that calculates the forward-and-backward-direction force by using a measured value of a strain gauge for measuring the forward-and-backward-direction force, for example. Note that thedata acquisition unit 502 does not acquire time-series data of a measured value of the acceleration of thevehicle body 11 in the right and left direction, time-series data of measured values of the accelerations of thebogies - The first
frequency adjustment unit 503 reduces (preferably removes) the signal strength of the low-frequency component included in the time-series data of the measured value of the forward-and-backward-direction force acquired by thedata acquisition unit 502. A signal of this low-frequency component is a signal that is not measured when the railway vehicle is traveling on the linear track whose curvature is 0 (zero), 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 whose curvature is 0 (zero). - The present inventors devised a model in which an AR (Auto-regressive) model is corrected. Further, the present inventors came up with an idea of reducing the signal strength of the low-frequency component included 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.
- A value of time-series data of a physical quantity y at a time k (1 ≦ k ≦ M) is set to yk. M is a number indicating, as the time-series data of the physical quantity y, data until when is 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 yk 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 yk-l of the physical quantity at a time k-l (1 ≦ l ≦ m) prior to the time k in the data y. Note that y^k is expressed by adding ^ above yk in (45) Equation.
[Mathematical equation 18] - 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 yk of the data y at the time k in the AR model, and is a number among values yk-1 to yk-m of the data y at continuous times k-1 to k-m prior to the time k. m is an integer of less than M. As m, for example, 1500 can be used.
- Subsequently, 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 yk 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 yk, 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 yk, 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 yk. (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 yk.
[Mathematical equation 19] -
-
-
- 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 yk 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] -
- In general, when deriving the coefficient of the AR model, a method of solving (50) Equation regarding a coefficient α 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 yk 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 included in the value yk of the data y at each time.
- Accordingly, 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 of this, the present inventors found out that it is possible to reduce the influence of a high-frequency component included 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.
- There will be explained a concrete example of the above below.
- The autocorrelation matrix R is subjected to singular value decomposition. Elements of the autocorrelation matrix R are symmetric. Therefore, 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] - 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 σ11, σ22, · · ·, σ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 orthogonal matrix U is set to u1, u2, ···, um. There is a correspondence relation in which the eigenvalue of the autocorrelation matrix R responsive to an eigenvector uj 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] - The values of σ11, σ22, ···, σ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 the present 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] - A matrix Us in (54) Equation is a matrix of m ×s composed of s pieces of the column component vectors, which are chosen from the left of the orthogonal matrix U of (52) Equation (eigenvectors corresponding to the eigenvalues to be used). That is, the matrix Us is a submatrix composed of the left elements of m×s cut out from the orthogonal matrix U. Further, Us T in (54) Equation is a transposed matrix of Us. Us T is a matrix of s×m composed of s pieces of row component vectors, which are chosen from the top of the matrix UT 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 Σ.
-
-
- (56) Equation is modified, and thereby (57) Equation below is obtained as the equation for 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] - The case where the values of σ11, σ22, ···, σ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 that case, the matrix Us 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.
- (57) Equation is an equation to be used for determining the coefficient of the corrected AR model. The matrix Us 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 UsΣsUs T in (57) Equation is a matrix (a first matrix) derived from the matrix Σs and the matrix Us.
- 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 that 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. Therefore, the autocorrelation of the data y in the present 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, R22 to Rmm are autocorrelation with a time lag of 0 (zero), but they may be replaced with R11.
- 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.
- 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 and the like of the railway vehicle.
- Accordingly, 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. 8 is a view illustrating one example of the distribution of the eigenvalues of the autocorrelation matrix R. InFig. 8 , eigenvalues σ11 to σmm, which are obtained by the autocorrelation matrix R in each of the time-series data of the measured value y of the forward-and-backward-direction force T1 in thewheel set 13a being subjected to singular value decomposition, are aligned in ascending order and are plotted. InFig. 8 , the horizontal axis is an index of the eigenvalue and the vertical axis is the value of the eigenvalue. - In the example illustrated in
Fig. 8 , there is one eigenvalue having significantly high value when compared to the other values. Further, there are two eigenvalues whose values are not so high when compared to the aforementioned eigenvalue having significantly high value, but are relatively large when compared to the other values and thus cannot be regarded as 0 (zero). Based on this, it is possible to employ two or three, for example, as the number s of eigenvalues to be extracted from the autocorrelation matrix R illustrated in (53) Equation. Even if either of them is employed, there is generated no significant difference in the result. Note that the number of eigenvalue having a value which is significantly higher than the other values may change depending on the configuration of the railway vehicle, the configuration of the track, and so on. Therefore, the number s of eigenvalues to be extracted from the autocorrelation matrix R is not limited to these values as long as it is one or more. - The first
frequency adjustment unit 503 performs the following processing every time thedata acquisition unit 502 acquires the value yk of the time-series data of the measured value y of the forward-and-backward-direction force at the time k at a predetermined sampling period. - First, the first
frequency adjustment unit 503 generates the autocorrelation matrix R by using (49) Equation and (51) Equation based on the time-series data of the measured value y of the forward-and-backward-direction force and preset numbers M, m. - Next, the first
frequency adjustment unit 503 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 σ11 to σmm of the autocorrelation matrix R from the diagonal matrix Σ. - Next, the first
frequency adjustment unit 503 chooses s pieces of the eigenvalues σ11 to σss from the largest from among the plural eigenvalues σ11 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. - Next, the first
frequency adjustment unit 503 determines the coefficient α of the corrected AR model by using (57) Equation based on the time-series data of the measured value y of the forward-and-backward-direction force, the eigenvalues σ11 to σss, and the orthogonal matrix U obtained by the singular value decomposition of the autocorrelation matrix R. - Subsequently, the first
frequency adjustment unit 503 derives the predicted value y^k of the time-series data of the measured value y 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 time-series data of the measured value y of the forward-and-backward-direction force. Time-series data of the predicted value y^k of the forward-and-backward-direction force results in the time-series data from which the low-frequency component included in the time-series data of the measured value y of the forward-and-backward-direction force has been extracted. -
Fig. 9 is a view illustrating one example of time-series data of measured values of forward-and-backward-direction forces (measured values) and time-series data of predicted values of the forward-and-backward-direction forces (calculated values). Note that in the present embodiment, the measured values of the four forward-and-backward-direction forces T1 to T4 are obtained. That is, four pieces of the data y of the forward-and-backward-direction force are obtained. InFig. 9 , the measured value and the calculated value of each of the four pieces of data y are illustrated. The horizontal axis inFig. 9 indicates a measuring time and a calculating time of the forward-and-backward-direction forces T1 to T4, 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 T1 to T4 (Nm). - In
Fig. 9 , the calculated value of the forward-and-backward-direction force T1 in thewheel set 13a is biased at about 15 seconds to 35 seconds. Specifically, the calculated value of the forward-and-backward-direction force T1 in thewheel set 13a exhibits, at about 15 seconds to 35 seconds, a value larger than that at another time. This period corresponds to the period when thewheel set 13a passes through the curved track. The calculated value of the forward-and-backward-direction force T2 in thewheel set 13b, the calculated value of the forward-and-backward-direction force T3 in the wheel set 13c, and the calculated value of the forward-and-backward-direction force T4 in thewheel set 13d are also biased during the period when the wheel sets 13b, 13c, 13d pass through the curved track, similarly to the calculated value of the forward-and-backward-direction force T1 in thewheel set 13a. - Accordingly, in
Fig. 9 , removal of the calculated values from the measured values of the forward-and-backward-direction forces T1 to T4 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 T1 to T4. That is, inFig. 9 , when the calculated values are removed from the measured values of the forward-and-backward-direction forces T1 to T4 in the wheel sets 13a to 13d, as the forward-and-backward-direction forces T1 to T4 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. - Accordingly, the first
frequency adjustment unit 503 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 yk 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 yk of the forward-and-backward-direction force is 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. -
Fig. 10 is a view illustrating one example of the time-series data of the high-frequency components of the forward-and-backward-direction forces. The vertical axis inFig. 10 indicates the time-series data of the high-frequency components of the forward-and-backward-direction forces T1, T2, T3, T4. That is, the high-frequency components of the forward-and-backward-direction forces T1, T2, T3, T4 illustrated on the vertical axis inFig. 10 are ones obtained by subtracting the calculated values from the measured values of the forward-and-backward-direction forces T1, T2, T3, T4 in thewheel sets Fig. 9 respectively. Further, the horizontal axis inFig. 10 indicates a measuring time and a calculating time of the forward-and-backward-direction forces T1 to T4, 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 inFig. 9 . - The first
frequency adjustment unit 503 derives the time-series data of the high-frequency components of the forward-and-backward-direction forces T1 to T4 as above. - The state
variable derivation unit 504 sets the observation equation as the observation equation stored by thestorage unit 501, sets the state equation as the state equation stored by thestorage unit 501, and determines estimated values of the state variables illustrated in (44) Equation by using the Kalman filter. At this time, the statevariable derivation unit 504 uses time-series data of the high-frequency components of the forward-and-backward-direction forces T1 to T4 generated by the firstfrequency adjustment unit 503. In the present embodiment, out of the time-series data of the measured value of the acceleration of thevehicle body 11 in the right and left direction, the time-series data of the measured values of the accelerations of thebogies - 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 reduce (minimize) the difference between, of an observable variable (observation variable), a measured value and an estimated value. The state
variable derivation unit 504 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. - 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. Note that the Kalman filter itself can be fabricated by a well-known technique, and thus its detailed explanation is omitted.
- In the technique described in
Patent Literature 1, as the values to be given as the measured values of the observation variables, the measured values (the measured value of the acceleration of thevehicle body 11 in the right and left direction, the measured values of the accelerations of thebogies variable derivation unit 504 derives the estimated value of the state variable so as to minimize an error of the estimated value of the observation variable with respect to the fixed value (0 (zero) in this case) or minimize an expected value of the error. - The state
variable derivation unit 504 determines the estimated values of the state variables illustrated in (44) Equation at a predetermined sampling period, to thereby generate time-series data of the estimated values of the state variables illustrated in (44) Equation. - When the signal strength of the low-frequency component included in the time-series data of the measured value of the forward-and-backward-direction force is not removed sufficiently by the first
frequency adjustment unit 503, the signal of the low-frequency component due to the railway vehicle traveling on the curved track may remain in the time-series data of the estimated values of the state variables generated by the statevariable derivation unit 504. Accordingly, the secondfrequency adjustment unit 505 reduces (preferably removes) the signal strength of the low-frequency component included in the time-series data of the estimated values of the state variables generated by the statevariable derivation unit 504. Note that 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 included in the time-series data of the measured value of the forward-and-backward-direction force by the firstfrequency adjustment unit 503, the processing of the secondfrequency adjustment unit 505 is no longer required. - In the present embodiment, the second
frequency adjustment unit 505 uses the corrected AR model to reduce the signal strength of the low-frequency component included in the time-series data of the estimated values of the state variables similarly to the firstfrequency adjustment unit 503. - The second
frequency adjustment unit 505 performs the following processing for each state variable at a predetermined sampling period. - 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 results in the time-series data of the estimated values of the state variables generated by the state
variable derivation unit 504. Each of the estimated values of the state variables varies according to the state of the railway vehicle. - First, the second
frequency adjustment unit 505 generates the autocorrelation matrix R by 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. - Next, the second
frequency adjustment unit 505 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 σ11 to σmm of the autocorrelation matrix R from the diagonal matrix Σ. - Next, the second
frequency adjustment unit 505 chooses s pieces of the eigenvalues σ11 to σss from the largest from among the plural eigenvalues σ11 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. The data y of the estimated value of each state variable can be obtained in a manner as explained so far in a state of causing the railway vehicle to travel on thetrack 16, for example. Subsequently, a distribution of the eigenvalues of the autocorrelation matrix R is made by the secondfrequency adjustment unit 505 individually for each state variable. From the distributions of the eigenvalues of the autocorrelation matrix R, the secondfrequency adjustment unit 505 determines the number s of the eigenvalues to be extracted from the autocorrelation matrix R illustrated in (53) Equation for each of the state variables. - Next, the second
frequency adjustment unit 505 determines the coefficient α of the corrected AR model by using (57) Equation based on the data y of the estimated value of the state variable, the eigenvalues σ11 to σss, and the orthogonal matrix U obtained by the singular value decomposition of the autocorrelation matrix R. - Subsequently, the second
frequency adjustment unit 505 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 results in the time-series data from which the low-frequency component included in the data y of the estimated value of the state variable has been extracted. - Subsequently, the second
frequency adjustment unit 505 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, a value of 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 at each sampling time is referred to as a value of a high-frequency component of the state variable as necessary. -
- In the present embodiment, as illustrated in (60) Equation to (63) Equation, relational expressions representing the relations between the forward-and-backward-direction forces T1 to T4 and the alignment irregularity amounts yR1 to yR4 at the positions of the wheel sets 13a to 13d are set.
- The track
state derivation unit 506 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. Subsequently, the trackstate derivation unit 506 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 by the secondfrequency adjustment unit 505, and the values of the high-frequency components of the forward-and-backward-direction forces T1 to T4 generated by the firstfrequency adjustment unit 503 to (60) Equation to (63) Equation, to thereby calculate the alignment irregularity amounts yR1 to yR4 at the positions of the wheel sets 13a to 13d. The state variables to be used here are the displacements yt1 and yt2 of thebogies bogies state derivation unit 506 performs such a calculation of the alignment irregularity amounts yR1 to yR4 as above at a predetermined sampling period, to thereby obtain the time-series data of the alignment irregularity amounts yR1 to yR4. - Subsequently, the track
state derivation unit 506 calculates a final alignment irregularity amount yR from the alignment irregularity amounts yR1 to yR4. For example, the trackstate derivation unit 506 matches phases of the time-series data of the alignment irregularity amounts yR2 to yR4 to a phase of the time-series data of the alignment irregularity amount yR1. That is, the trackstate derivation unit 506 calculates, from the distance between thewheel 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 thewheel set 13a passes through a certain position and the time when the wheel sets 13b to 13d pass through the certain position. The trackstate derivation unit 506 displaces the phases of the time-series data of the alignment irregularity amounts yR2 to yR4 by this delay time. - The track
state derivation unit 506 calculates an arithmetic mean value of the sum of the values of the alignment irregularity amounts yR1 to yR4 whose phases are matched at the same sampling time as the final alignment irregularity amount yR at this sampling time. The trackstate derivation unit 506 performs such a calculation at each sampling time, to thereby obtain time-series data of the final alignment irregularity amount yR. The phases of the alignment irregularity amounts yR2 to yR4 are matched to the phase of the alignment irregularity amount yR1, thereby making it possible to cancel disturbance factors existing in common in the time-series data of the alignment irregularity amounts yR1 to yR4. - Note that the track
state derivation unit 506 may find a moving average of each of the alignment irregularity amounts yR1 to yR4 whose phases are matched (namely, pass each of the alignment irregularity amounts yR1 to yR4 through a low-pass filter) and calculate the final alignment irregularity amount yR from the alignment irregularity amounts yR1 to yR4 whose moving averages have been found. - Further, the track
state derivation unit 506 may calculate, as the final alignment irregularity amount yR, an arithmetic mean value of two of the values of the alignment irregularity amounts yR1 to yR4 whose phases are matched at the same sampling time from which the maximum value and the minimum value are removed. - The
inspection apparatus 500 uses the time-series data of the measured value of the forward-and-backward-direction force at each sampling time acquired by thedata acquisition unit 502 while the railway vehicle is traveling in the traveling section being the target of deriving the alignment irregularity amount, to execute pieces of the processing of the firstfrequency adjustment unit 503, the statevariable derivation unit 504, the secondfrequency adjustment unit 505, and the trackstate derivation unit 506. - In this manner, the track
state derivation unit 506 can obtain the alignment irregularity amount yR at each sampling time while the railway vehicle is traveling in the traveling section being the target of deriving the alignment irregularity amount. The trackstate derivation unit 506 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. The traveling position of the railway vehicle can be set to the position of thewheel set 13a, for example. The trackstate derivation unit 506 derives the final alignment irregularity amount yR at each traveling position of the railway vehicle based on the alignment irregularity amount yR at each sampling time and the traveling position of the railway vehicle at each sampling time. - Note that the track
state derivation unit 506 does not always need to calculate the traveling position of the railway vehicle at each sampling time as described previously. For example, the trackstate derivation unit 506 may derive the traveling position of the railway vehicle at each sampling time by using a GPS (Global Positioning System). - The
output unit 507 outputs information of the final alignment irregularity amount yR that is calculated by the trackstate derivation unit 506. At this time, theoutput unit 507 may output information indicating that thetrack 16 is abnormal in the case where the final alignment irregularity amount yR 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. - As described above, in the present embodiment, the
inspection apparatus 500 gives the measured values of the forward-and-backward-direction forces T1 to T4 and the actual values of the transformation variables e1 to e4 to the Kalman filter, to derive the state variables (yw1· to yw4·, yw1 to yw4, yt1· and yt2·, yt1 and yt2, φ t1· and φ t2·, φ t1 and φ t2, φ t1· and φ t2· , φ t1 and φ t2, yb·, yb, φ b·, φ b, φ b·, φ b, φ y1, φ y2, φ a1, φ a2). At this time, the preset fixed value (0 (zero), for example) is used as the values to be normally given as the measured values of the observation variables (the accelerations of thevehicle body 11, thebogies inspection apparatus 500 uses the pivot amounts (angular displacements) φ t1 and φ t2 of thebogies inspection apparatus 500 substitutes the pivot amounts (angular displacements) φ w1 to φ w4 of the wheel sets 13a to 13d in the yawing direction, the state variables, and the measured values of the forward-and-backward-direction forces T1 to T4 into the motion equations that describe the yawings of the wheel sets 13a to 13d, to derive the alignment irregularity amounts yR1 to yR4 at the positions of the wheel sets 13a to 13d. Subsequently, theinspection apparatus 500 derives the final alignment irregularity amount yR from the alignment irregularity amounts yR1 to yR4. Therefore, it is possible to derive the alignment irregularity amounts yR1 to yR4 (the final alignment irregularity amount yR) without greatly decreasing accuracy, with no use of the measured values of the accelerations of thevehicle body 11, thebogies - Further, in the present embodiment, the
inspection apparatus 500 generates the autocorrelation matrix R from the time-series data of the measured value y of the forward-and-backward-direction force, and determines, 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, the coefficient α of the corrected AR model approximating the time-series data of the measured value y of the forward-and-backward-direction force. Therefore, it is possible to determine the coefficient α so as to make the signal of the low-frequency component included in the time-series data of the measured value y of the forward-and-backward-direction force remain and prevent the high-frequency component from remaining. Theinspection apparatus 500 calculates the predicted value y^k of the forward-and-backward-direction force at the time k by giving the time-series data of the measured value y of the forward-and-backward-direction force at the time k-l (1 ≦ 1 ≦ m), which is prior to the time k, to the corrected AR model whose coefficient α is determined in this manner. Therefore, 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 time-series data of the measured value y of the forward-and-backward-direction force without estimating a cutoff frequency beforehand. Subsequently, theinspection apparatus 500 reduces the signal strength of the low-frequency components included in the time-series data of the measured values of the forward-and-backward-direction forces T1 to T4 in this manner and generates the time-series data of the high-frequency components of the forward-and-backward-direction forces T1 to T4. Theinspection apparatus 500 gives the time-series data of the high-frequency components of the forward-and-backward-direction forces T1 to T4 to the relational expression between the forward-and-backward-direction forces T1 to T4 and the alignment irregularity amounts yR1 to yR4 at the positions of the wheel sets 13a to 13d, to thereby calculate the alignment irregularity amounts yR1 to yR4 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)). Therefore, it is possible to detect the alignment irregularity amounts yR1 to yR4 (the final alignment irregularity amount yR) in the curved track based on the motion equations that describe the motions of the railway vehicle when traveling on the linear track without using a special measuring apparatus. - In the present embodiment, the preset fixed value is given as the value to be normally given as the measured value of the observation variable when performing the data assimilation. This fixed value is not limited to 0 (zero). For example, it is also possible that the time-series data of the measured value of the acceleration of the
vehicle body 11 in the right and left direction, the time-series data of the measured values of the accelerations of thebogies inspection apparatus 500 mounted thereon or a railway vehicle equivalent to this railway vehicle (a railway vehicle whose structure is the same as that of this railway vehicle) travels on thetrack 16 being the target of deriving the alignment irregularity amounts yR1 to yR4 (the final alignment irregularity amount yR) are obtained, and an average value of the respective pieces of time-series data is used as the fixed value. Further, by using these measured values, time-series data of a predicted value of acceleration y^k of thevehicle body 11 in the right and left direction, time-series data of predicted values of accelerations y^k of thebogies track 16, and measured values of accelerations during a period in which the measured values of the forward-and-backward-direction forces T1 to T4 in the wheel sets 13a to 13d have been obtained are not used when deriving the state variables. - In the present embodiment, the case of using the corrected AR model has been explained as an example. 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 time-series data of the measured value y 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 time-series data of the measured value y of the forward-and-backward-direction force. Further, when the track on which the railway vehicle travels is a linear track (being an ideal track whose curvature is 0 (zero)) or a track which is designed as a linear track but has a curvature which is small enough to exert no influence on estimation accuracy of the alignment irregularity amount, the processing of the first
frequency adjustment unit 503 and the processing of the secondfrequency adjustment unit 505 are no longer required. - Further, in the present embodiment, the case where the wheel set to be a standard when performing the phase matching is the
wheel set 13a has been explained as an example. However, the wheel set to be a standard may be thewheel set wheel set 13a. - Further, in the present 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 of the estimated value of the observation variable with respect to the fixed 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. Note that as the error of the estimated value of the observation variable with respect to the fixed value, for example, a square error between the estimated value of the observation variable and the fixed value is cited.
- Further, in the present 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 information that reflects the track irregularity (appearance failure of the track 16) is derived as information reflecting the state of the
track 16. For example, in addition to or in place of the alignment irregularity amount, the following calculations of (64) Equation to (67) Equation may be performed, to thereby derive a lateral force that occurs when the railway vehicle travels on the linear track (stress in the right and left direction between the wheel and the rail). Note that Q1, Q2, Q3, Q4 are lateral forces in thewheels
[Mathematical equation 32] - Further, in the present embodiment, the case of including the state variables that represent the state of the
vehicle body 11 has been explained as an example. However, thevehicle body 11 is a part into which vibrations by acting forces between thewheels 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 thevehicle body 11 in the case where the influence by the propagation in thevehicle body 11 is judged to be small, for example. In such a case, out of the motion equations of (1) Equation to (21) Equation, (15) Equation to (17) Equation (the motion equations that describe the transversal vibration, the yawing, and the rolling of the vehicle body 11) and (18) Equation and (19) Equation (the motion equations that describe the yawing of the yaw damper disposed on thebogie 12a and the yawing of the yaw damper disposed on thebogie 12b) are no longer required. Further, in the motion equations of (1) Equation to (21) Equation, the state amount relating to the vehicle body (state amount including the subscript of b) and the value inside {} that includes 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)) are set to 0 (zero) . - Further, in the present embodiment, the case of the
bogies bogies - Next, a second embodiment will be described. In the first embodiment, the case where the state variables are derived by using the filter (Kalman filter) that performs the data assimilation by setting the values to be normally given as the measured values of the observation variables (the acceleration of the
vehicle body 11 in the right and left direction, the accelerations of thebogies Fig. 1 to Fig. 10 are added to the same parts as those in the first embodiment, or the like, and their detailed explanations are omitted. -
- (68) Equation is one example of equation obtained by changing an equation representing the motion equations of (9) Equation, (10) Equation, (13) Equation to (21) Equation, and (34) Equation to (39) Equation by using the state variables illustrated in (44) Equation (an equation corresponding to (68) Equation in which c is set to 1), so that a temporal change in the state variables becomes smaller than that of this equation. Concretely, (68) Equation is one obtained in a manner that, in an equation representing the motion equations of (9) Equation, (10) Equation, (13) Equation to (21) Equation, and (34) Equation to (39) Equation by using the state variables illustrated in (44) Equation, a term which is connected to a term of first-order time differential (X·) of the state variable by an equal sign is multiplied by a forgetting factor c. Specifically, (68) Equation corresponds to the state equation of (59) Equation in which the forgetting factor c is introduced, to thereby set the system noise W to 0 (zero).
- The forgetting factor c is a preset value, and is (theoretically) a value of greater than 0 and 1 or less (0 < c ≦ 1). The forgetting factor c functions in a manner that the smaller the value thereof, the greater the degree of forgetfulness of a past observation value. In (68) Equation, the smaller the value of the forgetting factor c, the smaller the influence of the measured value of the forward-and-backward-direction force on the estimated value (solution) of the state variable. For this reason, from a viewpoint of correctly obtaining the estimated value (solution) of the state variable, the value of the forgetting factor c is desirably close to 1. On the other hand, when the value of the forgetting factor c is excessively large, there is a high possibility of divergence of the estimated value (solution) of the state variable. In the present embodiment, (68) Equation is directly solved without using the filter that performs the data assimilation. Accordingly, there is a need to suppress the divergence of the estimated value (solution) of the state variable. The value of the forgetting factor c is determined based on the viewpoint as described above. As the forgetting factor c, for example, a value of greater than 0.0 and 1.0 or less (0.0 < c ≦ 1.0), preferably a value of 0.90 or more and 1.0 or less (0.90 ≦ c ≦ 1.0), more preferably a value of 0.95 or more and 1.0 or less (0.95 ≦ c ≦ 1.0), still more preferably a value of 0.99 or more and 1.0 or less (0.99 ≦ c ≦ 1.0), and most preferably 1.0 is selected among them.
- However, it is essential that the forgetting factor c is selected so as to prevent the divergence of the estimated value (solution) of the state variable obtained by solving (68) Equation. If the divergence of the estimated value (solution) of the state variable obtained by solving (68) Equation does not occur, the estimated value (solution) of this state variable when the value of the forgetting factor is 1.0 becomes the most accurate solution. However, when the value of the forgetting factor c is 1.0, it is highly likely that the estimated value (solution) of the state variable obtained by solving (68) Equation diverges (the solution is not determined).
- From such a viewpoint, the forgetting factor c may be selected by setting an upper limit value of the forgetting factor c to less than 1.0. Specifically, the value of the forgetting factor c can be selected from, for example, a value of greater than 0.0 and less than 1.0 (0.0 < c < 1.0), preferably a value of 0.90 or more and less than 1.0 (0.90 ≦ c < 1.0), more preferably a value of 0.95 or more and less than 1.0 (0.95 ≦ c < 1.0), and still more preferably a value of 0.99 or more and less than 1.0 (0.99 ≦ c < 1.0).
- Note that when the value of the forgetting factor c is 1.0, (68) Equation corresponds to the motion equations themselves of (9) Equation, (10) Equation, (13) Equation to (21) Equation, and (34) Equation to (39) Equation (an equation only representing the motion equations by using the state variables).
- The state
variable derivation unit 504 uses the time-series data of the high-frequency components of the forward-and-backward-direction forces T1 to T4 generated by the firstfrequency adjustment unit 503 to derive the actual values of the transformation variables e1 to e4 and substitutes the actual values into (34) Equation to (37) Equation, and substitutes the time-series data of the high-frequency components of the forward-and-backward-direction forces T1 to T4 generated by the firstfrequency adjustment unit 503 into (38) Equation and (39) Equation as the measured values of the forward-and-backward-direction forces T1 to T4 to solve the equation of (68) Equation, thereby determining the estimated values of the state variables illustrated in (44) Equation. The method of solving the equation of (68) Equation can be realized by a well-known numerical solution (Euler method or the like), for example. Therefore, when deriving the estimated values of the state variables, the statevariable derivation unit 504 does not use the time-series data of the measured value of the acceleration of thevehicle body 11 in the right and left direction, the time-series data of the measured values of the accelerations of thebogies - As described above, in the present embodiment, the
inspection apparatus 500 gives the measured values of the forward-and-backward-direction forces T1 to T4 and the actual values of the transformation variables e1 to e4 to the equation being the state equation with the system noise W set to 0 (zero) in which the term other than the time differential term X· of the state variable is multiplied by the forgetting factor c, to thereby derive the state variables (yw1· to yw4·, yw1 to yw4, yt1· and yt2·, yt1 and yt2, φ t1· and φ t2·, φ t1 and φ t2, φ t1· and φ t2·, φ t1 and φ t2, yb·, yb, φ b·, φ b, φ b·, φ b, φ y1, φ y2, φ a1, φ a2). Therefore, it is possible to derive the alignment irregularity amounts yR1 to yR4 (the final alignment irregularity amount yR) without greatly decreasing accuracy, with no use of the measured values of the accelerations of thevehicle body 11, thebogies -
- G is a vector in which the term that does not depend on the state equation is stored. f is a matrix corresponding to the vector G.
- Next, there will be explained a third embodiment.
- In the first and second embodiments, the case where the
inspection apparatus 500 mounted on the railway vehicle calculates the final alignment irregularity amount yR has been explained as an example. In contrast to this, in the present embodiment, a data processing device in which some functions of theinspection apparatus 500 are mounted is disposed in an operation center. The data processing device receives measured data transmitted from the railway vehicle and calculates the final alignment irregularity amount yR by using the received measured data. In this manner, in the present embodiment, the functions possessed by theinspection apparatus 500 in the first and second embodiments are shared and executed by the railway vehicle and the operation center. Constitutions and processing due to this are mainly different between the present embodiment and the first and second embodiments. Accordingly, in the explanation of the present embodiment, the same reference numerals and symbols as those added toFig. 1 to Fig. 10 are added to the same parts as those in the first and second embodiments, or the like, and their detailed explanations are omitted. Note that the present embodiment can be applied to any of the first and second embodiments. -
Fig. 11 is a view illustrating one example of a configuration of an inspection system. InFig. 11 , the inspection system includesdata collecting devices data processing device 1120. InFig. 11 , one example of functional configurations of thedata collecting devices data processing device 1120 is also illustrated. Note that each hardware of thedata collecting devices data processing device 1120 can be fabricated by the one illustrated inFig. 6 , for example. Accordingly, detailed explanations of the hardware configurations of thedata collecting devices data processing device 1120 are omitted. - The
data collecting devices data processing device 1120 is disposed at the operation center. The operation center centrally manages operations of a plurality of railway vehicles, for example. - The
data collecting devices data collecting devices data acquisition units data transmission units - The
data acquisition units data acquisition unit 502. That is, thedata acquisition units data acquisition unit 502. The configuration for obtaining the measured value of the forward-and-backward-direction force is the same as that explained in the first embodiment. - The
data transmission units data acquisition units data processing device 1120. In the present embodiment, thedata transmission units data acquisition units data processing device 1120 by radio. At this time, thedata transmission units data collecting devices data acquisition units data transmission units - A
data reception unit 1121 receives the time-series data of the measured values of the forward-and-backward-direction forces transmitted by thedata transmission units - A
data storage unit 1122 stores the time-series data of the measured values of the forward-and-backward-direction forces received by thedata reception unit 1121. Thedata storage unit 1122 stores the time-series data of the measured value of the forward-and-backward-direction force for each identification number of the railway vehicle. Thedata storage unit 1122 specifies the position of the railway vehicle at the time of receipt of the time-series data of the measured value of the forward-and-backward-direction force based on the current operation situation of the railway vehicle and the time of receipt of the time-series data of the measured value of the forward-and-backward-direction force, and stores information of the specified position and the time-series data of the measured value of the forward-and-backward-direction force in association with each other. Note that thedata collecting devices - A
data reading unit 1123 reads the time-series data of the measured value of the forward-and-backward-direction force stored in thedata storage unit 1122. Thedata reading unit 1123 can read, out of the time-series data of the measured values of the forward-and-backward-direction forces stored in thedata storage unit 1122, data designated by an operator. Further, thedata reading unit 1123 can also read the time-series data of the measured value of the forward-and-backward-direction force matching a preset condition at a preset timing. In the present embodiment, the time-series data of the measured value of the forward-and-backward-direction force read by thedata reading unit 1123 is determined based on at least any one of the identification number and the position of the railway vehicle, for example. - A
storage unit 501, a firstfrequency adjustment unit 503, a statevariable derivation unit 504, a secondfrequency adjustment unit 505, a trackstate derivation unit 506, and anoutput unit 507 are the same as those explained in the first embodiment. Accordingly, their detailed explanations are omitted here. Note that the firstfrequency adjustment unit 503 uses the time-series data of the measured value of the forward-and-backward-direction force read by thedata reading unit 1123 in place of using the time-series data of the measured value of the forward-and-backward-direction force acquired by thedata acquisition unit 502, and generates the time-series data of the high-frequency components of the forward-and-backward-direction forces T1 to T4. - As described above, in the present embodiment, the
data collecting devices data processing device 1120. Thedata processing device 1120 disposed at the operation center stores the time-series data of the measured values of the forward-and-backward-direction forces received from thedata collecting devices data processing device 1120 can calculate the final alignment irregularity amount yR at an arbitrary timing by reading the measured data at an arbitrary timing. Further, thedata processing device 1120 can output time-series variation of the final alignment irregularity amount yR at the same position. Further, thedata processing device 1120 can output the final alignment irregularity amounts yR in a plurality of routes for each route. - In the present embodiment, the case where the
data collecting devices data processing device 1120 has been explained as an example. However, it is not always necessary to design as above. An inspection system may be built by using cloud computing, for example. - Besides, also in the present embodiment, it is possible to employ the various modified examples explained in the first and second embodiments.
- Further, in the first and second embodiments, the case where the
storage unit 501, thedata acquisition unit 502, the firstfrequency adjustment unit 503, the statevariable derivation unit 504, the secondfrequency adjustment unit 505, the trackstate derivation unit 506, and theoutput unit 507 are included in one apparatus has been explained as an example. However, it is not always necessary to design as above. Functions of thestorage unit 501, thedata acquisition unit 502, the firstfrequency adjustment unit 503, the statevariable derivation unit 504, the secondfrequency adjustment unit 505, the trackstate derivation unit 506, and theoutput unit 507 may be fabricated by a plurality of apparatuses. In this case, the inspection system is constituted by using these plural apparatuses. - Next, calculation examples will be explained. In the present calculation example, derivation of a final alignment irregularity amount yR based on the method of the first embodiment and derivation of a final alignment irregularity amount yR based on the method of the second embodiment were performed. In the method of the first embodiment, when performing the data assimilation, the value to be normally given as the measured value of the observation variable (fixed value) was set to 0 (zero). Further, in the method of the second embodiment, the forgetting factor c was set to 0.9987.
- Further, in contrast to the method of the first embodiment, derivation of a final alignment irregularity amount yR was performed based on a method in which as the values to be given as the measured values of the observation variables (the measured values of the accelerations of the
vehicle body 11, thebogies -
Fig. 12 illustrates the present calculation example, and is a view illustrating thecurvature 1/R of thetrack 16 being the target of deriving the alignment irregularity amount and the traveling velocity v of the railway vehicle. InFig. 12 , agraph 1201 indicates the traveling velocity of the railway vehicle, and agraph 1202 indicates thecurvature 1/R of thetrack 16. Note that the horizontal axis inFig. 12 indicates an elapsed time (second) from a reference time when the reference time is set to 0 (zero). - Each of
Fig. 13A and Fig. 13B illustrates the present calculation example, and is a view illustrating a distribution of eigenvalues of an autocorrelation matrix R.Fig. 13A illustrates a distribution of eigenvalues of an autocorrelation matrix R with respect to the forward-and-backward-direction force T1 in thewheel set 13a, andFig. 13B illustrates a distribution of eigenvalues of an autocorrelation matrix R with respect to the forward-and-backward-direction force T2 in thewheel set 13b. -
Fig. 14 illustrates the present calculation example, and is a view illustrating time-series data of the measured values y of the forward-and-backward-direction forces T1, T2 and time-series data of predicted values y^k of the forward-and-backward-direction forces T1, T2 (time-series data obtained by extracting low-frequency components included in the time-series data of the measured values y of the forward-and-backward-direction forces). InFig. 14 , the measured value indicates the time-series data of the measured value y of the forward-and-backward-direction force, and bias indicates the time-series data of the predicted value y^k of the forward-and-backward-direction force. Note that the horizontal axis inFig. 14 indicates a measuring time and a calculating time of the forward-and-backward-direction forces T1 to T4, each of which is an elapsed time (second) from a reference time when the reference time is set to 0 (zero). -
Fig. 15 illustrates the present calculation example, and is a view illustrating the time-series data of the high-frequency components of the forward-and-backward-direction forces T1, T2. The time-series data of the high-frequency components of the forward-and-backward-direction forces T1, T2 can be obtained by subtracting the time-series data of the predicted values y^k of the forward-and-backward-direction forces T1, T2 from the time-series data of the measured values y of the forward-and-backward-direction forces T1, T2 illustrated inFig. 14 . Note that the horizontal axis inFig. 15 indicates an elapsed time (second) from a reference time when the reference time is set to 0 (zero) and indicates a calculating time of the time-series data of the high-frequency components of the forward-and-backward-direction forces T1, T2. -
Fig. 16A andFig. 16B are views illustrating the alignment irregularity amounts yR derived by the method of the first embodiment and the method described inPatent Literature 1 by using the time-series data of the high-frequency components of the forward-and-backward-direction forces T1, T2 illustrated inFig. 15 . InFig. 16A , the calculated value indicates the alignment irregularity amount yR derived by the method described inPatent Literature 1, and the measured value indicates the measured value of the alignment irregularity amount yR. InFig. 16B , the calculated value indicates the alignment irregularity amount yR derived by the method of the first embodiment, and the measured value indicates the measured value of the alignment irregularity amount yR. Here, as the calculated value of the alignment irregularity amount yR, an average value of the alignment irregularity amount yR at the position of thewheel set 13a and the alignment irregularity amount yR2 at the position of thewheel set 13b was used. Further, the measured value illustrated inFig. 16A and the measured value illustrated inFig. 16B are the same. Note that the horizontal axis in each ofFig. 16A andFig. 16B indicates an elapsed time (second) from a reference time when the reference
time is set to 0 (zero) and indicates a time corresponding to the position where the alignment irregularity amount yR exists. Further, inFig. 16A , the illustration of data of a part at which a distance from a starting point of the railway vehicle is small is omitted for convenience of illustration. -
Fig. 17A is a view illustrating the alignment irregularity amount yR derived by the method of the second embodiment by using the time-series data of the high-frequency components of the forward-and-backward-direction forces T1, T2 illustrated inFig. 15 .Fig. 17B is a view illustrating the alignment irregularity amount yR derived by the method described inPatent Literature 1 by using the time-series data of the high-frequency components of the forward-and-backward-direction forces T1, T2 illustrated inFig. 15 . InFig. 17A , the calculated value indicates the alignment irregularity amount yR derived by the method described inPatent Literature 1, and the measured value indicates the measured value of the alignment irregularity amount yR. InFig. 17B , the calculated value indicates the alignment irregularity amount yR derived by the method of the second embodiment, and the measured value indicates the measured value of the alignment irregularity amount yR. Here, as the calculated value of the alignment irregularity amount yR, an average value of the alignment irregularity amount yR1 at the position of thewheel set 13a and the alignment irregularity amount yR2 at the position of thewheel set 13b was used. Further, the measured value illustrated inFig. 17A and the measured value illustrated inFig. 17B are the same (these measured values are also the same as the measured values illustrated inFig. 16A andFig. 16B ). Note that the horizontal axis in each ofFig. 17A andFig. 17B indicates an elapsed time (second) from a reference time when the reference time is set to 0 (zero) and indicates a time corresponding to the position where the alignment irregularity amount yR exists. Further, inFig. 17A andFig. 17B , the illustration of data of a part at which a distance from a starting point of the railway vehicle is small is omitted for convenience of illustration. - When the calculated value in
Fig. 16A and the calculated value inFig. 16B are compared, it can be understood that the alignment irregularity amount yR derived by the method of the first embodiment matches the alignment irregularity amount yR derived by the method described inPatent Literature 1 with good accuracy. Further, it can be understood that the calculated value and the measured value also match with good accuracy. In like manner, when the calculated value inFig. 17A and the calculated value inFig. 17B are compared, it can be understood that the alignment irregularity amount yR derived by the method of the second embodiment matches the alignment irregularity amount yR derived by the method described inPatent Literature 1 with good accuracy. Further, it can be understood that the calculated value and the measured value also match with good accuracy. Further, when the calculated value inFig. 16B and the calculated value inFig. 17B are compared, it can be understood that the both values are almost the same, and thus the equivalent alignment irregularity amount yR can be derived by both the method of the first embodiment and the method of the second embodiment. - Note that 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.
- 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 spirit or main features thereof.
- Note that the entire contents of the description and drawings of
Patent Literature 1 can be incorporated herein by reference. - The present invention can be utilized for inspecting railway vehicles, for example.
Claims (16)
- An inspection system, comprising:a data acquisition means that acquires data of a measured value of a forward-and-backward-direction force as data of a measured value to be measured by causing a railway vehicle including a vehicle body, a bogie, and a wheel set to travel on a track;a state variable derivation means that derives state variables being variables to be determined in a state equation constituted by using motion equations that describe motions of the railway vehicle, by using the measured value of the forward-and-backward-direction force; anda track state derivation means that derives information reflecting a state of the track, whereinthe forward-and-backward-direction force is a force in a forward and backward direction that occurs in a member disposed between the wheel set and the bogie on which the wheel set is provided, and is a force to be determined according to a difference between an angular displacement of the wheel set in a yawing direction and an angular displacement of the bogie on which the wheel set is provided in the yawing direction,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 yawing direction is a pivoting direction with an up and down direction being a direction vertical to the track set as a pivot axis,the state equation is an equation described by using the state variables, the forward-and-backward-direction force, and a transformation variable,the state variables include a displacement and a velocity of the bogie in a right and left direction, an angular displacement and an angular velocity of the bogie in the yawing direction, an angular displacement and an angular velocity of the bogie in a rolling direction, a displacement and a velocity of the wheel set in the right and left direction, and an angular displacement of an air spring attached to the railway vehicle in the rolling direction, and do not include an angular displacement and an angular velocity of the wheel set in the yawing direction,the rolling direction is a pivoting direction with the forward and backward direction set as a pivot axis,the transformation variable is a variable that performs mutual transformation between the angular displacement of the wheel set in the yawing direction and the angular displacement of the bogie in the yawing direction,the track state derivation means uses the angular displacement of the bogie in the yawing direction, which is one of the state variables derived by the state variable derivation means, and an actual value of the transformation variable to derive an estimated value of the angular displacement of the wheel set in the yawing direction, and uses the derived estimated value of the angular displacement of the wheel set in the yawing direction to derive the information reflecting the state of the track,the actual value of the transformation variable is derived by using the measured value of the forward-and-backward-direction force, andthe state variable derivation means derives the state variables without using measured values of accelerations of the bogie, the wheel set, and the vehicle body in the right and left direction during a period in which the measured value of the forward-and-backward-direction force has been obtained.
- The inspection system according to claim 1, whereinthe state variable derivation means performs an operation using a filter performing data assimilation by using the state equation and an observation equation, to derive the state variables,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 and the wheel set in the right and left direction, andthe state variable derivation means uses the state equation into which the measured value of the forward-and-backward-direction force and the actual value of the transformation variable are substituted, and the observation equation into which the actual value of the transformation variable is substituted, by setting a value to be normally given as a measured value of the observation variable to a preset fixed value when performing the data assimilation, and derives the state variables when an error of a calculated value of the observation variable with respect to the fixed value or an expected value of the error becomes minimum.
- The inspection system according to claim 2, wherein
the fixed value is 0. - The inspection system according to claim 2 or 3, whereinthe state equation is constituted by using a motion equation that describes motion of the wheel set in the right and left direction, a motion equation that describes motion of the bogie in the right and left direction, a motion equation that describes motion of the bogie in the yawing direction, a motion equation that describes motion of the bogie 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 bogie in the yawing direction is a motion equation described by using the forward-and-backward-direction force in place of using the angular displacement and the angular velocity of the wheel set in the yawing direction,the observation equation is constituted by using the motion equation that describes the motion of the wheel set in the right and left direction and the motion equation that describes the motion of the bogie in the right and left direction,the motion equation that describes the motion of the wheel set in the right and left direction is a motion equation described by using the transformation variable in place of using the angular displacement of the wheel set in the yawing direction, andthe transformation variable is represented by a difference between the angular displacement of the bogie in the yawing direction and the angular displacement of the wheel set in the yawing direction.
- The inspection system according to claim 4, whereinthe state equation is constituted by further using a motion equation that describes motion of the vehicle body in the right and left direction, a motion equation that describes motion of the vehicle body in the yawing direction, a motion equation that describes motion of the vehicle body in the rolling direction, and a motion equation that describes motion of a yaw damper attached to the railway vehicle in the yawing direction,the observation equation is constituted by further using the motion equation that describes the motion of the vehicle body in the right and left direction,the observation variable further includes the acceleration of the vehicle body in the right and left direction, andthe state variables further include a displacement and a velocity of the vehicle body in the right and left direction, an angular displacement and an angular velocity of the vehicle body in the yawing direction, an angular displacement and an angular velocity of the vehicle body in the rolling direction, and an angular displacement of the yaw damper in the yawing direction.
- The inspection system according to any one of claims 2 to 5, wherein
the filter is a Kalman filter. - The inspection system according to claim 1, wherein
the state variable derivation means does not solve the state equation but solves an equation being an equation representing the motion equations that describe the motions of the railway vehicle by using the state variables, the forward-and-backward-direction force, and the transformation variable, into which the measured value of the forward-and-backward-direction force and the actual value of the transformation variable are substituted, to derive the state variables. - The inspection system according to claim 7, wherein
the state variable derivation means does not use the state equation but uses an equation being an equation as a result of changing the equation representing the motion equations that describe the motions of the railway vehicle by using the state variables, the forward-and-backward-direction force, and the transformation variable so that a temporal change in the state variables becomes smaller than that of the equation, into which the measured value of the forward-and-backward-direction force and the actual value of the transformation variable are substituted, to derive the state variables. - The inspection system according to claim 8, whereinthe state variable derivation means solves an equation being the equation representing the motion equations that describe the motions of the railway vehicle by using the state variables, the forward-and-backward-direction force, and the transformation variable in which each term which is connected to a term of first-order time differential of the state variable by an equal sign is multiplied by a forgetting factor, to derive the state variables, andthe forgetting factor is a preset value of 0.95 or more and less than 1.
- The inspection system according to any one of claims 7 to 9, whereinthe motion equations that describe the motions of the railway include a motion equation that describes motion of the wheel set in the right and left direction, a motion equation that describes motion of the bogie in the right and left direction, a motion equation that describes motion of the bogie in the yawing direction, a motion equation that describes motion of the bogie 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 in the right and left direction is a motion equation described by using the transformation variable in place of using the angular displacement of the wheel set in the yawing direction,the motion equation that describes the motion of the bogie in the yawing direction is a motion equation described by using the forward-and-backward-direction force in place of using the angular displacement and the angular velocity of the wheel set in the yawing direction, andthe transformation variable is represented by a difference between the angular displacement of the bogie in the yawing direction and the angular displacement of the wheel set in the yawing direction.
- The inspection system according to claim 10, whereinthe motion equations that describe the motions of the railway further include a motion equation that describes motion of the vehicle body in the right and left direction, a motion equation that describes motion of the vehicle body in the yawing direction, a motion equation that describes motion of the vehicle body in the rolling direction, and a motion equation that describes motion of a yaw damper attached to the railway vehicle in the yawing direction, andthe state variables further include a displacement and a velocity of the vehicle body in the right and left direction, an angular displacement and an angular velocity of the vehicle body in the yawing direction, an angular displacement and an angular velocity of the vehicle body in the rolling direction, and an angular displacement of the yaw damper in the yawing direction.
- The inspection system according to any one of claims 1 to 11, whereinthe track state derivation means derives an alignment irregularity amount of the track as information reflecting the state of the track, based on the displacement and the velocity of the bogie in the right and left direction, each of which is one of the state variables derived by the state variable derivation means, the displacement and the velocity of the wheel set in the right and left direction, each of which is one of the state variables derived by the state variable derivation means, the estimated value of the angular displacement of the wheel set in the yawing direction, the measured value of the forward-and-backward-direction force, and the motion equation that describes motion of the wheel set in the yawing direction, andthe motion equation that describes the motion of the wheel set in the yawing direction includes the forward-and-backward-direction force and the alignment irregularity amount of the track as variables.
- The inspection system according to any one of claims 1 to 11, wherein
the track state derivation means derives a lateral force being a stress in the right and left direction between a wheel provided on the wheel set and the track, as the information reflecting the state of the track, based on the angular displacement of the wheel set in the yawing direction and the velocity of the wheel set in the right and left direction which is one of the state variables. - The inspection system according to any one of claims 1 to 13, further comprisinga frequency adjustment means that 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 from time-series data of a physical quantity whose value varies according to a state of the railway vehicle, whereinthe frequency adjustment means includes a first frequency adjustment means that 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 from time-series data of a measured value of the forward-and-backward-direction force being one of the physical quantity, andthe state variable derivation means uses the value of the forward-and-backward-direction force from which the signal strength of the low-frequency component has been reduced by the first frequency adjustment means to derive the state variables.
- An inspection method, comprising:a data acquisition step of acquiring data of a measured value of a forward-and-backward-direction force as data of a measured value to be measured by causing a railway vehicle including a vehicle body, a bogie, and a wheel set to travel on a track;a state variable derivation step of deriving state variables being variables to be determined in a state equation constituted by using motion equations that describe motions of the railway vehicle, by using the measured value of the forward-and-backward-direction force; anda track state derivation step of deriving information reflecting a state of the track, whereinthe forward-and-backward-direction force is a force in a forward and backward direction that occurs in a member disposed between the wheel set and the bogie on which the wheel set is provided, and is a force to be determined according to a difference between an angular displacement of the wheel set in a yawing direction and an angular displacement of the bogie on which the wheel set is provided in the yawing direction,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 yawing direction is a pivoting direction with an up and down direction being a direction vertical to the track set as a pivot axis,the state equation is an equation described by using the state variables, the forward-and-backward-direction force, and a transformation variable,the state variables include a displacement and a velocity of the bogie in a right and left direction, an angular displacement and an angular velocity of the bogie in the yawing direction, an angular displacement and an angular velocity of the bogie in a rolling direction, a displacement and a velocity of the wheel set in the right and left direction, and an angular displacement of an air spring attached to the railway vehicle in the rolling direction, and do not include an angular displacement and an angular velocity of the wheel set in the yawing direction,the rolling direction is a pivoting direction with the forward and backward direction set as a pivot axis,the transformation variable is a variable that performs mutual transformation between the angular displacement of the wheel set in the yawing direction and the angular displacement of the bogie in the yawing direction,the track state derivation step uses the angular displacement of the bogie in the yawing direction, which is one of the state variables derived by the state variable derivation step, and an actual value of the transformation variable to derive an estimated value of the angular displacement of the wheel set in the yawing direction, and uses the derived estimated value of the angular displacement of the wheel set in the yawing direction to derive the information reflecting the state of the track,the actual value of the transformation variable is derived by using the measured value of the forward-and-backward-direction force, andthe state variable derivation step derives the state variables without using measured values of accelerations of the bogie, the wheel set, and the vehicle body in the right and left direction during a period in which the measured value of the forward-and-backward-direction force has been obtained.
- A computer program causing a computer to execute:a data acquisition step of acquiring data of a measured value of a forward-and-backward-direction force as data of a measured value to be measured by causing a railway vehicle including a vehicle body, a bogie, and a wheel set to travel on a track;a state variable derivation step of deriving state variables being variables to be determined in a state equation constituted by using motion equations that describe motions of the railway vehicle, by using the measured value of the forward-and-backward-direction force; anda track state derivation step of deriving information reflecting a state of the track, whereinthe forward-and-backward-direction force is a force in a forward and backward direction that occurs in a member disposed between the wheel set and the bogie on which the wheel set is provided, and is a force to be determined according to a difference between an angular displacement of the wheel set in a yawing direction and an angular displacement of the bogie on which the wheel set is provided in the yawing direction,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 yawing direction is a pivoting direction with an up and down direction being a direction vertical to the track set as a pivot axis,the state equation is an equation described by using the state variables, the forward-and-backward-direction force, and a transformation variable,the state variables include a displacement and a velocity of the bogie in a right and left direction, an angular displacement and an angular velocity of the bogie in the yawing direction, an angular displacement and an angular velocity of the bogie in a rolling direction, a displacement and a velocity of the wheel set in the right and left direction, and an angular displacement of an air spring attached to the railway vehicle in the rolling direction, and do not include an angular displacement and an angular velocity of the wheel set in the yawing direction,the rolling direction is a pivoting direction with the forward and backward direction set as a pivot axis,the transformation variable is a variable that performs mutual transformation between the angular displacement of the wheel set in the yawing direction and the angular displacement of the bogie in the yawing direction,the track state derivation step uses the angular displacement of the bogie in the yawing direction, which is one of the state variables derived by the state variable derivation step, and an actual value of the transformation variable to derive an estimated value of the angular displacement of the wheel set in the yawing direction, and uses the derived estimated value of the angular displacement of the wheel set in the yawing direction to derive the information reflecting the state of the track,the actual value of the transformation variable is derived by using the measured value of the forward-and-backward-direction force, andthe state variable derivation step derives the state variables without using measured values of accelerations of the bogie, the wheel set, and the vehicle body in the right and left direction during a period in which the measured value of the forward-and-backward-direction force has been obtained.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018230834 | 2018-12-10 | ||
PCT/JP2019/047673 WO2020121943A1 (en) | 2018-12-10 | 2019-12-05 | Inspection system, inspecting method, and program |
Publications (2)
Publication Number | Publication Date |
---|---|
EP3895955A1 true EP3895955A1 (en) | 2021-10-20 |
EP3895955A4 EP3895955A4 (en) | 2022-10-12 |
Family
ID=71076963
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP19894582.6A Pending EP3895955A4 (en) | 2018-12-10 | 2019-12-05 | Inspection system, inspecting method, and program |
Country Status (4)
Country | Link |
---|---|
EP (1) | EP3895955A4 (en) |
JP (1) | JP7088315B2 (en) |
CN (1) | CN113165676B (en) |
WO (1) | WO2020121943A1 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116258040B (en) * | 2022-12-30 | 2024-01-23 | 武汉理工大学 | Track irregularity detection method |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005231427A (en) * | 2004-02-18 | 2005-09-02 | B M C:Kk | Track monitoring device |
JP6096424B2 (en) * | 2012-05-25 | 2017-03-15 | 東海旅客鉄道株式会社 | Orbital condition monitoring device |
JP5959378B2 (en) * | 2012-09-11 | 2016-08-02 | 川崎重工業株式会社 | Load measuring method and apparatus, railway vehicle equipped with load measuring apparatus, and load management system |
EP3040251B1 (en) * | 2013-08-28 | 2018-10-17 | Nippon Steel & Sumitomo Metal Corporation | Method of decreasing lateral pressure in railroad vehicle |
EP3434552B1 (en) * | 2016-03-23 | 2021-05-05 | Nippon Steel Corporation | Inspection system, inspection method and program |
JP6747598B2 (en) * | 2017-08-31 | 2020-08-26 | 日本製鉄株式会社 | Inspection system, inspection method, and program |
JP7017179B2 (en) | 2018-07-03 | 2022-02-08 | 日本製鉄株式会社 | Inspection system, inspection method, and program |
JP7024874B6 (en) * | 2018-07-31 | 2023-12-15 | 日本製鉄株式会社 | Inspection systems, inspection methods, and programs |
-
2019
- 2019-12-05 WO PCT/JP2019/047673 patent/WO2020121943A1/en unknown
- 2019-12-05 JP JP2020560017A patent/JP7088315B2/en active Active
- 2019-12-05 EP EP19894582.6A patent/EP3895955A4/en active Pending
- 2019-12-05 CN CN201980079590.2A patent/CN113165676B/en active Active
Also Published As
Publication number | Publication date |
---|---|
WO2020121943A1 (en) | 2020-06-18 |
CN113165676A (en) | 2021-07-23 |
JPWO2020121943A1 (en) | 2021-10-28 |
JP7088315B2 (en) | 2022-06-21 |
CN113165676B (en) | 2024-04-02 |
EP3895955A4 (en) | 2022-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3434552B1 (en) | Inspection system, inspection method and program | |
EP3677485A1 (en) | Inspection system, inspection method, and program | |
Doumiati et al. | Estimation of road profile for vehicle dynamics motion: experimental validation | |
CN110874450B (en) | Railway bridge track irregularity calculation method based on vehicle-mounted monitoring | |
Muñoz et al. | Estimation of lateral track irregularity through Kalman filtering techniques | |
EP3819186B1 (en) | Inspection system, inspection method, and program | |
EP3895955A1 (en) | Inspection system, inspecting method, and program | |
EP3832284A1 (en) | Inspection system, inspecting method, and program | |
JP6939540B2 (en) | Contact angle estimation system, contact angle estimation method, and program | |
EP3992052A1 (en) | Estimation device, estimation method and program | |
WO2021014854A1 (en) | Inspection system, inspection method, and program | |
JP2022028374A (en) | Estimation device, estimation method, and program | |
Erviti et al. | Error analysis of the application of combined subspace identification to the modal analysis of railway vehicles | |
Muñoz Moreno et al. | Estimation of lateral track irregularity through Kalman filtering techniques | |
WO2024047741A1 (en) | Inference method, inference device, and inference program | |
Best et al. | Full vehicle and tyre identification using unscented and extended identifying Kalman filters | |
Ruga | Vehicle dynamics performance assessment by development of model-supported testing | |
JP2022042798A (en) | Processing device, processing method, and program | |
Saber et al. | Unknown Road Profile Input Estimation with Two Levels of Complexity Models | |
Barbosa | Identification of railway track quality and safety through the dynamic inertial response of railway carbody and truck | |
Baek et al. | Robust Estimation of Vehicle Dynamic Behavior by Using Fusion Sensors |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE |
|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
|
17P | Request for examination filed |
Effective date: 20210608 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
DAV | Request for validation of the european patent (deleted) | ||
DAX | Request for extension of the european patent (deleted) | ||
A4 | Supplementary search report drawn up and despatched |
Effective date: 20220908 |
|
RIC1 | Information provided on ipc code assigned before grant |
Ipc: B61K 9/08 20060101ALI20220902BHEP Ipc: B61L 23/04 20060101AFI20220902BHEP |