CN116295804B - Method for constructing vibration map - Google Patents
Method for constructing vibration map Download PDFInfo
- Publication number
- CN116295804B CN116295804B CN202310122131.6A CN202310122131A CN116295804B CN 116295804 B CN116295804 B CN 116295804B CN 202310122131 A CN202310122131 A CN 202310122131A CN 116295804 B CN116295804 B CN 116295804B
- Authority
- CN
- China
- Prior art keywords
- vibration
- map
- section
- vlzmax
- load
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000004458 analytical method Methods 0.000 claims abstract description 33
- 238000011156 evaluation Methods 0.000 claims abstract description 30
- 238000010586 diagram Methods 0.000 claims abstract description 6
- 239000002689 soil Substances 0.000 claims description 24
- 238000012360 testing method Methods 0.000 claims description 15
- 230000004044 response Effects 0.000 claims description 13
- 238000005516 engineering process Methods 0.000 claims description 10
- 238000009826 distribution Methods 0.000 claims description 9
- 239000003086 colorant Substances 0.000 claims description 8
- 230000004927 fusion Effects 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 7
- 241001669679 Eleotris Species 0.000 claims description 6
- 238000004364 calculation method 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
- 230000005284 excitation Effects 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000013075 data extraction Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000010183 spectrum analysis Methods 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 4
- 239000000725 suspension Substances 0.000 description 4
- 230000007613 environmental effect Effects 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000009931 harmful effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method for constructing a vibration map, which divides urban rail transit into a plurality of sections and then collects vibration data of each section; then, carrying out 1/3 frequency multiplication analysis on the vibration data to obtain a vibration influence evaluation index VLzmax, and comparing the vibration influence evaluation index VLzmax with a limit value index in a standard to obtain a vibration superscript point; and then drawing the vibration influence evaluation index VLzmax in a two-dimensional coordinate formed by the longitudinal mileage of the line and the transverse distance of the vertical line in the form of cloud pictures, thermodynamic diagrams and the like to form a vibration map. According to the invention, the vibration data along the urban rail transit line are collected and drawn in the two-dimensional coordinates consisting of the longitudinal mileage of the line and the transverse distance of the vertical line to form the vibration map, so that the urban vibration condition can be conveniently and intuitively observed.
Description
Technical Field
The invention relates to the technical field of urban area environment vibration evaluation, in particular to a method for constructing a vibration map.
Background
Along with the development of urban modern traffic, the vibration caused by the development of urban modern traffic is increasingly frequent, and the vibration is transmitted to two sides of a road through a ground terrace, a building and the like, so that the vibration has a great influence on precise instruments, equipment, buildings and residents sensitive to the adjacent vibration. When the vibration frequency is matched with the natural frequency of certain organs of the human body, the vibration frequency has harmful effects on the human body. According to the statistics of the related countries, traffic vibration is most intense in public reflection except factories, enterprises and constructional engineering, so that urban vibration conditions are monitored in the process of urban management so as to carry out vibration analysis and later management operation of urban areas.
The vibration and the noise are accompanied, the larger the vibration is, the larger the noise is correspondingly, the noise generated by the vibration can directly influence people living around, and the track traffic is continuously increased, so that the urban residents can travel conveniently, and meanwhile, the distance between the residents and the vibration is shortened. Therefore, the problem of vibration and noise pollution generated by the vibration is an unavoidable problem for improving urban living quality, and the urban vibration condition and the influence of the vibration on surrounding residents can be checked by monitoring the vibration condition of urban rail transit.
The prior art mainly monitors the urban monitoring through a noise map, and has no vibration map which is convenient for urban vibration monitoring and convenient to use, so that how to conveniently acquire environmental vibration data and intuitively display the environmental vibration data, and convenience for vibration condition inquiry and urban management become an important subject to be solved currently and urgently.
Disclosure of Invention
Aiming at the defects that the rail transit vibration is difficult to monitor and the urban vibration condition cannot be visually observed, the invention provides a method for constructing the vibration map, which is convenient for urban vibration monitoring, and has convenient use and good effect.
In order to achieve the above purpose, the technical scheme of the invention is as follows:
a method for constructing a vibration map comprises the following steps:
s1: dividing urban rail transit into a plurality of sections along a line area in the longitudinal direction of the line;
s2: obtaining ground vibration data of each section through simulation calculation or field test;
the operation steps of obtaining ground vibration data through simulation calculation comprise:
s2.11: establishing a dynamic model of a vehicle-track-tunnel subsystem based on a vehicle-track coupling dynamic theory, and carrying out frequency domain solving by adopting a random dynamic theory to obtain fulcrum counter force acting on a fastener, namely a vibration source excitation load;
s2.12: carrying out longitudinal linear equivalence on the excitation load of the vibration source by adopting a formula 2-1 to obtain an equivalent longitudinal uniform load;
L F (f i )=[P 1 (f i )+P 2 (f i )×(1-d)/d]×N b ×N w /L car 2-1
wherein: l (L) F For equivalent longitudinal uniform load, the unit is kN/m; p (P) 1 、P 2 Is the fulcrum reaction load of two fasteners adjacent to the lower part of the wheel set, and the unit is kN, N b Is the number of bogies of each section of vehicle, N b =2，N w Is the number of bogies of each section of vehicle, N w =2；L car Is 1 section of vehicle length load, the unit is m, and d is sleeper spacing，f i The unit is Hz for frequency load;
s2.13: defining the transverse influence range of each section;
s2.14: extracting the number of layers of each section soil body layering, the soil quality type and thickness information of each layer of soil body, and obtaining physical parameters of each layer of soil body according to engineering geological survey data, wherein the physical parameters comprise density, elastic modulus, dynamic elastic modulus, poisson ratio and damping parameters;
s2.15: based on the harmonic response analysis technology, a two-dimensional frequency domain analysis model of a single section track bed-tunnel-soil body coupling system is established;
s2.16: selecting at least 3 sections, establishing a two-dimensional frequency domain analysis model of each section ballast bed-tunnel-earth subsystem based on a harmonic response analysis technology, applying equivalent longitudinal uniform load obtained by a formula 2-1 to the model, and solving and extracting ground vibration A 1 (f) Wherein f represents frequency; performing field test on the selected section, wherein data extraction points in a two-dimensional frequency domain analysis model of measuring point arrangement correspond to each other one by one to obtain actual ground vibration A 2 (t), wherein t represents time; for actual measurement of ground vibration A 2 (t) to be equal to L F Spectrum analysis is carried out at the same frequency interval to obtain A 2 (f) Comparing the frequency spectrums of the ground vibration responses obtained by solving the field test and the two-dimensional frequency domain analysis model, correcting the equivalent longitudinal uniform load by adopting the field test result, and correcting the coefficient alpha (f) on each frequency i ) Calculating by adopting a formula 2-2;
α(f i )=A 2 (f i )/A 1 (f i ) 2-2
s2.17: adopting a formula 2-2 to correct the equivalent longitudinal uniform load;
L' F (f i )=α(f i )[P 1 (f i )+P 2 (f i )×(1-d)/d]×N b ×N w /L car 2-3
wherein: l'. F For the corrected equivalent longitudinal uniform load, the unit is kN/m; p (P) 1 、P 2 Is adjacent to the underside of the wheel pairThe fulcrum counterforces of the two fasteners are in the units of kN and N b Is the number of bogies of each section of vehicle, N b =2，N w Is the number of bogies of each section of vehicle, N w =2；L car 1 section of vehicle length, the unit is m, d is sleeper spacing, f i Frequency in Hz;
s2.18: based on the harmonic response analysis technology, the establishment and the solution of the two-dimensional frequency domain analysis model of each section ballast bed-tunnel-earth subsystem are sequentially completed in sequence, so that the prediction of all section ground vibration data can be rapidly completed; the load applied in the model is the equivalent longitudinal uniform load obtained by the correction of the formula 2-3;
s3: performing 1/3 frequency multiplication analysis on the obtained ground vibration data to obtain a vibration influence evaluation index VLzmax of each section;
s4: comparing the calculated vibration influence evaluation index VLzmax with a limit value in the current standard, judging whether a vibration exceeding standard point exists or not, and if so, recording the exceeding standard point;
s5: taking the interval distribution direction of the sections as an ordinate and the direction of the ground vertical ordinate as an abscissa, interpolating and fitting all obtained vibration influence evaluation indexes VLzmax, representing the vibration influence evaluation index VLzmax values by different colors and color depths, drawing a distribution diagram of the vibration influence evaluation index VLzmax, and marking out overproof points to obtain the vibration map.
Further, in the step S2.18, the order is from small to large or from large to small.
Further, the operation step of obtaining ground vibration data through field test comprises the following steps:
s2.21: defining an influence range along the rail transit;
s2.22: dividing an influence range into grids with smaller lengths and widths, arranging measuring points on corner points of the grids, and installing a vibration acceleration sensor;
s2.23: and (5) installing vibration signal acquisition and storage equipment and acquiring vibration data.
Further, the vibration map obtained in the step S5 is imported with google and hundred degree map information, and road and building information in the area corresponding to the vibration map are displayed in a fusion mode.
Further, the method further comprises the step S6: drawing a vibration standard exceeding map, wherein the drawing process is as follows:
s61: comparing the calculated vibration influence evaluation index VLzmax with a limit value in the current standard;
s62: and displaying the vibration influence evaluation index VLzmax which is lower than the limit value and exceeds the limit value by adopting two different colors, namely the vibration superscalar map.
Further, map information such as google and hundred degrees is imported into the vibration standard exceeding map obtained in the step S6, and road and building information in the area corresponding to the vibration standard exceeding map are displayed in a fusion mode.
Further, google and hundred-degree map information is imported into the vibration standard exceeding map obtained in the step S6, and road and building information in the area corresponding to the vibration standard exceeding map are displayed in a fusion mode.
Compared with the prior art, the invention has the following advantages and beneficial effects:
1. dividing an urban rail transit line area into a plurality of sections along a line longitudinal direction, acquiring ground vibration data of each section, analyzing the vibration data to obtain vibration influence evaluation indexes, and finally drawing a vibration distribution cloud picture (namely a vibration map); the vibration map is simple in overall manufacturing process, convenient to use, convenient to carry out visual display of vibration data and the phenomenon that the vibration data exceeds standard, and then convenient to monitor influence of rail transit on urban vibration and convenient to carry out urban management.
2. The invention sequentially completes the establishment and the solution of the two-dimensional frequency domain analysis model of each section ballast bed-tunnel-earth subsystem based on the harmonic response analysis technology, and can rapidly complete the prediction of the ground vibration data of all sections; the vibration data are obtained by using the method, the operation is simple, the accuracy of the vibration data obtained by prediction is high, the vibration data acquisition of each section is not required by staff, and the method is efficient and trouble-saving.
3. The invention can adopt the vibration data obtained by simulation to carry out vibration map drawing, and adopts the vibration data actually measured on site to carry out vibration map drawing.
4. According to the invention, the vibration map can be imported with map information such as google, hundred degrees and the like, so that the vibration data and the information such as regional roads, buildings and the like can be fused and displayed, and the vibration distribution condition can be checked conveniently; and the vibration exceeding map obtained by representing the vibration influence evaluation indexes VLzmax exceeding the limit value and not exceeding the limit value by adopting different colors is more convenient for observing the exceeding area and the exceeding condition.
Drawings
Fig. 1 is a construction flow chart of a construction method of a vibration map of the present invention.
Fig. 2 is a schematic diagram of soil distribution and section division of urban rail transit according to an embodiment of a method for constructing a vibration map of the present invention.
Fig. 3 shows each of the soil type numbers and the corresponding soil type in fig. 2.
Fig. 4 is a vehicle-track-tunnel subsystem dynamics model of a vibration map construction method of the present invention.
FIG. 5 is a two-dimensional frequency domain analysis model of a ballast-tunnel-earth subsystem of a vibration map construction method of the present invention.
Fig. 6 is a vibration map of rail transit in the form of a cloud image along the line drawn by a vibration map construction method of the present invention.
Fig. 7 is a vibration superscalar map of rail transit in the form of a cloud image along the line drawn by the method for constructing a vibration map of the present invention.
Description of the embodiments
The invention is further described below with reference to the accompanying drawings. In the figure 2 of the accompanying drawings, gray vertical lines represent calculated section division, gray broken lines represent soil layering interfaces, numerals represent soil type numbers corresponding to soil layers, distances between two adjacent intersection points of the vertical lines and the broken lines represent soil layer thicknesses, and in the figure, soil types corresponding to the soil type numbers are shown in figure 2. Mc in the diagram of fig. 4 is the body mass; jc is the inertia of the vehicle body point head; βc is the movement of the vehicle body nodding; zc is the vertical movement of the vehicle body; mt is bogie mass; jt is bogie nodding inertia; βt1 is the front bogie nodding motion; βt2 is the rear bogie nodding motion; zt1 is the vertical movement of the front bogie; zt2 is the vertical movement of the rear bogie; mw is the wheel set mass; zw1 is the vertical movement of wheel set 1; zw2 is the vertical movement of wheel set 2; zw3 is the vertical movement of wheel set 3; zw4 is the vertical movement of wheel set 4; ks1 is the primary suspension stiffness; cs1 is a series of suspension damping; ks2 is the secondary suspension stiffness; cs2 is a secondary suspension damper; kh is the wheel track Hertz contact equivalent spring rate; p1 is the wheeltrack force at wheelset 1; p2 is the wheeltrack force at wheelset 2; p3 is the wheeltrack force at wheelset 3; p4 is the wheeltrack force at wheelset 4; mr is the mass density of the steel rail; EI is a rail bending parameter; zr is the vertical movement of the steel rail; kp is the fastener spring rate; cp is fastener spring damping; v is the driving speed.
Example 1: a method for constructing a vibration map comprises the following steps:
s1: dividing urban rail transit into a plurality of sections along a line area at intervals of 5-10 m along the longitudinal direction of the line;
s2: obtaining ground vibration data of each section through simulation calculation;
the operation steps of obtaining ground vibration data through simulation calculation comprise:
s2.11: establishing a dynamic model of a vehicle-track-tunnel subsystem based on a vehicle-track coupling dynamic theory, and carrying out frequency domain solving by adopting a random dynamic theory to obtain fulcrum counter force acting on a fastener, namely a vibration source excitation load;
s2.12: carrying out longitudinal linear equivalence on the excitation load of the vibration source by adopting a formula 2-1 to obtain an equivalent longitudinal uniform load;
L F (f i )=[P 1 (f i )+P 2 (f i )×(1-d)/d]×N b ×N w /L car 2-1
wherein: l (L) F For equivalent longitudinal uniform load, the unit is kN/m; p (P) 1 、P 2 Is the fulcrum reaction load of two fasteners adjacent to the lower part of the wheel set, and the unit is kN, N b Is the number of bogies of each section of vehicle, N b =2，N w Is the number of bogies per vehicle section,N w =2；L car is 1 section of vehicle length load, the unit is m, d is sleeper spacing, f i The unit is Hz for frequency load;
s2.13: defining the transverse influence range of each section;
s2.14: extracting the number of layers of each section soil body layering, the soil quality type and thickness information of each layer of soil body, and obtaining physical parameters of each layer of soil body according to engineering geological survey data, wherein the physical parameters comprise density, elastic modulus, dynamic elastic modulus, poisson ratio and damping parameters;
s2.15: based on the harmonic response analysis technology, a two-dimensional frequency domain analysis model of a single section track bed-tunnel-soil body coupling system is established;
s2.16: selecting at least 3 sections, establishing a two-dimensional frequency domain analysis model of each section ballast bed-tunnel-earth subsystem based on a harmonic response analysis technology, applying equivalent longitudinal uniform load obtained by a formula 2-1 to the model, and solving and extracting ground vibration A 1 (f) Wherein f represents frequency; performing field test on the selected section, wherein data extraction points in a two-dimensional frequency domain analysis model of measuring point arrangement correspond to each other one by one to obtain actual ground vibration A 2 (t), wherein t represents time; for actual measurement of ground vibration A 2 (t) to be equal to L F Spectrum analysis is carried out at the same frequency interval to obtain A 2 (f) Comparing the frequency spectrums of the ground vibration responses obtained by solving the field test and the two-dimensional frequency domain analysis model, correcting the equivalent longitudinal uniform load by adopting the field test result, and correcting the coefficient alpha (f) on each frequency i ) Calculating by adopting a formula 2-2;
α(f i )=A 2 (f i )/A 1 (f i ) 2-2
s2.17: adopting a formula 2-2 to correct the equivalent longitudinal uniform load;
L' F (f i )=α(f i )[P 1 (f i )+P 2 (f i )×(1-d)/d]×N b ×N w /L car 2-3
wherein: l'. F For the corrected equivalent longitudinal directionUniformly distributing load with the unit of kN/m; p (P) 1 、P 2 Is the fulcrum reaction load of two fasteners adjacent to the lower part of the wheel set, and the unit is kN, N b Is the number of bogies of each section of vehicle, N b =2，N w Is the number of bogies of each section of vehicle, N w =2；L car Is 1 section of vehicle length load, the unit is m, d is sleeper spacing, f i The unit is Hz for frequency load;
s2.18: based on the harmonic response analysis technology, the establishment and the solution of the two-dimensional frequency domain analysis model of each section ballast bed-tunnel-earth subsystem are sequentially completed according to the sequence, and the prediction of the ground vibration data of all sections can be rapidly completed; the load applied in the model is the equivalent longitudinally uniform load obtained by the correction of the formula 2-3.
S3: performing 1/3 frequency multiplication analysis on the obtained ground vibration data to obtain a vibration influence evaluation index VLzmax of each section;
s4: comparing the calculated vibration influence evaluation index VLzmax with a limit value in the current standard, judging whether a vibration exceeding standard point exists or not, and if so, recording the exceeding standard point;
s5: taking the interval distribution direction of the sections as an ordinate and the direction of the ground vertical ordinate as an abscissa, interpolating and fitting all obtained vibration influence evaluation indexes VLzmax, representing the vibration influence evaluation index VLzmax values by different colors and color depths, drawing a distribution diagram of the vibration influence evaluation index VLzmax, and marking out overproof points to obtain the vibration map.
Wherein, 1/3 frequency multiplication analysis executes relevant standards of the current urban area environmental vibration influence evaluation in China.
In the step S2.18, the order is from small to large or from large to small.
Example 2: the difference from embodiment 1 is that the vibration map obtained in step S5 is imported with map information such as google and hundred degrees, and road and building information in the region corresponding to the vibration map is displayed in a fusion manner.
S6: drawing a map with vibration exceeding standard, wherein the drawing steps are as follows
S61: the obtained vibration influence VLzmax is compared with a standard limit value of the vibration influence evaluation index VLzmax.
S62: and displaying the vibration influence evaluation index VLzmax which is lower than the limit value and exceeds the limit value by adopting two different colors, namely the vibration superscalar map.
And (3) importing map information such as google, hundred degrees and the like into the vibration standard exceeding map obtained in the step (S6), and carrying out fusion display on road and building information in the area corresponding to the vibration standard exceeding map.
The working principle of this embodiment is the same as that of embodiment 1.
Example 3: the difference from the embodiment 1 is that the step S2 is to acquire vibration data of each section by using a field test;
the operation steps of acquiring ground vibration data through field test comprise:
s2.21: defining an influence range along the rail transit;
s2.22: dividing an influence range into grids with smaller lengths and widths, arranging measuring points on corner points of the grids, and installing a vibration acceleration sensor;
s2.23: and (5) installing vibration signal acquisition and storage equipment and acquiring vibration data.
And (5) importing google and hundred-degree map information into the vibration map obtained in the step (S5), and fusing and displaying road and building information in the area corresponding to the vibration map.
S6: drawing a map with vibration exceeding standard, wherein the drawing steps are as follows
S61: the obtained vibration influence VLzmax is compared with a standard limit value of the vibration influence evaluation index VLzmax.
S62: and displaying the vibration influence evaluation index VLzmax which is lower than the limit value and exceeds the limit value by adopting two different colors, namely the vibration superscalar map.
And (3) importing google and hundred-degree map information into the vibration standard exceeding map obtained in the step (S6), and carrying out fusion display on road and building information in the area corresponding to the vibration standard exceeding map.
The working principle of this embodiment is the same as that of embodiment 1.
In the description of the present invention, numerous specific details are set forth. However, it is understood that embodiments of the invention may be practiced without these specific details. In some instances, well-known methods, structures and techniques have not been shown in detail in order not to obscure an understanding of this description.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some or all of the technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit of the invention, and are intended to be included within the scope of the appended claims and description.
Claims (6)
1. A method for constructing a vibration map is characterized in that: the method comprises the following steps:
s1: dividing urban rail transit into a plurality of sections along a line area in the longitudinal direction of the line;
s2: obtaining ground vibration data of each section through simulation calculation or field test;
the operation steps of obtaining ground vibration data through simulation calculation comprise:
s2.11: establishing a dynamic model of a vehicle-track-tunnel subsystem based on a vehicle-track coupling dynamic theory, and carrying out frequency domain solving by adopting a random dynamic theory to obtain fulcrum counter force acting on a fastener, namely a vibration source excitation load;
s2.12: carrying out longitudinal linear equivalence on the excitation load of the vibration source by adopting a formula 2-1 to obtain an equivalent longitudinal uniform load;
L F (f i )=[P 1 (f i )+P 2 (f i )×(1-d)/d]×N b ×N w /L car 2-1
wherein: l (L) F Is equal toThe load is uniformly distributed longitudinally, and the unit is kN/m; p (P) 1 、P 2 Is the fulcrum reaction load of two fasteners adjacent to the lower part of the wheel set, and the unit is kN, N b Is the number of bogies of each section of vehicle, N b =2，N w Is the number of bogies of each section of vehicle, N w =2；L car Is 1 section of vehicle length load, the unit is m, d is sleeper spacing, f i The unit is Hz for frequency load;
s2.13: defining the transverse influence range of each section;
s2.14: extracting the number of layers of each section soil body layering, the soil quality type and thickness information of each layer of soil body, and obtaining physical parameters of each layer of soil body according to engineering geological survey data, wherein the physical parameters comprise density, elastic modulus, dynamic elastic modulus, poisson ratio and damping parameters;
s2.15: based on the harmonic response analysis technology, a two-dimensional frequency domain analysis model of a single section track bed-tunnel-soil body coupling system is established;
s2.16: selecting at least 3 sections, establishing a two-dimensional frequency domain analysis model of each section ballast bed-tunnel-earth subsystem based on a harmonic response analysis technology, applying equivalent longitudinal uniform load obtained by a formula 2-1 to the model, and solving and extracting ground vibration A 1 (f) Wherein f represents frequency; performing field test on the selected section, wherein data extraction points in a two-dimensional frequency domain analysis model of measuring point arrangement correspond to each other one by one to obtain actual ground vibration A 2 (t), wherein t represents time; for actual measurement of ground vibration A 2 (t) to be equal to L F Spectrum analysis is carried out at the same frequency interval to obtain A 2 (f) Comparing the frequency spectrums of the ground vibration responses obtained by solving the field test and the two-dimensional frequency domain analysis model, correcting the equivalent longitudinal uniform load by adopting the field test result, and correcting the coefficient alpha (f) on each frequency i ) Calculating by adopting a formula 2-2;
α(f i )=A 2 (f i )/A 1 (f i ) 2-2
s2.17: adopting a formula 2-2 to correct the equivalent longitudinal uniform load;
L' F (f i )=α(f i )[P 1 (f i )+P 2 (f i )×(1-d)/d]×N b ×N w /L car 2-3
wherein: l'. F For the corrected equivalent longitudinal uniform load, the unit is kN/m; p (P) 1 、P 2 Is the fulcrum counterforce of two fasteners adjacent to the lower part of the wheel set, and the unit is kN, N b Is the number of bogies of each section of vehicle, N b =2，N w Is the number of bogies of each section of vehicle, N w =2；L car 1 section of vehicle length, the unit is m, d is sleeper spacing, f i Frequency in Hz;
s2.18: based on the harmonic response analysis technology, the establishment and the solution of the two-dimensional frequency domain analysis model of each section ballast bed-tunnel-earth subsystem are sequentially completed in sequence, so that the prediction of all section ground vibration data can be rapidly completed; the load applied in the model is the equivalent longitudinal uniform load obtained by the correction of the formula 2-3;
s3: performing 1/3 frequency multiplication analysis on the obtained ground vibration data to obtain a vibration influence evaluation index VLzmax of each section;
s4: comparing the calculated vibration influence evaluation index VLzmax with a limit value in the current standard, judging whether a vibration exceeding standard point exists or not, and if so, recording the exceeding standard point;
s5: taking the interval distribution direction of the sections as an ordinate and the direction of the ground vertical ordinate as an abscissa, interpolating and fitting all obtained vibration influence evaluation indexes VLzmax, representing the vibration influence evaluation index VLzmax values by different colors and color depths, drawing a distribution diagram of the vibration influence evaluation index VLzmax, and marking out overproof points to obtain the vibration map.
2. The method for constructing a vibration map according to claim 1, wherein: in the step S2.18, the order is from small to large or from large to small.
3. The method for constructing a vibration map according to claim 1, wherein: the operation steps of acquiring ground vibration data through field test comprise:
s2.21: defining an influence range along the rail transit;
s2.22: dividing an influence range into grids with smaller lengths and widths, arranging measuring points on corner points of the grids, and installing a vibration acceleration sensor;
s2.23: and (5) installing vibration signal acquisition and storage equipment and acquiring vibration data.
4. The method for constructing a vibration map according to claim 1, wherein: and (5) importing google and hundred-degree map information into the vibration map obtained in the step (S5), and fusing and displaying road and building information in the area corresponding to the vibration map.
5. The method for constructing a vibration map according to claim 1, wherein: further comprising step S6: drawing a vibration standard exceeding map, wherein the drawing process is as follows:
s61: comparing the calculated vibration influence evaluation index VLzmax with a limit value in the current standard;
s62: and displaying the vibration influence evaluation index VLzmax which is lower than the limit value and exceeds the limit value by adopting two different colors, namely the vibration superscalar map.
6. The method for constructing a vibration map according to claim 5, wherein: and (3) importing google and hundred-degree map information into the vibration standard exceeding map obtained in the step (S6), and carrying out fusion display on road and building information in the area corresponding to the vibration standard exceeding map.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310122131.6A CN116295804B (en) | 2023-02-16 | 2023-02-16 | Method for constructing vibration map |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310122131.6A CN116295804B (en) | 2023-02-16 | 2023-02-16 | Method for constructing vibration map |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116295804A CN116295804A (en) | 2023-06-23 |
| CN116295804B true CN116295804B (en) | 2023-11-03 |
Family
ID=86821428
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310122131.6A Active CN116295804B (en) | 2023-02-16 | 2023-02-16 | Method for constructing vibration map |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116295804B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117115374B (en) * | 2023-08-28 | 2024-04-16 | 中海环境科技(上海）股份有限公司 | Vibration map drawing method based on track traffic on-line monitoring |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1054754A (en) * | 1996-08-12 | 1998-02-24 | Nishimatsu Constr Co Ltd | Noise and vibration measuring method |
| JP2003294850A (en) * | 2002-04-04 | 2003-10-15 | Ohbayashi Corp | Ground response-analyzing method and system, program for making computer to execute the ground response- analyzing method, and record medium with recorded program |
| JP2004219237A (en) * | 2003-01-15 | 2004-08-05 | Taisei Corp | Method and system for predicting vibration |
| CN104598757A (en) * | 2015-02-12 | 2015-05-06 | 西南交通大学 | Method for predicting traffic noise of rail regions |
| CN104715138A (en) * | 2015-01-28 | 2015-06-17 | 西南交通大学 | Rail fastener counter-force computing method and subway vibration frequency domain response computing method |
| CN107589686A (en) * | 2016-07-09 | 2018-01-16 | 韩宇龙 | A kind of vision and vibration array perceive the intelligent safety defense monitoring system of integration |
| CN112767575A (en) * | 2021-01-19 | 2021-05-07 | 同济大学 | Vehicle-mounted adding instrument and method for rail transit operation period based on smart phone |
| CN114086967A (en) * | 2021-11-22 | 2022-02-25 | 中国地震局地震研究所 | Anti-excavation vibration monitoring method for subway tunnel |
| CN114519286A (en) * | 2022-02-17 | 2022-05-20 | 中铁二院工程集团有限责任公司 | Rapid calculation method for secondary structure noise caused by vibration of rail transit vehicle |
| CN115510711A (en) * | 2022-09-29 | 2022-12-23 | 东南大学 | Subway river-crossing tunnel vibration noise characteristic analysis method |
-
2023
- 2023-02-16 CN CN202310122131.6A patent/CN116295804B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1054754A (en) * | 1996-08-12 | 1998-02-24 | Nishimatsu Constr Co Ltd | Noise and vibration measuring method |
| JP2003294850A (en) * | 2002-04-04 | 2003-10-15 | Ohbayashi Corp | Ground response-analyzing method and system, program for making computer to execute the ground response- analyzing method, and record medium with recorded program |
| JP2004219237A (en) * | 2003-01-15 | 2004-08-05 | Taisei Corp | Method and system for predicting vibration |
| CN104715138A (en) * | 2015-01-28 | 2015-06-17 | 西南交通大学 | Rail fastener counter-force computing method and subway vibration frequency domain response computing method |
| CN104598757A (en) * | 2015-02-12 | 2015-05-06 | 西南交通大学 | Method for predicting traffic noise of rail regions |
| CN107589686A (en) * | 2016-07-09 | 2018-01-16 | 韩宇龙 | A kind of vision and vibration array perceive the intelligent safety defense monitoring system of integration |
| CN112767575A (en) * | 2021-01-19 | 2021-05-07 | 同济大学 | Vehicle-mounted adding instrument and method for rail transit operation period based on smart phone |
| CN114086967A (en) * | 2021-11-22 | 2022-02-25 | 中国地震局地震研究所 | Anti-excavation vibration monitoring method for subway tunnel |
| CN114519286A (en) * | 2022-02-17 | 2022-05-20 | 中铁二院工程集团有限责任公司 | Rapid calculation method for secondary structure noise caused by vibration of rail transit vehicle |
| CN115510711A (en) * | 2022-09-29 | 2022-12-23 | 东南大学 | Subway river-crossing tunnel vibration noise characteristic analysis method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116295804A (en) | 2023-06-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Abulizi et al. | Measuring and evaluating of road roughness conditions with a compact road profiler and ArcGIS | |
| CN104792937B (en) | Bridge head bump detection evaluation method based on vehicle-mounted gravitational acceleration sensor | |
| Sun et al. | Modeling carbon emissions from urban traffic system using mobile monitoring | |
| Wang et al. | Identification of moving vehicle parameters using bridge responses and estimated bridge pavement roughness | |
| Yang et al. | Impact of accessibility on housing prices in Dalian city of China based on a geographically weighted regression model | |
| CN116295804B (en) | Method for constructing vibration map | |
| WO2019148950A1 (en) | Method and device for analyzing nonlinear seismic response history of urban building group | |
| CN115034697B (en) | Multi-domain analysis-based multi-element intelligent compaction index grading optimization method and system | |
| CN111259605A (en) | Soil landslide monitoring, early warning and evaluating method | |
| JP2003287574A (en) | System, method and program for predicting earthquake damage | |
| CN115808711B (en) | Earthquake disaster loss evaluation method and system based on intensity quick report | |
| CN112926110A (en) | Real-time visual early warning method for risks in construction process of subway station | |
| CN111290295A (en) | Decision support system for wheel-rail interface lubrication and friction control | |
| CN102636364A (en) | Vehicular safety monitoring system for shapes and structures of bridge floors and detection method | |
| CN116307379A (en) | Urban rail transit planning method | |
| CN104929018B (en) | Method for calculating construction times of road roller | |
| CN116303855A (en) | Vibration map management system | |
| CN116414933A (en) | Urban rail transit vibration map management system | |
| CN113447220B (en) | Analogy prediction method and system for indoor vibration of subway vehicle section upper cover building | |
| CN110070234B (en) | Earthquake landslide personnel death number prediction method and application thereof | |
| Zhou et al. | Dynamic interactions between spatial change and travel behavior variation in old town fringe | |
| Abulizi et al. | Measurement and Comparison Study of Road Roughness Conditions by using Two Compact Road Profiler for Evaluating Ride Quality | |
| CN114673051B (en) | Roadbed compaction quality continuous detection system, roadbed compaction quality continuous detection method and roadbed compaction quality continuous detection application | |
| Li et al. | An experimental study on dynamic response of cement concrete pavement under vehicle load using IoT MEMS acceleration sensor network | |
| NUERAIHEMAITIJIANG | Efficient Road Roughness Monitoring Using Mobile Profiler and GIS in Establishing Pavement Management System |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |