Vibration Analysis of Thin Plate Structures Subjected to aMoving Force Using Frequency-Domain Spectral Element Method

A frequency-domain spectral element method (SEM) is proposed for the vibration analysis of thin plate structures subjected to a moving point force. )e thin plate structures may consist of multiple rectangular thin plates with arbitrary boundary conditions that form multispan thin plate structures, such as bridges. )e time-domain point force moving on a thin rectangular plate with arbitrary trajectory is transformed into a series of stationary point forces in the frequency domain. )e vibration responses induced by the moving point force are then obtained by superposing all vibration responses excited by each stationary point force. For the vibration response of a specific stationary point force, the plate subjected to the specific stationary point force is represented by four spectral finite plate elements, which were developed in the authors’ previous work. )e SEM-based vibration analysis technique is first presented for single-span thin plate structures and then extended to the multispan thin plate structures. )e high accuracy and computational efficiency of the proposed SEM-based vibration analysis technique are verified by comparison with other well-known solution methods, such as the exact theory, integral transformmethod, finite element method, and the commercial finite element analysis package ANSYS.


Introduction
In the structural engineering field, moving load problems, which deal with vibrations of structures under moving forces or masses, have been important research subjects over the last few decades. Typical examples of moving load problems are overhead cranes, ballistic systems, cableways, railway tracks, magnetic disks, building slabs, and bridges. ese structures under moving loads have been modeled either as one-dimensional (1D) or two-dimensional (2D) structures. For the reliable design and maintenance of such structures, it is essential to predict the vibration characteristics of the structures accurately, using an efficient analysis method.
To investigate the vibration characteristics of complex structures subjected to moving loads, the FEM has been the most widely used. e shape functions used in the FEM are independent of the vibration frequencies of a structure. us, the FEM generally requires a very fine discretization of the structure to acquire reliable solutions, especially in the high-frequency range. is may result in a significant increase in computational cost.
us, the frequency-domain SEM based on the fast Fourier transform (FFT) theory [26,27] has been suggested as an alternative to the FEM. e spectral element matrix (or dynamic stiffness matrix) used in the SEM is formulated from the frequency-dependent dynamic shape functions derived from free-wave solutions that satisfy the governing equations of motion in the frequency domain.
us, the SEM can provide highly accurate solutions by representing a uniform structure member as a single finite element, regardless of its size, and thus significantly reduce the computational cost compared to the FEM.
Despite the outstanding features of the SEM, it has been applied to very few moving load problems, for instance, in [21][22][23][24] for 1D structures and in [25] for 2D structures. For 1D structures, Azizi et al. [21] seem to be the first to have applied the SEM to continuous beams and bridges under a moving point force. ey discretized a uniform beam structure into many finite elements (more than two) and represented a moving point force using effective nodal forces and moments acting on each of the nodes of a finite beam element, as commonly done using the FEM. Later, Sarvestan et al. [22] adopted the SEM technique developed by Azizi et al. [21] for the vibration analysis of a cracked Bernoulli-Euler beam subjected to a moving point force.
Recently, for the vibration analysis of a Timoshenko beam subjected to a moving point force, Song et al. [23] presented a new SEM technique. ey represented the time-domain moving point force by a series of stationary point forces in the frequency domain. ereafter, by using the superposition principle, they obtained the time-domain vibration responses to the moving point force by summing and taking the inverse Fourier transform of all individual frequencydomain vibration responses to each stationary point force. For the vibration responses excited by a specific stationary point force, they used the so-called two-element method, where the Timoshenko beam was discretized into two finite elements so that the stationary point force was located at their joint. To avoid the structural discretization required for the two-element method used in Song et al. [23], Kim and Lee [28] introduced a spectral element modeling technique, called the modified one-element method, where exact solutions can be obtained by simply adding some correction terms, formulated in analytical forms, to the solutions by the one-element model. Kim and Lee [24] applied the modified one-element method to the vibration analysis of multispan beams subjected to a moving point force. For the 2D structures, Shirmohammadi et al. [25] seem to be the first to have applied the SEM to the vibration analysis of rectangular plates subjected to a moving point force. ey assumed that the point force was moving at a constant speed and the trajectory of the moving point force was parallel to the one edge of a rectangular plate. For the SEM analysis, they used the spectral finite strip element model developed for a Levytype plate, and, in the frequency domain, they transformed the moving point force into the modal form by using the separation of variable method. us, the application of their SEM technique is limited to the Levy-type plates. rough a thorough literature survey, we found that there have been no studies yet in which the SEM has been applied to moving load-induced vibrations of a thin plate structure with arbitrary boundary conditions.
For the SEM-based vibration analysis of a thin plate structure with arbitrary boundary conditions, the formulation of a generic type of spectral plate element model that can be applied to any boundary conditions has been required. Recently, Park et al. [29,30] developed such a generic type of spectral plate element model for the thin isotropic and composite rectangular plates with arbitrary boundary conditions by using the boundary splitting method [31] and the spectral super element method [32]. However, to the   authors' best knowledge, there are no reports in the literature  on the application of the generic type of spectral plate element model to the vibration analysis of a thin plate  structure with arbitrary boundary conditions, subjected to  a moving point force with arbitrary trajectory. us, we propose a new SEM-based vibration analysis method for isotropic single-span or multispan thin plate structures with arbitrary boundary conditions, subjected to a moving point force by using the generic type of spectral plate element model developed by Park et al. [29]. By extending the frequency-domain representation of a moving point force moving on a beam structure [23], a point force that moves along an arbitrary trajectory on a thin plate is transformed into the frequency domain as a series of stationary point forces distributed along the trajectory of this moving point force. e vibration responses induced by the moving point force are then obtained by superposing all individual vibration responses excited by each stationary point force. An individual vibration response excited by a specific stationary point force is obtained by representing the rectangular plate span on which the specific stationary point force is located using the four-element model, where the specific stationary point force is located at the junction of four spectral plate elements. All the other plate spans, without the stationary point force, are represented by the two-element model. e high accuracy of the solution and computational efficiency of the newly proposed SEM-based vibration analysis method are verified by comparison with other solution techniques, such as the exact theory, ITM, FEM, and the well-known commercial finite element analysis package ANSYS. In addition, the effects of the boundary conditions, moving speed, and the trajectory of a moving force on the vibration responses of the single-span square plates and three-span rectangular plates are numerically investigated.

