# Hydrodynamic Pressure and Velocity Distributions in the Interlayer Crack of Ballastless Track under High-Speed Train Load: A Theoretical Analysis

Authors: Shihao Cao; Shufang Zhai; Feng Dai; Shijie Deng; Hui Wang
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

In the abundant rain or poor drainage areas, there will be serious water damage existing in the interlayer of ballastless track. The essence of water damage of ballastless track is a dynamic development process of damage shape under the combined action of train load, hydrodynamic pressure, and flow velocity. In view of the distribution characteristics of hydrodynamic pressure and flow velocity in the interlayer crack of ballastless track under the high-speed train load, the simplified mechanics model for water and composite slab with interfacial crack is proposed in accordance with the water damage characteristics of ballastless track. Based on the conservation of mass and momentum theorem, the analytical expressions of water pressure and velocity in the saturated water crack were deduced. Similarly, the analytical expressions of water pressure and velocity in the unsaturated water crack were deduced by adding the state equation of ideal air. Considering that the water does not flow back timely, the analytical expressions of water pressure and velocity in another kind of unsaturated water crack were deduced. The results show that the hydrodynamic pressure and flow velocity in the interlayer crack are synthetically determined by multiple influencing factors such as fluid viscosity, load characteristic, crack shape, absolute pressure at the crack mouth, and initial volume of the air in the crack, and there is an intersecting phenomenon between influencing factors. When there is a small amount of air at the crack tip, the pressure and velocity distribution in the crack can be divided into three parts in terms of air-water interface and stagnation point. When the crack is filled with water, the hydrodynamic pressure tends to decrease along the direction of the crack mouth, and the distribution along the whole crack approximates to a cubic polynomial curve. Similarly, the flow velocity increases along the direction of crack mouth, and the distribution along the whole crack approximates to a quadratic polynomial curve.


---

## Body

## 1. Introduction

Ballastless track has become the main track structure of rapid transit railway because of its advantages of high comfort, high stability, and less maintenance. It has been developed rapidly and used widely. The design life of ballastless track is six decades which is estimated on the basis of the force characteristics acting on the track and the number of loads [1, 2]. However, the ballastless track, as the driving foundation, is exposed in the complex natural environment and bears the repeated action of train load, temperature load, rain erosion, and so on. These influencing factors make the component materials of the ballastless track in a continuous deterioration process, and the damage of ballastless track is inevitable [3–5]. The field investigation found that the damage rate in abundant rain or poor drainage areas is much faster than in dry areas; that is, the existence of water obviously reduces the life of the track.Due to the differences of structural types, environmental factors, and operational conditions of ballastless track in various countries, their damage forms and related research priorities are also different [6]. The Japanese slab track adopts unit structure and takes the train load as the main design load without considering the influence of temperature. During the long-term operation, the track is in good condition, and the main damage forms are resin and CA mortar damage. Through analyzing the formation mechanism of damage, it is found that the longitudinal and transverse deformation of the track slab caused by the train load and the temperature load is the main cause of damage [7, 8]. But the damage mechanism of the other ballastless track components has not been deeply studied.The ballastless track in Germany is a longitudinal continuous structure. Temperature load and concrete shrinkage are the main considerations. The longitudinal reinforcement of the continuous concrete track bed is set in consideration of the temperature variation and the concrete shrinkage [9–11]. Where there is a possibility of interlayer separation, the pins are added to ensure enough interlayer bonding. The construction quality is strictly controlled under construction. Furthermore, the cracks are blocked timely and attention is paid to the drainage design. So far, the research on the water damage of ballastless track in Germany has not been reported.In the early years, the integrated ballastless track in China was mainly used in tunnels, and the water damage of ballastless track caused by poor drainage of tunnels has ever happened [12]. Since the first comprehensive test section of ballastless track was built on Suining-Chongqing line in 2006, the ballastless track has developed rapidly and widely applied [13]. Because of the vast territory and large differences in geological and environmental conditions in China, there are many diseases on operational lines. These diseases include obvious damage caused by foundation deformation or extreme climate and also the water damage. Attention has been paid to the obvious damage caused by foundation deformation and extreme climate by the maintenance department and the researchers because of its clear cause and severe harm. A large number of related studies have been carried out at present [14–17]. The disease of ballastless track caused by water is a gradual process, and the initial damage will not cause serious consequences. Hence, not enough attention has been paid to this problem by the management department, and the related research has just started.Aiming at the water damage of ballastless track, the CRTS series ballastless tracks were investigated nationally. On the basis of the field investigation data, the apparent characteristics and main causes of the water damage were clearly defined. The simplified mechanics model for water and composite slab with interfacial crack was presented according to the water damage characteristics, and the analytical expressions of hydrodynamic pressure and flow velocity in interlayer crack under the dynamic load were derived. The research results can provide hydrodynamic pressure and flow velocity boundary conditions for the further research on the dynamic evolution process of water damage shape of ballastless track.


## 2. Characteristics and Control Factors of Water Damage of Ballastless Track

The drainage problem is considered at the design stage of track structure. In general, when the rain falls onto the surface of ballastless track, most of water will drain through drainage facilities, and the remaining small amount of water will not affect the track structure. Through the nationwide investigation of the service situation of ballastless track, it is found that cracks are unavoidable on the surface or layer of the ballastless track due to the improper construction and maintenance and the long-term effect of train load and environmental factors, as shown in Figures1(a) and 1(b). In the areas where the rain is rich or the drainage is poor, the water will seep into the cracks and become the precondition for the formation and development of the water diseases. Under the long-term soaking, the material mechanics properties of ballastless track will gradually decrease. A certain hydrodynamic pressure and flow velocity will be generated in the crack under the action of train load. Then the hydrodynamic pressure reacts on the crack surface and accelerates the crack propagation. Meanwhile, the flow velocity acting on the surface of the cracks will scour the surface particles repeatedly. The surface particles are exfoliated and result in a gradual increase of crack opening and the oozing slurry phenomenon at the crack mouth, as shown in Figures 1(c) and 1(d). Under the combined action of material softening, crack tip expansion, and crack surface scouring, the water damage of ballastless track will be further developed into more obvious diseases, such as sleeper loosening, void under track slab, or track bed, as shown in Figures 1(e) and 1(f).Figure 1
Development process of water damage: (a) surface crack; (b) interlayer crack; (c) oozing slurry at sleeper; (d) oozing slurry in interlayer crack; (e) sleeper loosening; (f) void under track bed.
(a) (b) (c) (d) (e) (f)


## 3. Simplified Mechanical Model of Fluid-Structure Interaction

Hydrodynamic pressure and velocity are the driving forces for the development of ballastless track water diseases, and clarifying their distribution characteristics is the precondition for the study of the development mechanism of water diseases. However, due to the complexity of fluid boundary conditions and the difficulty of numerical simulation, it is difficult to obtain the analytical expressions of three-dimensional hydrodynamic pressure and flow velocity distribution in the crack. In addition, the hydrodynamic pressure and flow velocity in the crack of ballastless track are synthetically determined by multiple factors, and the factors will have cross effects on each other [18], which means that the formation and variation of velocity are too complex to be analyzed with complete numeric models. Therefore, the complete numeric model usually needs to be simplified in order to analyze conveniently.Through the numerical analysis, it is found that the flow velocity distribution in the interlaminar crack of ballastless track is a three-dimensional problem, but the longitudinal and vertical components of flow velocity are significantly less than lateral component of flow velocity, which means that the effect of longitudinal and vertical components on velocity distribution is very limited. Therefore, it is also feasible to simplify the mechanical model into two dimensions when the small longitudinal flow of water (out-of-plane velocity) is ignored in this paper, as shown in Figure2, where L is the crack depth anda is the crack opening.Figure 2
Simplified mechanical model of fluid-structure interaction.In the model shown in Figure2, the crack is assumed to be located at the interface of the composite slab. The crack is filled with water, and the compressibility of water is neglected. Without considering the influence of the elastic foundation, the fixed constraint is applied at the bottom of the model.


## 4. Analytical Expression for Hydrodynamic Pressure and Velocity in Saturated Water Crack

As a common liquid in nature, the movement of water under the compression of crack surfaces must obey the most basic and universal motion laws, such as the mass conservation and the momentum theorem. For the incompressible fluid, the integral expression of the mass conservation law is [19, 20](1)dmdt=ρddt∫AdA+ρ∮Γu·ndΓ=0The integral expression of the momentum theorem is(2)dMdt=ρddt∫AudA+ρ∮Γuu·ndΓ=∑Fwhere m is the mass of liquid, t is the time, ρ is the density of liquid, A is the area of the control volume, u is the velocity vector of liquid, n is the outer normal direction of the control volume boundary, Γ is the control volume boundary, M is the momentum of liquid, and F is the force vector acting on the control volume.The concrete strength grade of ballastless track slab can reach up to C55 and belongs to high strength concrete. It can be known from the uniaxial tensile test that [21] when the stress is less than 70 percent of concrete tensile strength, the constitutive relation is basically linear. Actually, the design of ballastless concrete track slab adopts the allowable stress method, and its safety is too conservative [22]. Under the action of train load, the mechanical behavior of concrete track slab can still be approximately linear. In the elastic range, the upper surface displacement δ of the crack in Figure 3 is linearly related to the load F. That is,(3)δx=F·βxFigure 3
Schematic diagram of crack shape before and after loading.When the crack shape and load position are unchanged,β(x) is a determined distribution function.Before loading, the crack openingh at the position x is(4)hx=b+a-bxLAfter loading, the crack openingω at the position x is(5)ωx=b+a-bLx+F·βxWhen considering the effect of fluid viscosity, the flow state of water in the interlaminar crack of ballastless track can be determined by Reynolds number (Re). The flow state is laminar flow when the Reynolds number Re < 2300, and the flow state is turbulence flow when the Reynolds number Re > 2300. For the laminar flow, the velocityv of arbitrary section in crack is parabolic [23]; thus(6)vx,y=vmx1-2y-ωx2ω2xThe relation between the maximum flow velocityvm and the average velocity v¯ of arbitrary section is as follows:(7)v¯x=23vmxSelect control volume 1 (fluid domain:0~x) shown in Figure 3, combined with the formulas (5), (6), and (7). Then,(8)ρddt∫AdA=ρddt∫0xωxdx=ρF′∫0xβxdx(9)ρ∮Γu·ndΓ=ρ∫0ωxvx,ydy=23ρvmxωxSubstitute formulas (8) and (9) into formula (1), the following can be obtained: (10)vmx=-3F′2ωx∫0xβxdxEquation (10) is the analytical expression of water velocity distribution in crack under the dynamic load. According to the analytical expression, the water velocity in the crack is not only affected by the geometric shape of the crack, but also sensitive to the first derivative of the loading function.Control volume 2 shown in Figure3 (fluid domain: x~L) is selected as an object of study, and the forces acting on the control volume are shown in Figure 4.Figure 4
Schematic diagram of forces acting on control volume 2.HereτFSI is the shear stress in the fluid-structure interaction surface; pL is the pressure at the crack mouth; px is the water pressure at the abscissa of x.For Newton fluid, the relationship between shear stress and velocity gradient is linear [23], as shown in the following formula:(11)τx,y=μdvx,ydywhere μ is the dynamic viscosity coefficient.The shear stressτFSI along the boundary can be obtained by formulas (6) and (11):(12)τFSIx=4μvmxωxThe component of the shear force acting on control volume 2 along thex direction is(13)∫xLτFSIxdx+∫xLτFSIxdxcos⁡θcos⁡θ=12μF′∫xLdxω2x∫0xβxdxWhenδ<<h, the angle between two crack faces is approximately(14)tan⁡θ≈a-bLAccording to Newton’s third law, the normal pressure acting on the control volume of the solid boundary is equal to the fluid pressurep, and the component of the normal pressures acting on control volume 2 along the x direction is(15)∫xLpxdxcos⁡θsin⁡θ=a-bL∫xLpxdxThe component of the resultant force acting on control volume 2 along thex direction is(16)∑Fx=pxωx+a-bL∫xLpxdx-pLωL-12μF′∫xLdxω2x∫0xβxdxIt can be obtained according to formulas (5) and (6) that(17)ρddt∫AudA=ρddt∫xLv¯xωxdx(18)ρ∮Γuu·ndΓ=65ρv¯2LωL-v¯2xωxSubstitute formulas (16), (17), and (18) into formula (2), and obtain(19)a-bL∫xLpxdx+pxωx=-ρF′′Ax︸A+1.2ρF′2Bx︸B+12μF′Cx︸C+pLωL︸Dwhere(20)Ax=∫xLdx∫0xβxdxBx=∫0Lβx2ωL-∫0xβx2ωxCx=∫xLdxω2x∫0xβxdxFormula (19) is the analytical formula of the hydrodynamic pressure in interlayer crack of laminate under dynamic load. From the analytical expression, the hydrodynamic pressure not only is influenced by the geometric shape of the crack, the fluid viscosity, the fluid density, and the atmosphere, but also is sensitive to the first and second derivative of the loading function.


