CN113446939B - Wheel-rail contact geometric parameter measuring method - Google Patents
Wheel-rail contact geometric parameter measuring method Download PDFInfo
- Publication number
- CN113446939B CN113446939B CN202110752175.8A CN202110752175A CN113446939B CN 113446939 B CN113446939 B CN 113446939B CN 202110752175 A CN202110752175 A CN 202110752175A CN 113446939 B CN113446939 B CN 113446939B
- Authority
- CN
- China
- Prior art keywords
- wheel
- point
- point cloud
- rail
- center
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
-
- 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
Abstract
The invention provides a wheel rail contact geometric parameter measuring method, which comprises the following steps: detecting plane elements and cylinder elements of the wheel-track point cloud; extracting elements where the inner side surfaces of the wheel rims are located, and roughly dividing point clouds of the wheel rims; detecting circular ring elements, extracting corresponding elements of the wheel rim part and the rail bottom part to obtain corresponding element parameters; calculating the central point of the wheel set and the central point of the rail bottom to obtain the direction of the wheel shaft; measuring the wheel set transverse displacement, the wheel set side rolling angle, the wheel set pan head angle and the wheel set sinking and floating amount, and positioning points belonging to the rail tread in the point cloud; slicing the point cloud, and removing points belonging to the rail tread in each point cloud slice; performing structural information fitting on the removed points on the slicing circle based on an RANSAC algorithm to obtain the circle center and the radius of the slicing circle; calculating the distance from a point on the rail tread in each slice to the center of a slice circle, and determining a contact patch; the long and short axis lengths of the contact patches were measured. The method can accurately measure the wheel track parameters and provide technical support for safe operation of the train.
Description
Technical Field
The invention relates to the technical field of three-dimensional vision measurement, in particular to a point cloud-based wheel-rail contact geometric parameter measurement method.
Background
The rapid development of high-speed trains enables the running speed and the density of the trains to be continuously improved, further the wheel-rail contact mechanism becomes more and more complex, the interaction between the wheel rails is stronger, particularly the wheel-rail rolling contact fatigue caused by abnormal contact postures between the wheel rails, so that the damages such as wheel-rail surface peeling, cracks and the like are caused, the service life of the track is shortened, the economy, the comfort and the running safety of the trains are seriously influenced, and the maintenance and replacement cost of the track is remarkably increased. The wheel-rail relationship is one of the leading subjects in the railway engineering research field, and the research of the wheel-rail contact relationship is the key of the wheel-rail contact relationship. The attitude and position of the wheel set relative to the rail is a very important part of the wheel-rail contact geometry, and the wheel-rail contact geometry parameters reflect the three-dimensional space geometrical relationship formed between the wheel rails under different contact attitudes and different stresses. But at the wheel rail contact area is located at the wheel rail interface area. The traditional contact measurement method cannot directly measure the wheel rail contact geometric parameters.
Therefore, a method for measuring the geometrical parameters of wheel rail contact is needed.
Disclosure of Invention
The invention provides a wheel-rail contact geometric parameter measuring method, which aims to overcome the defects in the prior art.
In order to achieve the purpose, the invention adopts the following technical scheme.
The embodiment provides a wheel-rail contact geometric parameter measuring method, which comprises the following steps:
s1, detecting the point cloud of the wheel track by plane elements and cylinder elements;
s2 extracting the elements where the inner side surface of the rim is located from the detection result of the step S1, and roughly dividing the point cloud of the rim part according to the information of the elements where the inner side surface of the rim is located;
s3, performing circular ring element detection on the point cloud of the roughly divided rim part, and extracting corresponding elements of the rim part and the rail bottom part to obtain corresponding element parameters;
s4, calculating the wheel set central point and the rail bottom central point according to the elementary parameters to obtain the wheel axle direction;
s5 measuring wheel set transverse displacement, wheel set side rolling angle, wheel set rolling angle and wheel set sinking and floating amount according to the calculation result of S4, and positioning points belonging to the rail tread in the point cloud;
s6, slicing the point clouds, and removing points belonging to the rail tread in each point cloud slice;
s7, performing slicing circle structure information fitting on the point cloud points with the tread points removed based on a RANSAC algorithm to obtain the circle center and the radius of the slicing circle;
S8, calculating the distance from the point on the rail tread to the center of the slice circle in each slice, wherein the set of points smaller than a certain threshold value is the contact patch;
s9, the length of the long axis and the length of the short axis of the contact patch are measured according to the coordinates of the center point of the contact patch.
Preferably, the wheel-track point cloud data is processed before S1 to establish wheel sets and tracksTwo different coordinate systems, when the wheel set is at the normal operation position, the middle point of the axle center of the wheel set is used as the original point to establish O w -X w Y w Z w A wheelset coordinate system that moves with the wheelset relative to the track; establishing O by using the central point of the centers of the two side rail surfaces as the origin r -X r Y r Z r A traveling coordinate system, a plane of the rail bottom and a traveling coordinate system O r -X r Y r Z r X in (1) r O r Y r Plane parallel, wherein, X r In the direction of horizontal extension of the rail surface, Y r On the plane of the rail foot and in contact with X r Perpendicular, Z r Perpendicular to the plane of the rail bottom. X w 、Y w 、Z w Are each independently of X r 、Y r 、Z r And establishing a coordinate system under a normal condition, wherein each coordinate axis of the wheel set coordinate system and the following coordinate system is pairwise parallel, and only Zr and Zw have a difference value in the Z-axis direction.
Preferably, S2 includes:
step S21: extracting the inner side surfaces omega of the two-side wheel rims from the element detection result in S1 according to the constraint conditions that the inner side distance of the wheel pair and the inner side surfaces of the wheel rims are parallel to each other 1 As the inner side of the flange and omega of the left wheel pair 2 The inner side surface of the wheel rim is used as a right wheel set;
step S22: calculating the omega of the plane according to the thickness of the standard wheel rim 1 And Ω 2 Plane omega separated from the thickness of the standard rim 1 ' and omega 2 The plane parameters of' correspond to the front view of the wheel-track point cloud data;
step S23: according to
And
roughly dividing the rim part in the wheel track original point cloud, wherein
Is a wheelFrom any point p in the orbit point cloud to plane omega i The Euclidean distance of (a) is,
from any point p to plane omega in the wheel track point cloud i The euclidean distance of' i ═ 1, 2.
Preferably, extracting the rim part and the rail bottom part corresponding cells in S3 and acquiring corresponding cell parameters includes:
step S31: the extraction of the rail bottom element is carried out according to the following three conditions: 1) the type of the model element to be extracted is a plane; 2) extracting normal vector and Z of plane r The directions are parallel; 3) extracting the minimum z coordinate value of the plane;
step S32: the elements of the rim portion are extracted according to the following three conditions: 1) the extracted primitive type is a circular ring primitive, comprising
Eight parameters, wherein,
the method comprises the steps of setting a circular ring axis direction as a circular ring axis direction, setting Center as a circular ring axis Center and setting three parameters, setting radiuMax as a circular ring large radius and setting a parameter, and setting randiuMin as a circular ring small radius and setting a parameter; 2) a wheel rim radius; 3) normal vector of extracted circular ring element To 
Normal to the inner side of the rim
Are parallel and satisfy
Wherein, the first and the second end of the pipe are connected with each other,
to set the threshold.
Preferably, calculating a wheel set center point and a rail bottom center point according to the element parameters to obtain a wheel axle direction, and the method comprises the following steps:
s41: calculating the center of the wheel pair according to the center point coordinate of the intersection point of the straight line passing through the center point of the circular ring where the wheel rim is located and the inner side faces of the two wheel rims along the direction of the wheel shaft;
s42: calculating the center of the rail bottom according to the center point coordinate of the center point of the plane bounding box of the rail plane;
s43: the normal vector of the plane element of the inner side surface of the wheel rim is the wheel axle direction.
Preferably, S5 includes:
s51: measuring the transverse displacement of the wheel set, wherein the center of the wheel set and the center point of the rail bottom are in Y r Spacing in the direction;
s52: measuring the sinking and floating amount of the wheel set, wherein the center of the wheel set and the center point of the rail bottom are in Z r Spacing and initial position Z in axial direction r The difference in the amount of shaft traverse;
s53: measuring the side rolling angle of the wheel set, and taking the axis of the wheel set as O in a following coordinate system r -X r Y r Z r X in (1) r O r Y r The included angle of the plane;
s54: measuring the pan angle of the wheel pair, as the perpendicular to the wheel axis relative to X r The angle of the axes.
Preferably, in S5, the points belonging to the rail tread in the point cloud are located according to the structural prior information of the rail, which is determined according to the railroad industry standard.
Preferably, S6 includes:
s61: estimating the point cloud density rho by adopting a nuclear density estimation method, wherein a slice thickness delta calculation formula is shown as the following formula (1);
δ＝M·ρ (1)
wherein M is an empirical value of 4-8;
s62: the coordinate system of the point cloud slice is that the X axis is along the central line direction of the point cloud at the bottom, the Y axis is parallel to the central line direction of the point cloud at the upper part of the combined point cloud, the Z axis is vertical to the bottom surface of the point cloud, and the Y axis direction is selected as the point cloud slice direction;
s63: and projecting the point set in each point cloud slice to the central plane to realize the dimension reduction processing of local data.
Preferably, S7 includes:
s71: projecting points in each point cloud piece to a point cloud piece central plane to obtain a slice projection point set, then removing points on a bottom tread in the point cloud slice by using structure prior information, and only keeping point clouds of a wheel to perform parameter fitting of a circle;
s72: randomly selecting three non-collinear points in the slice projection point set as three points capable of forming a circle model, calculating the circle center and the radius of the circle model, calculating the absolute value of the difference between the distances from the other points except the three points of the circle model in the projection point set to the circle center of the circle model and the radius of the circle model, taking the points with the absolute values of the differences smaller than a set threshold value as local points of the current circle model, taking the circle model obtained when the number of the local points reaches a set proportion as a slice circle, and taking the circle center O and the radius R as the circle center and the radius of the slice circle.
Preferably, S8 includes: calculating the distance D from the point on the tread of the bottom rail in each slice to the center of the slice i Setting a certain threshold value as R + T, T as elastic threshold value, R as radius of slice circle, if D i Satisfy D i And if the radius is less than or equal to R + T, the contact points are obtained by extraction, and the set of all the contact points is the contact spots, wherein the elastic threshold value T is 0.001-0.005 times of the radius R.
The technical scheme provided by the wheel-rail contact geometric parameter measuring method is that based on the machine vision theory, the wheel-rail point cloud is generated by using a structured light non-contact three-dimensional reconstruction technology, and on the basis, the wheel-rail contact geometric parameters are measured, so that the transverse displacement, the sinking displacement, the pan angle and the side rolling angle of a wheel pair relative to a rail can be accurately measured, a wheel-rail contact spot area is extracted, and the long axis and short axis parameters of the wheel-rail contact spot are obtained through calculation, so that the rapid, accurate and reliable technical support is provided for the safe operation of a train, an experimental means is provided for the theoretical research of the wheel-rail contact relationship, and the method has important practical significance for the theoretical research of the safe operation and the wheel-rail contact relationship.
Additional aspects and advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a schematic flow chart of a method for measuring geometric parameters of wheel-rail contact according to this embodiment;
FIG. 2 is a schematic structural diagram of two coordinate systems of the present embodiment;
FIG. 3 is a diagram illustrating the rough segmentation result;
FIG. 4 is a schematic diagram of positions of variables of a wheel set lateral shift amount, a wheel set side roll angle, a wheel set rock angle, and a wheel set heave amount, which are designed by the wheel-rail contact geometric parameter measurement method of the embodiment;
FIG. 5 is a schematic diagram of removing tread points from a slice point cloud and fitting circle information;
FIG. 6 is a schematic diagram of computing contact points in a slice using structural information;
fig. 7 is a schematic view of a wheel-rail contact patch measurement scenario.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention.
As used herein, the singular forms "a", "an", "the" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof. It will be understood that when an element is referred to as being "connected" or "coupled" to another element, it can be directly connected or coupled to the other element or intervening elements may also be present. Further, "connected" or "coupled" as used herein may include wirelessly connected or coupled. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
It will be understood by those skilled in the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.
For the convenience of understanding the embodiments of the present invention, the following description will be further explained by taking several specific embodiments as examples in conjunction with the drawings, and the embodiments are not to be construed as limiting the embodiments of the present invention.
Examples
Fig. 1 is a schematic flow chart of a method for measuring geometric parameters of wheel-rail contact according to this embodiment, and referring to fig. 1, the method includes the following steps:
processing the point cloud data of the wheel track to establish two different coordinate systems for the wheel set and the track, wherein fig. 2 is a schematic structural diagram of the two coordinate systems of the embodiment, and referring to fig. 2, the point cloud data of the wheel set is used as an origin to establish an O w -X w Y w Z w A wheelset coordinate system that moves with the wheelset relative to the track; establishing O by using the central point of the centers of the two side rail surfaces as the origin r -X r Y r Z r A traveling coordinate system, which is formed by the plane of the rail bottom and a traveling coordinate system O r -X r Y r Z r X in (1) r O r Y r Plane parallel, wherein, X r In the direction of horizontal extension of the rail surface, Y r On the plane of the rail foot and in contact with X r Perpendicular, Z r Perpendicular to the plane of the rail bottom. X w 、Y w 、Z w Are each independently of X r 、Y r 、Z r Parallel. The coordinate system is established under normal conditions, coordinate axes of the wheel set coordinate system and the following coordinate system are pairwise parallel, and only Zr and Zw have difference values in the Z-axis direction. At any position where the attitude of the wheelset changes, if the wheelset deviates from the normal position, the wheelset coordinate system is relative to Y r And Z r Deviation DeltaY r And Δ Z r Rotated by ψ about the X axis and θ about the Z axis.
S1 performs plane element and cylinder element detection on the wheeltrack point cloud.
S2, extracting the element of the inner side of the rim from the detection result of the step S1, and roughly dividing the point cloud of the rim part according to the information of the element of the inner side of the rim.
Step S21: according to the constraint condition that the inner side distance of the wheel set and the inner side surfaces of the wheel rims are parallel to each other, the inner side distance of the wheel set in the embodiment is 1353 +/-2 mm, and the inner side surfaces omega of the double-side wheel rims are extracted from the element detection result in S1 1 As the inner side of the rim and omega of the left wheel set 2 The inner side of the wheel rim of the right wheel pair is used as the inner side of the wheel rim of the right wheel pair;
step S22: calculating the omega of the plane according to the thickness of the standard wheel rim 1 And Ω 2 Plane omega separated from the thickness of the standard rim 1 ' and omega 2 ' the plane parameters, in this embodiment, the standard rim thickness is 32 + -2 mm, corresponding to the front view of the point cloud data of the wheel track, specifically, omega 1 ' is on the left side of the plane and on the plane omega 1 Planes spaced apart by 32mm, omega 2 Is on the right side of the plane and omega on the plane 2 Planes 32mm apart;
step S23: according to
And
roughly dividing the rim part in the wheel track original point cloud, wherein
From any point p to plane omega in the wheel track point cloud i The Euclidean distance of (a) is,
from any point p in the wheel-track point cloud to plane omega i The euclidean distance of' i ═ 1, 2. Fig. 3 is a schematic diagram of the rough segmentation result.
S3, performing circular ring element detection on the point cloud of the roughly divided rim part, and extracting corresponding elements of the rim part and the rail bottom part to obtain corresponding element parameters.
The method specifically comprises the following steps:
step S31: the extraction of the rail bottom element is carried out according to the following three conditions: 1) the type of the model element to be extracted is a plane; 2) extracting normal vector and Z of plane r The directions are parallel; 3) extracting the minimum z coordinate value of the plane;
step S32: the elements of the rim portion are extracted according to the following three conditions: 1) the extracted primitive type is a circular ring primitive, comprising
Eight parameters, wherein,
the method comprises the steps that three parameters are included in the direction of a circular ring shaft, Center is the Center of the circular ring shaft and comprises the three parameters, radiuMax is the large radius of the circular ring and comprises one parameter, and radiuMax is the small radius of the circular ring and comprises one parameter; 2) the radius of the wheel rim is 458 +/-2 mm; 3) normal vector of extracted circular ring element To
Normal to the inner side of the rim
Parallel and satisfy
Wherein, the first and the second end of the pipe are connected with each other,
to set the threshold, take
S4, calculating the wheel set central point and the rail bottom central point according to the elementary parameters to obtain the wheel axle direction.
The method specifically comprises the following steps:
s41: calculating the center of the wheel pair according to the center point coordinate of the intersection point of the straight line passing through the center point of the circular ring where the wheel rim is located and the inner side faces of the two wheel rims along the direction of the wheel shaft;
s42: calculating the center of the rail bottom according to the center point coordinate of the center point of the plane bounding box of the rail plane;
s43: the normal vector of the plane element of the inner side surface of the wheel rim is the wheel axle direction.
S5, measuring wheel pair transverse displacement, wheel pair side rolling angle, wheel pair rolling angle and wheel pair sinking and floating amount according to the calculation result of S4, and positioning points belonging to the rail tread in the point cloud according to the structure prior information of the rail.
Fig. 4 is a schematic diagram of positions of variables of a wheel set transverse displacement amount, a wheel set side rolling angle, a wheel set rock angle and a wheel set sinking and floating amount, which are designed by the wheel-rail contact geometric parameter measurement method of the embodiment, and referring to fig. 4, the wheel set transverse displacement amount is Δ y; the sinking and floating amount of the wheel set is delta z; the side rolling angle of the wheel pair is psi; the wheelset yaw angle is θ.
S51: measuring the transverse displacement of the wheel set, wherein the center of the wheel set and the center point of the rail bottom are in Y r Spacing in the direction;
s52: measuring the sinking and floating amount of the wheel set, wherein the center of the wheel set and the center point of the rail bottom are in Z r Spacing and initial position Z in axial direction r The difference in the amount of shaft traverse;
s53: measuring the side roll angle of the wheel set as O in the wheel axis relative to the following coordinate system r -X r Y r Z r X in (1) r O r Y r The included angle of the plane;
s54: measuring the yaw angle of the wheel pair, with the perpendicular to the wheel axis being relative to X r The angle of the axes.
And positioning points belonging to the rail tread in the point cloud according to the structure prior information of the rail, wherein the structure prior information is determined according to the rail industry standard, and specifically the TB/T2341.3-9360 kg/m steel rail type standard.
S6, slicing the point clouds, and removing points belonging to the rail treads in each point cloud slice.
S61: estimating the point cloud density rho by adopting a nuclear density estimation method, wherein a slice thickness delta calculation formula is shown as the following formula (1);
δ＝M·ρ (1)
wherein M is an empirical value of 4-8;
s62: the coordinate system of the point cloud slice is as follows: the X axis is along the central line direction of the bottom point cloud, the Y axis is parallel to the central line direction of the upper point cloud of the combined point cloud, and the Z axis is perpendicular to the bottom surface of the point cloud. Selecting a Y-axis direction from the point cloud slicing direction;
s63: and projecting the point set in each point cloud slice to the central plane to realize the dimension reduction processing of local data.
S7, performing slicing circle structure information fitting on the point cloud points with the tread points removed based on the RANSAC algorithm to obtain the circle center and the radius of the slicing circle.
S71: and projecting the points in the cloud slices to the central plane of the cloud slices to obtain a slice projection point set, removing the points on the bottom tread in the cloud slices by using the structure prior information, and only keeping the point cloud of the wheel to perform the parameter fitting of the circle.
S72: randomly selecting three non-collinear points from the projection point set of the slice as three points capable of forming a circular model, calculating the center and radius of the circular model, and calculating the distances D from the other points except the three points of the circular model in the projection point set to the center of the circular model i Absolute value of difference abs (D) from radius R of circular model i -R) And when the absolute value of the difference is smaller than a set threshold and the number of the points meeting the set threshold reaches a set proportion, taking the obtained circular model as a slicing circle, wherein the circle center O and the radius R are the circle center and the radius of the slicing circle.
To ensure that the fitting results in a suitable model, the difference threshold T is set to ensure that most points are on the fitted circle 1 When abs (D) i -R)<T 1 Considering the current point as the local point, obtaining the total number of the local points under the current model parameter after completely traversing, if the proportion of the total number of the local points is not less than the threshold value T 2 Then, a better model is found, and the center O and the radius R are recorded, and fig. 5 is a schematic diagram of removing tread points from the slice point cloud and fitting circle information. The difference threshold value T 1 Generally, the radius of the circular model R can be set to be about one thousandth of the radius R of the circular model, and in order to obtain a more accurate model proportion threshold value T 2 Not less than 85%.
S8, calculating the distance between the point on the rail tread surface in each slice and the center of the circle of the slice, wherein the set of the points smaller than a certain threshold value is the contact patch.
Calculating the distance D from the point in the bottom rail surface determined by the structure prior information to the center of the slice in each slice i Setting a certain threshold value as R + T, T as elastic threshold value, R as radius of slice circle, if D i Satisfy D i And if the value is less than or equal to R + T, the contact point is obtained through extraction, the schematic diagram of the contact point is calculated by utilizing structural information in the slice in fig. 6, the set of all the contact points is the contact spot, and the elastic threshold value T is 0.001-0.005 times of the radius R.
S9, the length of the long axis and the length of the short axis of the contact patch are measured according to the coordinates of the center point of the contact patch.
The maximum value of the component (i.e., abscissa) of the coordinates of all the points in the contact patch in the X-axis direction is X max Minimum value of X min And then the long axis length XLEN of the contact spot is calculated as XLEN ═ X max -X min . Similarly, the minor axis length YLen can be obtained by the component (i.e. ordinate) of the coordinates of all points in the contact patch in the Y-axis direction, and the calculation formula is that YLen is Y max -Y min . Fig. 7 is a schematic view of a wheel track contact patch measurement scene, and length of a long axis and a short axis in the X, Y direction of a contact patch is measured according to coordinates of a midpoint of the contact patch to obtain XLen and YLen, which are lengths of the long axis and the short axis of the contact patch.
From the above description of the embodiments, it is clear to those skilled in the art that the present invention can be implemented by software plus necessary general hardware platform. Based on such understanding, the technical solutions of the present invention may be embodied in the form of a software product, which may be stored in a storage medium, such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the method according to the embodiments or some parts of the embodiments.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention should be subject to the protection scope of the claims.
Claims (10)
1. A wheeltrack contact geometric parameter measuring method is characterized by comprising the following steps:
S1, detecting the point cloud of the wheel track by plane elements and cylinder elements;
s2 extracting the elements where the inner side surface of the rim is located from the detection result of the step S1, and roughly dividing the point cloud of the rim part according to the information of the elements where the inner side surface of the rim is located;
s3, performing circular ring element detection on the point cloud of the roughly divided rim part, and extracting corresponding elements of the rim part and the rail bottom part to obtain corresponding element parameters;
s4, calculating the wheel set central point and the rail bottom central point according to the elementary parameters to obtain the wheel axle direction;
s5 measuring wheel set transverse displacement, wheel set side rolling angle, wheel set rolling angle and wheel set sinking and floating amount according to the calculation result of S4, and positioning points belonging to the rail tread in the point cloud;
s6, slicing the point clouds, and removing points belonging to the rail tread in each point cloud slice;
s7, performing slicing circle structure information fitting on the point cloud points with the tread points removed based on a RANSAC algorithm to obtain the circle center and the radius of the slicing circle;
s8, calculating the distance from the point on the rail tread to the center of the slice circle in each slice, wherein the set of points smaller than a certain threshold value is the contact patch;
s9 measures the long and short axial lengths of the contact patch from the coordinates of all points in the contact patch.
2. The method of claim 1, wherein said step S1 is preceded by processing the point cloud data of the wheel-track to establish two different coordinate systems for the wheel set and the track, and establishing O with the midpoint of the axle center of the wheel set as the origin when the wheel set is in the normal operating position w -X w Y w Z w A wheelset coordinate system that moves with the wheelset relative to the track; establishing O by using the central point of the centers of the two side rail surfaces as the origin r -X r Y r Z r A traveling coordinate system, a plane of the rail bottom and a traveling coordinate system O r -X r Y r Z r X in (1) r O r Y r Plane parallel, wherein, X r In the direction of horizontal extension of the rail surface, Y r On the plane of the rail foot and in contact with X r Perpendicular, Z r Perpendicular to the plane of the rail foot, X w 、Y w 、Z w Are each independently of X r 、Y r 、Z r And establishing a coordinate system under a normal condition, wherein each coordinate axis of the wheel set coordinate system and the following coordinate system is pairwise parallel, and only Zr and Zw have a difference value in the Z-axis direction.
3. The method according to claim 2, wherein the S2 includes:
step S21: according to the parallel constraint strips of the inner side distance of the wheel set and the inner side surface of the wheel rimExtracting the inner side surface omega of the double-side rim from the element detection result in S1 1 As the inner side of the rim and omega of the left wheel set 2 As the inner side surface of the wheel rim of the right wheel pair;
step S22: calculating the omega of the plane according to the thickness of the standard wheel rim 1 And Ω 2 Plane omega separated from the thickness of the standard rim 1 ' and omega 2 The plane parameters of' correspond to the front view of the wheel-track point cloud data;
step S23: according to
And
roughly dividing the rim part in the wheel track original point cloud, wherein 
From any point p to plane omega in the wheel track point cloud i The Euclidean distance of (a) is,
from any point p to plane omega in the wheel track point cloud i The euclidean distance of' i ═ 1, 2.
4. The method according to claim 2, wherein the extracting of the rim portion and the rail bottom portion corresponding cells in S3 and obtaining corresponding cell parameters comprises:
step S31: the extraction of the rail bottom element is carried out according to the following three conditions: 1) the type of the model element to be extracted is a plane; 2) extracting normal vector and Z of plane r The directions are parallel; 3) extracting the minimum z coordinate value of the plane;
step S32: the elements of the rim portion are extracted according to the following three conditions: 1) the extracted primitive type is a circular ring primitive, comprising
Eight parameters, wherein,
the method comprises the steps that three parameters are included in the direction of a circular ring shaft, Center is the Center of the circular ring shaft and comprises the three parameters, radiuMax is the large radius of the circular ring and comprises one parameter, and radiuMin is the small radius of the circular ring and comprises one parameter; 2) a wheel rim radius; 3) normal vector of extracted circular ring element To
Normal to the inner side of the rim
Are parallel and satisfy
Where θ is a set threshold.
5. The method according to claim 2, wherein the calculating wheel set center point and rail base center point according to the primitive parameters to obtain wheel axle direction comprises:
S41: calculating the center of the wheel pair according to the center point coordinate of the intersection point of the straight line passing through the center point of the circular ring where the wheel rim is located and the inner side faces of the two wheel rims along the direction of the wheel shaft;
s42: calculating the center of the rail bottom according to the center point coordinate of the center point of the plane bounding box of the rail plane;
s43: the normal vector of the plane element of the inner side surface of the wheel rim is the wheel axle direction.
6. The method according to claim 2, wherein the S5 includes:
s51: measuring the transverse displacement of the wheel set, wherein the center of the wheel set and the center point of the rail bottom are in Y r Spacing in the direction;
s52: measuring the sinking and floating amount of the wheel set, wherein the center of the wheel set and the center point of the rail bottom are in Z r Spacing and initial position Z in axial direction r Transverse shaft movementA difference in amount;
s53: measuring the side rolling angle of the wheel set, and taking the axis of the wheel set as O in a following coordinate system r -X r Y r Z r X in (1) r O r Y r The included angle of the plane;
s54: measuring the pan angle of the wheel pair, as the perpendicular to the wheel axis relative to X r The angle of the axes.
7. The method of claim 2, wherein in step S5, the points belonging to the rail tread in the point cloud are located according to a priori structural information of the rail, the priori structural information being determined according to railroad industry standards.
8. The method according to claim 2, wherein the S6 includes:
S61: estimating the point cloud density rho by adopting a nuclear density estimation method, wherein a slice thickness delta calculation formula is shown as the following formula (1);
δ＝M·ρ (1)
wherein M is an empirical value of 4-8;
s62: the coordinate system of the point cloud slice is that the X axis is along the central line direction of the point cloud at the bottom, the Y axis is parallel to the central line direction of the point cloud at the upper part of the combined point cloud, the Z axis is vertical to the bottom surface of the point cloud, and the Y axis direction is selected as the point cloud slice direction;
s63: and projecting the point set in each point cloud slice to the central plane to realize the dimension reduction processing of local data.
9. The method according to claim 2, wherein the S7 includes:
s71: projecting points in each point cloud piece to a point cloud piece central plane to obtain a slice projection point set, then removing points on a bottom tread in the point cloud slice by using structure prior information, and only keeping point clouds of a wheel to perform parameter fitting of a circle;
s72: randomly selecting three non-collinear points in the slice projection point set as three points capable of forming a circle model, calculating the circle center and the radius of the circle model, calculating the absolute value of the difference between the distances from the other points except the three points of the circle model in the projection point set to the circle center of the circle model and the radius of the circle model, taking the points with the absolute values of the differences smaller than a set threshold value as local points of the current circle model, taking the circle model obtained when the number of the local points reaches a set proportion as a slice circle, and taking the circle center O and the radius R as the circle center and the radius of the slice circle.
10. The method according to claim 2, wherein the S8 includes: calculating the distance D from a point in the bottom rail tread in each slice to the center of the slice i Setting a certain threshold value as R + T, T as elastic threshold value, R as radius of slice circle, if D i Satisfy D i And if the radius is less than or equal to R + T, the contact points are obtained by extraction, and the set of all the contact points is the contact spots, wherein the elastic threshold value T is 0.001-0.005 times of the radius R.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110752175.8A CN113446939B (en) | 2021-07-02 | 2021-07-02 | Wheel-rail contact geometric parameter measuring method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110752175.8A CN113446939B (en) | 2021-07-02 | 2021-07-02 | Wheel-rail contact geometric parameter measuring method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113446939A CN113446939A (en) | 2021-09-28 |
| CN113446939B true CN113446939B (en) | 2022-07-29 |
Family
ID=77815026
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110752175.8A Active CN113446939B (en) | 2021-07-02 | 2021-07-02 | Wheel-rail contact geometric parameter measuring method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113446939B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114087989B (en) * | 2021-11-19 | 2023-09-22 | 江苏理工学院 | Method and system for measuring three-dimensional coordinates of circle center of positioning hole of automobile cylinder workpiece |
| CN114426039B (en) * | 2022-03-14 | 2023-06-30 | 浙江师范大学 | Point cloud data processing method based on urban rail train wheel tread |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108657222A (en) * | 2018-05-03 | 2018-10-16 | 西南交通大学 | Railroad track gauge and horizontal parameters measurement method based on vehicle-mounted Lidar points cloud |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6909514B2 (en) * | 2002-11-18 | 2005-06-21 | Beena Vision Systems, Inc. | Wheel profile inspection apparatus and method |
| CN101117127A (en) * | 2007-08-21 | 2008-02-06 | 西南交通大学 | Method for describing railway wheeltrack spacing dynamic connecting state |
| CN102607439A (en) * | 2012-02-17 | 2012-07-25 | 上海交通大学 | System and method for carrying out on-line monitoring on railway wheel-rail contact relationship on basis of structured light |
| JP2016097868A (en) * | 2014-11-25 | 2016-05-30 | 公益財団法人鉄道総合技術研究所 | Method of numerically identifying contact position between wheel of running railway vehicle and rail, and automatic identifying device of contact position between wheel and rail |
| CN107097807A (en) * | 2017-03-27 | 2017-08-29 | 北京交通大学 | A kind of measurement system of high-speed railway wheel track dynamic Contact state |
| JP6932099B2 (en) * | 2018-03-28 | 2021-09-08 | 公益財団法人鉄道総合技術研究所 | Rail unevenness measurement method and its equipment |
| EP3575179A1 (en) * | 2018-06-01 | 2019-12-04 | ALSTOM Transport Technologies | Train wheel measurement process, and associated system |
| CN110425977B (en) * | 2019-07-10 | 2021-03-26 | 上海工程技术大学 | Contact rail geometric parameter detection system and detection method |
| CN110667643B (en) * | 2019-09-25 | 2020-10-23 | 北京交通大学 | Laser detection system and method for wheel-rail contact state and wheel tread fault |
-
2021
- 2021-07-02 CN CN202110752175.8A patent/CN113446939B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108657222A (en) * | 2018-05-03 | 2018-10-16 | 西南交通大学 | Railroad track gauge and horizontal parameters measurement method based on vehicle-mounted Lidar points cloud |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113446939A (en) | 2021-09-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113446939B (en) | Wheel-rail contact geometric parameter measuring method | |
| CN102866672B (en) | Online detecting method of numerical control machining middle state of plane structural member | |
| CN109653045B (en) | Track gauge measuring method and device | |
| JPWO2016076307A1 (en) | Flange wear measurement method for railway vehicle wheels | |
| CN107101594B (en) | Method for extracting lowest point of wheel flange of wheel track wheel space | |
| CN104878667B (en) | The steel rail grinding method in changing points area | |
| CN107685748B (en) | Train wheel dimension parameter online detection method based on laser displacement sensor | |
| CN110298853B (en) | Visual inspection method for surface difference | |
| CN111391881B (en) | Method and device for detecting fastening state of steel rail fastener | |
| CN102901457A (en) | Dynamic measurement method and system for train wheel diameter | |
| AU2022349296A1 (en) | Road surface technical condition detection method and device based on three-dimensional contour | |
| CN114332212B (en) | Track superelevation and front-back height detection method based on vehicle-mounted mobile laser point cloud | |
| CN112164081B (en) | Vehicle-mounted LiDAR point cloud railway cross section contour extraction method | |
| CN109353372A (en) | A kind of tramcar wheelset profile on-line monitoring system and method | |
| CN106112698A (en) | The method of non-pulling wheel lathe wheel processing front and rear profile circumference accuracy detection | |
| CN110232388A (en) | A method of identifying honeycomb side from honeycomb core surface measurement data | |
| CN115112044A (en) | Wheel set size measurement method based on light spot cloud data of multi-line structure | |
| CN104590314A (en) | Device and method for measuring diameter of urban rail vehicle wheel based on multiple sensors | |
| CN113504545A (en) | Contact network data detection method based on laser radar | |
| CN107451095B (en) | Urban rail vehicle wheel pair curve fitting method | |
| CN113091835A (en) | Network rail tunnel comprehensive detection device and method based on edge calculation and big data | |
| Li et al. | Robust registration of rail profile and complete detection of outliers in complex field environment | |
| Chen et al. | Vision-Based Position Deviation Measurement of Rail Corrugation Chord Measuring Points Under Bi-Linear Laser Assistance | |
| CN112581521B (en) | Method for extracting central line of magnetic suspension track | |
| Ma et al. | Metro wheel set size measurement method based on CDD image recognition |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |