Abstract
This paper introduces a threedimensional model for the steadystate response of a pavementsubgradesoft ground system subjected to moving traffic load. A semianalytical wave propagation model is introduced which is subjected to four rectangular moving loads and based on a calculation method of the dynamic stiffness matrix of the ground. In order to model a complete road system, the effect of a simple road model is taken into account including pavement, subgrade and soft subsoil. The pavement and the subgrade are regarded as two elastic layers resting on a poroelastic halfspace soil medium. The priority has been given to a simple formulation based on the principle of spatial Fourier transforms compatible with good numerical efficiency and yet providing quick solutions. The frequency wavenumber domain solution of the road system is obtained by the compatibility condition at the interface of the structural layers. By introducing FFT (Fast Fourier Transform) algorithm, the numerical results are derived and the influences of the observation coordinates, the load speed and excitation frequency, the permeability of the soft subsoil, and the rigidity of the subgrade on the response of the pavementsubgradesoft ground system are investigated. The numerical results show that the influences of these parameters on the dynamic response of the road system are significant.
1. Introduction
A growing attention on highway development has raised many questions about the comfort of transportation and the economic maintenance of systems, which consist of vehicles, road construction proper and subsoil. Rational design and maintenance procedures for road structures, such as pavement and subgrade, require elaborate predictive methods to assess the dynamic performance induced by moving traffic loads. On the other hand, with the rapid development of mechanical technology, vehicles on the highway can reach speeds of over 200 km/h. It is known that the Rayleigh wave velocity is slow in soft soils; thus, vehicles can reach or exceed the Rayleigh wave velocity. The subject of the dynamic response of road systems to moving loads is of great importance in the field of geotechnical and transportationrelated engineering.
A multilayered soil model is usually employed to simulate the road system. There are many kinds of analytical and numerical methods that can be used to calculate the dynamic responses of the layered road structure. For example, Kausel and Roesset [1] and Kausel and Peek [2] presented a simple technique based on layer stiffnesses to determine the displacements and stresses due to any dynamic loads on or within the layered elastic media. Chow et al. [3] applied the aforementioned methodology to predict the response of a layered elastic medium to a stationary strip footing subjected to a vertical or horizontal harmonic load. To consider a load in actual motion, De Barros and Luco [4, 5] developed a semianalytical model for a point load moving on a viscoelastic multilayered halfspace by determining the stiffness matrices and combining the stiffness matrices of each layer. This model was based on a double Fourier transform and showed the presence of Mach cones in the first layer at superRayleigh speeds. This method was then used by Gunaratine and Sanders in analysing a road [6] and by Auersh [7], Picoux et al. [8], Picoux and Houedec [9] in analysing track excitation. It has also been used by Jones and Petyt for rectangular loads or strip loads with different ground conditions: half space [10], elastic layers on a half space [11] or on a rigid foundation [12]. In order to assess the effect of road unevenness on the pavement responses, the lateral and vertical vibrations of moving vehicles were investigated by Skvireckas [13] and Lu [14].
All the studies mentioned above treated the soil medium as a singlephase elastic or viscoelastic solid, but ground water may exist in the soil medium and affect wave propagation during the passage of the car. Therefore, a fully saturated poroelastic soil model is superior to linear elastic or viscoelastic models for analyzing the dynamic response of pavement systems. Biot [15, 16] pioneered the development of an elastodynamic theory for a fluidfilled elastic porous medium. This theory, termed Biot’s theory, has wide applications in the geotechnical profession for analyzing wave propagation characteristics under dynamic loads. By applying Biot’s theory, Theodorakopoulos [17, 18] obtained the steadystate displacements and stresses of the poroelastic soil subjected to a line load under plane strain conditions with the relative motion between the solid and fluid phases. Jin [19] studied the dynamic response of a poroelastic halfspace generated by a highspeed point load under plane strain conditions using a model of an infinite Euler beam resting on a poroelastic halfspace. Zeng and Rajapakse [20], Jin [21], Senjuntichai et al. [22] and Cai et al. [23] examined the dynamic response of a plate on a poroelastic halfspace to a vertical harmonic point load and rectangular load. More recently, LefeuveMesgouez and Mesgouez [24, 25] developed a number of semianalytical tools to study ground vibrations induced by dynamic loads on multilayered poroviscoelastic media. Though these previous works have made great contributions to the development of the elastodynamic theory for a fluidfilled elastic porous medium, they did not study the problems of a complete road system including pavement, subgrade and soft ground subjected to moving traffic load.
A review of the literature shows that few studies have investigated the dynamic response of a pavementsubgradesoft ground system subjected to moving traffic loads, using both a singlephase elastic soil model and a saturated poroelastic soil model. Due to its importance to engineering practice, further investigation of the dynamic response of pavement systems subjected to moving traffic loads is needed.
In this paper, a semianalytical approach is used to investigate the dynamic response of a pavementsubgradesoft ground system subjected to moving traffic load. The traffic loads are simulated by four rectangular load pressures, and the pavement and the subgrade are regarded as two elastic layers resting on a poroelastic halfspace soil medium. The exact analytical solutions of the displacements, stresses, and accelerations of the road structures are obtained in the form of integral based on a calculation method of the dynamic stiffness matrix of the ground. The governing equations are solved by Fourier transform techniques, and, to calculate the inverse Fourier transform, the fast Fourier transform procedure is used. The numerical results show that a moving load with a high speed will generate a larger response in soft subsoil than a lower speed load, and the load excitation frequency, the permeability of the soft subsoil, and the rigidity of the subgrade affect the dynamic response significantly.
2. Mathematical model
Two dimensional models are not able to reproduce mechanisms of wave propagation in the ground correctly since the loading zone represents a reduced area compared to infinite ground surface, thus excluding the hypothesis of plane deformation. This limitation has motivated researchers to obtain solutions using 3D models (see Fig. 1).
2.1. Description of the model
The calculation method uses the Fourier transform formalism for a semianalytical solution in the wave number domain. The dynamic stiffness matrix for pavementsubgradesoft ground system is written and accompanied by a fitted phase argument in Helmoltz functions. This provides us a fast numerical approach to the problem. Equations are written in the wave number domain as well as for steady state solutions. In this paper, the steadystate response means all the components of response are assumed to be timeharmonic form, which is consistent with dynamic loading. An inverse numerical Fourier transform is then applied when the matrix equation of the whole system (pavementsubgradesoft ground system) has been solved.
Fig. 13D model of pavementsubgradesoft ground system and the tandemaxle loads
(a)
(b)
2.2. Analytical solution of singlephase elastic soil
For an elastic, homogeneous and isotropic layer, the equation of motion is:
where ${\mathbf{u}}^{e}$ is the vector of displacements for elastic soil, ${\rho}_{\text{e}}$ is the mass density of the solid and ${G}_{\text{e}}$, ${\lambda}_{e}$ are the Lame’s constants, and for the behaviour law:
where ${\mathbf{\tau}}^{e}$, ${\mathbf{\epsilon}}^{e}$ are the vector of stresses and strains for elastic soil, ${\delta}_{ij}$ is the Kronecker delta, and ${\theta}_{e}$ is the volumetric strain of the soil. To define elastic body waves ($P$ – and $S$ – waves), Helmholtz’s potentials $\varphi $ and $\psi $ are usually used, so that:
In this paper, the Fourier transforms with respect to $x$ and $y$ coordinate are defined as:
The inverse relationship is given by the following form:
where $\beta $, $\gamma $ are the wave number in the $x$ and $y$ direction, $f$ is a variable in space domain, and $\overline{f}$ is the corresponding variable in the transformed domain.
By using the Fourier transforms given in Eq. (4), Eqs. (1)(3) and the solid strain equation are transformed from partial differential equations to ordinary differential equations. The typical vehicle excitation is composed of a set of vertical components that symbolize the effect of each wheel on the pavement. Each of these forces is taken as harmonic. The “fitted phase angle” is defined by Helmholtz’s potentials [7], so that:
where ${\alpha}_{i}=\sqrt{{\beta}^{2}+{\gamma}^{2}{k}_{i}^{2}}$, $i=1,\mathrm{}2\text{;}$${k}_{i}=({\mathrm{\Omega}}_{0}\beta c)/{c}_{i}$,$\hspace{0.17em}i=1,\mathrm{}2\text{;}$${c}_{\text{1}}=\sqrt{({\lambda}_{e}+2{G}_{e})/{\rho}_{e}}$ and ${c}_{2}=\sqrt{{G}_{e}/{\rho}_{e}}$ are the compression ($P$waves) and shear wave ($S$waves) speed of singlephase elastic soil, respectively; ${\mathrm{\Omega}}_{0}$ is the excitation frequency, $c$ is the vehicle speed and $h$ is the thickness of the soil layer. Because of numerical difficulties due to exponential terms in the rigidity matrix, layers must be divided into several sublayers increasing the number of iterations. In fact, the role of the fitted phase translation is to solve this problem by excluding exponential terms (growing exponential) in the dynamic stiffness matrix thus providing good condition.
A change of variables (translation in the mobile reference linked to the load) and functions (evaluation of the steady state solution) as well as the use of spatial Fourier transform lead to the calculation of displacements and stresses in the wave number domain, as in Lu [26]:
where $({\overline{u}}_{x}^{e},{\overline{u}}_{y}^{e},{\overline{u}}_{z}^{e})$ is the transformed displacement vector and $({\overline{\tau}}_{xz}^{e},{\overline{\tau}}_{yz}^{e},{\overline{\tau}}_{zz}^{e})$ is the stress vector; ${A}_{1}~{A}_{6}$ are arbitrary constants.
2.3. Analytical solution of poroelastic halfspace
The soft subsoil medium is modeled as a uniform poroelastic halfspace, fully saturated by the viscous fluid. The fluid is free to flow throughout the upper surface of the soil. The governing equations of a fully fluidsaturated poroelastic soil, based on the assumption of incompressible solid grains and neglecting the apparent mass density, are given by Biot [16]:
The constitutive relationships of the saturated soil medium are given as follows:
where ${u}_{i}^{s}$ and ${w}_{i}^{s}$ are the displacements of the solid skeleton and the infiltration displacement of the pore fluid, respectively, in the $x$, $y$ and $z$ directions; ${G}_{s}$ and ${\lambda}_{s}$ are Lame constants; $\alpha $ and $M$ denote of the compressibility of the soil and fluid (under the assumption of incompressible solid grains, $\alpha =$1, $M={K}_{f}/n$); and ${K}_{f}$ is the bulk modulus of the fully saturated fluid. The parameters ${\eta}_{w}$, ${k}_{d}$ represent the fluid viscosity and the intrinsic permeability of the soil medium, respectively; $\rho =n{\rho}_{f}+(1n){\rho}_{s}$, where ${\rho}_{f}$ and ${\rho}_{s}$ are the mass densities of the fluid and solid and $n$ is the porosity. ${\tau}_{ij}^{s}$ is the total stress components of the bulk material, and ${p}_{s}$ is the pore water pressure.
By applying the Fourier transforms given in Eq. (4), Eqs. (13)(16) and the solid strain equation are transformed from partial differential equations to ordinary differential equations. By resolving the ordinary differential equations, one can obtain the responses of the displacements and the stresses of the poroelastic halfspace, as given by Lu [26]:
where ${A}_{7}~{A}_{10}$ are arbitrary constants determined by boundary conditions, other symbols in Eqs. (17) and (18) are detailed as follows:
It should be noted that, as $z\to \infty $, the dynamic responses become zero; therefore, the real part of ${\alpha}_{3}$, ${\alpha}_{4}$ and ${\alpha}_{5}$ must be positive.
2.4. Dynamic stiffness matrix of pavementsubgradesoft ground system
In this part, the previous model (Fig. 1) is used with two singlephase elastic layers on a poroelastic halfspace. For the pavement with thickness ${h}_{p}$, the displacements at the surface and bottom of the layer are written in terms of the constants ${A}_{1}$, ${A}_{2}$, ${A}_{3}$, ${A}_{4}$, ${A}_{5}$ and ${A}_{6}$ according to Eqs. (7)(9):
where the superscripts ‘(0)’ and ‘(1)’ denote the surface and bottom of the layer, respectively. $\left[{D}^{p}\right]$ is displacement matrix for pavement.
Similarly, for the same layer, stresses can be expressed in terms of the previous constants:
where $\left[{S}^{p}\right]$ is stress matrix for pavement.
According to Eqs. (24) and (25), a relation between displacements and stresses can be written into a matrix formulation:
where $\left[{K}^{p}\right]=\left[{S}^{p}\right]{\left[{D}^{p}\right]}^{1}$ is dynamic stiffness matrix for pavement.
Similarly, for the second layer of elastic medium, a relation between displacements and stresses can be written into a matrix formulation:
whereis $\left[{K}^{e}\right]=\left[{S}^{e}\right]{\left[{D}^{e}\right]}^{1}$ dynamic stiffness matrix for subgrade.
For the poroelastic halfspace, the displacements and stresses at the surface of the halfspace are written in terms of the constants ${A}_{7}$, ${A}_{8}$, ${A}_{9}$ and ${A}_{10}$ according to Eqs. (17)(18):
where the superscript ‘(2)’ denotes the surface of the halfspace. $\left[{D}^{s}\right]$, $\left[{S}^{s}\right]$ are displacement matrix and stress matrix for the poroelastic halfspace.
From Eqs. (28) and (29), the following relation between displacements and stresses can be derived:
where $\left[{K}^{s}\right]=\left[{S}^{s}\right]{\left[{D}^{s}\right]}^{1}$ is dynamic stiffness matrix for poroelastic halfspace.
Assuming that the contact conditions among the layers are smooth, the displacements and the stresses are continuous at the interfaces. Thus, the following continuity conditions hold at each interface:
As the contact surface between the subgrade and the poroelastic halfspace is assumed to be fully permeable, the following expression can be derived:
From Eqs. (26)(27) and Eqs. (30)(33), the 3D stiffness matrix is the combination of the stiffness matrices of each layer, the following relationship is given in a matrix form as:
Consequently, the displacement vector at the interface of each layer is given by:
where $\left[B\right]$ is the inverse stiffness matrix whose size is 10×10 of pavementsubgradesoft ground system.
2.5. Description of the traffic load
The shape of the tirepavement contact area is assumed to be square, as shown in Fig. 1, and any change in the shape during load variation is neglected. The variations of the load amplitude resulting from the pavement surface roughness and the mechanical system of the vehicle are neglected. The traffic load is modeled as four rectangular load pressures expressed as in Eq. (36):
where ${q}_{0}$ and ${Q}_{0}$ are the magnitude of the load pressure for the front and the rear wheel, respectively. $2{l}_{1}$, $2{l}_{2}$ are the loaded lengths of a tire print in the ${x}_{1}$ and ${y}_{1}$ directions, respectively, for the front wheel, and $2{L}_{1}$, $2{L}_{2}$ are the loaded lengths of a tire print in the ${x}_{1}$ and ${y}_{1}$ directions, respectively, for the rear wheel. $2{d}_{a}$ is the distance between the front and rear wheels, and $2{d}_{w}$ is the distance between the left and right wheels.
A change of variables (translation in the mobile reference linked to the load) and functions (evaluation of the steady state solution) as well as the use of spatial Fourier transform lead to the expression of traffic load in the wave number domain:
As no lateral and longitudinal stress is considered in the present problem, the following boundary conditions, omitting the common time factor, ${e}^{{i\mathrm{\Omega}}_{0}t}$, can be obtained at the surface of pavement:
Substituting Eqs. (37) and (38) into Eq. (35), the displacements in the transformed domain can be obtained as follows:
After the displacements at the interface of layers are determined, it is rather straightforward to obtain all the stresses at an arbitrary interface. Furthermore, the constants of each layer can be determined, and the stresses and displacements in an arbitrary depth of the road structures can also be derived.
Once the transformed stress and displacement components are determined, the actual medium response can be established by means of the Fourier inversion theorem from Eq. (5). It should be noted that all the exponential functions involved in the present formulation contain only nonpositive exponents, which are critical for numerical stability. Taking the surface of halfspace for example, the displacements in the time domain can be expressed as follows:
2.6. Comparison with existing work
To compute the inverse transform accurately with a discrete transform, the mesh of the calculated functions must be fine enough to accurately represent the details of the functions, and the integrals must be truncated at values high enough to avoid distorting the results. To satisfy the above requirements, a FFT algorithm is employed with a grid of 2048×2048 points with ranges of 16 m^{1}$<\beta <$16 m^{1} and 16 m^{1}$<\gamma <$16 m^{1} is adopted.
Table 1Parameters for fully watersaturated poroelastic soil medium [27]
Lame constant ${\lambda}_{s}/\times {10}^{7}Pa$  1.0 
Lame constant ${G}_{s}/\times {10}^{7}Pa$  3.0 
Solid density ${\rho}_{s}/$ (g/cm^{3})  2.5 
Compresswave velocity ${c}_{SS}/$ (m/s)  140 
Shearwave velocity ${c}_{SH}/$ (m/s)  121 
Porosity $n$  0.3 
Compressibility parameter $\alpha $  0.95 
Compressibility parameter $M/\times {10}^{9}Pa$  5.0 
Ratio between fluid viscosity and the permeability $b/\left({10}^{10}{m}^{1}s\right)$  1.0 
Water density ${\rho}_{w}/$ (g/cm^{3})  1.0 
To verify the accuracy of this work, Fig. 2 compares the vertical stresses and displacements obtained from our model with that given by Lu et al. [27]. Their work studied the threedimensional steadystate response of a poroelastic halfspace soil medium subjected to a moving point load. By selecting the parameters of pavement and subgrade (${E}_{p}={E}_{e}=$0, ${h}_{p}={h}_{e}=$0) and the same parameters of the soil medium (as shown in Table. 1) and load as those in Lu et al. [27], one can calculate the stresses and displacements. Then, the stresses and displacements are nondimensionalized as ${\sigma}_{zz}^{*}={\sigma}_{zz}{a}_{R}^{2}/T$ and ${u}_{i}^{*}=({u}_{i}{G}_{R}\cdot {a}_{R})/T$$(i=x,y,z)$, respectively, and the coordinates are nondimensionalized as ${x}^{*}=x/{a}_{R}$, ${y}^{*}=y/{a}_{R}$, ${z}^{*}=z/{a}_{R}$, where ${a}_{R}=$1 m, is the reference length; ${G}_{R}=$3.0×10^{9} Pa, is the reference shear modulus; ${c}_{SH}$ is the shearwave velocity which is defined as ${c}_{SH}=\sqrt{{G}_{s}/\rho}$, and $T$ is the load pressure. The dynamic responses of the observation with coordinates $y=z=$1.0 are calculated in Fig. 2. The agreement between the two solutions, as far as the vertical stress and displacement is concerned, is excellent.
Fig. 2Comparison between present work and Lu et al.’s [27] work
(a) Vertical stress
(b) Vertical displacement
3. Numerical results and discussion
In order to investigate the response characteristics of the pavementsubgradesoft ground system, Saga airport road [28] is selected and analyzed by the proposed method. Saga airport road was constructed from 1990 to 1992 with a total length of 9 km. At the site, there is about 20mthick highly compressible and highly sensitive Ariake clay. The typical cross section of Saga airport road is composed of three layers: asphalt surface layer, subgrade with decomposed granite and soft subsoil with Ariake clay. The physical and mechanical properties of the road structure layers are given in Table 2, and they correspond to those given in Chai [28], and the parameters of the tandemaxle loads are given in Table 3, corresponding to the heavy truck adopted by Chai [28]. The main physical quantities of interest in this investigation are the vertical displacements, stresses, and accelerations of the soft subsoil. The variations of these quantities with respect to the coordinates, the load velocity and frequency, the thickness and the modulus of the sugrade, and the permeability coefficient of the soft subsoil are studied.
In Fig. 3, the vertical stresses of the points with coordinates –$\text{15}\text{m}\le x\le \text{15}\text{}\text{m}$, $z=\text{1.1}\text{}\text{m}$ and $y=$0.0 m, 1.0 m, 2.0 m, respectively, are calculated for two kinds of load velocities, $c=0.1{c}_{ss}$ and $c=0.9{c}_{ss}$, where ${c}_{ss}$ is the shear wave velocity of the saturated poroelastic halfspace. It can be seen from Fig. 3(a) that, when $c=0.1{c}_{\text{ss}}$, two stress peaks appear at the points where the load is applied. The values of the peaks for rear wheel are larger than that for front wheel, and this is because the rearaxle weight of the heavy trucks is larger than frontaxle weight. In terms of different points $y=$0.0 m, 1.0 m, 2.0 m, the stress for $y=$1.0 m is the largest, and the stress for $y=$0.0 m is larger than that for $y=$2.0 m, which indicates that the maximum response occurs beneath the point where the load is applied. When $c=0.9{c}_{ss}$, as shown in Fig. 3(b), the magnitude of the stresses is much larger than that in Fig. 3(a), and significant ground vibration is observed, which is obviously due to the fact that velocity $\text{0.9}{c}_{ss}$ is very close to the Rayleigh wave velocity of the porous medium. The vertical stresses for different points $y=$0.0 m, 1.0 m, 2.0 m, are also presented in Fig. 3(b). The stress for $y=$1.0 m is the largest, which is the same of the lowloadvelocity condition. In order to observe the response more clearly, the vertical stresses for two kinds of load velocities, $c=0.1{c}_{ss}$ and $c=0.9{c}_{ss}$ are calculated in plane –$\text{1}\text{0}\text{}\text{m}\le x\le \text{10}\text{}\text{m}$, –$\text{1}\text{0}\text{}\text{m}\le y\le \text{10}\text{}\text{m}$, as shown in Fig. 4. It is clear in Fig. 4 that four peaks can be observed in the plane, which corresponds to four wheels of the vehicle. The magnitude of the stresses for the high velocity of $c=0.9{c}_{ss}$ is much larger than that for the low velocity of $c=0.1{c}_{ss}$, and the stresses apparently vibrate to a wide extent.
Table 2Physical and mechanical properties of each layer (Saga Airport Road)
Layers  Asphalt concrete  Subgrade  Soft subsoil 
Thickness / m  0.1  1.0  – 
Elastic modulus $E$/ MPa  1000  35  10 
Poisson’s ratio $\nu $  0.3  0.35  0.4 
Solid density ${\rho}_{s}$/ kg m^{3}  –  –  1700 
Water density ${\rho}_{f}$/ kg m^{3}  –  –  1000 
Soil density $\rho $/ kg m^{3}  2200  2000  1490 
Porosity $n$  –  –  0.3 
Compressibility parameter $\alpha $  –  –  0.95 
Compressibility parameter $M$/ MPa  –  –  5.0×10^{9} 
Fluid viscosity $\eta $  –  –  0.1 
Intrinsic permeability $k$/ (m^{3}s/kg)  –  –  1.0×10^{11} 
Table 3Parameters of tandem axle loads for heavy truck
Weight of front wheel ${q}_{0}$ (kN)  30 
Weight of rear wheel ${Q}_{0}$ (kN)  70 
Loaded length of front wheel $2{l}_{1}$ (m)  0.32 
Loaded width of front wheel $2{l}_{2}$ (m)  0.32 
Loaded length of rear wheel $2{L}_{1}$ (m)  0.44 
Loaded width of rear wheel $2{L}_{2}$ (m)  0.44 
Distance between front and rear wheels $2{d}_{a}$ (m)  6 
Distance between right and left wheels $2{d}_{w}$ (m)  2 
Fig. 3Distribution of the vertical stresses for different observation points
(a) $c=0.1{c}_{ss}$
(b) $c=0.9{c}_{ss}$
Fig. 4Distribution of the vertical stresses in (x,y) plane
(a) $c=0.1{c}_{ss}$
(b) $c=0.9{c}_{ss}$
Fig. 5 shows the vertical displacements for different load excitation frequency, where $f=$0 Hz, $f=$2 Hz, $f=4$ Hz and $f=$8 Hz are considered. For the sake of comparison, the displacements are presented for $z=$1.1 m. In Fig. 5(a), for $c=0.1{c}_{ss}$, the displacements responses decrease as load excitation frequency increases. Two peaks appear where the load is applied when $f=$0 Hz and $f=$2 Hz, but, as load excitation frequency increases to $4$ Hz or more, a negative peak is observed near the center of the tandem axle. As the load velocity is close to the Rayleigh wave speed, as shown in Fig. 5(b), the magnitude of the displacements is much larger than that in Fig. 5(a), and the displacements response decreases rapidly with increasing load frequency. The displacements fluctuate before the load reaches the observation point.
Fig. 5Distribution of the vertical displacements for different excitation frequency
(a) $c=0.1{c}_{ss}$
(b) $c=0.9{c}_{ss}$
The variations of the vertical accelerations with distance $y$ for the road system at $z=$1.1 m are investigated in Fig. 6, where $f=$0 Hz, $f=$4 Hz and $f=$8 Hz are considered. It is obvious that the vertical accelerations increase rapidly with increasing load frequency. For a low load velocity of $c=0.1{c}_{ss}$, the acceleration has two peaks at the point where the load is applied when $f=$0 Hz and $f=$4 Hz. As load excitation frequency increases to $f=$8 Hz, a negative peak stress is caused near the center of the tandem axle. In Fig. 6(b), the vertical accelerations for high load velocity $c=0.9{c}_{ss}$ is presented. The vertical accelerations of the soft subsoil fluctuate significantly in front of the load, and the magnitude of the vertical accelerations increases with increasing load frequency. A closer inspection of Fig. 6(b) reveals that the frequency of the acceleration fluctuations in front of the loads becomes larger as load frequency increases. Fig. 7 shows the vertical accelerations in plane –$\text{1}\text{0}\text{}\text{m}\le x\le \text{10}\text{}\text{m}$, –$\text{1}\text{0}\text{}\text{m}\le y\le \text{10}\text{}\text{m}$ for load excitation frequency $f=$4 Hz. When $c=0.1{c}_{ss}$, as shown in Fig. 7(a), the peak value of the vertical acceleration occurs where the load is applied, and no positive acceleration is caused. For the high loadvelocity condition, the magnitude of the vertical acceleration is significantly larger than that for $c=0.1{c}_{ss}$. The positive peak stresses occur near the center of the tandem axle, and the vertical acceleration fluctuates in the area in front of the load.
Fig. 8 presents the variations of the vertical accelerations with distance $x$ for different permeable coefficients when $f=$4 Hz, where ${k}_{d}={\text{10}}^{\text{6}}\text{m/s}$, ${k}_{d}={\text{10}}^{\text{}\text{8}}\text{m/s}$ and ${k}_{d}={\text{10}}^{\text{}\text{10}}\text{m/s}$ are considered. From Fig. 8(a), it can be seen that the acceleration has two peaks at the point where the load is applied. The permeability of the subsoil has little effect on soil acceleration. When the load velocity is close to the Rayleigh wave speed, as shown in Fig. 8(b), the accelerations apparently vibrate to a wide extent, the positive acceleration peak appears near the center of the tandem axle and the magnitude of the accelerations is much larger than that for the low loadvelocity condition. The permeability affects the acceleration response slightly when the velocity is close to the critical speed. It also can be seen from Fig. 8 that the effect of the permeability of the soft subsoil on the magnitude of the vertical accelerations at the halfspace surface is not pronounced.
Fig. 6Distribution of the vertical accelerations for different excitation frequency
(a) $c=0.1{c}_{ss}$
(b) $c=0.9{c}_{ss}$
Fig. 7Distribution of the vertical accelerations in (x,y) plane
(a) $c=0.1{c}_{ss}$
(b) $c=0.9{c}_{ss}$
Fig. 8Distribution of the vertical accelerations for different permeable coefficients
(a) $c=0.1{c}_{ss}$
(b) $c=0.9{c}_{ss}$
The effects of the properties of the road structure layers on its mechanical response are important in the design and maintenance of the road infrastructure. Therefore, further investigation of the dynamic response of road systems with different material properties is needed. In this section, the vertical stresses of different subgrade heights and Young’s modulus, which represent subgrade rigidity, are discussed. The observation point is located at the surface of the soft subsoil with coordinate –$\text{15}\text{}\text{m}\le x\le \text{15}\text{}\text{m}$ and $y=$1 m, as shown in Fig. 1. In this section, the load velocity is $\text{0.1}{c}_{ss}$ and frequency is 4 Hz.
In Fig. 9, the calculation parameters of the road structure layers are fixed and shown in Table 2, where the height and modulus of the subgrade varied from 1.0 to 3.0 m, and from 15 to 60 MPa, respectively. It can be seen from Fig. 9(a) that the magnitude of the vertical stresses increased sharply at the wheel action point and then gradually decreased to zero far from the wheel action point. As the subgrade height increases, the vertical stress clearly decreases, and the variations of the maximum vertical stress at different subgrade height are extremely important near the loading area but less so far from this zone. The effect of subgrade modulus on the vertical stress is presented in Fig. 9(b). The magnitude of the vertical stress decreases apparently when subgrade modulus increases from 15 to 30 MPa, but, as subgrade modulus reaches 45 MPa or more, the effect of subgrade modulus on the vertical stress is not pronounced.
Fig. 9Effect of the subgrade rigidity on the vertical stress (c=0.1css)
(a) Effect of the subgrade height on the stress
(b) Effect of the subgrade modulus on the stress
The influence of the subgrade height and Young’s modulus on the vertical stress and acceleration is further investigated, as shown in Fig. 10. In this figure the vertical stress ratio and the vertical acceleration ratio are presented as follows:
$\mathrm{}\frac{\text{Vertical\hspace{0.17em}stresses\hspace{0.17em}for\hspace{0.17em}different\hspace{0.17em}subgrade\hspace{0.17em}heights\hspace{0.17em}and\hspace{0.17em}modulus}}{\text{Vertical\hspace{0.17em}stresses\hspace{0.17em}for\hspace{0.17em}fixed\hspace{0.17em}subgrade\hspace{0.17em}height\hspace{0.17em}and\hspace{0.17em}}{\text{modulus}(h}_{e}=\text{1.0 m}{,E}_{e}=\text{10}\text{}\text{MPa)}},$
$\frac{\text{Vertical accelerations\hspace{0.17em}for\hspace{0.17em}different\hspace{0.17em}subgrade\hspace{0.17em}heights\hspace{0.17em}and\hspace{0.17em}modulus}}{\text{Vertical accelerations\hspace{0.17em}for\hspace{0.17em}fixed\hspace{0.17em}subgrade\hspace{0.17em}height\hspace{0.17em}and\hspace{0.17em}}{\mathrm{m}\mathrm{o}\mathrm{d}\mathrm{u}\mathrm{l}\mathrm{u}\mathrm{s}\mathrm{}(h}_{e}\text{= 1.0 m}{,E}_{e}\text{= 10 MPa)}}.$
Fig. 10Effect of the subgrade rigidity on the response components (c=0.1css)
(a) Effect of the subgrade height on the responses
(b) Effect of the subgrade modulus on the responses
Fig. 10(a) shows that for a fixed subgrade modulus of 10 MPa and for the subgrade height varying from 1.0 to 3.0 m, the dynamic responses are sensitive to the variation in subgrade height. As the subgrade height increases, the vertical stress ratio and acceleration ratio decrease at small subgrade heights and gradually become constant at larger subgrade heights. This result explains why this factor is so important in the design of subgrade. It suggests disastrous effects on the subgrade in the case of an underestimated height and predicts inefficiency in the case of an overestimated height. Finally, in Fig. 10(b), the height of the subgrade is fixed at 1.0 m and its Young’s modulus varies from 10 to 80 MPa. For low subgrade modulus, the vertical stress ratio decreases until the trend become constant as a larger modulus, and the vertical acceleration ratio trends become linear with the increasing subgrade modulus. This result implies that a reasonable subgrade modulus should be chosen in the design of road.
4. Conclusions
In this paper, the steadystate response of a pavementsubgradesoft ground system subjected to moving traffic load is studied using a semianalytical approach. The pavement and the subgrade are regarded as two elastic layers resting on a poroelastic halfspace soil medium. Four rectangular loads with a constant magnitude are used to represent the traffic load. Elastic dynamic theory for the singlephase elastic soil and Biot’s fully dynamic poroelastic theory for the saturated halfspace are used to account for the model of threelayer road structure. Using the double Fourier transform, the governing equations of motion in the frequency wavenumber domain are analytically solved. The fitted phase angle is employed to exclude the growing exponential terms in the dynamic stiffness matrix, which are critical for numerical stability. The expressions of the vertical stress, displacement, and acceleration in the time domain are evaluated by the inverse fast Fourier transform. The main conclusions of this study can be summarized as follows:
1. The responses in the condition of high velocity shows obvious difference from those of low velocity case, which proves that the influence of moving speed on the dynamic responses of the road system is significant.
2. The distributions of the stress, displacement, and acceleration of the soft subsoil are different for various values of loading frequency. For small values of frequency, the maximum responses occur at the point where the load is applied. As frequency increases, the positive response peaks appear and the maximum responses occur near the center of the tandem axle.
3. The acceleration has two peaks at the point where the load is applied, and the permeability of the subsoil has little effect on soil acceleration when the velocity is low. When the load velocity is close to the Rayleigh wave speed, the acceleration is much larger than that for the low loadvelocity condition and apparently vibrates to a wide extent. Generally speaking, the effect of the permeability of the soft subsoil on the magnitude of the vertical accelerations at the halfspace surface is not pronounced.
4. The dynamic responses of the road structure layers are sensitive to the variations in the subgrade height and Young’s modulus. The numerical calculations shows that these two parameters are extremely important in the design of subgrade and suggests potentially disastrous effects on the subgrade in the case of an underestimated height and suggests inefficiency in the case of an overestimated height.
References