## 5. Analytical Expression for Hydrodynamic Pressure and Velocity in Unsaturated Water Crack

In fact, there may be a certain amount of air in the crack before loading; that is, the crack will not be full of water, which is called unsaturated water crack. Under the action of dynamic load, the flow and pressure distribution of fluids including water and air in the crack are shown in Figure5.Figure 5
Flow condition and pressure distribution of fluids including water and air in the crack.HereL is the crack depth, Lap is the air depth in the crack after loading, Lsta is the distance from stagnation point to the crack tip, pai is the initial pressure of air in the crack, pap is the pressure at the interface between air and water, and psta is the water pressure at the stagnation point.Under the action of dynamic load, the distribution of hydrodynamic pressure in saturated water crack is a typical two-phase coupling problem of liquid and solid, and the pressure distribution in the unsaturated water crack is a typical three-phase coupling problem of liquid, solid, and gas. Under the action of dynamic load, there will be a point with the flow velocityv=0 in the unsaturated water crack, which is called stagnation point, as shown in Figure 5. While loading, the liquid flows from stagnation point to two opposite directions and tends to increase along two opposite directions. The location of stagnation point is not fixed but will change with the initial volume of air in the crack. The shape change of the unsaturated water crack before and after loading is shown in Figure 6.Figure 6
Schematic diagram of crack deformation.Herea is the initial crack opening at the crack mouth, ω(x) is the crack opening at the position x after loading, δ(x) is the upper surface displacement of the crack at the position x after loading, Lai is the air depth in the unsaturated water crack at the initial time, and F is the force acting on the solid structure.


### 5.1. Analytical Expression of Hydrodynamic Pressure and Velocity at the Stage0~Lsta

At the initial time, the interface between air and water in unsaturated water crack is atx=Lai, and the initial volume Vai of the air field in crack can be obtained by(21)Vai=2bL-a-bxLai2LBefore loading, the initial pressure in the two regions (air and water) is equal to the atmospheric pressurepL because the atmosphere, water domain, and air domain exiting in the crack are in the static equilibrium state.After loading, the interface between air and water moves fromLai to Lap. On the basis of the incompressibility of water, it is known that volume 1 is equal to volume 2; that is,(22)∫LapLaiωxdx=∫LaiLstaδxdxThe volume change of the air region,ΔV, is as follows:(23)ΔV=∫0Lstaδxdx=F∫0LstaβxdxAccording to the state equation of ideal gas, the product of the pressure and volume of the gas is a constant when the mass and temperature of the gas remain unchanged [23]; that is,(24)pV=CCombined with formulas (21), (23), and (24), the pressure pair of the air region under dynamic load is(25)pair=pai·VaiVai+ΔV=pai·2bL+a-bxLai2bL+a-bxLai+2LF∫0LstaβxdxThe pressurepap at the interface between air and water is equal to the pressure in the air region.(26)pap=pair


