Dynamic Response of a Rigid Pavement Plate Based on an Inertial Soil
This work presents the dynamic response of a pavement plate resting on a soil whose inertia is taken into account in the design of pavements by rational methods. Thus, the pavement is modeled as a thin plate with finite dimensions, supported longitudinally by dowels and laterally by tie bars. The subgrade is modeled via Pasternak-Vlasov type (three-parameter type) foundation models and the moving traffic load is expressed as a concentrated dynamic load of harmonically varying magnitude, moving straight along the plate with a constant acceleration. The governing equation of the problem is solved using the modified Bolotin method for determining the natural frequencies and the wavenumbers of the system. The orthogonal properties of eigenfunctions are used to find the general solution of the problem. Considering the load over the center of the plate, the results showed that the deflections of the plate are maximum about the middle of the plate but are not null at its edges. It is therefore observed that the deflection decreased 18.33 percent when the inertia of the soil is taken into account. This result shows the possible economic gain when taking into account the inertia of soil in pavement dynamic design.
Pavements are an essential feature of the urban communication system and provide an efficient means of transportation of goods and services. Depending on its rigidity compared to the subsoil, pavements are classified as flexible, rigid, and semiflexible . To design these classes of pavements, the rational methods models are often used. The backbone of a method of roadway design is the mechanical model used to define the structure .
In the case of rigid pavements, the most used models are the multilayer elastic model of Burmister  and the Westergaard model, assuming pavement as a plate resting on the Winkler soil type [2, 4–6]. The differences between the pavement and the soil rigidities were conducted by Ullidtz , to deduce that Burmister model is generally not considered as an appropriate tool for the analysis of a rigid pavement response. Then, the large used design model of existing rigid pavement is Westergaard’s, using the Winkler soil type. Although Winkler model leads to relatively simplified results, it has serious limitations . Firstly, there is the deflection discontinuity between the charged and the uncharged part of pavement plate. Secondly, both models consider loads usually as a static one applied on a plate [2–4, 8]. According to Sun and Greenberg (2000) cited by St-Laurent , traffic loads on the pavement induce inertial effects that must be supported by foresaid pavement. Thus, the static load model does not reflect accurately the actual conditions of load on the pavement . During the last decade, many researchers have examined the problem assuming the loads as a dynamic one [1, 8, 10]. In the studies mentioned above, the soil used in the structure modeling are Winkler soil type. But, according to the design guide of United States National Cooperative Highway Research Program (NCHRP) , the two-parameter Pasternak model is designated in 1998 as the best pavement option to model the foundation of pavements. From that, many researchers base theirs works on using this type of soil [12, 13]. Alisjahbana and Wangsadinata  looked at the dynamic analysis of a rigid pavement under mobile load resting on Pasternak soil type. Based on this model, the determination of soil parameters is only based on the elasticity modulus and Poisson’s ratio. On the other hand, the Pasternak-Vlasov model that takes into account the logarithmic decrement of soil was used by Rahman and Anam  using the finite element method. The study of Rahman also showed that the Pasternak-Vlasov model is more economical than that of Winkler and cannot arbitrarily set the values of the intrinsic characteristics of the soil.
In most models used previously, the dynamic effect is taken into account only by the inertia of the plate [11, 15]. Concerning the soil, inertia is neglected in dynamic modeling of pavement structure. However, Civalek  took into account the inertia of the soil but had defined as constant independents values the intrinsic parameters of soil. He concluded that the effect of foundation inertia on the central deflection of the finite plate is not considerable. But according to the results of Pan and Atluri’s work , in engineering practice, this is still not always the case, and this factor may have significant effects on the dynamic response of the plate modeling the pavement. For these reasons, several studies have tried to modify the Pasternak-Vlasov soil introducing the inertia of the soil to a depth of soil susceptible to dynamic forces applied to the structure. Gibigaye in his work used an inertia soil model foundation for the study of the behavior of shells modeling underground shells . He concluded that it is important to take into account the inertia of the foundation soil on the dynamic response of civil engineering structures. Besides, Dimitrovová  noted that the classical formula which predicts a critical velocity of load significantly is overestimated compared to the one experienced in reality. She indicated that this formula should be revised by introducing two important notions: the effective finite depth of the foundation that is dynamically activated and the inertial effect of this activated foundation layer.
This work investigates dynamic response of a rigid pavement resting on an inertial soil. For that, the rigid pavement is modeled as a thin plate with dowels and tie bars in its edges.
In order to take into account its inertia, the soil is modeled as a three-parameter type (, , and , resp., integral characteristics in compression and in shearing and reduced mass of the foundation soil). The boundary conditions of these plates are modeled by two linear relationships between strains and stresses at the plate edges. The homogeneous solution of the problem is achieved by the method of separation of variables, so that the superposition gives a solution satisfying the boundary conditions. Since the deformation is expressed as eigenfunctions products, the solution of the dynamic problem is obtained based on the orthogonality properties of eigenfunctions. The general solution of the problem is obtained from the specific properties of the Dirac delta function. This paper provides an overview of the dynamic analysis of rigid pavements response as described above.
2. Materials and Methods
2.1. Governing Equation
In this research work, an isotropic homogeneous elastic rectangular plate resting on an elastic three-parameter soil is considered to model a pavement. The adjacent plates are supposed to be joined by dowels and tie bars. Based on the work of Asik  and Gibigaye work in the cylindrical axis system , the soil response according to the deflection at a given point inside the soil layer is equivalent towhere , , , and are space-time coordinates of the soil studying point; is the deflection inside the soil layer defined as ; and is a vertical decay function of soil that must verify and ; , , and are, respectively, integral characteristics in compression and in shearing and linear reduced mass of the foundation soil, supposed to be homogeneous and monolayer. They are expressed as follows :where is the effective finite depth of the foundation that is dynamically activated; is the density of the subgrade; is a constant, named logarithmic decrement of the soil, which determinates the rate of decrease of the deflections depending on the depth; is Young’s Modulus; is Poisson’s ratio.
According to the classic theory of thin plates and if taking into account the reduced mass of soil, the transverse deflection of the Kirchhoff plate satisfies the following partial differential equation: is the deflection of Kirchhoff plate which is equal to the deflection of the plate/soil interface.
is the load transmitted to the pavement . Here ; are the geometrical position of load at the time ; is the magnitude of the moving wheel load; is the acceleration of the load; is the angular frequency of the applied load; is the Dirac function; , , and are the dimensions of the finite plate and is the flexural stiffness of the plate.
The boundary conditions (Figure 1) are modeled as follows:
(i) The restriction of the elastic vertical translation is characterized by the four equations :where , , , and are the elastic vertical translation stiffness and , are the vertical shear forces of the plate.
(ii) The restriction of the elastic rotation is characterized by the following four equations :where , , , and are the elastic rotational stiffness and , are the bending moments of the finite plate.
The initial conditions ( s) are
2.2. Resolution of the Problem
2.2.1. Determination of the Eigenfrequencies
In order to solve governing equation (3) of the problem, it is assumed that the principal elastic axes of the plate are parallel to its edges. The free vibrations solution of the problem is set as [14, 20]
is the circular frequency of plate and is the function of position coordinates determined for the mode numbers and in - and -directions. This form satisfies the initial conditions and the undamped free vibrations equation. Therefore natural modes satisfy the equation below :
Here, and are mode numbers of the plate. They are real numbers because of the boundary conditions of the problem [14, 22, 25] and , are their respective roundness to the nearest integer number. The eigenfrequencies solutions of (9) are for the first auxiliary problem, those which satisfy boundary conditions (5a), (5b), (5c), and (5d):
2.2.2. Determination of Eigenmodes of the Plate
To obtain the mode numbers and and the eigenfunctions, modified Bolotin method is used. This consists of the resolution of the two auxiliary Levy’s problems. That can permit determining the eigenmode of the plate.
(i) First Auxiliary Levy Problem. The solution of (9) for the first auxiliary problem that satisfies the boundary conditions of (4a); (4b); (5a); and (5b) can be expressed as where is the eigenmode of the plate in the -direction.
The solutions of the characteristic equation of (13) are
For , becomes
Equation (15) gives the general form of the eigenmode of the plate in the -direction.
Boundary conditions along -axis permit determining the coefficients:where coefficients are given by
In order to obtain no trivial solution, it is necessary to propose that the determinant of (16) is zero, so
(ii) Second Auxiliary Levy Problem. The solution of (9) for the second auxiliary problem that satisfies the boundary conditions of (4c); (4d); (5c); and (5d) can be expressed aswhere is the eigenmode of the plate in the -direction of the plate.
The solutions of the characteristic equation of (20) arewhere .
Equation (21) gives the general form of the eigenmode of the plate in the -direction.
Boundary conditions along -axis permit determining the coefficients:Coefficients were determined analogously to .
In order to obtain no trivial solution, it is necessary to propose that the determinant of (22) is zero, so
2.2.3. Determination of Mode Numbers
To obtain the couples that permit having no trivial solutions, the transcendental equation system formed by (18) and (23) is solved. The solution cannot be determined analytically  so we used the Wolfram Mathematica software version 8.0.1 to get the numerical solutions. The triangulation of the systems of two equations permits determining the and coefficients after normalizing and .
The natural mode of the plate is therefore given by
2.2.4. Determination of the Time Function
Suppose the solution of governing equation (3) is in the form likewhere is a function of the spatial coordinates named modal function or natural mode of the plate and a function of time.
The corresponding homogeneous solutions of (26) can be written:According to the initial conditions defined in (6), . A particular and the general solution of (26) are both given by where is the normalizing function defined byFinally, the deflection solution of governing equation (3) is in the form of
3. Numerical Applications, Results, and Discussion
Using the procedure described above, a rigid roadway pavement subjected to a dynamic traffic load is analyzed. In this work, a finite rectangular plate doweled along its edges is considered as shown in Figure 1. The structural properties of the plate include the size of 5 m × 3.5 m, the thickness of 0.25 m, and physical characteristics of the plate like the density of m−3, Poisson’s ratio , and the longitudinal elastic modulus Pa. The density of the subgrade is taken equal to kgm−3, with Poisson’s ratio and a longitudinal elastic modulus Pa. Finally the moving load magnitude is supposed to be N and circular frequency rads−1, acceleration ms−2, and a speed of ms−1 . These parameters are typical material and structural properties of highway and airport pavements according to French Central Laboratory of Bridges and Pavements . It is also assumed that a damping ratio of the system equals % . For comparison purposes, three types of soil, (i) Pasternak soil type, (ii) Pasternak-Vlasov soil type, and (iii) Pasternak-Vlasov accounting the soil inertia type (three-parameter soil type), are considered.
3.1. Variation of Deflection as a Function of Time
Figure 2 shows, for the three-parameter soil, the variation of deflection under load as a function of time. This figure shows the time-displacement curve for two types of load: the uniform step load and the harmonically load. These loads evolve with the same constant acceleration and initial velocity. It is found that the deflection of the plate initially greatly increases with rapid oscillations and high amplitude, up to the moment s. This observation characterizes the transient domain for both types of load. After this phase, a stabilization of the oscillations is noted and the plate enters the stationary domain. In the stationary domain, the amplitude of deflection does not vary depending on the time. For harmonic load, the deflection varies harmonically with the time, since this load type is harmonic; but for step load, the deflection does not vary considerably in stationary domain. Yet, the deflection seemed affected by the boundary condition, since the deflection varies considerably for both types of load when the load approaches boundary of the plate. The maximum of deflection in transient domain appears at time s and this value is 66 percent more than those in stationary domain which appears at time s for harmonic load. This shows how it is important to take into account the transient domain response in the rigid pavement structures resistance analysis. In the specific case, the part of the plate concerned by this deflection is the segment m. Beyond these points, the value of maximum deflection is m. Besides, it is noticed that the response of uniform concentrated time step load is less than the harmonic load. That occurs in both transient domain and steady-state domain of plate response.
Figure 3 shows the variation of the deflection as a function of time at a fixed point with coordinates and m, when the load is moving along the central axis of the plate ( m and m). It is noted that the maximum deflection obtained neighborhood observation point along the traveling direction of the load and decreases in magnitude as well, as the load is removing. This shows that the maximum deflection occurs near the application point of the load.
3.2. Dynamically Activated Soil Depth () Effect on Dynamic Response
Based on the data listed above, the first five mode numbers of the plate modeling the pavement were determined in the -direction and the first five mode numbers were determined in the -direction. We plotted the variation of the deflection curves according to different variables: , , and .
Figure 4 shows the variations of the deflection at the center of the plate ( m; m) depending on the dynamically activated soil depth for different types of soil when the load is located at the center of the plate. It can be seen that, for the Pasternak soil, the value of the deflection of the plate is constant ( m) whatever the considered depth of the foundation. For Pasternak-Vlasov soil type and three-parameter soil type, the deflection increases to reach a maximum value at a given depth ( m for Vlasov soil type with m and m for three-parameter soil type with m). For lower values of up to m, the deflection of the two types of soil is the same. After this depth, the value of the deflection for three-parameter soil type is less than those of Vlasov soil type. The difference between the two responses increases with the depth of the dynamically activated soil up to 52 percent at m.
It is deduced from these observations that the dynamically activated soil depth greatly influences the response of the plate. We have chosen in this study the depth maximizing the three-parameter soil type; m.
3.3. Influence of the Inertia of the Soil on the Dynamic Response
Figure 5 shows the changes in the deflection along the central axis of the pavement plate ( m; m), for different types of soil, at time s when the moving load arrived at the center of the plate (m and m). It is noticed that the deflection is larger throughout the plate considering the Pasternak soil, compared to the deflection values for the three-parameter soil. The deflections of the plate when the soil is Pasternak-Vlasov type have values between those obtained for the three-parameter soil and Pasternak soil types. Taking the values of the deflections of the plate on the Pasternak-Vlasov soil type as a reference, it can be seen that the deflections at the center of the plate are reduced up to 18.33 percent. The more the depth is, the higher the gap is. In conclusion, inertial soil greatly reduces (up to 18.33%) the dynamic response of the pavement plate when the moving load is over the center of the plate.
3.4. Effect of the Variation of Load Magnitude on the Displacement for Different Values of Dynamically Activated Depth
Figure 6 shows the deflection of the plate versus the load magnitude for dynamically activated soil depths: m; 2.5 m; 5 m; 7.5 m; 10 m. It is noticed that the displacement increases linearly with the increase in load magnitude. Besides, the increasing of the dynamically activated depth parameter increases the deflection of the plate for low values of . Nevertheless, for great values of dynamically activated load, the deflection of the plate is not affected by the variation of the dynamically activated depth.
3.5. Variation of Displacement at Optimal Time versus Soil Parameters for Different Values of
Table 1 presents the values of displacement under load versus soil parameters for different values of (; ; ) at time s. This table is obtained for plates covering the same area. The table shows that deflection increases with the increase in . So the deflection is great for square plate compared to the rectangular plate. The size of pavement plate considerably affected the displacement of the pavement plate. It is more useful for engineer to design the rigid pavement plate as a rectangular one.
This paper dealt with some significant results from a study of the dynamic analysis of rigid pavements. The soil models used in this work are the well-known Pasternak model, the Pasternak-Vlasov model which takes into account the interaction between soil layers, and the improved three-parameter model considering the inertia of the soil. The main conclusions of this study are the following:(i)The soil inertia influences the pavement response at the middle of the plate when the load evolving along its centerline arrived at the center. This indicates a possible overdesign of pavements when using the two-parameter soil model.(ii)The effect of dynamically activated depth of Pasternak-Vlasov soil and three-parameter soil on the response is found to be significant for both soils types but more for three-parameter type than the Pasternak-Vlasov type.(iii)Before the stationary domain of oscillations, a transient response of the plate during seconds is noted. The transient response is bigger than the stationary one and it will be necessary to investigate deeply in further study.(iv)This study only covers plates pavement interconnected by dowels and tie bars. So, we could extend it to continuous pavements plates.(v)This study does not take into account the cyclic effect of the load. So, the fatigue response of the studied system could be further analyzed.(vi)The resonance is not studied in this work despite the effect that soil inertia can have in that, so we intend to study it in further work.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
S. O. Rahman and I. Anam, “Dynamic analysis of concrete pavement under moving loads,” Journal of Civil and Environmental Engineering, vol. 1, no. 1, pp. 1–6, 2005.View at: Google Scholar
M. S. Wei Tu, “Response modelling of pavement subjected to dynamic surface loading based on stress-based multi-layered plate theory,” 2007.View at: Google Scholar
LCPC, Conception et Dimensionnement des Chaussées Neuves: Guide Technique, LCPC, Ed., Ministère de l'Équipement des Transports et du Tourisme, Paris, France, 1994.
R. L. Baus and N. R. Stires, Mechanistic-Empirical Pavement Design Guide, University of South Carolina, Columbia, South Carolina, 2010.
P. Ullidtz, Modelling Flexible Pavement Response and Performance, Technical University of Denmark Polytekn, Kongens Lyngby, Denmark, 1998.
D. St-Laurent, Synthèse des Outils de Modélisation de Chausssées Actuellement Disponibles au LCPC, LCPC, Paris, France, 2008, http://www.lcpc.fr/.
ARA and ERES Division, “Appendix QQ: structural response models for rigid pavements,” in Guide for Mechanistic-Empirical Design of New and Rehabilitated Pavement Structures, pp. 1–91, NCHRP, Champaign, Ill, USA, 2003.View at: Google Scholar
T. Nguyen-Thoi, H. Luong-Van, P. Phung-Van, T. Rabczuk, and D. Tran-Trung, “Dynamic responses of composite plates on the pasternak foundation subjected to a moving mass by a cell-based smoothed discrete shear gap (CS-FEM-DSG3) method,” International Journal of Composite Materials, vol. 3, no. 6, pp. 19–27, 2013.View at: Publisher Site | Google Scholar
C. Aron and E. Jonas, Structural element approaches for soil-structure interaction [M.S. thesis], Chalmers University of Technology, Göteborg, Sweden, 2012.
M. Gibigaye, Problems of Dynamic Response of Shells Contacting with an Inertial Medium, at a Transient Domain of Vibrations, Industrial Institut of Zaporozhye, Zaporozhye, Ukraine, 1992.
Z. Dimitrovová, “Enhanced formula for a critical velocity of a uniformly moving load,” in Proceedings of the 6th European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS '12), pp. 1–10, Vienna, Austria, September 2012.View at: Google Scholar
M. S. Asik, “Vertical vibration analysis of rigid footings on a soil layer with a rigid base,” 1993, https://repositories.tdl.org/ttu-ir/bitstream/handle/2346/20107/31295007715294.pdf?sequence.View at: Google Scholar
A. Turhan, A Consistent Vlasov Model for Analysis of Plates on Elastic Foundations Using the Finite Element Method, Graduate Faculty, Texas Tech University, Lubbock, Tex, USA, 1992.
V. I. Samul, The Basis of Elasticity and Plasticity Theories: Pedagogic Book for High Schools Students, Moscow High School, Moscow, Idaho, USA, 2nd edition, 1982.
R. Harberman, Elementary Applied Partial Differential Equations with Fourier Serie and Boundary Value Problem, Prentice Hall, Englewood Cliffs, NJ, USA, 2nd edition, 1987.
S. Alisjahbana, “Dynamic response of clamped orthotropic plates to dynamic moving loads,” in Proceedings of the 13th World Conference on Earthquake Engineering, Vancouver, Canada, November 2004.View at: Google Scholar