CN101906758A - Method for determining intermittent force of sound barrier in rapid transit railway and application - Google Patents
Method for determining intermittent force of sound barrier in rapid transit railway and application Download PDFInfo
- Publication number
- CN101906758A CN101906758A CN2010102256688A CN201010225668A CN101906758A CN 101906758 A CN101906758 A CN 101906758A CN 2010102256688 A CN2010102256688 A CN 2010102256688A CN 201010225668 A CN201010225668 A CN 201010225668A CN 101906758 A CN101906758 A CN 101906758A
- Authority
- CN
- China
- Prior art keywords
- sound barrier
- pulsating force
- speed
- speed railway
- train
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
The invention relates to a method for determining an intermittent force of a sound barrier in a rapid transit railway, which comprises the following steps of: adopting relevant parameters: acquiring the distance from the center line of rapid transit railway rails to the sound barrier, acquiring train speed; and determining the intermittent force of the sound barrier: in the intermittent force of the sound barrier in the rapid transit railway, P represents the intermittent force of the sound barrier, v represents the train speed, and D represents the distance from the center line of the rapid transit railway rails to the sound barrier. In the method for determining the intermittent force of the sound barrier in the rapid transit railway and the application, the motion boundary is calculated by adopting a flow field analysis method, constructing the rapid transit railway model and carrying out mixed grid optimization with a movable grid technology for the operating train in the rapid transit railway. The invention can accurately determine the intermittent force of the sound barrier in the rapid transit railway by using the fitting intermittent force of the sound barrier simply and practically.
Description
Technical field
The present invention relates to a kind of definite method and application of sound barrier pulsating force size, relate in particular to the definite method and the application of sound barrier pulsating force size in a kind of high-speed railway.
Background technology
On October 1st, 1964, the first in the world row bullet train slowly starts at the station, Tokyo, and accelerate to 210km/h gradually, go at express speed towards the terminal and go, the pessimistic view of " railway is in the setting sun epoch " has thoroughly been pulverized in the successful operation of the Shinkansen train, and from then on world's railway steps into the high speed epoch.From then on after, countries in the world are used for reference Japanese experience one after another, build up the high-speed railway of oneself successively, adopt and created many advanced persons' technology.The process of railway high speedization has been accelerated in the variation that China's five speed raisings since nineteen ninety-six make railway produce essence.On April 18th, 2007, Chinese Railway has been carried out the 6th speed-raising and has been changed the line map, and " Harmony " CRH series EMUs begin the speed operation with 200km/h, and part section highest running speed reaches 250km/h, and from then on Chinese Railway enters the high speed epoch.High-speed railway has following technical advantage: road speed height, capacity are big, safe, round-the-clock running, energy consumption are low, save the area, environment pollution is light, comfort level is high, these advantages make high-speed railway obtain immense success in competition, have obtained tangible economic and social benefit.
For the measurement by experiment of determining mainly of sound barrier pulsating force size in the high-speed railway, there is not unified theoretical formula to solve the size of determining sound barrier pulsating force in the high-speed railway at present.
Summary of the invention
The technical problem that the present invention solves is: definite method of sound barrier pulsating force size in a kind of high-speed railway is provided, overcomes the technical problem that prior art can not accurately be determined sound barrier pulsating force size in the high-speed railway.
Technical scheme of the present invention is: definite method of sound barrier pulsating force size in a kind of high-speed railway is provided, comprises the steps:
Adopt relevant parameter: gather the distance of high speed railway track center line, gather train speed to sound barrier; Determine the sound barrier pulsating force: sound barrier pulsating force in the high-speed railway
Wherein: P represents the sound barrier pulsating force; V represents train speed; D represents the distance of high speed railway track center line to sound barrier.
Further technical scheme of the present invention is: in definite sound barrier pulsating force step, also comprise and adopt the high-speed motion of dynamic mesh modeling train in the flow field to calculate moving boundaries.
Further technical scheme of the present invention is: in definite sound barrier pulsating force step, in the process that adopts the high-speed motion of dynamic mesh modeling train in the flow field, the motor area of the non-structural network of tetrahedron is set.
Further technical scheme of the present invention is: in definite sound barrier pulsating force step, in the process that adopts the high-speed motion of dynamic mesh modeling train in the flow field, generate near wall boundary layer grid for the surface of the bullet train in the motor area.
Further technical scheme of the present invention is: in definite sound barrier pulsating force step, computational fields is carried out the subregion grid dividing, adopt the hybrid grid division methods.
Further technical scheme of the present invention is: in definite sound barrier pulsating force step, comprise and adopt the hybrid grid model to realize the bullet train motion analysis.
Further technical scheme of the present invention is: in definite sound barrier pulsating force step, comprise the aerodynamic force curve that obtains the sound barrier surface.
Further technical scheme of the present invention is: in definite sound barrier pulsating force step, obtain the high speed railway track center line to the influence to the sound barrier pulsating force of the distance of sound barrier.
Technical scheme of the present invention is: definite method of sound barrier pulsating force size in the high-speed railway is applied to high-speed railway.
Technique effect of the present invention is: by the train that moves in the high-speed railway being adopted the method for flow field analysis, and make up the bullet train model, carry out hybrid grid optimization with the dynamic mesh technology and calculate moving boundaries.The sound barrier pulsating force of match of the present invention can accurately be determined the size of sound barrier pulsating force in the high-speed railway, and is simple and practical.
Description of drawings
Fig. 1 is a flow chart of the present invention.
The hybrid grid model computational fields that Fig. 2 optimizes for the present invention.
The hybrid grid model computational fields grid dividing figure that Fig. 3 optimizes for the present invention.
Fig. 4 is train dimensional drawing in the model of the present invention.
Fig. 5 is a grid sectional drawing of the present invention.
Fig. 6 is a bullet train of the present invention border arrow diagram layer by layer.
Fig. 7 divides for bullet train surface mesh of the present invention.
Fig. 8 is that the present invention is by crossing pneumatic the trying hard to of train to building or member generation.
Fig. 9 is the pulsating force curve map of 200km/h for the speed of a motor vehicle of the present invention.
Figure 10 is the pulsating force curve map of 250km/h for the speed of a motor vehicle of the present invention.
Figure 11 is the pulsating force curve map of 300km/h for the speed of a motor vehicle of the present invention.
Figure 12 is the pulsating force curve map of 350km/h for the speed of a motor vehicle of the present invention.
Figure 13 is the different speed of a motor vehicle pulsating force of the present invention amplitude comparison diagram.
The pulsating force positive peak comparison diagram of different distance when Figure 14 is 200km/h for the speed of a motor vehicle of the present invention.
The pulsating force positive peak comparison diagram of different distance when Figure 15 is 250km/h for the speed of a motor vehicle of the present invention.
The pulsating force positive peak comparison diagram of different distance when Figure 16 is 300km/h for the speed of a motor vehicle of the present invention.
The pulsating force positive peak comparison diagram of different distance when Figure 17 is 350km/h for the speed of a motor vehicle of the present invention.
Figure 18 is the present invention pulsating force formula calculated value and the temporary figure that compares that advises when D=6m.
Figure 19 is the present invention pulsating force formula calculated value and the temporary figure that compares that advises when D=3m.
The specific embodiment
Below in conjunction with specific embodiment, technical solution of the present invention is further specified.
As shown in Figure 1, the specific embodiment of the present invention is: definite method of sound barrier pulsating force size in a kind of high-speed railway is provided, comprises the steps:
Step 100: adopt relevant parameter, that is, gather the distance of high speed railway track center line, gather train speed to sound barrier.
Step 200: determine the sound barrier pulsating force: sound barrier pulsating force in the high-speed railway
Wherein: P represents the sound barrier pulsating force; V represents train speed; D represents the distance of high speed railway track center line to sound barrier.
Among the present invention, sound barrier pulsating force in the high-speed railway
Wherein: P represents the sound barrier pulsating force; V represents train speed; D represents the distance of high speed railway track center line to sound barrier.
Detailed process is as follows:
One, train exterior flow field characteristic
The flowing gas state in the outer flow field of train mainly contains two kinds of forms: laminar flow and turbulent flow, laminar flow are meant that the mobile of fluid is layering, do not mix mutually between promptly two-layer, do not disturb mutually; Turbulent flow then in contrast, in turbulent flow, velocity component also is superimposed with turbulent flow pulsation at random on its average.The boundary layer can be differentiated by critical Reynolds number by the condition that laminar flow is converted into turbulent flow, and the definition of critical Reynolds number is as follows:
In the formula: x
CrBe the laminar flow distance in the boundary layer; u
∞Be speed of incoming flow; υ is the kinematic viscosity coefficient of fluid, and kinematic viscosity coefficient equals the dynamic viscosity coefficient of fluid and the ratio of density, that is:
Wherein, μ is the dynamic viscosity coefficient of fluid, and ρ is the density of fluid.For the flat-plate flow field, Re ≈ 5 * 10
5Therefore, judge the flow behavior in flow field usually with the size of Reynolds number, for train, the definition of reynolds number Re is as follows:
In the formula: u
∞Be train speed; L is a characteristic size.
When the minimum speed of a motor vehicle of train is 200km/h, be characteristic size with the train width, the train width is 3.2m in the embodiment of the invention, kinematic coefficient of viscosity is 15 * 10
-6m
2/ s, calculating as can be known, the outer flow field of bullet train Reynolds number is about 1.18 * 10
7, therefore, the outer flow field of bullet train is in turbulence state, and computation model should be selected turbulence model for use.
Two, the fundamental equation of fluid motion
The fundamental equation that hydrodynamics control flow field fluid flows is continuity equation, Navier-Stokes equation and the energy conservation equation that obtains by convection cell micro unit application quality conservation principle, principle of conservation of momentum and conservation of energy principle.Do not consider the flow field variations in temperature in the bullet train model calculates, its fluid governing equation is continuity equation and Navier-Stokes equation.
(1) continuity equation
If incompressible fluid, density are constant, continuity equation is:
In the formula, v
x, v
y, v
zFor flow field velocity at three coordinate x, y, the component of z direction.
(2) Navier-Stokes equation
For incompressible fluid, fluid dynamic coefficient of viscosity μ promptly can be used as constant, obtains the Navier-Stokes equation that incompressible viscous fluid flows:
In the formula, v
x, v
y, v
zFor flow field velocity at three coordinate x, y, the component of z direction; f
x, f
y, f
zBe x, y, unit mass fluid volume power on the z direction; ρ is a fluid density; P is a fluid pressure; υ is the kinematic viscosity coefficient of fluid.
(3) state equation:
Have according to perfect gas:
p＝ρRT (7)
In the formula, p is a fluid pressure; ρ is a fluid density; R is a mol gas constant; T is a gas temperature.
In theory, by above equation and specified fringe conditions and primary condition, can determine train velocity field, pressure field and thermal field on every side.But, since turbulent flow be a kind of at random, non-permanent, Three-dimensional Flow, the scrambling of fluid particle track brings very big difficulty for the research of turbulent motion rule in the turbulent flow.It is generally acknowledged that the Navier-Stokes equation can be described the instantaneous rule of turbulent flow, therefore, Fluid Mechanics Computation worker often can utilize the numerical simulation of Navier-Stokes equation to solve engineering problem.
Three, Fluid Mechanics Computation
Fluid Mechanics Computation (Computational Fluid Dynamics is called for short CFD) is the special tool that is used for carrying out flow field analysis, calculating, prediction.It is instrument with the computer, and with the mathematical method of various discretizations, all kinds of problems of convection cell mechanics are analyzed and researched.Briefly, CFD is equivalent to the virtual experiment of doing in computer, come the mobility status of analog simulation real fluid with this, and then finds the solution the differential equation of control fluid, draws the Discrete Distribution of flow field on the continuum.
Four, calculate moving boundaries.
In definite sound barrier pulsating force step, also comprise and adopt the high-speed motion of dynamic mesh modeling train in the flow field to calculate moving boundaries.
The dynamic mesh technology is used to calculate the moving boundaries problem, and the mobile problem that certain object is arranged in border or the computational fields, therefore, should adopt the dynamic mesh model to carry out analog computation for the high-speed motion of train in the flow field herein.Because FLUENT carries out the reconstruction of computational fields automatically according to moving of border or object, all will remove to divide the grid of whole computational fields in other words again at each time step, this for the requirement of computer resource than higher.In the dynamic mesh technology, the generation of grid and motion are to re-generate the change that adapts to the zoning by stretching, compression or increase, minimizing and part.
(1) the dynamic mesh setting of bullet train motion analysis.
1, the composition of computational fields.
Because bullet train is long, volume ratio is bigger, the corresponding calculated zone also can be very big, number of grid can increase greatly so virtually, and, the velocity ratio that train moves is very fast, the distortion meeting of grid is very big, if when upgrading grid with elasticity fairing method, mesh quality can be very poor, not only computational accuracy is low, might calculate not convergence, therefore, must adopt dynamic layered technology to upgrade grid, but dynamic layered method is relatively harsher to the requirement of trellis-type, therefore needs to adopt a kind of new grid subregion and division methods.The present invention adopts the hybrid grid method, has optimized the simulate effect of dynamic mesh to the nonstationary flow field, shown in Fig. 4-3.Hybrid grid comprises following a few part:
As shown in Figure 2, the motor area of the non-structured grid of tetrahedron is set, it is set in software does rectilinear motion with train.Because the airflow design on train surface, make the moving region that surrounds train form irregular shape (promptly no longer becoming a cube), therefore can only divide with tetrahedron, if upgrade grid with dynamic layered technology, the grid adjacent with moving boundaries is necessary for structured grid, in order to solve this contradiction, contacted former and later two faces in motor area and deformed area are divided with quadrangle, whole moving region is divided with tetrahedron, before and after the motor area, respectively form one deck pyramid grid like this, it is the pentahedron grid, inside, motor area is tetrahedral grid, the deformed region of front and back is divided with hexahedral mesh, so, the interface of motor area and deformed area is a quadrilateral mesh, and the grid of deformed area changes along with moving of motor area, selects for use dynamic layered method to define the distortion of this area grid and re-generates.In addition, generate near wall boundary layer grid for the surface of the bullet train in the motor area with Boundary Layer, Boundary Layer order is used for direction and the size in nearly wall mandatory provision grid growth, it can control near wall create-rule, high-quality grid, to guarantee the computational accuracy of key area.
In order to improve computational efficiency, the distortion of the mesh workload that minimizing is less to the result of calculation meaning, suitably be set to the quiescent centre in distance zone in addition apart from train in the computational fields, the hexahedron structure grid is adopted in the quiescent centre, dynamic mesh utilizes the slip grid function of Grid interface to be connected with Jing Wanggejiaojiequ, promptly allows fluid to pass through in the mesh overlay district.
2, computational fields grid dividing.
Carry out subregion for computational fields and divide grid, adopt the hybrid grid division methods of optimizing.Its grid dividing schematic diagram as shown in Figure 3, wherein: body V1 is the motor area, body V2, V3 are the deformed area, body V4, V5 are the quiescent centre, face F1, F2 are the interfaces of motor area and deformed area, face F3, F4 are the deformed areas.Among Fig. 3, face is represented body, and line is represented face.
(2) dynamic mesh analytical procedure
Adopt the hybrid grid model of FLUENT software and optimization to realize that the detailed process of bullet train motion analysis is as follows:
1, in Gambit, set up the motor area, the deformed area, the geometrical model of quiescent centre, and divide each regional grid.
2, each regional fringe conditions is set respectively, the condition of continuity is exported each regional grid file, promptly different mesh files.
3, the grid file that will represent zones of different in Tgrid software calls, and carries out grid and merge, and merges into a grid file, is about to each zone and is merged into as a wholely, generates a mesh file.
4, in Fluent, read in grid file, computational methods are set, turbulence model, the physical property of fluid, information such as fringe conditions; This paper adopts three-dimensional incompressible non-permanent transient state N-S method, Realizablek-ε turbulence model, and fluid is an air, selects the self-defining air characteristics of FLUENT for use.
5, load UDF (User-Defined Function, User-Defined Functions) file, the dynamic mesh parameter is set; Promptly select for use dynamic layered method to upgrade grid.
6, parameter is found the solution in setting, selects the SIMPLE algorithm for use, and the second order upstreame scheme to reduce numerical error, improves computational accuracy; The convergence residual error is not less than 10
-3, utilize many CPU concurrent technique to find the solution.
Five, sound barrier panel surface pulsating force analysis.
(1) computation model and grid control.
Among the present invention, the bullet train model adopts a car and trailer directly to connect extension, and the headstock curved surface adopts a complete curved surface, ignores headlight, window, ridges such as door handle.
The analog computation parameter: train chief 200m, wide 3.2m, high 3.5m, as shown in Figure 4, track centerline is 4m to the distance of sound barrier, the high 4m of sound barrier, long 100m.The high 10m in zoning, wide 50m, long 530m.
Adopt the grid optimization method, the motor area of the non-structured grid of tetrahedron is set, the deformed area of hexahedron structure grid and quiescent centre, owing to obtain the blast on sound barrier panel surface, therefore, zone with the sound barrier place, the grid dividing of train region thinner, thicker for the grid dividing of some outer peripheral areas is under the situation that has guaranteed computational accuracy, control mesh quantity reduces computing time.Fig. 5 is the grid sectional drawing that cuts downwards along short transverse, and Fig. 6 is the boundary layer grid on train surface, and Fig. 7 is the grid dividing on train surface.
(2) fringe conditions is set
All CFD problems all will have fringe conditions, and so-called fringe conditions is exactly the mathematical physics condition that the flow field variable should satisfy on computation bound, has only given appropriate boundary condition, just can calculate separating of flow field.When utilizing FLUENT simulation bullet train in air, to move, be provided with downstream condition:
The entrance boundary: speed inlet (VELOCITY_INLET), in the value of entrance boundary given speed and other scalar attributes.
The outlet border: pressure export (PRESSURE_OUTLET) is used in the exit, flow field to the value of constant static-pressure and other scalar variables.
Sound barrier, the train wall, the bottom surface: wall condition (WALL), because turbulent flow develops into laminar flow in the near wall district, therefore need adopt the wall function method at the near wall district, the given value on the wall is incorporated into the source item of the discrete equation of interior nodes.
Lateral surface, end face: symmetrical fringe conditions (SYMMETRY), do not need given any parameter, only need to determine rational symmetric position.
The interface of dynamic grid and static grid: interface (INTERFACE), when two adjacent area coincidences, but interface two faces for overlapping still, and the grid node on two faces do not overlap mutually, then two faces that overlap will be defined as interface.
In addition, be rigid motion with the movement definition of train, its motion mode is finished by loading the UDF function, and UDF has mainly defined the speed and the direction of train running.
(3) pulsating force is calculated standard.
Propose in " the Beijing-Shanghai High-Speed Railway design temporary provisions " that China railways the 3rd surveying and designing institute writes to run at high speed train to the structure gas pressure by following calculating:
By crossing Pneumatic pressure and the aerodynamic suction that train causes, should form by a 5m long all-moving surface load+q and the long all-moving surface load-q of 5m, aerodynamic force is divided into horizontal gas power q
hWith vertical aerodynamic force q
vThe horizontal gas dynamic pressure can be looked into by the curve of Fig. 8 and get, and promptly for the different speed of a motor vehicle, track is looked into and got pulsating wind pressure to the different distance of sound barrier.Vertical gas power q
vCan be calculated as follows:
In the formula: q
hBe horizontal gas power (kN/m
2); D is that position is to wire center distance (m).
(4) interpretation of result.
1, aerodynamic force curve.
In cloth monitoring point, sound barrier panel surface, and its surface wind compressed into the row monitoring, Fig. 4-10 is to Fig. 4-the 13rd, and bullet train is at speed of a motor vehicle 200km/h, 250km/h, and 300km/h, during 350km/h, the time dependent aerodynamic force curve of monitoring point wind pressure.Wherein, when the speed of a motor vehicle is 200km/h, the position of monitoring point is 50m from the distance of sound barrier panel starting end, and when the speed of a motor vehicle was 250km/h, the position of monitoring point was 70m from the distance of sound barrier panel starting end, when the speed of a motor vehicle is 300km/h, the position of monitoring point is 80m from the distance of sound barrier panel starting end, and when the speed of a motor vehicle was 350km/h, the position of monitoring point was 80m from the distance of sound barrier panel starting end, the height of monitoring point all in the middle of plate, promptly is along plate height direction 2m place.
Can obtain from above-mentioned pulsating force curve: when bullet train enters sound barrier, during away from the monitoring point, the monitoring point blast is zero substantially, when headstock near the monitoring point, blast sharply increases, reach maximum value, the malleation peak point occurs, when headstock crosses the monitoring point, blast sharply reduces, reach minimum value, the negative pressure peak point occurs.Explain that theoretically sail from afar when bullet train, headstock the place ahead is because of the retardance of headstock, the air that causes the train front is piled up, flow velocity reduces, and the normal pressure of monitoring point increases gradually, when train during near the monitoring point, the malleation peak point appears, then air is to the shunting of train head both sides, and flow velocity raises, and pressure reduces to form local depression, the negative pressure peak point occurs, positive peaked minimum value is equal substantially.
2, compare with standard.
Provided the train that crosses in " Beijing-Shanghai High-Speed Railway design temporary provisions " to the aerodynamic force reference value that building or member produce, will advise temporarily with FLUENT result of calculation and make comparisons, as table 1 and shown in Figure 13.
Table 1 FLUENT malleation value result of calculation compares with temporary rule
| The speed of a motor vehicle (km/h) | 200 | 250 | 300 | 350 |
| FLUENT result (Pa) | 266.72 | 421.44 | 592.37 | 830.88 |
| Advise result (Pa) temporarily | 250 | 400 | 570 | 790 |
| Error | 6.0％ | 4.9％ | 3.7％ | 4.8％ |
As can be seen from Figure 13, the fluctuating wind malleation value of FLUENT analog computation and temporarily rule the result who gets that looks into more approaching, illustrated that it is feasible calculating pulsating force with the method for numerical simulation, further verified the correctness of numerical simulation, and, also as can be seen, the pulsating force peak value is approximate to become the secondary relation to increase with train speed from Figure 13.
3, track centerline is to the influence to pulsating force of the distance D of sound barrier.
Table 2 is that FLUENT compares with the pulsating force peak value size of advising temporarily under different D (track centerline is to the distance of sound barrier) value.
The different D values of table 2 FLUENT malleation value result of calculation down compare with temporary rule value
| Gauge D (m) | The speed of a motor vehicle (km/h) | Normal value (Pa) | FLUENT calculated value (Pa) |
| 3.5 | 200 | 320 | 345.77 |
| 3.5 | 250 | 500 | 538.73 |
| 3.5 | 300 | 725 | 768.25 |
| 3.5 | 350 | 1000 | 1033.06 |
| 4 | 200 | 250 | 266.72 |
| 4 | 250 | 400 | 421.44 |
| 4 | 300 | 570 | 592.37 |
| 4 | 350 | 790 | 830.88 |
| 4.5 | 200 | 205 | 211.37 |
| 4.5 | 250 | 335 | 324.09 |
| 4.5 | 300 | 475 | 464.17 |
| 4.5 | 350 | 650 | 637.19 |
| 5 | 200 | 170 | 168.87 |
| 5 | 250 | 275 | 283.74 |
| 5 | 300 | 390 | 393.54 |
| 5 | 350 | 540 | 522.91 |
Figure 14 to Figure 17 is that the speed of a motor vehicle is respectively 200km/h, 250km/h, and 300km/h, during 350km/h, FLUENT calculated value under the different gauge D values and temporary rule value comparison diagram.
By table 2 and Figure 14 to-Figure 17 as can be known, FLUENT and temporarily the rule result be more or less the same, and along with the increase of track centerline to the sound barrier distance, the pulsating force peak value decreases, from figure as can be seen, the pulsating force size roughly becomes negative quadratic relationship with gauge D.
Five, the match of pulsating force peak value formula.
By the pulsating force comparison diagram as can be seen, raising along with train speed, aerodynamic force is approximate to become the secondary relation to increase with the speed of a motor vehicle, and, in the research of train car compression wave, also show, pressure intensity of wave and amplitude are approximate observes square being directly proportional of train speed sum with train speed and 1/8, promptly two travel relatively train crossing the time, the air pressure that can cause another train meeting one side surface is undergone mutation, between a few tens of milliseconds, just occurring in succession, the peak of negative pressure, the pressure amplitude on observation train surface and train speed and 1/8 are observed square being directly proportional of train speed sum.
Can obtain the pressure amplitude and gauge is approximated to negative quadratic relationship from Figure 14 to Figure 17, therefore, suppose train by the time act on pulsating wind pressure on the sound barrier panel on the occasion of become positive quadratic relationship with the speed of a motor vehicle, become to bear quadratic relationship with gauge, set up following relational expression:
In the formula: v is a train speed, km/h; D is the distance of track centerline to sound barrier, m; Since when v=0, P=0, so a=0 by changing the value of independent variable, can obtain the series of analysis data, utilize the Origin mathematical software to carry out match again and obtain following formula at last:
Calculating the speed of a motor vehicle with above formula is 160km/h, 180km/h, 200km/h, 220km/h, 250km/h, 300km/h, 350km/h, D=6m, the pulsating wind pressure during D=3m on the occasion of, and compare with temporary rule, as table 3, table 4, Figure 18, shown in Figure 19:
Table 3 formula result of calculation compares (D=6m) with temporary rule
| The speed of a motor vehicle (km/h) | 160 | 180 | 200 | 220 | 250 | 300 | 350 |
| Advise result (Pa) temporarily | 80 | 100 | 120 | 160 | 200 | 280 | 380 |
| Formulae results (Pa) | 75.11 | 95.07 | 117.37 | 142.01 | 183.39 | 264.08 | 359.44 |
Table 4 formula result of calculation compares (D=3m) with temporary rule
| The speed of a motor vehicle (km/h) | 160 | 180 | 200 | 220 | 250 | 300 | 350 |
| Advise result (Pa) temporarily | 280 | 340 | 420 | 520 | 680 | 960 | 1300 |
| Formulae results (Pa) | 300.46 | 380.27 | 469.47 | 568.05 | 733.54 | 1056.30 | 1437.74 |
From Figure 18, Figure 19 as can be seen, calculated value and temporarily the rule value is more approaching, therefore, and can be with this formula preresearch estimates in the different speed of a motor vehicle, under the different gauge effects, the pulsating force on sound barrier panel surface.
The specific embodiment of the present invention is: definite method of sound barrier pulsating force size in the high-speed railway is applied to high-speed railway.
Above content be in conjunction with concrete preferred implementation to further describing that the present invention did, can not assert that concrete enforcement of the present invention is confined to these explanations.For the general technical staff of the technical field of the invention, without departing from the inventive concept of the premise, can also make some simple deduction or replace, all should be considered as belonging to protection scope of the present invention.
Claims (9)
1. definite method of sound barrier pulsating force size in the high-speed railway comprises the steps: to adopt relevant parameter: gather the distance of high speed railway track center line to sound barrier, gather train speed; Determine the sound barrier pulsating force: sound barrier pulsating force in the high-speed railway
Wherein: P represents the sound barrier pulsating force; V represents train speed; D represents the distance of high speed railway track center line to sound barrier.
2. definite method of sound barrier pulsating force size is characterized in that in the high-speed railway according to claim 1, in definite sound barrier pulsating force step, also comprises and adopts the high-speed motion of dynamic mesh modeling train in the flow field to calculate moving boundaries.
3. definite method of sound barrier pulsating force size in the high-speed railway according to claim 2, it is characterized in that, in definite sound barrier pulsating force step, in the process that adopts the high-speed motion of dynamic mesh modeling train in the flow field, the motor area of the non-structural network of tetrahedron is set.
4. definite method of sound barrier pulsating force size in the high-speed railway according to claim 2, it is characterized in that, in definite sound barrier pulsating force step, in the process that adopts the high-speed motion of dynamic mesh modeling train in the flow field, generate near wall boundary layer grid for the surface of the bullet train in the motor area.
5. definite method of sound barrier pulsating force size is characterized in that in the high-speed railway according to claim 1, in definite sound barrier pulsating force step, computational fields is carried out the subregion grid dividing, adopts the hybrid grid division methods.
6. definite method of sound barrier pulsating force size is characterized in that in the high-speed railway according to claim 1, in definite sound barrier pulsating force step, comprises and adopts the hybrid grid model to realize the bullet train motion analysis.
7. definite method of sound barrier pulsating force size is characterized in that in the high-speed railway according to claim 1, in definite sound barrier pulsating force step, comprises the aerodynamic force curve that obtains the sound barrier surface.
8. definite method of sound barrier pulsating force size is characterized in that in the high-speed railway according to claim 1, in definite sound barrier pulsating force step, obtains the high speed railway track center line to the influence to the sound barrier pulsating force of the distance of sound barrier.
9. a high-speed railway of using definite method of sound barrier pulsating force size in the high-speed railway is characterized in that, definite method of sound barrier pulsating force size is applied to high-speed railway in the high-speed railway.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010102256688A CN101906758A (en) | 2010-07-12 | 2010-07-12 | Method for determining intermittent force of sound barrier in rapid transit railway and application |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010102256688A CN101906758A (en) | 2010-07-12 | 2010-07-12 | Method for determining intermittent force of sound barrier in rapid transit railway and application |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN101906758A true CN101906758A (en) | 2010-12-08 |
Family
ID=43262334
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2010102256688A Pending CN101906758A (en) | 2010-07-12 | 2010-07-12 | Method for determining intermittent force of sound barrier in rapid transit railway and application |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101906758A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105243224A (en) * | 2015-10-28 | 2016-01-13 | 唐山轨道客车有限责任公司 | Parameterized simulation method and apparatus for flow field in guest room, and train |
| CN107560657A (en) * | 2017-09-26 | 2018-01-09 | 北京铁路局保定工务段 | A kind of sound barrier monitoring system and method |
| CN108021760A (en) * | 2017-12-08 | 2018-05-11 | 中国电建集团中南勘测设计研究院有限公司 | A kind of fluid mean camber distribution of force computational methods based on numerical simulation |
| CN111104764A (en) * | 2020-01-21 | 2020-05-05 | 湖南科技大学 | Structured grid optimization division method for alternating current motor rotor conducting bar thermal analysis model |
| CN111222202A (en) * | 2020-01-19 | 2020-06-02 | 智慧航海（青岛）科技有限公司 | Automatic division method for ship body grids based on virtual test platform |
-
2010
- 2010-07-12 CN CN2010102256688A patent/CN101906758A/en active Pending
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105243224A (en) * | 2015-10-28 | 2016-01-13 | 唐山轨道客车有限责任公司 | Parameterized simulation method and apparatus for flow field in guest room, and train |
| CN107560657A (en) * | 2017-09-26 | 2018-01-09 | 北京铁路局保定工务段 | A kind of sound barrier monitoring system and method |
| CN108021760A (en) * | 2017-12-08 | 2018-05-11 | 中国电建集团中南勘测设计研究院有限公司 | A kind of fluid mean camber distribution of force computational methods based on numerical simulation |
| CN111222202A (en) * | 2020-01-19 | 2020-06-02 | 智慧航海（青岛）科技有限公司 | Automatic division method for ship body grids based on virtual test platform |
| CN111104764A (en) * | 2020-01-21 | 2020-05-05 | 湖南科技大学 | Structured grid optimization division method for alternating current motor rotor conducting bar thermal analysis model |
| CN111104764B (en) * | 2020-01-21 | 2020-08-25 | 湖南科技大学 | Structured grid optimization division method for alternating current motor rotor conducting bar thermal analysis model |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | The performance of different turbulence models (URANS, SAS and DES) for predicting high-speed train slipstream | |
| CN104866648B (en) | One kind is run at high speed traffic tool simulating model pre-treatment template and processing method | |
| Khayrullina et al. | CFD simulation of train aerodynamics: train-induced wind conditions at an underground railroad passenger platform | |
| Li et al. | Dynamic analysis of the flow fields around single-and double-unit trains | |
| Baker | A review of train aerodynamics Part 2–Applications | |
| Mohebbi et al. | Analysis of the effects of lateral wind on a high speed train on a double routed railway track with porous shelters | |
| Li et al. | Comparative study on the unsteady slipstream induced by a single train and two trains passing each other in a tunnel | |
| Liang et al. | Numerical simulation of pressure transients caused by high-speed train passage through a railway station | |
| CN101906758A (en) | Method for determining intermittent force of sound barrier in rapid transit railway and application | |
| CN102043871A (en) | Method for solving aerodynamic effect of long tunnel or tunnel group | |
| Yuan et al. | Study of the unsteady aerodynamic performance of an inter-city train passing through a station in a tunnel | |
| Yang et al. | Numerical study on transient aerodynamic behaviors in a subway tunnel caused by a metro train running between adjacent platforms | |
| CN101697173A (en) | Method for calculating micro-pressure waves generated in process of passing through tunnel by high-speed train | |
| Cui et al. | Investigation of train safety domain in cross wind in respect of attitude change | |
| Liang et al. | Effect of the ballast height on the slipstream and wake flow of high-speed train | |
| Wang et al. | Effect of surface roughness on the aerodynamics of a high-speed train subjected to crosswinds | |
| Zhou et al. | Effect of the wind speed on aerodynamic behaviours during the acceleration of a high-speed train under crosswinds | |
| Zhou et al. | Numerical investigation of the evolution of aerodynamic behaviour when a high-speed train accelerates under crosswind conditions | |
| Li et al. | Influence of blockage ratio on slipstreams in a high-speed railway tunnel | |
| Xiong et al. | Unsteady slipstream of a train passing through a high-speed railway tunnel with a cave | |
| Sun et al. | Study on the critical diameter of the subway tunnel based on the pressure variation | |
| Wang et al. | Study on the mitigation effect of a new type of connected structure on micro-pressure waves around 400 km/h railway tunnel exit | |
| Zhuo-ming et al. | Flow structure and far-field noise of high-speed train under ballast track | |
| CN101650757A (en) | Method for calculating human body aerodynamic characteristic values of side personnel under action of train wind | |
| Xue et al. | Influence of variable cross-section on pressure transients and unsteady slipstream in a long tunnel when high-speed train passes through |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C12 | Rejection of a patent application after its publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20101208 |