### 5.2. Analytical Expression of Hydrodynamic Pressure and Velocity at the StageLsta~L

Select the control volume 1 shown in Figure6 (fluid domain: Lsta~x), combined with formula (1), (5), (6) and (7). The following can be obtained:(27)vmx=-3F′2ωx∫LstaxβxdxFormula (27) is the analytical expression of flow velocity at the stage Lsta~L in unsaturated water crack under the dynamic load.Selecting the control volume 2 shown in Figure6 (fluid domain: x~L) as a research object, the force diagram is shown in Figure 4. The component of resultant force acting on the control volume along the x direction is(28)∑Fx=pxωx+a-bL∫xLpxdx-pLωL-12μF′∫xLdxω2x∫LstaxβxdxSubstitute formulas (5), (27), and (28) into (2), and obtain(29)a-bL∫xLpxdx+pxωx=-ρF′′A1x︸A+1.2ρF′2B1x︸B+12μF′C1x︸C+pLωL︸Dwhere A1(x)=∫xLdx∫Lstaxβ(x)dx.(30)B1x=∫LstaLβx2ωL-∫Lstaxβx2ωxC1x=∫xLdxω2x∫LstaxβxdxFormula (29) is the analytical expression of hydrodynamic pressure at the stage Lsta~L in unsaturated water crack under dynamic load.


### 5.3. Analytical Expression of hydrodynamic Pressure and Velocity at the StageLap~Lsta

Similarly, control volume 3 shown in Figure6 was selected to obtain the analytical expression of flow velocity at the stage Lap~Lsta:(31)vmx=-3F′2ωx∫xLstaβxdxSelecting control volume 4, the analytical expression of hydrodynamic pressure at this stage can be obtained as follows:(32)a-bL∫xLpxdx-pxωx=-ρF′′A2x︸A+1.2ρF′2B2x︸B-12μF′C2x︸C-papωLap︸Dwhere (33)A2x=∫Lapxdx∫xLstaβxdxB2x=∫LapLstaβx2ωLap-∫xLstaβx2ωxC2x=∫Lapxdxω2x∫xLstaβxdxFormula (32) is the analytical expression of hydrodynamic pressure at stage Lap~Lsta in the unsaturated water crack under dynamic load.According to formulas (27) and (31), the flow velocity in unsaturated water crack is influenced by not only the crack geometry and load characteristics, but also the position of the stagnation point in the crack. The stagnation point is mainly determined by the initial volume of the air inside the crack. In terms of the formulas (29) and (32), the hydrodynamic pressure in unsaturated water crack is influenced by not only the geometric shape of the crack, the fluid viscosity, the fluid density, and the absolute pressure at crack mouth, but also the position of the stagnation point in the crack.When the initial volume of the air region in the crack is determined, the flow velocity distribution in the crack can be obtained by formulas (6), (27), and (31). Similarly, the hydrodynamic pressure distribution in the whole crack can be obtained by formulas (25), (29), and (32).


## 5.1. Analytical Expression of Hydrodynamic Pressure and Velocity at the Stage0~Lsta

At the initial time, the interface between air and water in unsaturated water crack is atx=Lai, and the initial volume Vai of the air field in crack can be obtained by(21)Vai=2bL-a-bxLai2LBefore loading, the initial pressure in the two regions (air and water) is equal to the atmospheric pressurepL because the atmosphere, water domain, and air domain exiting in the crack are in the static equilibrium state.After loading, the interface between air and water moves fromLai to Lap. On the basis of the incompressibility of water, it is known that volume 1 is equal to volume 2; that is,(22)∫LapLaiωxdx=∫LaiLstaδxdxThe volume change of the air region,ΔV, is as follows:(23)ΔV=∫0Lstaδxdx=F∫0LstaβxdxAccording to the state equation of ideal gas, the product of the pressure and volume of the gas is a constant when the mass and temperature of the gas remain unchanged [23]; that is,(24)pV=CCombined with formulas (21), (23), and (24), the pressure pair of the air region under dynamic load is(25)pair=pai·VaiVai+ΔV=pai·2bL+a-bxLai2bL+a-bxLai+2LF∫0LstaβxdxThe pressurepap at the interface between air and water is equal to the pressure in the air region.(26)pap=pair


## 5.2. Analytical Expression of Hydrodynamic Pressure and Velocity at the StageLsta~L

Select the control volume 1 shown in Figure6 (fluid domain: Lsta~x), combined with formula (1), (5), (6) and (7). The following can be obtained:(27)vmx=-3F′2ωx∫LstaxβxdxFormula (27) is the analytical expression of flow velocity at the stage Lsta~L in unsaturated water crack under the dynamic load.Selecting the control volume 2 shown in Figure6 (fluid domain: x~L) as a research object, the force diagram is shown in Figure 4. The component of resultant force acting on the control volume along the x direction is(28)∑Fx=pxωx+a-bL∫xLpxdx-pLωL-12μF′∫xLdxω2x∫LstaxβxdxSubstitute formulas (5), (27), and (28) into (2), and obtain(29)a-bL∫xLpxdx+pxωx=-ρF′′A1x︸A+1.2ρF′2B1x︸B+12μF′C1x︸C+pLωL︸Dwhere A1(x)=∫xLdx∫Lstaxβ(x)dx.(30)B1x=∫LstaLβx2ωL-∫Lstaxβx2ωxC1x=∫xLdxω2x∫LstaxβxdxFormula (29) is the analytical expression of hydrodynamic pressure at the stage Lsta~L in unsaturated water crack under dynamic load.


## 5.3. Analytical Expression of hydrodynamic Pressure and Velocity at the StageLap~Lsta

Similarly, control volume 3 shown in Figure6 was selected to obtain the analytical expression of flow velocity at the stage Lap~Lsta:(31)vmx=-3F′2ωx∫xLstaβxdxSelecting control volume 4, the analytical expression of hydrodynamic pressure at this stage can be obtained as follows:(32)a-bL∫xLpxdx-pxωx=-ρF′′A2x︸A+1.2ρF′2B2x︸B-12μF′C2x︸C-papωLap︸Dwhere (33)A2x=∫Lapxdx∫xLstaβxdxB2x=∫LapLstaβx2ωLap-∫xLstaβx2ωxC2x=∫Lapxdxω2x∫xLstaβxdxFormula (32) is the analytical expression of hydrodynamic pressure at stage Lap~Lsta in the unsaturated water crack under dynamic load.According to formulas (27) and (31), the flow velocity in unsaturated water crack is influenced by not only the crack geometry and load characteristics, but also the position of the stagnation point in the crack. The stagnation point is mainly determined by the initial volume of the air inside the crack. In terms of the formulas (29) and (32), the hydrodynamic pressure in unsaturated water crack is influenced by not only the geometric shape of the crack, the fluid viscosity, the fluid density, and the absolute pressure at crack mouth, but also the position of the stagnation point in the crack.When the initial volume of the air region in the crack is determined, the flow velocity distribution in the crack can be obtained by formulas (6), (27), and (31). Similarly, the hydrodynamic pressure distribution in the whole crack can be obtained by formulas (25), (29), and (32).


## 6. Analytical Expression for Hydrodynamic Pressure and Velocity in Another Kind of Unsaturated Water Crack

