US20230324581A1 - Ballastless track roadbed damage forewarning method and system considering uncertainty - Google Patents
Ballastless track roadbed damage forewarning method and system considering uncertainty Download PDFInfo
- Publication number
- US20230324581A1 US20230324581A1 US18/297,455 US202318297455A US2023324581A1 US 20230324581 A1 US20230324581 A1 US 20230324581A1 US 202318297455 A US202318297455 A US 202318297455A US 2023324581 A1 US2023324581 A1 US 2023324581A1
- Authority
- US
- United States
- Prior art keywords
- roadbed
- surface layer
- damage
- spatial correlation
- worst
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000002344 surface layer Substances 0.000 claims abstract description 55
- 239000002245 particle Substances 0.000 claims abstract description 43
- 239000000463 material Substances 0.000 claims abstract description 24
- 238000001556 precipitation Methods 0.000 claims abstract description 23
- 238000005070 sampling Methods 0.000 claims abstract description 5
- 239000012530 fluid Substances 0.000 claims abstract description 4
- 230000004044 response Effects 0.000 claims description 27
- 239000010410 layer Substances 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 9
- 238000012544 monitoring process Methods 0.000 claims description 8
- 238000005314 correlation function Methods 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 238000012614 Monte-Carlo sampling Methods 0.000 claims description 3
- 238000012821 model calculation Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000007789 sealing Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000013013 elastic material Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- 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
-
- G01V99/005—
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L23/00—Control, warning or like safety means along the route or between vehicles or trains
- B61L23/04—Control, warning or like safety means along the route or between vehicles or trains for monitoring the mechanical state of the route
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L23/00—Control, warning or like safety means along the route or between vehicles or trains
- B61L23/04—Control, warning or like safety means along the route or between vehicles or trains for monitoring the mechanical state of the route
- B61L23/042—Track changes detection
- B61L23/048—Road bed changes, e.g. road bed erosion
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L27/00—Central railway traffic control systems; Trackside control; Communication systems specially adapted therefor
- B61L27/50—Trackside diagnosis or maintenance, e.g. software upgrades
- B61L27/53—Trackside diagnosis or maintenance, e.g. software upgrades for trackside elements or systems, e.g. trackside supervision of trackside control system conditions
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- E—FIXED CONSTRUCTIONS
- E01—CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
- E01B—PERMANENT WAY; PERMANENT-WAY TOOLS; MACHINES FOR MAKING RAILWAYS OF ALL KINDS
- E01B1/00—Ballastway; Other means for supporting the sleepers or the track; Drainage of the ballastway
- E01B1/002—Ballastless track, e.g. concrete slab trackway, or with asphalt layers
-
- 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
Definitions
- the disclosure relates to the filed of railway engineering service safe, in particularly to a ballastless track roadbed damage forewarning method and system considering uncertainty.
- the disclosure provides a ballastless track roadbed damage forewarning method and system considering uncertainty.
- a ballastless track roadbed damage forewarning method considering uncertainty includes the following steps:
- step 1 obtaining roadbed material parameters and precipitation over a period of time
- step 2 establishing a fast Lagrangian analysis of continuum (FLAC)-particle flow code (PFC) model of the ballastless track roadbed, calibrating particle materials of a roadbed surface layer according to the roadbed material parameters obtained in step 1 and obtaining mesoscopic parameters of the particle materials;
- FLAC Lagrangian analysis of continuum
- PFC particle flow code
- step 3 generating a lognormal random field of particle contact friction coefficients, and assigning the lognormal random field to particle contact nodes of the roadbed surface layer in the FLAC-PFC model obtained in step 2;
- step 4 performing Monte Carlo sampling based on a probability density function of the precipitation over the period of time obtained in step 1, and adjusting a fluid domain of the roadbed surface layer according to the samples (i.e., precipitation);
- step 5 increasing a spatial correlation distance ⁇ in a step by step manner, calculating a structure dynamic response, and determining a worst spatial correlation distance ⁇ worst ;
- step 6 calculating a damage probability pf of the roadbed surface layer under the worst spatial correlation distance ⁇ worst obtained in step 5; and outputting alarm information when the damage probability pf exceeds an alarm threshold p, otherwise, quitting.
- the roadbed material parameters in step 1 are obtained according to a field test, and include: base plate parameters, and material parameters of the roadbed surface layer and a roadbed bottom layer; and the base plate parameters include: a density and modulus of elasticity; and the material parameters of the roadbed surface layer and the roadbed bottom layer include: a density, a gradation, modulus of elasticity, a cohesion and a friction angle.
- step 2 particle contacts of the roadbed surface layer adopt a Hertz contact model, and the mesoscopic parameters include: shear modulus, a Poisson's ratio and a friction coefficient; FLAC is a kind of continuous medium mechanics analysis software, which is used to establish a finite difference model of the base plate and a roadbed bottom layer; PFC is a particle flow program built in FLAC, which is used to establish a discrete element model of the roadbed surface layer; and the FLAC-PFC model is to call PFC program through a built-in interface of FLAC, so as to perform a coupling analysis between the finite difference model of the base plate and the roadbed bottom layer established by FLAC and the discrete element model of the roadbed surface layer established by PFC.
- FLAC is a kind of continuous medium mechanics analysis software, which is used to establish a finite difference model of the base plate and a roadbed bottom layer
- PFC is a particle flow program built in FLAC, which is used to establish a discrete element model of the roadbed surface layer
- the generating a lognormal random field of particle contact friction coefficients in step 3 includes:
- a formula of the Markov spatial correlation function in the MATLAB script is as follows:
- ⁇ ⁇ ( ⁇ x , ⁇ y ) exp ( - 2 ⁇ ⁇ x 2 + ⁇ y 2 ⁇ )
- ⁇ represents the correlation coefficient between the two coordinate points of the roadbed surface layer
- ⁇ x and ⁇ y respectively correspond to absolute distances between the two coordinate points in horizontal and vertical directions
- ⁇ represents a spatial correlation length
- obtaining the worst spatial correlation distance ⁇ worst in step 5 includes:
- obtaining the damage probability pf of the roadbed surface layer in step 6, includes:
- a forewarning system by using the above ballastless track roadbed damage forewarning method considering uncertainty includes: an information collection module, a processing module, and a forewarning module;
- the information collection module is configured to collect material parameters of the ballastless track roadbed and real-time precipitation information
- the processing module is configured to load the FLAC-PFC model and calibrate the mesoscopic parameters based on the information collected by the information collection module; load the MATLAB script, generate a random field, and assign the random field to the particle contact nodes; calculate the worst spatial correlation distance ⁇ worst and calculate the damage probability pf; and
- the forewarning module is configured to determine whether the damage probability exceeds the alarm threshold; if damage probability exceeds the alarm threshold, output the alarm information; if the damage probability does not exceed the alarm threshold, exit.
- each of the information collection module, the processing module and the forewarning module is embodied by at least one processor and at least one memory coupled to the at least one processor, and the at least one memory stores programs executable by the at least one processor.
- the randomness of the precipitation in the roadbed area and the uncertainty of the roadbed surface layer material parameters changing with space are considered, and the damage probability of the roadbed surface layer under the action of rain and dynamic load is quantitatively predicted by combining a FLAC-PFC coupling numerical analysis method.
- the disclosure provides a powerful means for damage monitoring of the ballastless track roadbed, and the analysis method has clear flow and strong reliability.
- FIG. 1 illustrates a flowchart of a forewarning method according to an embodiment of the disclosure.
- FIG. 2 illustrates a schematic diagram of a FLAC-PFC coupling model of a roadbed of a high-speed rail ballastless track according to an embodiment of the disclosure; in which, 1 : base plate, 2 : roadbed surface layer, 3 : roadbed bottom layer.
- FIG. 1 a ballastless track roadbed damage forewarning method considering uncertainty provided by an embodiment of the disclosure is illustrated, including the following steps.
- Step 101 obtaining precipitation in a monitoring area for a period of time, and calculating a probability density function of the precipitation using a maximum entropy principle; obtaining material parameters of a base plate 1 , a roadbed surface layer 2 , and a roadbed bottom layer 3 of the ballastless track.
- Step 102 establishing a finite difference model of the base plate 1 and the roadbed bottom layer 3 through fast Lagrangian analysis of continuum (FLAC), assigning values based on the material parameters obtained in step 101, establishing a discrete element model of the roadbed surface layer 2 through particle flow code (PFC), and calibrating mesoscopic parameters of particle materials of the roadbed surface layer 2 based on the material parameters obtained in step 101.
- FLAC Lagrangian analysis of continuum
- PFC particle flow code
- Step 103 extracting spatial positions of particle contact nodes and friction coefficients of particle contacts of the discrete element model of the roadbed surface layer 2 , and importing them into a MATLAB script; setting a spatial correlation length ⁇ of the friction coefficients of the particle contacts of the roadbed surface layer 2 in the MATLAB script, calculating a correlation coefficient ⁇ between the friction coefficients at any two coordinate points by using a Markov spatial correlation function, generating a lognormal random field of the friction coefficients of the particle contacts and giving it to the corresponding particle contact nodes in the FLAC-PFC model, and the Markov spatial correlation function is as follows:
- ⁇ ⁇ ( ⁇ x , ⁇ y ) exp ( - 2 ⁇ ⁇ x 2 + ⁇ y 2 ⁇ )
- ⁇ represents the correlation coefficient between two coordinate points on the roadbed surface layer
- ⁇ x and ⁇ y respectively correspond to absolute distances between two points in horizontal and vertical directions
- ⁇ represents the spatial correlation length
- Step 104 performing Monte Carlo sampling on the precipitation according to its probability density function, and adjusting a corresponding fluid domain of the roadbed surface layer according to the samples (i.e., precipitation).
- Step 105 applying a high-speed train load, and calculating a dynamic response of the base plate; repeating step 103, and increasing ⁇ in a step by step manner. Under each ⁇ i , the random field generation of the friction coefficients of the particle contacts of the roadbed surface layer is repeated for N times, and the precipitation is sampled for N times.
- the high-speed train load is applied and brought into the model calculation, and an average value of N groups of dynamic response of the base plate is taken as a representative value of the dynamic response of the base plate under ⁇ i .
- a worst spatial correlation distance ⁇ worst is determined. Where i is a positive integer.
- Step 106 under ⁇ worst , the random field generation of the friction coefficients of the particle contacts of the roadbed surface layer is repeated for M times, and the precipitation is sampled for M times to obtain M groups of samples, and the dynamic response results of the base plate in each group of samples are calculated by applying the high-speed train load; a safety threshold of the dynamic response according to the standard and a damage alarm threshold p of the roadbed surface layer according to the operation conditions are input.
- the dynamic response of the base plate includes the maximum vibration amplitude, maximum vibration acceleration, maximum dynamic stress and maximum dynamic strain of the base plate. When any dynamic response of the base plate exceeds the specified macro-safety threshold, it is determined that roadbed surface layer is damaged and endangers the superstructure.
- the damage probability pf of the roadbed surface layer is calculated, when pf exceeds the alarm threshold p, the alarm information is sent to the monitoring personnel, otherwise quitting.
- the FLAC-PFC coupling numerical model is established, which includes the base plate 1 of the ballastless track based on an elastic material model, with the parameters including density and modulus of elasticity; the roadbed bottom layer 3 based on a Mohr-Coulomb model, with parameters including density, modulus of elasticity, cohesion, and friction angle; and the roadbed surface layer 2 based on a Hertz contact model, with the parameters including gradation, shear modulus, Poisson's ratio, and friction coefficient.
- the scheme provided by the embodiment of the disclosure takes into account the randomness of precipitation in the subgrade area and the spatial uncertainty of the material parameters of the surface layer of the roadbed, and combines the FLAC-PFC coupling numerical analysis method to predict the damage of the track roadbed under dynamic loads and rain, which helps to ensure the safety of high-speed railway traffic.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Structural Engineering (AREA)
- Architecture (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A ballastless track roadbed damage forewarning method considering uncertainty includes the following steps. Step 1: counting roadbed material parameters and precipitation; Step 2: establishing a FLAC-PFC model of a ballastless track roadbed and calibrating mesoscopic parameters of a roadbed surface layer; Step 3: generating a lognormal random field of particle contact friction coefficients and assigning it to particle contact nodes of the roadbed surface layer; Step 4: perform sampling on the precipitation and adjusting a fluid domain of the roadbed surface layer; Step 5: determining a worst spatial correlation distance in the random field; and Step 6: calculating a damage probability pf of the roadbed surface layer under the worst spatial correlation distance; outputting alarm information when pf exceeds an alarm threshold, otherwise, quitting. The method monitors and gives early warning of damage to the roadbed surface layer under dynamic loads, ensuring driving safety.
Description
- The disclosure relates to the filed of railway engineering service safe, in particularly to a ballastless track roadbed damage forewarning method and system considering uncertainty.
- Chinese high-speed railway already turned from large-scale construction to conserve operation stage. Under the action of traffic cycle load and long-term rainfall, and the durability of sealing layer materials having a certain time limit, it is inevitable that the rain penetrates into the roadbed through the failure of the sealing layer and expansion joint, thus inducing the damage of ballastless track roadbed. The surface layer damage of the roadbed will cause the void between the roadbed and the track structure, and the structure dynamic response will increase significantly and then affect the safety of train operation and passenger comfort, which has become an urgent problem to be solved in the safety field of railway engineering roadbed service.
- Due to the uncertainty of environmental precipitation, the uncertainty of filler parameters and the error brought by construction, the accuracy of traditional deterministic analysis method is low, and it is difficult to achieve the advanced monitoring and early warning of roadbed damage of ballastless track, which may lead to miss the best opportunity for roadbed maintenance and cause serious economic losses. Therefore, it is urgent to develop a method and system for damage forewarning of the ballastless track roadbed considering uncertainty.
- In order to solve the technical problem of the disclosure and overcome the deficiency of the prior art, the disclosure provides a ballastless track roadbed damage forewarning method and system considering uncertainty.
- The technical solutions of the disclosure are as follows:
- A ballastless track roadbed damage forewarning method considering uncertainty includes the following steps:
- step 1: obtaining roadbed material parameters and precipitation over a period of time;
- step 2: establishing a fast Lagrangian analysis of continuum (FLAC)-particle flow code (PFC) model of the ballastless track roadbed, calibrating particle materials of a roadbed surface layer according to the roadbed material parameters obtained in step 1 and obtaining mesoscopic parameters of the particle materials;
- step 3: generating a lognormal random field of particle contact friction coefficients, and assigning the lognormal random field to particle contact nodes of the roadbed surface layer in the FLAC-PFC model obtained in
step 2; - step 4: performing Monte Carlo sampling based on a probability density function of the precipitation over the period of time obtained in step 1, and adjusting a fluid domain of the roadbed surface layer according to the samples (i.e., precipitation);
- step 5, increasing a spatial correlation distance θ in a step by step manner, calculating a structure dynamic response, and determining a worst spatial correlation distance θworst; and
- step 6, calculating a damage probability pf of the roadbed surface layer under the worst spatial correlation distance θworst obtained in step 5; and outputting alarm information when the damage probability pf exceeds an alarm threshold p, otherwise, quitting.
- In an embodiment, the roadbed material parameters in step 1 are obtained according to a field test, and include: base plate parameters, and material parameters of the roadbed surface layer and a roadbed bottom layer; and the base plate parameters include: a density and modulus of elasticity; and the material parameters of the roadbed surface layer and the roadbed bottom layer include: a density, a gradation, modulus of elasticity, a cohesion and a friction angle.
- In an embodiment, in
step 2, particle contacts of the roadbed surface layer adopt a Hertz contact model, and the mesoscopic parameters include: shear modulus, a Poisson's ratio and a friction coefficient; FLAC is a kind of continuous medium mechanics analysis software, which is used to establish a finite difference model of the base plate and a roadbed bottom layer; PFC is a particle flow program built in FLAC, which is used to establish a discrete element model of the roadbed surface layer; and the FLAC-PFC model is to call PFC program through a built-in interface of FLAC, so as to perform a coupling analysis between the finite difference model of the base plate and the roadbed bottom layer established by FLAC and the discrete element model of the roadbed surface layer established by PFC. - In an embodiment, the generating a lognormal random field of particle contact friction coefficients in step 3, includes:
- extracting spatial positions of the particle contact nodes and the particle contact friction coefficients in the discrete element model of the roadbed surface layer, and importing the spatial positions and the particle contact friction coefficients into a MATLAB script; and
- setting the spatial correlation distance θ of the particle contacts of the roadbed surface layer in the MATLAB script, calculating a correlation coefficient ρ between friction coefficients at any two coordinate points by using a Markov spatial correlation function, and generating the lognormal random field of the particle contact friction coefficients.
- In an embodiment, a formula of the Markov spatial correlation function in the MATLAB script is as follows:
-
- where ρ represents the correlation coefficient between the two coordinate points of the roadbed surface layer, τx and τy respectively correspond to absolute distances between the two coordinate points in horizontal and vertical directions, and θ represents a spatial correlation length.
- In an embodiment, obtaining the worst spatial correlation distance θworst in step 5, includes:
- increasing the spatial correlation distance θ in the step by step manner, generating N random fields of the particle contact friction coefficients of the roadbed surface layer (i.e., repeatedly performing the random field generation on the particle contact friction coefficients of the roadbed surface layer for N times) and performing N sampling on the precipitation under a i-th spatial correlation distance θi, applying a high-speed train load and bringing the high-speed train load into model calculation to obtain N groups of dynamic response results of the base plate, and taking an average value of the N groups of dynamic response results as a representative value of the dynamic response of the base plate under θi; and
- comparing the representative values of the dynamic response of the base plate under different θi, taking a spatial correlation distance corresponding to the representative value of the maximum dynamic response as the worst spatial correlation distance θworst.
- In an embodiment, obtaining the damage probability pf of the roadbed surface layer in step 6, includes:
- generating M random fields of the particle contact friction coefficients of the roadbed surface layer and performing M sampling on the precipitation under the worst spatial correlation distance θworst to obtain M groups of samples, and applying the high-speed train load to calculate a dynamic response of the base plate in each of the M groups of samples; and
- when any dynamic response of the base plate exceeds a specified safety threshold, determining that the roadbed surface layer is damaged and endangers a superstructure, and counting the number of damaged samples i, and calculating the damage probability using a formula: pf=i/M.
- A forewarning system by using the above ballastless track roadbed damage forewarning method considering uncertainty, includes: an information collection module, a processing module, and a forewarning module;
- the information collection module is configured to collect material parameters of the ballastless track roadbed and real-time precipitation information;
- the processing module is configured to load the FLAC-PFC model and calibrate the mesoscopic parameters based on the information collected by the information collection module; load the MATLAB script, generate a random field, and assign the random field to the particle contact nodes; calculate the worst spatial correlation distance θworst and calculate the damage probability pf; and
- the forewarning module is configured to determine whether the damage probability exceeds the alarm threshold; if damage probability exceeds the alarm threshold, output the alarm information; if the damage probability does not exceed the alarm threshold, exit.
- In an embodiment, each of the information collection module, the processing module and the forewarning module is embodied by at least one processor and at least one memory coupled to the at least one processor, and the at least one memory stores programs executable by the at least one processor.
- Compared with the prior art, the beneficial effects of the disclosure are as follows.
- According to the scheme provided by the embodiment of the disclosure, the randomness of the precipitation in the roadbed area and the uncertainty of the roadbed surface layer material parameters changing with space are considered, and the damage probability of the roadbed surface layer under the action of rain and dynamic load is quantitatively predicted by combining a FLAC-PFC coupling numerical analysis method. The disclosure provides a powerful means for damage monitoring of the ballastless track roadbed, and the analysis method has clear flow and strong reliability.
-
FIG. 1 illustrates a flowchart of a forewarning method according to an embodiment of the disclosure. -
FIG. 2 illustrates a schematic diagram of a FLAC-PFC coupling model of a roadbed of a high-speed rail ballastless track according to an embodiment of the disclosure; in which, 1: base plate, 2: roadbed surface layer, 3: roadbed bottom layer. - The content not described in detail in the specification of the disclosure belongs to the publicly known technique of those skilled in the art.
- The following will provide a clear and complete description of the technical solution in embodiments of the disclosure in conjunction with the accompanying drawings. Apparently, the described embodiments are some of the embodiments of the disclosure, not all of them. Based on the embodiments in the disclosure, all other embodiments obtained by those skilled in the art without creative work fall within the scope of protection of the disclosure.
- Referring to
FIG. 1 , a ballastless track roadbed damage forewarning method considering uncertainty provided by an embodiment of the disclosure is illustrated, including the following steps. - Step 101: obtaining precipitation in a monitoring area for a period of time, and calculating a probability density function of the precipitation using a maximum entropy principle; obtaining material parameters of a base plate 1, a
roadbed surface layer 2, and a roadbed bottom layer 3 of the ballastless track. - Step 102: establishing a finite difference model of the base plate 1 and the roadbed bottom layer 3 through fast Lagrangian analysis of continuum (FLAC), assigning values based on the material parameters obtained in
step 101, establishing a discrete element model of theroadbed surface layer 2 through particle flow code (PFC), and calibrating mesoscopic parameters of particle materials of theroadbed surface layer 2 based on the material parameters obtained instep 101. - Step 103: extracting spatial positions of particle contact nodes and friction coefficients of particle contacts of the discrete element model of the
roadbed surface layer 2, and importing them into a MATLAB script; setting a spatial correlation length θ of the friction coefficients of the particle contacts of theroadbed surface layer 2 in the MATLAB script, calculating a correlation coefficient ρ between the friction coefficients at any two coordinate points by using a Markov spatial correlation function, generating a lognormal random field of the friction coefficients of the particle contacts and giving it to the corresponding particle contact nodes in the FLAC-PFC model, and the Markov spatial correlation function is as follows: -
- where ρ represents the correlation coefficient between two coordinate points on the roadbed surface layer, τx and τy respectively correspond to absolute distances between two points in horizontal and vertical directions, and θ represents the spatial correlation length.
- Step 104: performing Monte Carlo sampling on the precipitation according to its probability density function, and adjusting a corresponding fluid domain of the roadbed surface layer according to the samples (i.e., precipitation).
- Step 105: applying a high-speed train load, and calculating a dynamic response of the base plate; repeating
step 103, and increasing θ in a step by step manner. Under each θi, the random field generation of the friction coefficients of the particle contacts of the roadbed surface layer is repeated for N times, and the precipitation is sampled for N times. The high-speed train load is applied and brought into the model calculation, and an average value of N groups of dynamic response of the base plate is taken as a representative value of the dynamic response of the base plate under θi. By comparing the representative values of the dynamic response of the base plate under different θi, a worst spatial correlation distance θworst is determined. Where i is a positive integer. - Step 106: under θworst, the random field generation of the friction coefficients of the particle contacts of the roadbed surface layer is repeated for M times, and the precipitation is sampled for M times to obtain M groups of samples, and the dynamic response results of the base plate in each group of samples are calculated by applying the high-speed train load; a safety threshold of the dynamic response according to the standard and a damage alarm threshold p of the roadbed surface layer according to the operation conditions are input. The dynamic response of the base plate includes the maximum vibration amplitude, maximum vibration acceleration, maximum dynamic stress and maximum dynamic strain of the base plate. When any dynamic response of the base plate exceeds the specified macro-safety threshold, it is determined that roadbed surface layer is damaged and endangers the superstructure. The damage probability pf of the roadbed surface layer is calculated, when pf exceeds the alarm threshold p, the alarm information is sent to the monitoring personnel, otherwise quitting. The calculation method of the damage probability pf includes: counting the number of damaged samples i in the M groups of samples, pf=i/M.
- Referring to
FIG. 2 , the FLAC-PFC coupling numerical model is established, which includes the base plate 1 of the ballastless track based on an elastic material model, with the parameters including density and modulus of elasticity; the roadbed bottom layer 3 based on a Mohr-Coulomb model, with parameters including density, modulus of elasticity, cohesion, and friction angle; and theroadbed surface layer 2 based on a Hertz contact model, with the parameters including gradation, shear modulus, Poisson's ratio, and friction coefficient. - The scheme provided by the embodiment of the disclosure takes into account the randomness of precipitation in the subgrade area and the spatial uncertainty of the material parameters of the surface layer of the roadbed, and combines the FLAC-PFC coupling numerical analysis method to predict the damage of the track roadbed under dynamic loads and rain, which helps to ensure the safety of high-speed railway traffic.
- It should be noted that in this article, the terms “include”, “contain”, or any other variation thereof are intended to cover non-exclusive inclusion, so that a process, method, article, or device that includes a series of elements includes not only those elements, but also other elements that are not explicitly listed, or also elements inherent in such the process, method, article, or device. Without further restrictions, the element defined by the statement “includes a . . . ” does not preclude the existence of another identical element in the process, method, article, or device that includes the element.
- The embodiments of the disclosure have been described above with reference to the accompanying drawings, but the disclosure is not limited to the specific embodiments described above, which are merely illustrative and not restrictive. With the inspiration of the disclosure, those skilled in the art can also make many forms that fall within the protection of the disclosure, without departing from the scope protected by the purpose and claims of the disclosure.
Claims (7)
1. A ballastless track roadbed damage forewarning method, comprising the following steps:
step 1: obtaining roadbed material parameters and precipitation over a period of time;
step 2: establishing a fast Lagrangian analysis of continuum (FLAC)-particle flow code (PFC) model of a ballastless track roadbed, calibrating particle materials of a roadbed surface layer according to the roadbed material parameters obtained in step 1, and obtaining mesoscopic parameters of the particle materials; wherein FLAC is a kind of continuous medium mechanics analysis software, which is used to establish a finite difference model of a base plate and a roadbed bottom layer; PFC is a particle flow program built in FLAC, which is used to establish a discrete element model of the roadbed surface layer; and the FLAC-PFC model is to call PFC program through a built-in interface of FLAC, so as to perform a coupling analysis between the finite difference model of the base plate and the roadbed bottom layer established by FLAC and the discrete element model of the roadbed surface layer established by PFC;
step 3: generating a lognormal random field of particle contact friction coefficients, and assigning the lognormal random field to particle contact nodes of the roadbed surface layer in the FLAC-PFC model obtained in step 2;
step 4: performing Monte Carlo sampling based on a probability density function of the precipitation over the period of time obtained in step 1, and adjusting a fluid domain of the roadbed surface layer according to the precipitation;
step 5: increasing a spatial correlation distance θ in a step by step manner, calculating a dynamic response of the base plate, and determining a worst spatial correlation distance θworst; and
step 6: calculating a damage probability pf of the roadbed surface layer under the worst spatial correlation distance θworst obtained in step 5; outputting alarm information to monitoring personnel when the damage probability pf exceeds an alarm threshold p, thereby maintaining, by the monitoring personnel, the ballastless track roadbed in response to the alarm information; otherwise, quitting.
2. The ballastless track roadbed damage forewarning method according to claim 1 , wherein the roadbed material parameters in step 1 comprise: base plate parameters, and material parameters of the roadbed surface layer and the roadbed bottom layer; and
wherein the base plate parameters comprise: a density and modulus of elasticity; and the material parameters of the roadbed surface layer and the roadbed bottom layer comprise: a density, a gradation, modulus of elasticity, a cohesion and a friction angle.
3. The ballastless track roadbed damage forewarning method according to claim 1 , wherein generating the lognormal random field of the particle contact friction coefficients in step 3, comprises:
extracting spatial positions of the particle contact nodes and the particle contact friction coefficients in the discrete element model of the roadbed surface layer, and importing the spatial positions and the particle contact friction coefficients into a MATLAB script; and
setting the spatial correlation distance θ of the particle contact friction coefficients of the roadbed surface layer in the MATLAB script, calculating a correlation coefficient ρ between friction coefficients at any two coordinate points by using a Markov spatial correlation function, and generating the lognormal random field of the particle contact friction coefficients.
4. The ballastless track roadbed damage forewarning method according to claim 3 , wherein a formula of the Markov spatial correlation function in the MATLAB script is as follows:
wherein ρ represents the correlation coefficient between two coordinate points of the roadbed surface layer, τx and τy respectively correspond to absolute distances between the two coordinate points in horizontal and vertical directions, and θ represents a spatial correlation length.
5. The ballastless track roadbed damage forewarning method according to claim 1 , wherein determining the worst spatial correlation distance θworst in step 5, comprises:
increasing the spatial correlation distance θ in the step by step manner, generating N random fields of the particle contact friction coefficients of the roadbed surface layer and performing N sampling on the precipitation under a i-th spatial correlation distance θi, applying a high-speed train load and bringing the high-speed train load into model calculation to obtain N groups of dynamic response results of the base plate, and taking an average value of the N groups of dynamic response results as a representative value of the dynamic response of the base plate under θi; and
comparing representative values of the dynamic response of the base plate under different θi, taking the spatial correlation distance corresponding to the representative value of the maximum dynamic response as the worst spatial correlation distance θworst.
6. The ballastless track roadbed damage forewarning method according to claim 1 , wherein calculating the damage probability pf of the roadbed surface layer in step 6 comprises:
generating M random fields of the particle contact friction coefficients of the roadbed surface layer and performing M sampling on the precipitation under the worst spatial correlation distance θworst to obtain M groups of samples, and applying a high-speed train load to calculate a dynamic response of the base plate in each of the M groups of samples; and
when any dynamic response of the base plate exceeds a specified safety threshold, determining that the roadbed surface layer is damaged and endangers a superstructure, counting a number of damaged samples i, and calculating the damage probability using a formula: pf=i/M.
7. A forewarning system by using the ballastless track roadbed damage forewarning method according to claim 1 , comprising: an information collection module, a processing module, and a forewarning module; wherein,
the information collection module is configured to collect material parameters of the ballastless track roadbed and real-time precipitation information;
the processing module is configured to load the FLAC-PFC model and calibrate the mesoscopic parameters based on information collected by the information collection module; load the MATLAB script, generate a random field, and assign the random field to the particle contact nodes; calculate the worst spatial correlation distance θworst and calculate the damage probability pf, and
the forewarning module is configured to determine whether the damage probability pf exceeds the alarm threshold; if the damage probability pf exceeds the alarm threshold, output the alarm information to monitoring personnel, and the monitoring personnel maintain, in response to the alarm information, the ballastless track roadbed; if the damage probability pf does not exceed the alarm threshold, exit.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210366148.1 | 2022-04-08 | ||
CN2022103661481 | 2022-04-08 | ||
CN202210366148.1A CN114491781B (en) | 2022-04-08 | 2022-04-08 | Ballastless track foundation bed damage early warning method and system considering uncertainty |
Publications (2)
Publication Number | Publication Date |
---|---|
US11782185B1 US11782185B1 (en) | 2023-10-10 |
US20230324581A1 true US20230324581A1 (en) | 2023-10-12 |
Family
ID=81487618
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US18/297,455 Active US11782185B1 (en) | 2022-04-08 | 2023-04-07 | Ballastless track roadbed damage forewarning method and system considering uncertainty |
Country Status (2)
Country | Link |
---|---|
US (1) | US11782185B1 (en) |
CN (1) | CN114491781B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115266929B (en) * | 2022-08-09 | 2023-05-23 | 成都工业职业技术学院(成都市工业职业技术学校) | High-speed railway base resonance early warning method and system based on deep learning |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105507085B (en) * | 2015-12-01 | 2017-05-10 | 北京交通大学 | Ballast track bed structure |
CN209447406U (en) * | 2018-05-21 | 2019-09-27 | 西南交通大学 | The simulator of non-fragment orbit surface layer of subgrade bed frost boiling disease |
CN108447377A (en) * | 2018-05-21 | 2018-08-24 | 西南交通大学 | The simulator and analogy method of non-fragment orbit surface layer of subgrade bed frost boiling disease |
CN109610253B (en) * | 2018-11-19 | 2020-04-17 | 西南交通大学 | Method for evaluating slurry leakage hazard grade of ballastless track subgrade |
CN109629347B (en) * | 2018-11-19 | 2020-04-17 | 西南交通大学 | Method for evaluating slurry leakage hazard grade of ballastless track subgrade |
CN110188442B (en) * | 2019-05-23 | 2023-04-21 | 南京金蓝智慧城市规划设计有限公司 | Finite element simulation analysis method for coupling power of roadbed foundation of high-speed railway ballastless track |
CN110390176B (en) * | 2019-07-31 | 2020-05-19 | 西南交通大学 | Ballastless track freezing and damage behavior calculation method |
CN110728000B (en) * | 2019-09-29 | 2021-06-29 | 华中科技大学 | Ballastless track elastic fastener damage identification method |
CN110781628B (en) * | 2019-11-08 | 2020-11-06 | 西南交通大学 | Method for identifying damage to connecting surface of double-block ballastless track bed plate and sleeper block |
CN112214919B (en) * | 2020-10-19 | 2021-08-17 | 西南交通大学 | Method for calculating high-cycle fatigue damage behavior of CA mortar of plate-type ballastless track |
CN112364553B (en) * | 2021-01-13 | 2021-06-25 | 西南交通大学 | Method for evaluating coupling of finite element-discrete element of surface layer seepage corrosion of ballastless track foundation bed of high-speed rail |
CN113657041A (en) * | 2021-08-25 | 2021-11-16 | 哈尔滨工业大学 | Intelligent sensing and forecasting system for physical and mechanical states of roadbed in alpine region |
CN114228786B (en) * | 2021-12-29 | 2024-03-26 | 中铁第四勘察设计院集团有限公司 | Plate-type ballastless track and plate seam dislocation monitoring system and health monitoring method thereof |
-
2022
- 2022-04-08 CN CN202210366148.1A patent/CN114491781B/en active Active
-
2023
- 2023-04-07 US US18/297,455 patent/US11782185B1/en active Active
Also Published As
Publication number | Publication date |
---|---|
US11782185B1 (en) | 2023-10-10 |
CN114491781A (en) | 2022-05-13 |
CN114491781B (en) | 2022-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111859724B (en) | Fatigue life prediction method and system for mixed-drive ballastless track | |
Wang et al. | A novel hybrid-Garch model based on ARIMA and SVM for PM2. 5 concentrations forecasting | |
US11782185B1 (en) | Ballastless track roadbed damage forewarning method and system considering uncertainty | |
Lam et al. | Identification of rail-sleeper-ballast system through time-domain Markov chain Monte Carlo-based Bayesian approach | |
Bashor et al. | Full-scale performance evaluation of tall buildings under wind | |
CN111368476A (en) | Method for evaluating applicability of high-speed rail bridge-track system under action of near-fault earthquake | |
Zhan et al. | Wireless rail fastener looseness detection based on MEMS accelerometer and vibration entropy | |
Xu et al. | Reliability assessment on earthquake early warning: A case study from Taiwan | |
CN103913512B (en) | The damage reason location system of suspension cable periodic detection | |
CN112945490A (en) | Method for testing bearing capacity of bridge based on deflection influence line | |
Muin et al. | Human-machine collaboration framework for structural health monitoring and resiliency | |
Wei et al. | Dynamic performance of a deep buried pile-plank structure transition section for a high-speed railway—Field tests and numerical analyses | |
Hu et al. | A feasibility study on void detection of cement-emulsified asphalt mortar for slab track system utilizing measured vibration data | |
CN104698936A (en) | Big data concept-based bridge structure health monitoring system | |
Jiang et al. | Seismic-induced track spectrum characteristics of high-speed railway bridges | |
Zhang et al. | Response analysis and effect evaluation of dynamic stabilization for ballasted track | |
CN106926755A (en) | Dynamic locator gradient defect inspection method and system | |
Aktan et al. | Opportunities and challenges in health monitoring of constructed systems by modal analysis | |
Lin et al. | Development of on-site earthquake early warning system for Taiwan | |
CN103911958B (en) | The damage reason location system of suspension bridge and arch bridge suspender periodic detection and method thereof | |
Zhu et al. | Investigation on the pattern for train-induced strains of the long-span steel truss railway bridge | |
Chase | Smarter bridges, why and how? | |
CN108385551A (en) | A kind of metal sound barrier with monitoring function | |
Chase | High-tech inspection | |
Zhangabylova et al. | Experimental determination of dynamic characteristics of a railway track |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: SOUTHWEST JIAOTONG UNIVERSITY, CHINA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:LIU, KAIWEN;QIU, RUIZHE;SU, QIAN;AND OTHERS;REEL/FRAME:063262/0790 Effective date: 20230407 |
|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO SMALL (ORIGINAL EVENT CODE: SMAL); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |