CN117150612A - Rail grinding reference profile design method based on multi-objective optimization - Google Patents
Rail grinding reference profile design method based on multi-objective optimization Download PDFInfo
- Publication number
- CN117150612A CN117150612A CN202311060358.9A CN202311060358A CN117150612A CN 117150612 A CN117150612 A CN 117150612A CN 202311060358 A CN202311060358 A CN 202311060358A CN 117150612 A CN117150612 A CN 117150612A
- Authority
- CN
- China
- Prior art keywords
- rail
- wheel
- profile
- abrasion
- vehicle
- 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
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000013461 design Methods 0.000 title claims abstract description 36
- 238000005457 optimization Methods 0.000 title claims abstract description 19
- 229910000831 Steel Inorganic materials 0.000 claims abstract description 83
- 239000010959 steel Substances 0.000 claims abstract description 83
- 238000005299 abrasion Methods 0.000 claims abstract description 77
- 238000001514 detection method Methods 0.000 claims abstract description 13
- 238000005498 polishing Methods 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims description 44
- 239000000463 material Substances 0.000 claims description 16
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 230000002068 genetic effect Effects 0.000 claims description 13
- 230000001133 acceleration Effects 0.000 claims description 7
- 238000012216 screening Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 230000004927 fusion Effects 0.000 claims description 5
- 238000011160 research Methods 0.000 claims description 5
- 230000003137 locomotive effect Effects 0.000 claims description 4
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000035939 shock Effects 0.000 claims description 3
- 238000012938 design process Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 10
- 230000002829 reductive effect Effects 0.000 description 5
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000003993 interaction Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 101100346171 Arabidopsis thaliana MORC3 gene Proteins 0.000 description 1
- 101100168604 Candida albicans (strain SC5314 / ATCC MYA-2876) CRH12 gene Proteins 0.000 description 1
- 101100168607 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) UTR2 gene Proteins 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Structural Engineering (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Machines For Laying And Maintaining Railways (AREA)
Abstract
The invention discloses a rail grinding reference profile design method based on multi-objective optimization, which relates to the technical field of rail profile design. The rail polishing reference profile design method provided by the invention is based on the multi-objective optimized reference profile design, can comprehensively consider factors such as vehicle body dynamics performance, vehicle body stability, wheel-rail contact force, rail abrasion and the like, ensures the safety and comfort of a line, simultaneously reduces the field polishing workload as much as possible, and prolongs the service life of the rail. The whole steel rail profile detection and design process realizes full automation without a large amount of field manual work.
Description
Technical Field
The invention relates to the technical field of steel rail profile design, in particular to a steel rail polishing reference profile design method based on multi-objective optimization.
Background
Abnormal abrasion of the wheel rail is one of important diseases of a railway curve rail, and causes the service life of the rail to be reduced, the maintenance cost of the rail and the wheel is increased, and even the safety and the comfort of the running of a vehicle are affected. The good profile matching of the wheel rail is obtained, so that the contact geometric relationship of the wheel rail can be effectively improved, the adverse abrasion between the wheel rails is reduced, and the contact stress of the wheel rail is reduced.
Many numerical studies on optimization of the contour of the vehicle are being carried out in China, and with the progress of time, students are continually striving to find reasonable wheel treads so as to design the treads to achieve the set targets. In 2000, the design method of the high-order curve abrasion tread is provided by the Liuzhujun on the basis of Chebyshev value approximation method. The new tread optimized by the method can improve the curve passing performance and ensure the stability in a safe range. 2002. Shen Gang and She Zhisen have proposed a wheel tread optimization method based on the reverse of the contact angle curve. And obtaining a desired contact angle curve through a minimum distance searching method, and expanding the partially optimized new profile according to the original profile, so as to obtain a new tread meeting the contact characteristic. After that, two persons develop the method to the optimal design of the independent wheel set of the low-floor vehicle, and the influence of different contact angle curves on the profile of the vehicle is analyzed. After years Shen Gang and Zhong Xiaobo have also proposed a method of optimizing the tread of a wheel that targets the wheel diameter difference and the contact angle difference as a function of the amount of wheel set traverse. The result shows that the algorithm based on the wheel diameter difference function is suitable for tread appearance design of the rigid wheel set, and the algorithm based on the contact angle difference is suitable for tread appearance design of the independent wheel set. 2008. In the years, zhang Jian and the like try to select a proper region of the profile of the steel rail to expand to optimize the tread of the wheel, and the tread of the wheel after optimization can effectively reduce the contact stress.
The optimization research of the vehicle profile should involve various factors such as vehicle dynamics, wheel track contact geometry, wheel track material characteristics and the like, so that partial scholars in the last ten years adopt a multi-target optimization design method. 2011. In the year, the secondary-secondary set proposes a multi-objective numerical solution method based on Gaussian radial basis function response surfaces by taking a plurality of circular arc center coordinates and radiuses as design variables, determining a plurality of optimized dynamic indexes as a plurality of optimized targets according to different vehicle types and lines, taking the aspects of wheel geometric shapes, vehicle running safety and stability, material strength requirements and the like into consideration as constraint conditions, optimizing the wheel treads of CRH2 high-speed motor train units and HXD2 electric locomotives, and obtaining expected target effects. In 2016, cui Meng established a tread optimization model using a quadratic rational Bezier curve description method. The nonlinear critical speed and the abrasion index are taken as optimization targets, derailment coefficients, lateral forces of wheel tracks and vertical forces of wheel tracks are taken as constraint conditions, the generated control point coordinates are taken as design variables, and the optimal wheel tread is obtained based on a NSGA-II genetic algorithm and an Isight Pareto multi-target optimization solving method.
In summary, many studies on optimization of the profile of the wheel are focused in China, and methods proposed for different vehicle types and line conditions are different, wherein the methods comprise a reverse solution method based on contact characteristics of the wheel, such as a contact angle curve reverse-push method, a rolling circle radius difference reverse-push method and an equivalent taper curve reverse-push method. And a rail profile expansion method, a rail normal clearance inverse method and the like for achieving the desired rail contact conformal. In addition, in recent years, there are many wheel profile design methods using different dynamic indexes of a vehicle as a plurality of optimization targets, but the existing multi-target optimization lacks knowledge of mutual influence between sub-targets, so that it is difficult to ensure that the design result is globally optimal, and further intensive research is needed.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention discloses a rail grinding reference profile design method based on multi-objective optimization.
In order to achieve the above purpose, the present invention adopts the technical scheme that:
a rail grinding reference profile design method based on multi-objective optimization comprises the following steps:
s1, researching the basic information of a line and a vehicle, and establishing a vehicle dynamics model according to the research result;
preferably, in the step S1, the investigated line information includes a rail type, a curve radius, and a rail bottom slope breaking; the basic information of the investigated vehicles comprises the types of wheel treads, locomotive and carriage models, a vehicle body structure, a bogie structure and a shock absorption structure.
In the steps, the purpose of establishing the vehicle dynamics model is to simulate the running conditions of the vehicle under different profile conditions, so as to calculate the derailment coefficient, the comfort, the wheel track contact stress and the like, and the vehicle dynamics model is established by adopting SIMPACK multi-body dynamics simulation software.
S2, detecting the profile of the steel rail of the investigated line;
in the above steps, the steel rail profile detection means that the full section profile data of the steel rail is detected, the steel rail profile detection is carried out by adopting a minirof steel rail profile instrument, the purpose of carrying out steel rail profile detection on an investigated line is to obtain the current situation of the steel rail profile, and the abrasion condition of the steel rail profile is known so as to design the steel rail profile more in line with site conditions.
S3, screening and fusing the detected rail profile to obtain a representative profile of the section of line;
preferably, the step S3 includes the steps of:
s31, comparing the detected steel rail profile with a standard steel rail profile, and screening out a profile with profile points of 400-1000 and steel rail abrasion less than 2 mm;
s32, carrying out average fusion on all the screened profiles to obtain the representative profile of the section of line.
In the steps, firstly, the detected steel rail profile is compared with the standard steel rail profile, and only the profile with the number of profile points being less than 1000 and more than 400 and the steel rail abrasion being less than 2mm is reserved. And then carrying out average fusion on all the retained profiles, thereby obtaining the representative profile of the section.
For the representative profile of the steel rail, the existing mode at present collects one to two steel rail profiles to represent the profile of the whole line according to manual experience, and the mode has high randomness and cannot represent the whole line condition. Because of certain difference of the rail profiles of each line, in order to obtain a more representative rail profile, the invention needs to acquire the line profile in a high density, then screen the profile to acquire abnormal data, and fuse the screened and qualified profile, wherein the fused profile can represent most of the whole conditions.
S4, taking the representative profile as an initial reference profile, substituting the initial reference profile into a vehicle dynamics model to calculate vehicle dynamics performance, wheel rail contact performance and steel rail abrasion, and obtaining a plurality of parameters;
preferably, in the step S4, the vehicle dynamics calculation includes derailment coefficient calculation and vehicle running stability calculation, and the vehicle running stability calculation is:
wherein A is vehicle body vibration acceleration, F is vibration frequency, and F (F) is frequency correction coefficient.
In the present invention, the derailment coefficient is a ratio (Q/P) of a lateral force Q and a vertical force P acting on the wheel at a certain time, and is required to be 1 or less.
The vehicle running stability mainly evaluates that various vibrations generated during the running of the train affect the riding comfort of passengers and the integrity of goods. For passenger vehicles, the stability index is required to be less than 3.5, and for freight lines, the stability index is required to be less than 4.
Preferably, in the step S4, the wheel-rail contact performance calculation uses a hertz contact theory to calculate the wheel-rail contact stress.
Preferably, in the step S4, the rail abrasion calculation adopts an Archard abrasion calculation model, and the rail abrasion is:
wherein V is wear For the material wear volume, K, of the contact area of the wheel rail w Is a dimensionless wear coefficient, F n S is the relative sliding distance of the contact area of the wheel rail, and H is the hardness of the material.
In the invention, in a wheel and rail abrasion calculation model, an Archard abrasion model is a sliding abrasion model for researching the abrasion mechanism of materials, the model considers that abrasion among materials only occurs in a contact spot sliding area, the abrasion of the materials is in direct proportion to the pressure and the sliding distance in the contact area, and an abrasion calculation formula is shown in the formula, wherein the hardness H of the materials is generally consistent with the hardness of a default wheel and a steel rail.
Preferably, in the step S4, the plurality of parameters include derailment coefficient, wheel load shedding rate, wheel set attack angle, vehicle body lateral stability index, vehicle body vertical stability index, wheel rail contact stress and rail total abrasion.
In the present invention, the derailment coefficient is a ratio (Q/P) of a lateral force Q and a vertical force P acting on a wheel at a certain moment, which is used to evaluate the running safety of a vehicle, and when the derailment coefficient is too large, the derailment phenomenon of the vehicle is liable to occur.
The wheel weight load shedding rate is delta P/P, wherein delta P is the wheel weight load shedding amount of the wheel at the load shedding side, P is the average static wheel weight of the wheel at the load shedding side and the load increasing side, and the wheel weight load shedding rate and the derailment coefficient are used for evaluating the running safety of the vehicle and preventing the derailment phenomenon of the vehicle.
The wheel set attack angle refers to an included angle formed between the tangential direction of a steel rail at the contact point of a wheel rail and the axis of the wheel shaft, and is used for reflecting the curve passing capability and difficulty of the vehicle.
The vehicle transverse stability index refers to the magnitude of acceleration in the vehicle transverse direction, which is used to judge the magnitude of transverse vibration during the running process of the vehicle.
The vehicle vertical stability index refers to the vertical acceleration of the vehicle and is used for judging the vertical vibration in the running process of the vehicle.
The wheel rail contact stress refers to a vertical force generated in a contact area between a rail and a wheel, and the wheel rail contact stress not only affects the running safety of a vehicle, but also can cause the abrasion of a steel rail to be increased.
The total abrasion of the steel rail refers to the phenomenon that the abrasion of the steel rail head is generated by the friction between the wheels and the steel rail, the total abrasion = vertical abrasion + side abrasion is 0.5, the vertical abrasion is the abrasion of the steel rail in the vertical direction, the side abrasion is the abrasion of the steel rail side surface caused by the interaction with the wheels, the profile is changed, and the side abrasion is used for reflecting the serious abrasion condition of the steel rail.
S5, carrying out normalization processing on a plurality of parameters, and substituting the parameters into an objective function according to weights of different parameters to obtain a target value;
preferably, in the step S5, the weights of the different parameters are obtained according to analysis and statistics of the historical detection data and the rail profile.
Preferably, in the step S5, the objective function is:
wherein,is a target value; />Is the weight coefficient of the parameter, wherein W 1 -W 7 The weight coefficient of the derailment coefficient, the weight coefficient of the wheel load shedding rate, the weight coefficient of the wheel set attack angle, the weight coefficient of the vehicle body transverse stability index, the weight coefficient of the vehicle body vertical stability index, the weight coefficient of the wheel track contact stress and the weight coefficient of the total abrasion of the steel rail are respectively; />Is a parameter, wherein a 1 -a 7 Respectively the derailment coefficient, the wheel load relief rate, the wheel set attack angle, the vehicle body transverse stability index, the vehicle body vertical stability index, the wheel rail contact stress and the total abrasion of the steel rail。
The function of the above objective function is: the influences of different profile designs on derailment coefficients, wheel load relief rates, wheel set attack angles, vehicle transverse stability indexes, vehicle vertical stability indexes, wheel rail contact stress and total abrasion of steel rails can be comprehensively considered.
S6, adopting a genetic algorithm to select, intersect and mutate the steel rail profile, returning to the step S4, and iteratively calculating an objective function until the objective function reaches a convergence standard to obtain the steel rail polishing reference profile.
In the invention, each steel rail profile is subjected to vehicle dynamics performance calculation, wheel-rail contact performance calculation and steel rail abrasion calculation to obtain 7 parameter values, namely a1-a7, and the 7 parameters are substituted into an objective function to obtain an A value. However, the A value is not necessarily the optimal solution, so that a large number of steel rail profiles are generated through a genetic algorithm, and a large number of A values can be obtained through calculation of vehicle dynamics performance, wheel-rail contact performance and steel rail abrasion calculation, and the optimal solution of the objective function, namely the most suitable steel rail profile, namely the steel rail grinding reference profile, is obtained through comparison of the A values.
The convergence of genetic algorithms generally refers to the fact that a new population generated by genetic iteration reaches a certain steady state, and the evolution process is regarded as a random process and examined.
Preferably, in the step S6, the global convergence of the genetic algorithm is:
wherein S is the individual space, S N Is population space, P is S N Probability distribution on X t ={X 1 (t),X 2 (t),...,X M (t) } is a population of the t generation, X i (t) is the population of individuals of the t generation, i=1, 2,..m, M is the population size, Q is set t =max{f(X i (t) |i=1, 2,..m) } is the maximum value of the fitness function value of the individuals in the population, and y=max { f (X) |x e S } is the global optimal solution.
The invention has the beneficial effects that:
according to the rail grinding reference profile design method provided by the invention, before rail grinding, automatic statistical analysis can be performed on the rail profile of the detection line, and the most typical representative profile is calculated, so that the reference profile can be designed according to the latest line profile, automatic collection and calculation can be realized in the whole process, and the workload of manual collection and calculation is reduced.
The rail polishing reference profile design method provided by the invention is based on the multi-objective optimized reference profile design, can comprehensively consider factors such as vehicle body dynamics performance, vehicle body stability, wheel-rail contact force, rail abrasion and the like, ensures the safety and comfort of a line, simultaneously reduces the field polishing workload as much as possible, and prolongs the service life of the rail. The whole steel rail profile detection and design process realizes full automation without a large amount of field manual work.
Drawings
FIG. 1 is a flow chart of the present invention.
Detailed Description
The conception, specific structure, and technical effects produced by the present invention will be clearly and completely described below with reference to the embodiments and the drawings to fully understand the objects, features, and effects of the present invention.
Example 1
A rail grinding reference profile design method based on multi-objective optimization, as shown in fig. 1, comprises the following steps:
s1, researching the basic information of a line and a vehicle, and establishing a vehicle dynamics model according to the research result;
s2, detecting the profile of the steel rail of the investigated line;
s3, screening and fusing the detected rail profile to obtain a representative profile of the section of line;
s4, taking the representative profile as an initial reference profile, substituting the initial reference profile into a vehicle dynamics model to calculate vehicle dynamics performance, wheel rail contact performance and steel rail abrasion, and obtaining a plurality of parameters;
s5, carrying out normalization processing on a plurality of parameters, and substituting the parameters into an objective function according to weights of different parameters to obtain a target value;
s6, adopting a genetic algorithm to select, intersect and mutate the steel rail profile, returning to the step S4, and iteratively calculating an objective function until the objective function reaches a convergence standard to obtain the steel rail polishing reference profile.
In this embodiment, before each rail polishing, automatic statistical analysis can be performed on the rail profile of the detection line, and the most typical representative profile is calculated, so that the design of the reference profile according to the latest line profile is ensured, the automatic collection and calculation can be realized in the whole process, and the workload of manual collection and calculation is reduced.
Example 2
The embodiment further describes step S1 on the basis of embodiment 1, wherein in step S1, the investigated line information includes rail type, curve radius, rail bottom slope breaking; the basic information of the investigated vehicles comprises the types of wheel treads, locomotive and carriage models, a vehicle body structure, a bogie structure and a shock absorption structure.
In the S1 step, the purpose of establishing a vehicle dynamics model is to be used for simulating the running conditions of a vehicle under different profile conditions, so that derailment coefficients, comfortableness, wheel-rail contact stress and the like are calculated, and the vehicle dynamics model is established by adopting SIMPACK multi-body dynamics simulation software.
Example 3
In this embodiment, the step S2 is further described on the basis of embodiment 2, in the step S2, the rail profile detection means that the rail full section profile data is detected, the rail profile detection is performed by using a minitrof rail profile meter, and the purpose of performing rail profile detection on the investigated line is to obtain the current situation of the rail profile, and understand the rail profile abrasion condition so as to design the rail profile more conforming to the site conditions.
Example 4
This embodiment further describes step S3 on the basis of embodiment 3, where step S3 includes the following steps:
s31, comparing the detected steel rail profile with a standard steel rail profile, and screening out a profile with profile points of 400-1000 and steel rail abrasion less than 2 mm;
s32, carrying out average fusion on all the screened profiles to obtain the representative profile of the section of line.
In the steps, firstly, the detected steel rail profile is compared with the standard steel rail profile, and only the profile with the number of profile points being less than 1000 and more than 400 and the steel rail abrasion being less than 2mm is reserved. And then carrying out average fusion on all the retained profiles, thereby obtaining the representative profile of the section.
For the representative profile of the steel rail, the existing mode at present collects one to two steel rail profiles to represent the profile of the whole line according to manual experience, and the mode has high randomness and cannot represent the whole line condition. Because of certain difference of the rail profiles of each line, in order to obtain a more representative rail profile, the invention needs to acquire the line profile in a high density, then screen the profile to acquire abnormal data, and fuse the screened and qualified profile, wherein the fused profile can represent most of the whole conditions.
Example 5
The present embodiment further describes step S4 on the basis of embodiment 4, wherein in step S4, the vehicle dynamics performance calculation, the wheel-rail contact performance calculation, the rail abrasion calculation, and several parameters are specifically as follows:
1. vehicle dynamics calculation
In the step S4, the vehicle dynamics performance calculation includes derailment coefficient calculation and vehicle running stability calculation, where the vehicle running stability calculation is:
wherein A is vehicle body vibration acceleration, F is vibration frequency, and F (F) is frequency correction coefficient.
In the present invention, the derailment coefficient is a ratio (Q/P) of a lateral force Q and a vertical force P acting on the wheel at a certain time, and is required to be 1 or less.
The vehicle running stability mainly evaluates that various vibrations generated during the running of the train affect the riding comfort of passengers and the integrity of goods. For passenger vehicles, the stability index is required to be less than 3.5, and for freight lines, the stability index is required to be less than 4.
The frequency correction coefficient F (F) is specifically as follows:
TABLE 1 frequency correction coefficient
2. Wheel-rail contact performance calculation
In the step S4, the wheel-rail contact performance calculation adopts the hertz contact theory to calculate the wheel-rail contact stress, and specifically includes the following steps:
table 2 Hertz formula
3. Rail wear calculation
In the step S4, the steel rail abrasion calculation adopts an Archard abrasion calculation model, and the steel rail abrasion is as follows:
wherein V is wear For the material wear volume, K, of the contact area of the wheel rail w Is a dimensionless wear coefficient, F n S is the relative sliding distance of the contact area of the wheel rail, and H is the hardness of the material.
In the invention, in a wheel and rail abrasion calculation model, an Archard abrasion model is a sliding abrasion model for researching the abrasion mechanism of materials, the model considers that abrasion among materials only occurs in a contact spot sliding area, the abrasion of the materials is in direct proportion to the pressure and the sliding distance in the contact area, and an abrasion calculation formula is shown in the formula, wherein the hardness H of the materials is generally consistent with the hardness of a default wheel and a steel rail.
4. Several parameters
In the step S4, the parameters comprise derailment coefficient, wheel load shedding rate, wheel set attack angle, vehicle body transverse stability index, vehicle body vertical stability index, wheel rail contact stress and steel rail total abrasion.
In the present invention, the derailment coefficient is a ratio (Q/P) of a lateral force Q and a vertical force P acting on a wheel at a certain moment, which is used to evaluate the running safety of a vehicle, and when the derailment coefficient is too large, the derailment phenomenon of the vehicle is liable to occur.
The wheel weight load shedding rate is delta P/P, wherein delta P is the wheel weight load shedding amount of the wheel at the load shedding side, P is the average static wheel weight of the wheel at the load shedding side and the load increasing side, and the wheel weight load shedding rate and the derailment coefficient are used for evaluating the running safety of the vehicle and preventing the derailment phenomenon of the vehicle.
The wheel set attack angle refers to an included angle formed between the tangential direction of a steel rail at the contact point of a wheel rail and the axis of the wheel shaft, and is used for reflecting the curve passing capability and difficulty of the vehicle.
The vehicle transverse stability index refers to the magnitude of acceleration in the vehicle transverse direction, which is used to judge the magnitude of transverse vibration during the running process of the vehicle.
The vehicle vertical stability index refers to the vertical acceleration of the vehicle and is used for judging the vertical vibration in the running process of the vehicle.
The wheel rail contact stress refers to a vertical force generated in a contact area between a rail and a wheel, and the wheel rail contact stress not only affects the running safety of a vehicle, but also can cause the abrasion of a steel rail to be increased.
The total abrasion of the steel rail refers to the phenomenon that the abrasion of the steel rail head is generated by the friction between the wheels and the steel rail, the total abrasion = vertical abrasion + side abrasion is 0.5, the vertical abrasion is the abrasion of the steel rail in the vertical direction, the side abrasion is the abrasion of the steel rail side surface caused by the interaction with the wheels, the profile is changed, and the side abrasion is used for reflecting the serious abrasion condition of the steel rail.
Example 6
The embodiment further describes step S5 on the basis of embodiment 5, wherein in step S5, weights of different parameters are obtained according to analysis and statistics on historical detection data and rail profile. In this embodiment, the weights of the 7 parameters are as follows:
TABLE 3 parameter weights
In the step S5, the objective function is:
wherein,is a target value; />Is the weight coefficient of the parameter, wherein W 1 -W 7 The weight coefficient of the derailment coefficient, the weight coefficient of the wheel load shedding rate, the weight coefficient of the wheel set attack angle, the weight coefficient of the vehicle body transverse stability index, the weight coefficient of the vehicle body vertical stability index, the weight coefficient of the wheel track contact stress and the weight coefficient of the total abrasion of the steel rail are respectively; />Is a parameter, wherein a 1 -a 7 The derailment coefficient, the wheel load relief rate, the wheel set attack angle, the vehicle body transverse stability index, the vehicle body vertical stability index, the wheel rail contact stress and the total abrasion of the steel rail are respectively.
The function of the above objective function is: the influences of different profile designs on derailment coefficients, wheel load relief rates, wheel set attack angles, vehicle transverse stability indexes, vehicle vertical stability indexes, wheel rail contact stress and total abrasion of steel rails can be comprehensively considered.
Example 7
In this embodiment, the step S6 is further described on the basis of embodiment 6, in this embodiment, each rail profile performs calculation of vehicle dynamics performance, wheel-rail contact performance and rail wear calculation, so as to obtain 7 parameter values, namely a1-a7, and the 7 parameters are substituted into the objective function to obtain an a value. However, the A value is not necessarily the optimal solution, so that a large number of steel rail profiles are generated through a genetic algorithm, and a large number of A values can be obtained through calculation of vehicle dynamics performance, wheel-rail contact performance and steel rail abrasion calculation, and the optimal solution of the objective function, namely the most suitable steel rail profile, namely the steel rail grinding reference profile, is obtained through comparison of the A values.
The convergence of genetic algorithms generally refers to the fact that a new population generated by genetic iteration reaches a certain steady state, and the evolution process is regarded as a random process and examined.
In the step S6, the global convergence of the genetic algorithm is:
wherein S is the individual space, S N Is population space, P is S N Probability distribution on X t ={X 1 (t),X 2 (t),...,X M (t) } is a population of the t generation, X i (t) is the population of individuals of the t generation, i=1, 2,..m, M is the population size, Q is set t =max{f(X i (t) |i=1, 2,..m) } is the maximum value of the fitness function value of the individuals in the population, and y=max { f (X) |x e S } is the global optimal solution.
While the embodiments of the present invention have been described in detail, the present invention is not limited to the embodiments described above, and various equivalent modifications and substitutions can be made by those skilled in the art without departing from the spirit of the present invention, and these are intended to be included in the scope of the present invention as defined in the appended claims.
Claims (10)
1. A rail grinding reference profile design method based on multi-objective optimization is characterized by comprising the following steps:
s1, researching the basic information of a line and a vehicle, and establishing a vehicle dynamics model according to the research result;
s2, detecting the profile of the steel rail of the investigated line;
s3, screening and fusing the detected rail profile to obtain a representative profile of the section of line;
s4, taking the representative profile as an initial reference profile, substituting the initial reference profile into a vehicle dynamics model to calculate vehicle dynamics performance, wheel rail contact performance and steel rail abrasion, and obtaining a plurality of parameters;
s5, carrying out normalization processing on a plurality of parameters, and substituting the parameters into an objective function according to weights of different parameters to obtain a target value;
s6, adopting a genetic algorithm to select, intersect and mutate the steel rail profile, returning to the step S4, and iteratively calculating an objective function until the objective function reaches a convergence standard to obtain the steel rail polishing reference profile.
2. The method for designing a rail grinding reference profile according to claim 1, wherein in the step S1, the investigated line information includes rail type, curve radius, rail bottom slope breaking; the basic information of the investigated vehicles comprises the types of wheel treads, locomotive and carriage models, a vehicle body structure, a bogie structure and a shock absorption structure.
3. A rail grinding reference profile design method as set forth in claim 1, wherein said step S3 includes the steps of:
s31, comparing the detected steel rail profile with a standard steel rail profile, and screening out a profile with profile points of 400-1000 and steel rail abrasion less than 2 mm;
s32, carrying out average fusion on all the screened profiles to obtain the representative profile of the section of line.
4. The rail grinding reference profile design method as set forth in claim 1, wherein in the step S4, the vehicle dynamics calculation includes a derailment coefficient calculation and a vehicle running smoothness calculation, the vehicle running smoothness calculation being:
the method comprises the steps of carrying out a first treatment on the surface of the Wherein A is vehicle body vibration acceleration, f is vibration frequencyThe rate, F (F), is the frequency correction coefficient.
5. A rail grinding reference profile design method as set forth in claim 1, wherein in said step S4, said wheel-rail contact performance calculation uses a hertz contact theory to calculate wheel-rail contact stress.
6. The method for designing a rail grinding reference profile as claimed in claim 1, wherein in the step S4, the rail abrasion calculation is implemented by using an Archard abrasion calculation model, and the rail abrasion is:
the method comprises the steps of carrying out a first treatment on the surface of the Wherein V is wear For the material wear volume, K, of the contact area of the wheel rail w Is a dimensionless wear coefficient, F n S is the relative sliding distance of the contact area of the wheel rail, and H is the hardness of the material.
7. The method of claim 1, wherein in step S4, the plurality of parameters includes derailment coefficient, wheel load reduction rate, wheel set angle of attack, vehicle body lateral stability index, vehicle body vertical stability index, wheel rail contact stress, and rail total wear.
8. A rail grinding reference profile design method as set forth in claim 1, wherein in step S5, weights of the different parameters are obtained according to analysis statistics of historical detection data and rail profiles.
9. A rail grinding reference profile design method as set forth in claim 1, wherein in said step S5, said objective function is:
the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Is a target value; />Is the weight coefficient of the parameter, wherein W 1 -W 7 The weight coefficient of the derailment coefficient, the weight coefficient of the wheel load shedding rate, the weight coefficient of the wheel set attack angle, the weight coefficient of the vehicle body transverse stability index, the weight coefficient of the vehicle body vertical stability index, the weight coefficient of the wheel track contact stress and the weight coefficient of the total abrasion of the steel rail are respectively; />Is a parameter, wherein a 1 -a 7 The derailment coefficient, the wheel load relief rate, the wheel set attack angle, the vehicle body transverse stability index, the vehicle body vertical stability index, the wheel rail contact stress and the total abrasion of the steel rail are respectively.
10. A rail grinding reference profile design method as claimed in claim 1, wherein in step S6, the global convergence of the genetic algorithm is:
the method comprises the steps of carrying out a first treatment on the surface of the Wherein S is the individual space, S N Is population space, P is S N Probability distribution on X t ={X 1 (t),X 2 (t),...,X M (t) } is a population of the t generation, X i (t) is the population of individuals of the t generation, i=1, 2,..m, M is the population size, Q is set t =max{f(X i (t) |i=1, 2,..m) } is the maximum value of the fitness function value of the individuals in the population, and y=max { f (X) |x e S } is the global optimal solution.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311060358.9A CN117150612A (en) | 2023-08-22 | 2023-08-22 | Rail grinding reference profile design method based on multi-objective optimization |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311060358.9A CN117150612A (en) | 2023-08-22 | 2023-08-22 | Rail grinding reference profile design method based on multi-objective optimization |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117150612A true CN117150612A (en) | 2023-12-01 |
Family
ID=88883507
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311060358.9A Pending CN117150612A (en) | 2023-08-22 | 2023-08-22 | Rail grinding reference profile design method based on multi-objective optimization |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117150612A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117973086A (en) * | 2024-03-28 | 2024-05-03 | 中铁第四勘察设计院集团有限公司 | Multi-objective optimization-based whole-path steel rail polishing planning method and system |
-
2023
- 2023-08-22 CN CN202311060358.9A patent/CN117150612A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117973086A (en) * | 2024-03-28 | 2024-05-03 | 中铁第四勘察设计院集团有限公司 | Multi-objective optimization-based whole-path steel rail polishing planning method and system |
CN117973086B (en) * | 2024-03-28 | 2024-06-11 | 中铁第四勘察设计院集团有限公司 | Multi-objective optimization-based whole-path steel rail polishing planning method and system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112131678B (en) | Railway steel rail polishing profile design method | |
Persson et al. | Optimisation of railway wheel profiles using a genetic algorithm | |
CN108589449B (en) | Personalized railway steel rail grinding target profile design method | |
CN110222437B (en) | Method and device for evaluating health status of train, and storage medium | |
Jendel | Prediction of wheel profile wear—comparisons with field measurements | |
CN117150612A (en) | Rail grinding reference profile design method based on multi-objective optimization | |
Nia et al. | Prediction of RCF and wear evolution of iron-ore locomotive wheels | |
CN110377986B (en) | Method for predicting side grinding of outer rail of small-radius curve of subway | |
CN109214554A (en) | A kind of railway vehicle wheel under non-pulling wheel operating condition repairs strategic decision-making algorithm to vehicle rotation | |
CN111582656A (en) | Quantitative evaluation method for shallow layer state of high-speed railway steel rail | |
Ignesti et al. | Development of a wear model for the wheel profile optimisation on railway vehicles | |
Markine et al. | An inverse shape design method for railway wheel profiles | |
Lewis et al. | Integrating dynamics and wear modelling to predict railway wheel profile evolution | |
Kou et al. | Influence of crossing wear on rolling contact fatigue damage of frog rail | |
de Paula Pacheco et al. | Optimization of heavy haul railway wheel profile based on rolling contact fatigue and wear performance | |
Lu et al. | Optimization and design of a railway wheel profile based on interval uncertainty to reduce circular wear | |
Molatefi et al. | Advances in Iran railway wheel wear management: a practical approach for selection of wheel profile using numerical methods and comprehensive field tests | |
CN113806883A (en) | Stepless turning method for DIN5573 type subway wheel | |
Lingaitis et al. | The analysis of wear intensity of the locomotive wheel-sets | |
Ngigi et al. | Predictive control strategies used to solve challenges related to modern railway vehicles | |
CN113158327B (en) | Ultrahigh dynamic design method for high-speed railway curve with speed per hour of 400km and above | |
Staśkiewicz et al. | Verification of a tram wheel new profile dynamic behavior | |
CN108755285B (en) | Method for designing and manufacturing fixed frog and fixed frog | |
Staskiewicz et al. | Analysis of wheel/rail interaction for the development of a new tramway wheel profile | |
Nia et al. | An investigation of the iron-ore wheel damages using vehicle dynamics simulation |
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 |