For the saturated water crack as shown in Figure2, the water will periodically outflow and flow back due to the extrusion and pumping action of crack surfaces, respectively. When the water can not flow back timely, there is going to be another kind of unsaturated water crack, as shown in Figure 7.Figure 7
Schematic diagram of crack deformation.Herec is the liquid level, J(x) is the distance between liquid level and crack upper surface at the position x, Lib is the position of intersection point between liquid level and crack upper surface before loading, and Ldb is the position of intersection point between liquid level and crack upper surface after loading. In order to distinguish from the crack mouth L, Ldb is defined as fluid outlet in this paper.According to the crack geometry and liquid boundary parameter in the crack,J(x) can be obtained:(34)Jx=hx-c=ωx+δx-cIn fact, the position ofLdb is not fixed, but moving under the dynamic load. The position should satisfy(35)JLdb=δLdbCombined with formulas (3), (4), (34), and (35), the following is obtained:(36)a-bLdbL+b-c=FβLdbIn this case, the hydrodynamic pressure and flow velocity distributions in the crack can be divided into two regions in terms of fluid outletLdb: 0~Ldb and Ldb~L. For the water in region Ldb~L,the value of hydrodynamic pressure is approximately equal to the atmospheric pressure due to the direct contact with the atmosphere.For the water in region0~Ldb, the hydrodynamic pressure and flow velocity distribution expressions can be derived on the basis of the same method as saturated water crack, as shown in Section 4 in this paper. Selecting the control volume 1 shown in Figure 7 (fluid domain: 0~x), the analytical expression of flow velocity can be obtained as follows:(37)vmx=-3F′2ωx∫0xβxdxSimilarly, selecting control volume 2 (fluid domain:x~Ldb), the analytical expression of hydrodynamic pressure can be obtained as follows:(38)a-bL∫xLdbpxdx+pxωx=-ρF′′A3x︸A+1.2ρF′2B3x︸B+12μF′C3x︸C+pLdbωLdb︸Dwhere (39)A3x=∫xLdbdx∫0xβxdxB3x=∫0Ldbβx2ωLdb-∫0xβx2ωxC3x=∫xLdbdxω2x∫0xβxdx


## 7. Simplification of Analytical Expression

When the crack is filled with water andb=a (the crack shape is rectangular), the analytical expression (38) of the hydrodynamic pressure can be transformed into formula (41) under the sinusoidal load shown in the following formula:(40)F=Fm+Fosin⁡2πftwhere Fm is the average load, Fo is the loading amplitude, and f is the loading frequency.(41)px=4π2ρf2Fosin⁡2πft/aA4x︸A+24π2ρf2Fo2cos2⁡2πft/5a2B4x︸B+24πμfFocos⁡2πft/a3A4x︸C+pLwhere (42)A4x=∫xLdx∫0xβxdxB4x=∫0Lβx2-∫0xβx2The itemA of formula (41) is produced by the mutual extrusion of the fluid. The term B is produced by the shape change of the crack relative to the initial shape, and the term C is caused by the fluid viscosity. When δ<<h, the hydrodynamic pressure generated by the change of the crack shape is very small. Therefore, the effect of the term B on the hydrodynamic pressure can be negligible.In order to further study the primary and secondary relationship between termsA and C in formula (41), the ratio R of the maximum value of term C to the maximum value of term A is supposed; that is,(43)R=6μπρa2fIn view of formula (43), it is known that when the fluid properties (fluid viscosity µ and density ρ) are certain, R is affected by the loading frequency f and the crack opening a together. With the increase of the loading frequency and the crack opening, the influence of viscosity decreases gradually. Through analysis of formula (43), it is found that the influence of the water viscosity on the pressure is negligible basically when the crack opening a > 2 mm under the high frequency dynamic load. In this case, the hydrodynamic pressure is mainly determined by term A of analytical expression (41), and formula (41) can be simplified to the following formula:(44)px=4π2ρf2Fosin⁡2πftaAx+pLFormula (44) indicates that the pressure p in the interlayer crack is linear with the loading amplitude Fo under the dynamic load and has a square relationship with the loading frequency f and is inversely proportional to the crack opening a.Similarly, substitute formula (40) into formula (37). The following can be obtained: (45)vmx=-3πfFocos⁡2πftaDxwhere D(x)=∫0xβxdx.Formula (45) shows that the flow velocity v in the crack is linear with the loading amplitude Fo under the dynamic load and has a linear relationship with the loading frequency f and is inversely proportional to the crack opening a.


## 8. A Numerical Example