Problem Statement
Consider a thin isotropic uniform rectangular plate (simply, plate) having arbitrary boundary conditions, subjected to a moving point force with a constant magnitude P, as shown in Figure 1. e length and width of the plate are denoted by l x and l y , respectively. e red dotted line represents the arbitrary trajectory of the moving point force, and the current positions of the moving point force with respect to the x-and y-coordinates are denoted by s x (t) and s y (t), respectively.
Using the classical plate theory [33], the governing equation of the plate subjected to a moving point force P can be written as where w(x, y, t) is the transverse displacement; D � Eh 3 /[12 (1 − ] 2 )] is the bending stiffness of the plate, where E is Young's modulus, ] is Poisson's ratio, and h is the plate thickness; ρ is the mass density; and f (x, y, t) Pδ (x − s x (t)) δ (y − s y (t)) represents a moving point force of magnitude P, where δ(x) and δ(y) are delta functions [34].  [23], the 2D moving point force f (x, y, t) can be transformed into the frequency domain, using the discrete Fourier transform theory (DFT) [35], as follows:

Vibration
where f n (x, y) Pδ x − s x(n) δ y − s y(n) , where i −1 √ is the imaginary unit. In (2), f m (x, y, ω m ) are spectral components of f (x, y, t), N is the number of spectral components up to the Nyquist frequency [35], and t n n∆t, where Δt T/N is the time increment and T is the time window.
Using the delta function, a continuous function g(x, y) can be approximated in terms of its sampling as follows [36]: By using (4), (2) can be written as follows: Because the delta functions δ(s x(j) − s x(n) ) and δ(s y(k) − s y(n) ) di er from zero when n j and n k, respectively, (5) can be rewritten as (6) where K ≤ N − 1. By omitting the subscripts m for brevity, (6) can be written in a simple form as follows: where As shown in (7), the frequency-domain representation of a moving point force P is given by a series of stationary point forces F k (k 0, 1, 2, . . ., K) de ned by (8). As illustrated in Figure 2, all stationary point forces are distributed along the trajectory of the moving point force. e total number of stationary point forces on the plate is denoted by K + 1. K is dependent on the sizes of T and T A , where T A is the traveling time required for a point force to move from one end to the other end of its trajectory [23]: Figure 2, the position of a stationary point force F k is denoted by (s x(k) , s y(k) ), where  Figure 3. e rectangular nite plate element has dimensions l ζ and l η with respect to the local coordinates ζ and η, respectively. In [29], the spectral element model for a rectangular nite plate element was developed based on the Kantorovich-based nite strip element method where the rectangular nite plate element was divided into N ζ nite strip elements in the direction of coordinate ζ and N η nite strip elements in the direction of coordinate η. e spectral element model developed in [29] can be cited in the following form: where S (e) (ω) is the 8(N ζ + N η )-by-8(N ζ + N η ) symmetric spectral element matrix (or dynamic sti ness matrix) and d (e) represents the spectral nodal degrees of freedom (simply, nodal DOFs) vector de ned by x y s y (t) f(x,y,t)=Pδ(xs x (t))δ(ys y (t)) s x (t) l x l y Figure 1: A thin rectangular plate subjected to a moving point force P, where the red arrow indicates current position of the moving point force.

