Dynamic response of a rigid pavement plate based on an inertial soil.
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, Dimitrovova  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 ([k.sub.O], [C.sub.O], and [m.sub.O], 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 [bar.w](x, y, z, t) at a given point inside the soil layer is equivalent to
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)
where x, y, z, and t are space-time coordinates of the soil studying point; [bar.w](x, y, z, t) is the deflection inside the soil layer defined as [bar.w](x, y, z, t) = [bar.w](x, y, 0, t) x [phi](z); and [phi](z) = sinh [[gamma](1 - z/[H.sub.s])]/sinh [gamma] is a vertical decay function of soil that must verify [phi](0) = 1 and [phi]([H.sub.s]) = 0; [k.sub.O], [c.sub.O], and [m.sub.O] 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 :
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (2)
where [H.sub.s] is the effective finite depth of the foundation that is dynamically activated; m is the density of the subgrade; y is a constant, named logarithmic decrement of the soil, which determinates the rate of decrease of the deflections depending on the depth; [E.sub.O] is Young's Modulus; v 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:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (3)
w(x, y, t) = [bar.w](x, y, 0, t) is the deflection of Kirchhoff plate which is equal to the deflection of the plate/soil interface.
p(x,y,t) = [p.sub.o][1 + (1/2) cos([omega]-t)][delta][x-x(t)][delta][y-y(t)] is the load transmitted to the pavement . Here x(t) = [v.sub.o]t + (1/2)acc([t.sup.2]); y(t) = (1/2)b are the geometrical position of load at the time t; [p.sub.o] is the magnitude of the moving wheel load; acc is the acceleration of the load; [omega] is the angular frequency of the applied load; [delta] is the Dirac function; a, b, and h are the dimensions of the finite plate and D 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 :
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4a)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4b)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4c)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4d)
where [ks.sub.x1], [ks.sub.x2], [ks.sub.y1], and [ks.sub.y2] are the elastic vertical translation stiffness and [V.sub.x], [V.sub.y] are the vertical shear forces of the plate.
(ii) The restriction of the elastic rotation is characterized by the following four equations :
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5a)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5b)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5c)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5d)
where [kr.sub.x1], [kr.sub.x2], [kr.sub.y1], and [kr.sub.y2] are the elastic rotational stiffness and [M.sub.x], [M.sub.y] are the bending moments of the finite plate.
The initial conditions (t = 0 s) are
[partial derivative]w (x, y, 0)/[partial derivative]t = (x, y, 0) = 0. (6)
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]
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (7)
[[omega].sub.mn] is the circular frequency of plate and [W.sub.mn] is the function of position coordinates determined for the mode numbers m and n in x- and y-directions. This form satisfies the initial conditions and the undamped free vibrations equation. Therefore natural modes satisfy the equation below :
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (8)
Equation (7) in (8) gives
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (9)
This equation is independent of time as the function [W.sub.mn]. Such a problem has solutions which may be in Navier's form [23, 24]:
[W.sub.mn] (x,y) = [A.sub.mn] sin(p[pi]/a x) sin (q[pi]/b y). (10)
Here, p and q are mode numbers of the plate. They are real numbers because of the boundary conditions of the problem [14, 22, 25] and m, n 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):
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (11)
2.2.2. Determination of Eigenmodes of the Plate. To obtain the mode numbers p and q 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
[W.sub.mn] (x, y) = [X.sub.mn] (x) sin (q[pi]/b y), (12)
where [X.sub.mn](x) is the eigenmode of the plate in the x-direction.
Substituting (12) into (9), we obtain an ordinary differential equation for [X.sub.mn](y):
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (13)
The solutions of the characteristic equation of (13) are
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)
For [beta] = [square root of ([2q.sup.2][a.sup.2] + [p.sup.2][b.sup.2] + ([c.sub.O][a.sup.2][b.sup.2]/[D[pi].sup.2]), X(x) becomes
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (15)
Equation (15) gives the general form of the eigenmode of the plate in the x-direction.
Boundary conditions along x-axis permit determining the [A.sub.i] coefficients:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (16)
where [a.sub.ij] coefficients are given by
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (17)
In order to obtain no trivial solution, it is necessary to propose that the determinant of (16) is zero, so
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (18)
(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 as
[W.sub.mn] (x, y) = [Y.sub.mn] (y) sin (p[pi]/a x), (19)
where [Y.sub.mn](y) is the eigenmode of the plate in the y-direction of the plate.
Substituting (19) into (9) conducts to an ordinary differential equation for [Y.sub.mn](y):
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (20)
The solutions of the characteristic equation of (20) are
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (21)
where [theta] = [square root of ([2p.sup.2][b.sup.2] + [q.sup.2][a.sup.2] + ([c.sub.O][a.sup.2][b.sup.2]/[D[pi].sup.2]))].
Equation (21) gives the general form of the eigenmode of the plate in the y-direction.
Boundary conditions along y-axis permit determining the [B.sub.i] coefficients:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (22)
Coefficients [b.sub.ij] were determined analogously to .
In order to obtain no trivial solution, it is necessary to propose that the determinant of (22) is zero, so
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (23)
2.2.3. Determination of Mode Numbers. To obtain the couples [p, q] 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 At and Bt coefficients after normalizing [A.sub.1] and [B.sub.1].
The natural mode of the plate is therefore given by
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (24)
2.2.4. Determination of the Time Function [T.sub.mn]. Suppose the solution of governing equation (3) is in the form like
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (25)
where [W.sub.mn] is a function of the spatial coordinates named modal function or natural mode of the plate and [T.sub.mn] a function of time.
Thus, for [W.sub.mn] satisfying (9), the [T.sub.mn](t) function verifies the following equation [22, 24]:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (26)
with [2[alpha][omega].sub.mn] = [gamma]h/(ph + [m.sub.O]) and a being the damping ratio of system.
The corresponding homogeneous solutions of (26) can be written:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (27)
According to the initial conditions defined in (6), [T.sub.0mn](t) = [a.sub.mn] = [b.sub.mn] = 0. A particular and the general solution of (26) are both given by 
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (28)
where [Q.sub.mn] is the normalizing function defined by
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (29)
Finally, the deflection solution of governing equation (3) is in the form of
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (30)
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 x 3.5 m, the thickness of 0.25 m, and physical characteristics of the plate like the density of [rho] = 2500kg x [m.sup.-3], Poisson's ratio [upsdsilon] = 0.25, and the longitudinal elastic modulus Ep = [24.10.sup.9] Pa. The density of the subgrade is taken equal to m = 1800 kg-[m.sup.-3], with Poisson's ratio [[upsilon].sub.s] = 0.35 and a longitudinal elastic modulus E = [50.10.sup.6] Pa. Finally the moving load magnitude is supposed to be [P.sub.O] = [80.10.sup.3] N and circular frequency [omega] = 100 rad x [s.sup.-1], acceleration acc = 2 m x [s.sup.-2], and a speed of [[upsilon].sub.O] = 25 m x [s.sup.-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 a = 10% . 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 timedisplacement 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 t = 3.5 x [10.sup.-2] 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 t = 0.003 s and this value is 66 percent more than those in stationary domain which appears at time t = 4[pi]/[omega] = 0.1256 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 [0; 0.87] m. Beyond these points, the value of maximum deflection is w = 0.000551479 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 x = 0 and y = 1.75 m, when the load is moving along the central axis of the plate (0 [less than or equal to] x [less than or equal to] 5 m and y = 1.75 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 ([H.sub.s]) 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 x-direction and the first five mode numbers were determined in the y-direction. We plotted the variation of the deflection curves according to different variables: x, y, and t.
Figure 4 shows the variations of the deflection at the center of the plate (x = 2.5 m; y = 1.75 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 (w = 0.000197106 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 ([H.sub.s] = 4.0 m for Vlasov soil type with w = 0.000192939 m and [H.sub.s] = 2.5 m for three-parameter soil type with w = 0.000157565 m). For lower values of [H.sub.s] up to [H.sub.s] = 1.25 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 [H.sub.s] = 10 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; [H.sub.s] = 2.5 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 (y = 1.75 m; 0 [less than or equal to] x [less than or equal to] 5 m), for different types of soil, at time t = 0.099603 s when the moving load arrived at the center of the plate (x = 2.5 m and y = 1.75 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 PasternakVlasov type have values between those obtained for the threeparameter soil and Pasternaksoil 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 [H.sub.s] 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: [H.sub.s] = 0.5 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 [H.sub.s] increases the deflection of the plate for low values of [H.sub.s]. 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 b/a. Table 1 presents the values of displacement under load versus soil parameters for different values of b/a (b/a = 0.5; b/a = 0.7; b/a = 1.0) at time t = 0.1256 s. This table is obtained for plates covering the same area. The table shows that deflection increases with the increase in b/a. 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 threeparameter 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 3.5 x [10.sup.-2] 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.
 M. S. Wei Tu, "Response modelling of pavement subjected to dynamic surface loading based on stress-based multi-layered plate theory," 2007
 LCPC, Conception et Dimensionnement des Chaussees Neuves: Guide Technique, LCPC, Ed., Ministere de l'Equipement 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.
 M.-H. Huang and D. P. Thambiratnam, "Deflection response of plate on Winkler foundation to moving accelerated loads," Engineering Structures, vol. 23, no. 9, pp. 1134-1141, 2001.
 S.-M. Kirn and J. M. Roesset, "Movingloads onaplate on elastic foundation," Journal of Engineering Mechanics, vol. 124, no. 9, pp. 1010-1017, 1998.
 P. Ullidtz, Modelling Flexible Pavement Response and Performance, Technical University of Denmark Polytekn, Kongens Lyngby, Denmark, 1998.
 L. Sun, "Steady-state dynamic response of a Kirchhoff's slab on viscoelastic Kelvin's foundation to moving harmonic loads," Journal of Applied Mechanics, vol. 74, no. 6, pp. 1212-1224, 2007
 D. St-Laurent, Synthese des Outils de Modelisation de Chausssees Actuellement Disponibles au LCPC, LCPC, Paris, France, 2008, http://www.lcpc.fr/.
 S. Lu, "Analytical dynamic displacement response of rigid pavements to moving concentrated and line loads," International Journal of Solids and Structures, vol. 43, no. 14-15, pp. 4370-4383, 2006.
 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.
 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 cellbased smoothed discrete shear gap (CS-FEM-DSG3) method," International Journal of Composite Materials, vol. 3, no. 6, pp. 19-27, 2013.
 C. Aron and E. Jonas, Structural element approaches for soil-structure interaction [M.S. thesis], Chalmers University of Technology, Goteborg, Sweden, 2012.
 S. Alisjahbana and W. Wangsadinata, "Dynamic analysis of rigid roadway pavement under moving traffic loads with variable velocity," Interaction and Multiscale Mechanics, vol. 5, no. 2, pp. 105-114, 2012.
 S. Lu, "Dynamic displacement response of beam-type structures to moving line loads," International Journal of Solids and Structures, vol. 38, no. 48-49, pp. 8869-8878, 2001.
 O. Civalek, "Nonlinear analysis of thin rectangular plates on Winkler-Pasternak elastic foundations by DSC-HDQ methods," Applied Mathematical Modelling, vol. 31, no. 3, pp. 606624, 2007.
 G. Pan and S. N. Atluri, "Dynamic response of finite sized elastic runways subjected to movingloads: a coupled BEM/FEM approach," International Journal for Numerical Methods in Engineering, vol. 38, no. 18, pp. 3143-3166, 1995.
 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. Dimitrovova, "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.
 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.
 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.
 C. Xiang-sheng, "Dynamic response of plates on elastic foundations due to the moving loads," Applied Mathematics and Mechanics, vol. 8, no. 4, pp. 355-365, 1987
 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.
Mohamed Gibigaye, Crespin Prudence Yabi, and I. Ezechiel Alloba
University of Abomey-Calavi, 01 BP 2009 Cotonou, Benin
Correspondence should be addressed to Mohamed Gibigaye; email@example.com
Received 7 November 2015; Accepted 15 December 2015
Academic Editor: Omer Civalek
Caption: Figure 1: Modeling of doweled rigid pavement under moving load .
Caption: Figure 2: Variation of the deflection directly under load versus time for different types of load ([H.sub.s] = 2.5 m).
Caption: Figure 3: Variation of the deflection at the fixed point of coordinates (x = 0; y = 1.75 m), versus time for a dynamically activated soil depth: [H.sub.s] = 2.5 m.
Caption: Figure 4: Variation of the deflection in the center of the plate (x = 2.5 m, y = 1.75 m) as a function of the dynamically activated depth of soil at time t = 0.099603 s when the load is located at the center of the plate.
Caption: Figure 5: Variation of the deflection along the central axis of the plate (0 [less than or equal to] x [less than or equal to] 5 m; y = 1.75 m), for different types of soil, at time t = 0.099603 s, when the load is located at the center of the plate.
Caption: Figure 6: Variation of the deflection as a function of plate magnitude, for different values of dynamically activated soil depth, at time t = 0.1256 s.
Table 1: Displacement versus soils parameters for different values of b/a (b/a = 0.5; b/a = 0.7; b/a = 1.0) at time t = 0.1256 s. b/a b/a = 0.5 b/a = 0.7 [H.sub.s] = 0.5 [k.sub.0] = 1.662 x 3.77 x 3.78 x [10.sup.8] [10.sup.-4] m [10.sup.-4] m [C.sub.0] = 2.596 x [10.sup.6] [m.sub.0] = 252.326 [H.sub.s] = 2.5 [k.sub.0] = 3.325 x 5.27 x 5.50 x [10.sup.7] [10.sup.-4] m [10.sup.-4] m [C.sub.0] = 1.298 x [10.sup.7] [m.sub.0] = 1261.63 [H.sub.s] = 5 [k.sub.0] = 1.662 x 5.36 x 5.8 x [10.sup.]7 [10.sup.-4] m [10.sup.-4] m [C.sub.0] = 2.596 x [10.sup.7] [m.sub.0] = 2523.26 [H.sub.s] = 7.5 [k.sub.0] = 1.108 x 5.48 x 5.73 x [10.sup.7] [10.sup.-4] m [10.sup.-4] m [C.sub.0] = 3.894 x [10.sup.7] [m.sub.0] = 3784.89 [H.sub.s] 10 [k.sub.0] = 8.311 x 5.13 x 5.77 x [10.sup.]6 [10.sup.-4] m [10.sup.-4] m [C.sub.0] = 5.192 x [10.sup.7] [m.sub.0] = 5046.52 b/a b/a = 1.0 [H.sub.s] = 0.5 [k.sub.0] = 1.662 x 5.73 x [10.sup.8] [10.sup.-4] m [C.sub.0] = 2.596 x [10.sup.6] [m.sub.0] = 252.326 [H.sub.s] = 2.5 [k.sub.0] = 3.325 x 7.08 x [10.sup.7] [10.sup.-4] m [C.sub.0] = 1.298 x [10.sup.7] [m.sub.0] = 1261.63 [H.sub.s] = 5 [k.sub.0] = 1.662 x 7.81 x [10.sup.]7 [10.sup.-4] m [C.sub.0] = 2.596 x [10.sup.7] [m.sub.0] = 2523.26 [H.sub.s] = 7.5 [k.sub.0] = 1.108 x 8.48 x [10.sup.7] [10.sup.-4] m [C.sub.0] = 3.894 x [10.sup.7] [m.sub.0] = 3784.89 [H.sub.s] 10 [k.sub.0] = 8.311 x 9.37 x [10.sup.]6 [10.sup.-4] m [C.sub.0] = 5.192 x [10.sup.7] [m.sub.0] = 5046.52
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Gibigaye, Mohamed; Yabi, Crespin Prudence; Alloba, I. Ezechiel|
|Publication:||International Scholarly Research Notices|
|Date:||Jan 1, 2016|
|Previous Article:||A new class of almost Ricci solitons and their physical interpretation.|
|Next Article:||The impact of duodenal diverticuli and the execution of endoscopic retrograde cholangiopancreaticography.|