In order to verify the correctness of the analytical formula of hydrodynamic pressure, the model shown in Figure8 was taken as a numerical example. The calculation model of hydrodynamic pressure is composed of two parts: solid domain (concrete specimen with interlayer crack) and fluid domain (water in crack), and its corresponding ANSYS-CFX model is shown in Figure 9. The mesh size of the solid domain far from the crack is 10 mm, while near the crack it is 0.25 mm. The mesh size of the fluid domain along the depth direction is 0.25 mm, and that at the vertical opening direction is less than 0.2 mm.Figure 8
Calculation model of hydrodynamic pressure in the interlayer crack of concrete (mm).Figure 9
Finite calculation model: (a) solid domain; (b) fluid domain.
(a) (b)The solid material of calculation model is C30 concrete, and its elasticity modulus and Poisson’s ratio are 30 GPa and 0.2, respectively. The fluid in the crack is water, and its densityρ=998.2 kg/m3; the dynamic viscosity coefficient µ=1.002×10-3 pa.s. The reference pressure at the crack mouth is 0 Pa. The peak value of the hydrodynamic pressure at different positions was obtained by formula (41) and ANSYS-CFX under the sinusoidal load with loading frequency f=10 Hz and loading amplitude Fo=1 kN, as shown in Table 1, where x represents the distance between the calculated point and the crack tip.Table 1
Comparison of analytical and numerical results.
x/m analytical results numerical results A m a x/Pa B m a x/Pa C m a x/Pa P m a x/Pa P m a x/Pa 0 67.37 0.2154 1.435 67.37 62.56 0.05 65.75 0.215 1.401 65.75 61.82 0.1 61.55 0.2106 1.311 61.55 58.55 0.15 51.69 0.1912 1.101 51.69 49.61 0.2 32.57 0.134 0.694 32.57 31.53 0.25 0 0 0 0 0.001As shown in Table1, the peak value of hydrodynamic pressure calculated by the item A is much larger than those calculated by items B and C of hydrodynamic pressure analytical expression (41). This phenomenon indicates that the hydrodynamic pressure in the crack is mainly determined by itemA of equation (41) under the dynamic load. The primary and secondary influence relation of all the items is consistent with the theoretical analysis results. As can be seen from Figure 10(a), the water pressure in the interlayer crack of concrete decreases along the direction of crack mouth, and the peak value appears at the crack tip. The distribution law can be explained based on the conservation of mass and momentum theorem. When the water is pressed by the crack surfaces, it is known from the mass conservation of the fluid that the flow velocity along the crack mouth will increase gradually, as shown in Figure 10(b). The increase of flow velocity also represents the increase of momentum. According to the momentum theorem, the direction of the resultant force acting on the control volume is in accordance with the direction of momentum increase. This is the reason that water pressure will decrease along the direction of crack mouth. In addition, the peak value of hydrodynamic pressure calculated by analytical expression is 67.37 Pa, and the numerical value is 62.56 Pa. The similarity between analytical and numerical results can reach 92.86%. The consistency between analytical and numerical simulation results shows that the analytical formulas of hydrodynamic pressure were correctly derived.Figure 10
Distribution of water pressure and velocity along the direction of crack mouth: (a) pressure; (b) velocity.


## 9. Conclusions

In this paper, the simplified mechanics model for water and composite slab with interfacial crack was proposed in accordance with the water damage characteristics of ballastless track. Based on the conservation of mass and momentum theorem, the analytical expressions of water pressure and velocity in the saturated water crack were deduced. Similarly, the analytical expressions of water pressure and velocity in the unsaturated water crack were deduced by adding the state equation of ideal air. Considering the water not flowing back timely, the analytical expressions of water pressure and velocity in another kind of unsaturated water crack were deduced. In order to verify the correctness of the theoretical derivation, the ANSYS Workbench was used to conduct fluid-structure interaction analysis, and the similarity between analytical and numerical results can reach 92.86%. According to the theoretical results, the primary and secondary influence relationship between multiple factors and the mechanism of intersecting action are determined. The research results can provide hydrodynamic pressure and flow velocity boundary conditions for further research on the dynamic evolution process of water damage shape of ballastless track.



---

*Source: 3187979-2018-08-19.xml*