Four-Element Model for a Rectangular Plate Subjected to a Stationary Point Force.
A rectangular plate subjected to a stationary point force (say, the k-th stationary point force, F k ) can be represented by four nite plate elements as shown in Figure 4. e stationary point force F k is placed at the junction of four nite elements and its location is denoted by (s x , s y ). By assembling all spectral element equations for the four nite plate elements shown in Figure 4, and by using the same assembly process commonly used in the standard FEM, we can obtain an assembled spectral element equation as follows: where S(ω) and Figure 2: Frequency-domain representation of a moving point force P as a series of stationary point forces F k (k 0, 1, 2, . . ., K).  Shock and Vibration In (14), S (e) (ω) and d (e) (e 1, 2, 3, 4) are de ned by (10) and (11), respectively; additionally, T (e) are the 8(N x(e) + N y(e) )-by -12[(N x(1) + N x (2) ) + (N y(1) + N y(3) ) − 1] transformation matrices required for the assembly of four nite plate elements, and they satisfy the following relations: In (16) and (17), 0 denotes the zero vectors and f F {F k 0 0 0} T .

Spectral Element Vibration Analysis.
By imposing all geometric and natural boundary conditions on (13), and by using the same method commonly used in the standard FEM, we can obtain a reduced spectral element equation in the following form: where the symbol ( ) denotes the reduced matrices and vectors after all geometric and natural boundary conditions are imposed. To compute the natural frequencies, we consider the free vibration problem by removing all external force terms from (19). e natural frequencies are then determined by searching for the frequencies ω that satisfy the condition of det[ S(ω) ] 0. As the matrix S(ω) is a transcendental function of ω, a proper root nding method [27] can be used to compute the natural frequencies.
By solving (19) for d, using the relation between d and d, and nally using the relations between d and d (e) (e 1, 2, 3, 4), the frequency-domain dynamic responses due to a stationary point force F k can be obtained as follows: (2) , s x ≤ x ≤ l x , 0 ≤ y ≤ s y , where where N(x, y, ω) is the dynamic shape function matrix that is available from [29]. e above solution process is repeated to obtain the frequency-domain dynamic responses for all stationary point forces, F k (k 0, 1, 2, . . ., K). Based on the superposition principle, the total frequency-domain dynamic responses due to all stationary point forces can be obtained from Finally, the time-domain dynamic responses due to a moving point force can be obtained from (22), by using the inverse fast Fourier transform (IFFT) algorithm [35], as follows:

Frequency-Domain Representation of a Moving Point
Force. Figure 5 shows a multispan rectangular plate (simply, multispan plate) with arbitrary supports and boundary conditions, where (X, Y) denotes the global coordinates for the entire multispan plate. e total length and uniform width of the multispan plate are denoted by L X and L Y , respectively, and the length of each span is denoted by l x [q] , where q 1, 2, 3, . . ., Q. e distance from the left edge of the rst span (X 0) to the right edge of the q-th span is denoted by X [q] . e multispan plate is subjected to a moving point force with a constant magnitude P, and the red dotted line in Figure 5 represents the arbitrary trajectory of the moving point force. e current position of the moving point force is denoted by S X (t) and S Y (t), with respect to the global coordinates (X, Y), as follows: With the same approach used in Section 3.1, the moving point force represented by (24) can be transformed into the frequency domain as follows: Figure 4: Four-element model for a rectangular plate subjected to a stationary point force F k . ◇ the node at which the stationary point force F k is located, • the nodes on the boundary edges, and ○ the interim nodes without external loadings.
Shock and Vibration 5 where F k (k 0, 1, 2, . . ., K) is also de ned by (8) and K ≤ N − 1. As illustrated in Figure 6, all stationary point forces are distributed along the trajectory of the moving point force. As shown in Figure 7, the position of the stationary point force F k with respect to the global coordinates Similarly, the position of the stationary point force F k with respect to the local coordinates (x, y) for the p-th span can be represented by (s

Spectral Element Modeling for a Multispan Plate Subjected to a Stationary Point Force.
To obtain the frequency-domain vibration responses due to a stationary point force, F k , which is located on the p-th span as shown in Figure 7, we should assemble one nite plate element on which the stationary point force F k is located (i.e., the p-th span) and Q − 1 nite plate elements that are free from the stationary point force. By using the four-element model introduced in Section 3.2.2, the p-th span subjected to a stationary point force F k is represented by the four nite plate elements as shown in Figure  e spectral element equation for the p-th span subjected to a stationary point force can be readily obtained from (13) in the following form: e spectral element equations for all other spans which are free from the stationary point force can be derived by assembling the two nite plate elements shown in Figure 8(b) in the following form: Finally, by assembling all the spectral element equations given by (28) and (30), a global spectral element equation for the complete multi-span plate can be obtained as follows:

Spectral Element Vibration Analysis.
By imposing all geometric and natural boundary conditions on (30), we can obtain a reduced global spectral element equation in the following form: To compute natural frequencies, we consider free vibration problem by removing all external forcing terms from (31). As explained in Section 3.3, the natural frequencies are then computed from det[ S g (ω) ] 0 by using a proper root nding method [27].
To obtain the frequency-domain dynamic responses, we rst solve (31) for d g and then use the relations between d g and d g and between d g and d [q] to obtain the nodal DOFs for the q-th span (q 1, 2, 3, . . ., Q). By using the relations between d [p] and d [p(e)] (e 1, 2, 3, 4), the frequency-domain dynamic responses at the p-th span subjected to a stationary point force F k can be obtained, in the local coordinates (x, y), as follows: [2] where N [p(e)] (x,y,ω) denotes the dynamic shape functions de ned by (21) for the element [p(e)]. Similarly, by using the relation between d [q] and d [q(e)] (e 1, 2), the frequency-domain dynamic responses in the q-th span, which is free from the stationary point force, can be also obtained as follows: By using (32) and (33), the frequency-domain vibration responses of the entire multispan plate due to a stationary point force F k can be represented with respect to the global coordinates (X, Y) as follows: x y L X Figure 7: A multispan rectangular plate subjected to a stationary point force F k acting on the p-th span.

Shock and Vibration
where where h(X) denotes the Heaviside step function defined by [34] h(X) � 0, (X < 0), e above solution process is repeated to obtain the frequency-domain vibration responses for all stationary point forces. By using the superposition principle, the total frequency-domain dynamic responses can be obtained as follows: Finally, the time-domain dynamic responses can be obtained from (37), by using the IFFT algorithm, as follows:

Numerical Results and Discussions
For the numerical study, three single-span square plates and two three-span plates subjected to a moving point force P � 1 N are considered as example problems. All example plates have null initial conditions. As shown in Figure 9, For Examples 1, 2, and 3, three cases of different trajectories of a moving point force are considered, which are represented with respect to the coordinates (x, y) as follows: Here, v 0 denotes a constant moving speed. Case 1 represents the case where the moving point force travels along the straight centerline of the plate, Case 2 represents the case where the moving point force crosses the catercornered plate, and Case 3 represents the case where the moving point force travels along a sinusoidal trajectory with respect to the centerline of the plate. For Examples 4 and 5, the moving point force is assumed to travel along the straight centerline of the multispan plates, which can be represented with respect to the global coordinates (X, Y) as follows:

Evaluation of the Proposed Method.
It is necessary to verify the proposed method. us, the accuracy and computational efficiency of the proposed SEM were evaluated by comparing the natural frequencies and vibration responses with those obtained by other solution techniques: the exact theory [37], ITM [5], FEM [14,38], and the well-known commercial finite element analysis package ANSYS [39].  (40).
For the results obtained using the FEM, a four-node 12-DOF (nonconforming) finite rectangular thin plate element model (known as ACM element) [38] was used. e vibration responses by the FEM were obtained using the technique introduced by Wu et al. [14] for moving force problems. Table 1 compares the lowest first ten natural frequencies of the S-S-S-S square plate (Example 1) obtained by the exact theory [37], FEM, and SEM. e exact natural frequencies (in Hz) of a thin S-S-S-S classical (or Kirchhoff) plate are given by [37] f (m,n) � π m 2 + n 2 2L 2 where the subscripts (m, n) denote the mode number. In the FEM, the natural frequencies are determined by using MATLAB function "eigs" [40]. For the proposed SEM, a single finite element is used for the entire single-span plate, and the root-finding method in [27] is applied to obtain the 8 Shock and Vibration x y L (c)

Shock and Vibration
natural frequencies. Table 1 shows that the SEM results are almost identical to the exact solutions given by (41). Table 1 also shows that the FEM results converge to both the exact solutions and the SEM results as the number of finite elements used in FEM is increased to more than 100 × 100. Table 2 compares the lowest first ten natural frequencies of the three-span rectangular plate with S-S-S-S supports (Example 4) obtained by the FEM, ANSYS, and SEM. e finite element model created in ANSYS is shown in Figure 11. e ANSYS finite element model is obtained by using the SHELL63 element that is a four-node finite rectangular plate element model based on the classical (Kirchhoff) plate theory [39]. For the ANSYS results, the MODAL [39] is chosen as the type of analysis. e MODAL in ANSYS performs a modal analysis. For the SEM results, only three finite elements are used for the entire three-span plate (i.e., a single finite element for each span). Table 2 shows that the FEM results converge to the SEM results as more finite elements are employed in the FEM analysis. Table 2 also shows that the ANSYS results obtained by using a sufficient number of finite elements (e.g., 400 × 100 finite elements) are found to be identical to the SEM results. Figure 12 compares the vibration responses for 0 ≤ t ≤ T A � 0.2 s at the center (x � y � 0.5L) of the S-S-S-S square plate (Example 1) obtained using the ITM [5], FEM [14], and SEM, when the trajectory of the moving point force is given by Case 1 of (39) with v 0 � 5 m/s. e vibration responses obtained by the ITM [5] are given by · sin(nπ/2)sin(mπx/L)sin(nπy/L) (42) Figure 12(a) shows that the vibration responses obtained by the ITM converge gradually to the SEM results when the total number of vibration modes used in the ITM-based analysis exceeds 100. Figure 12(b) also shows that the vibration responses obtained by the FEM converge gradually to the SEM results as the total number of finite elements used in the FEM-based analysis is increased to more than 40 × 40. Figure 13 compares the vibration responses for 0 ≤ t ≤ T A � 0.2 s at the center (x � y � 0.5L) of the S-S-S-S square plate (Example 1) obtained by the present SEM and FEM for two different trajectories of a moving point force P � 1 N: (a) Case 2 of (39), where v 0 � 5 m/s, and (b) Case 3 of (39), where v 0 � 5 m/s. Figure 13 clearly shows that, for both Case 2 and Case 3, the vibration responses obtained by the FEM converge gradually to the SEM results as the total number of finite elements used in the FEM-based analysis is increased to more than 40 × 40. Figure 14 compares the vibration responses for 0 ≤ t ≤ T A � 0.2 s at the center (x � y � 0.5L) of the S-S-S-S square plate (Example 1) obtained by the present SEM and ANSYS [39] for the same moving point force P � 1 N as considered in Figure 13. e moving point force travels along two different trajectories: (a) Case 2 of (39), where v 0 � 5 m/s, and (b) Case 3 of (39), where v 0 � 5 m/s. For the ANSYS results, the finite element model for the S-S-S-S square plate (Example 1) is obtained by using the SHELL63 element in ANSYS [39]. e technique proposed by Wu et al. [14] is then used to represent a moving point force in the ANSYS finite element model. For the ANSYS simulations, the TRANS [39] is chosen as the type of analysis. e TRANS in ANSYS performs a transient dynamic analysis. Similar to Figure 13, Figure 14 also shows that the ANSYS results converge gradually to the SEM results as the total number of finite elements used in the ANSYS analysis is increased to more than 40 × 40. Figure 15 compares the vibration responses obtained by the FEM and SEM for 0 ≤ t ≤ T A � 0.8 s at three locations, that is, at the center of the first span (X � 0.5L, Y � 0.5L), at the center of the second span (X � 2L, Y � 0.5L), and at the center of the last span (X � 3.5L, Y � 0.5L) of the S-S-S-S three-span rectangular plate (Example 4) subjected to a moving point force whose trajectory is given by (40) where v 0 � 5 m/s. Figure 15 also clearly shows that the FEM results approach the SEM results as the number of finite elements used in the FEM-based analysis is increased to more than 160 × 40.
e CPU times required to compute the SEM and FEM results shown in Figure 15 are measured with a standard workstation PC equipped with two sockets of Intel Xeon E5-2630v3 processors clocked at 2.4 GHz, and 384 GB of DDR4 RAM memory. For the SEM analysis, the span subjected to a stationary point force defined by (8) is represented by the four finite elements as shown in Figure 8(a), and each finite element is represented by 20 finite strip elements in both xand y-directions. However, the span that is free from the stationary point force is divided into two finite elements as shown in Figure 8(b) and each finite element is represented by 40 finite strip elements in the x-direction and 20 finite strip elements in the y-direction. e total number of stationary point forces considered for the S-S-S-S three-span rectangular plate (Example 4) is 64. For the FEM analysis, the example S-S-S-S three-span rectangular plate is represented by 40 × 10, 80 × 20, and 160 × 40 finite elements, and the FEM results are obtained using the MATLAB function "ode45" [40]. e CPU time for the SEM results obtained by using 4 × 2 finite elements is approximately 1.5 h, while the CPU times required for the FEM results obtained by using 40 × 10 finite elements, 80 × 20 finite elements, and 160 × 40 finite elements are 0.1 h, 5.5 h, and 305 h, respectively.
Based on the above observations, we can conclude that the proposed SEM provides highly accurate solutions and demands a relatively smaller computational cost when compared with the standard FEM.

Effects of the Boundary Conditions and Moving Speed on the Vibration Responses.
It has been well known that the vibration responses of a structure are dependent on its boundary conditions and external loading. Since the high accuracy of the proposed SEM was veri ed in the previous section, we applied it to square plates with three di erent boundary conditions (Examples 1-3 in Figure 9) and threespan rectangular plates with two di erent boundary conditions (Examples 4-5 in Figure 10) to investigate the e ects of boundary conditions on the vibration responses of example plates depending on the speed of a moving point force. ree types of moving force trajectories given by (39) are considered for the example square plates, while a single   Figure 12: Comparison of the vibration responses obtained by the integral transform method (ITM) [5], FEM [14], and SEM at the center of the S-S-S-S square plate (Example 1) subjected to a moving point force whose trajectory is given by Case 1 of (39), where v 0 5 m/s. (a) ITM [5] versus SEM. (b) FEM [14] versus SEM. moving force trajectory given by (40) is considered for the example three-span rectangular plates. Figures 16-18 show the e ects of boundary conditions and moving speed v 0 on the vibration responses at the center (x y 0.5L) of the S-S-S-S (Example 1), S-F-S-F (Example 2), and C-F-F-F (Example 3) square plates when a moving point force travels along the trajectories given by Case 1, Case 2, and Case 3 of (39).
From Figures 16-18, the vibration amplitudes are found to be quite dependent on the boundary conditions and the moving speed v 0 as follows: (1) e vibration amplitude of the C-F-F-F plate (Example 3) tends to decrease as the moving speed v 0 decreases, and (2) the S-S-S-S plate (Example 1) seems to have the smaller vibration amplitude than the S-F-S-F (Example 2) for all assumed trajectories of a moving point force.
When a moving point force travels on the S-S-S-S plate (Example 1) at a constant speed v 0 along the trajectories given by Case 1 of (39), the lowest critical speed can be obtained by equating the period of the (1,1) mode to the time needed to pass through a length double that of the plate length, as follows [5]: v cr 2f 11 L 9.713 m/s.
(43) Figure 19 shows the e ects of moving speed on the vibration responses at three locations of the example S-S-S-S square plate: (x 0.25L, y 0.5L), (x y 0.5L), and (x 0.75L, y 0.5L). From Figure 19, we can nd the following: (1) When v 0 < v cr , the vibration amplitudes tend to return to the original undeformed state; (2) when v 0 ≥ v cr , the vibration amplitudes tend to keep increasing to have a maximum value until the moving force reaches the other edge of the plate; (3) the maximum de ection of the plate when v 0 < v cr is slightly larger than that when v 0 ≥ v cr ; and (4) the maximum de ection of the plate tends to decrease when v 0 becomes larger than v cr . ese ndings appear quite similar to those obtained for the simply supported beam subjected to a moving point force in [5]. Shock and Vibration Figure 20 shows the e ects of boundary conditions on the vibration responses at the center of the rst span (X Y 0.5L), second span (X 2L, Y 0.5L), and third span (X 3.5L, Y 0.5L) of the S-S-S-S (Example 4) and C-S-S-S (Example 5) supports of three-span rectangular plates subjected to the moving point force whose trajectory is given by (40). e vibration amplitudes of the S-S-S-S supports (Example 4) are found to be larger than those of the C-S-S-S  supports (Example 5), at all the locations and moving speeds considered. In general, the vibration amplitudes at the center of the rst and third spans tend to increase as the moving speed v 0 increases. Above observations con rm that the boundary conditions and moving speed of a moving point force are essential factors that determine the vibration responses of a multispan plate subjected to a moving point force. Especially, for a plate with preset boundary conditions, quite di erent dynamic behaviors are observed depending on whether the moving speed is larger than the critical speed or not.

E ects of the Trajectory of a Moving Point Force on the Vibration Responses.
e present SEM is developed to consider any types of moving force trajectories. us, we considered square plates with three di erent boundary conditions (Examples 1-3 in Figure 9) to investigate the e ects of the trajectories of a moving point force on the deformed shapes of the example plates with time when the moving point force travels at the speed of v 0 5 m/s. ree types of trajectories are considered: Case 1, Case 2, and Case 3 de ned by (39).

24
Shock and Vibration Example 1 shows such symmetry. However, for the trajectory given by Case 3 of (39), the symmetry of the deformed shapes cannot be observed for the example problems considered in this study. In addition to the observations made in the previous section, when example plates are subjected to the same loading conditions such as the moving speed and trajectory, the e ects of boundary conditions on the deformed shapes can be also clearly observed by comparing Figures 21(a)

Conclusions
In this study, a frequency-domain vibration analysis method was proposed for thin rectangular plates subjected to a moving point force, and it was applied to several single-span and multispan rectangular plates. e moving point force was transformed into the frequency domain as a series of stationary point forces distributed along the trajectory of the moving point force. Based on the spectral plate element model developed in the previous study [29], the plate span subjected to a stationary  e results are summarized as follows: (1) e proposed frequency-domain representation of a moving point force can be applied to any type of moving point forces with arbitrary trajectories (2) Based on the spectral plate element model developed in [29], the proposed SEM can be applied to any thin rectangular plates with arbitrary boundary conditions (3) e high accuracy and computational efficiency of the proposed SEM are verified by comparison with the exact theory, ITM, FEM, and ANSYS (4) Numerical results show that the time histories and deformed shapes of the single-span and multispan plates are dependent on the boundary conditions and on the speed and trajectory of the moving point force (5) e proposed SEM can be readily applied to plates subjected to multiple moving point forces or moving distributed forces based on the superposition principle e application of the spectral element model developed in this study will be limited to the thin plate structures represented by single or multiple rectangular finite plate elements.
us, ongoing studies will be focused on the development of a spectral element model for the thick plate structures represented by triangular finite plate elements.

Conflicts of Interest
e authors declare that they have no conflicts of interest.