Kausel E., Roesset J. M. Stiffness matrices for layered soil. Bulletin of the Seismological Society of America, Vol. 71, Issue 6, 1981, p. 17431761.

Kausel E., Peek R. Dynamic loads in the interior of a layered stratum: an explicit solution. Bulletin of the Seismological Society of America, Vol. 72, Issue 5, 1982, p. 14591481.

Chow Y. K., Swaddiwudhipong S., Lim S. A. Dynamic finite strip analysis of surface foundations. Earthquake Engineering and Structural Dynamics, Vol. 16, Issue 3, 1988, p. 457467.

De Barros F. C. P., Luco J. E. Response of a layered viscoelastic half space to a moving point load. Wave Motion, Vol. 19, 1994, p. 189210.

De Barros F. C. P., Luco J. E. Stresses and displacements in a layered half space for a moving line load. Applied Mathematics and Computation, Vol. 67, 1995, p. 103134.

Gunaratine M., Sanders O. Response of a layered elastic medium to a moving strip load. International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 20, 1996, p. 191208.

Auersh L. Wave propagation in layered soils: theoretical solution in wave number domain and experimental results of hammer and railway traffic excitation. Journal of Sound and Vibration, Vol. 173, Issue 2, 1994, p. 233264.

Picoux B., Rotinatb R., Regoina J. R., Houedec D. L. Prediction and measurements of vibrations from a railway track lying on a peaty ground. Journal of Sound and Vibration, Vol. 267, Issue 3, 2003, p. 575589.

Picoux B., Houedec D. L. Diagnosis and prediction of vibration from railway trains. Soil Dynamics and Earthquake Engineering, Vol. 25, 2005, p. 905921.

Jones D. V., Petyt M. Ground vibration in the vicinity of a rectangular load on a half space. Journal of Sound and Vibration, Vol. 166, Issue 1, 1993, p. 141159.

Jones D. V., Petyt M. Ground vibration in the vicinity of a strip load: an elastic layer on an elastic half space. Journal of Sound and Vibration, Vol. 161, Issue 1, 1993, p. 118.

Jones D. V., Petyt M. Ground vibration in the vicinity of a strip load: an elastic layer on a rigid foundation. Journal of Sound and Vibration, Vol. 152, Issue 3, 1992, p. 501515.

Skvireckas R., Kersys A., Kersys R., Lukosevicius V. Research of lateral vibrations of a passenger wagon running along the curved path. Journal of Vibroengineering, Vol. 14, Issue 2, 2012, p. 706714.

Lu Z., Yao H. L. Effects of the dynamic vehicleroad interaction on the pavement vibration due to road traffic. Journal of Vibroengineering, Vol. 15, Issue 3, 2013, p. 12911301.

Biot M. A. Theory of propagation of elastic waves in a fluidsaturated porous solid. 1. Lowfrequency range. Journal of the Acoustical Society of America, Vol. 28, 1956, p. 168178.

Biot M. A. Mechanics of deformation and acoustic propagation in porous media. Journal of Applied Physics, Vol. 33, 1962, p. 14821498.

Theodorakopoulos D. D. Dynamic analysis of a poroelastic halfspace soil medium under moving loads. Soil Dynamics and Earthquake Engineering, Vol. 23, 2003, p. 521533.

Theodorakopoulos D. D., Chassiakoos A. P. Dynamic effects of moving load on poroelastic soil medium by an approximate method. International Journal of Solids and Structures, Vol. 256, 2004, p. 18011822.

Jin B. Dynamic response of a poroelastic halfspace generated by high speed load. China Quarter of Mechanics, Vol. 25, 2004, p. 168174.

Zeng X., Rajapakse R. Vertical vibrations of a rigid disk embedded in a poroelastic medium. International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 23, 1999, p. 20752095.

Jin B. The vertical vibration of an elastic circular plate on a fluidsaturated poros half space. International Journal of Engineering Science, Vol. 37, 1999, p. 379393.

Senjuntichai T., Sapsathiarn Y. Forced vertical vibration of circular plate in multilayered poroelastic medium. Journal of Engineering Mechanics, Vol. 129, 2003, p. 13301341.

Cai Y., Cao Z., Sun H. Dynamic response of pavements on poroelastic halfspace soil medium to a moving traffic load. Computers and Geotechnics, Vol. 36, 2009, p. 5260.

LefeuveMesgouez G., Mesgouez A. Ground vibration due to a highspeed moving harmonic rectangular load on a poroviscoelastic halfspace. International Journal of Solids and Structures, Vol. 45, Issue 11, 2008, p. 3353–3374.

LefeuveMesgouez G., Mesgouez A. Threedimensional dynamic response of a porous multilayered ground under moving loads of various distributions. Advanced Engineering Software, Vol. 46, 2012, p. 7584.

Lu Z., Yao H. L., Cheng P. Ground vibration of soft subgrade subjected to a nonuniformly distributed train load. Rock and Soil Mechanics, Vol. 31, Issue 10, 2010, p. 32863294.

Lu J. F., Jeng D. S. A halfspace saturated poroelastic medium subjected to a moving point load. International Journal of Solids and Structures, Vol. 44, Issue 26, 2007, p. 573586.

Chai J. C., Miura N. Trafficloadinduced permanent deformation of road on soft subsoil. Journal of Geotechnical and Geoenvironmental Engineering, Vol. 128, Issue 11, 2002, p. 907916.
About this article
The authors gratefully acknowledge the financial support of the National 973 Project of China (No. 2013CB036405) and the Natural Science Foundation of China (No. 51209201 and 51279198).