Document Type : Research Paper
Authors
Faculty of Mechanical Engineering, K.N. Toosi University of Technology, Tehran, Iran
Abstract
Keywords

Mechanics of Advanced Composite Structures 2 (2015) 1730


Semnan University 
Mechanics of Advanced Composite Structures journal homepage: http://macs.journals.semnan.ac.ir 
The Effects of the Moving Load and the Attached MassSpringDamper System Interactions on the Dynamic Responses of the Composite Plates: An Analytical Approach
S. Fallahzadeh R., M. Shariyat^{*}
Faculty of Mechanical Engineering, K.N. Toosi University of Technology, Tehran, Iran
Paper INFO 

ABSTRACT 
Paper history: Received 30 October 2014 Received in revised form 21 Ocotber 2015 Accepted 27 October 2015 
In the current study, the effects of interactions of the moving loads and the attached massspringdamper systems of the composite plates on the resulting dynamic responses are investigated comprehensively, for the first time, using the classical plate theory. The solution of the coupled governing system of equations is accomplished through tracing the spatial variations using a Naviertype solution and the time variations by means of a Laplace transform. Therefore, the results are exact. The effects of various material, stiffness, and kinematic parameters of the system on the responses are investigated comprehensively and the results are illustrated graphically. Apart from the novelties presented in the modeling and solution stages, some practical conclusions have been drawn such as the fact that the amplitude of vibration increases for both the free and forced vibrations of the plate and the suspended mass, when the magnitude of suspended mass increases. 



Keywords: Composite plate Dynamic response Laplace transform Attached massspring system Moving load 


© 2015 Published by Semnan University Press. All rights reserved. 
The investigation of the forced vibration behavior of the structures carrying massspring systems is an important issue in many engineering applications, including the automotive, aerospace, marine, and civil structures. The control and measurement instruments attached to a control board and the navigation instrumentations, as well as the passengers and sprung masses (such as the differential) of a vehicle may be regarded as the attached springmass systems [1,2]. Moreover, it has been experimentally verified that a person or a crowd on a bridge or narrow platetype structure may be regarded as an attached massspring or massspringdamper system when the interaction between the human and the structure is of concern [3]. Sometimes, the attached massspring system may deliberately be included into the structure to serve as a dynamic vibration absorber [4,5].
The vibration of continuous systems with mass attachments has been the topic of some investigations recently. Turhan [6] analyzed the fundamental natural frequencies of an EulerBernoulli beam with an attached point mass. Kopmaz and Telli [7] investigated the free transverse vibration of a plate carrying a distributed mass through a simple model. Wong [8] analyzed the free bending vibration of a simply supported rectangular plate carrying distributed mass using the RayleighRitz method. Chiba and Sugimoto [9] studied the vibration characteristics of a cantilever plate with a springmass system attached upon it using RayleighRitz method. Li [10] represented an exact solution for free vibration analysis of rectangular plates with lineconcentrated mass and elastic linesupport. The vibratory characteristics of rectangular plates attached with continuously and uniformly distributed springmass which might represent free vibration of a humanstructure system were studied by Zhou and Tianjian [11]. Khodzhaev and Éshmatov [12,13] analyzed vibration and stability of viscoelastic plates with several attached masses using the numerical methods. Ciancio et al. [14] studied the free vibrations of a cantilever anisotropic plate carrying a concentrated mass approximately, using the energy method. Alibeigloo et al. [15] studied the free vibration of a simply supported laminated composite plate with a distributed patch mass using Hamilton’s principle and a double Fourier series solution. Watkins et al. [16] investigated vibration responses of an elastically pointsupported plate with attached masses, experimentally. Using a highorder plate theory, Malekzadeh et al. [17] investigated the free vibration characteristics of plates carrying distributed attached masses and the effects of the stiffness of the attached mass on the frequency of the system were included. In a more recent attempt, Amabili and Carra [18] used an experimental set up to investigate large amplitude vibrations of rectangular plates carrying concentrated masses.
The studies on the dynamic response of composite plates to moving loads and masses are of great importance from both the theoretical and practical points of view. The results of such investigations can be employed in various branches of transportation engineering, including designing road beds and bridges for trains, cars, trucks, etc., parking garages, aircraft runways, highspeed precision machining, magnetic disk drives, and so forth. Most of the previous studies on plates subjected to moving loads were accomplished based on the twodimensional theories such as the Classical Plate Theory (CPT) and the Firstorder Shear Deformation Theory (FSDT). Agrawal et al. [19] analyzed the dynamic responses of the orthotropic thin plates subjected to moving masses using Green’s function based on the CPT. Taheri and Ting [20] analyzed the moving load response of a thin plate using the finite element approach. Zaman et al. [21] presented a finite element algorithm to evaluate the dynamic response of a thick isotropic plate on viscoelastic foundation subjected to moving loads. Faria and Oguamanam [22] analyzed the vibration of Mindlin plates with moving concentrated loads, using the finite element method. Wu [23] investigated dynamic responses of a rectangular plate undergoing a transverse moving line load. Sun [24] studied the dynamic displacements of a thin plate induced by moving harmonic line and point loads. Malekzadeh et al. [25] accomplished threedimensional dynamic analysis of laminated composite plates subjected to moving loads using the analytical methods.
Several examples of plates with attached massspring systems and moving loads may be presented. Moving of the passengers/trolley across the airplane cabin that is carrying the passengers and seats (as mass and springs), moving of a lift truck on deck of a ship that is conveying goods and connexes, passage of the air flow (the induced aerodynamic forces) over the airplane wing with attached engine (that is usually attached trough flexible mounts), and passage of a high speed car on a bridge with some other stationary or low velocity commercial vehicles are some typical examples of a moving load on plates with attached massspring systems. In the present study, the dynamic responses of a rectangular composite plate with one or two attached massspringdamper systems to a moving load have been investigated analytically, for the first time. Since the method is fully analytical (including Laplace transform and Navier solutions) the results are exact. For the simple case of a symmetric crossply laminated plate, the present results are compared with the results of the available solutions and a good concordance is noticed. The effects of various parameters (i.e. magnitude of the suspended mass, coefficients of the spring and damper and velocity and magnitude of moving load) on the time history and the frequencies of vibrations of the plate and suspended mass have been investigated comprehensively and illustrated graphically.
According to the classical plate theory, the governing equations of motion of the composite plate are written as follows in the absence of the inplane forces [26].
(1) 
where , and
(2) 

(3) 
where
(4) 
where is measured for the midplane of the plate.
Assuming that: (i) the inplane displacement components u and v, have no direct effects on the suspended massspring system, i.e., the horizontal displacements of the mass are restrained and negligible, (ii) the deflections of the plate are small, and (iii) the equivalent spring exhibits a linear behavior, the mass –spring system can be modeled as a system of one degree of freedom whose motion can be described by the absolute vertical position of the suspended mass shown in Fig. 1. Therefore for the plate shown in Fig. 1, using Newton’s second law of motion and the Dirac delta function ( ), the distributed transverse load is expressed as follows:
(5) 
where can be described based on Fig. 1, as follows:
(6) 
On the other hand, the governing Equation of the suspended mass can be written as:
(7) 

Figure 1. A composite plate with an attached massspring system under moving load 
Substituting Eq. (6) with Eq. (5) and subsequently, substituting Eqs. (2)(5) with Eq. (1) lead to the governing equations of the plate in terms of the displacement components as the following:
(8) 

(9) 

(10) 
3.1. Composite Plate with Oneattached Massspring System
The simplysupported boundary conditions for the classical linear plate theory (CLPT) are as the following:
(11) 

(12) 
The boundary conditions appeared in Eq. (11) can be satisfied through the following Naviertype solutions
(13)

where and . can also be expressed in the following double Fourier series form:
(14) 
where
(15) 
Substituting from Eq. (5) with Eq. (15) leads to the following:
(16) 
Substituting Eq. (13) with Eqs. (2) and (3), and substituting the resulting equations with Eq. (12), reveal that the Navier solution (13) exists only if
It follows that plates with a single generally orthotropic layer, symmetrically laminated plates with multiple especially orthotropic layers, and antisymmetric crossply laminated plates, which satisfy Eq. (16), admit the Navier solutions for simplysupported boundary conditions. Substituting the Navier solution of Eq. (13) with governing Eqs. (8)(10) leads to the following equation in terms of unknown coefficients of the Navier solution and the vertical position of the suspended mass .
(17) 
where
(18) 
Equations (17) and (7) are the governing equations of the system shown in Fig. 1, which must be solved simultaneously, using the following initial conditions:
(19) 
Applying Laplace transform to Eqs. (7) and (17) leads to the following system of equations:
(20) 

(21) 
Substituting from Eq. (21) with Eq. (20) leads to the following equation:
(22) 
where
(23) 
Now suppose that
(24) 
Therefore,
(25) 
Rewriting the third equation of Eq. (25) considering Eq. (23) and using the Laplace transformed version of Eq. (13) lead to what follows:
(26) 
Replacing by in both sides of Eq. (26) results in the following equaition for
(27) 
where
(28) 

(29) 
Substituting from Eq. (27) with Eq. (25) and substituting the results into the Laplace transformed form of Eq. (13) lead to the Laplace transformed form of the displacement components ( , and . Eventually, using the inverse Laplace transform, Navier solution (13), i.e., , and may be obtained in the time domain. On the other hand, substituting from Eq. (27) with Eq. (21) results in the following equation for the vertical position of the suspended mass in the complex domain of Laplace transform:
(30) 
The inverse Laplace transform of Eq. (30) is determined by means of the residue theorem and by aid of Maple software. In order to successfully accomplish the process of applying the invers Laplace transform to any rational function by means of the residue theorem [27], the roots of the denominator must be found. In the case of Eq. (30) the most challenging part of denominator is , but the series of converges rapidly by growth of the number of its sentences, so the roots of can be found in any desired degree of accuracy and these roots can be used in implementing of the residue theorem.
3.2. A Composite Plate with Twoattached Massspring Systems
In the case of a composite plate with twoattached massspring systems, the total external mechanical load can be expressed as follows:
(31) 
where is the moving load defined by Eq. (6), and are magnitudes of the suspended masses, and are the vertical positions of the masses and , respectively and and are positions of the suspended masses, respectively. The governing equations of motion of the suspended masses are as the following:
(32) 

(33) 
According to Eqs. (15) and (31), the following expression may be used for a plate with two attached massspring systems:
(34) 
Thus, the Laplace transformed version of the relevant governing equation takes the following form.
(35) 

(36) 

(37) 
where
(38) 
Substituting and from Eqs. (36) and (37) with Eq. (38) results in:
(39) 
where
(40) 
Eq. (35) and the Laplace transformed version of Eq. (13) lead to the following:
(41) 

(42) 

(43) 
Replacing by and in both sides of Eq. (43) results in the following equations respectively:
(44) 

(45) 
where
(46) 

(46) 

(46) 

(46) 
So that
(47) 

(48) 
Substituting and from Eqs. (47) and (48) with Eqs. (41)(43) leads to the Laplace transformed version of the displacement components. Using the inverse Laplace transform, the displacement components may be obtained in the time domain. The most important conceptual consequence of this analysis is that the superposition principle is not allowed to be used in case of a plate with several massspringdamper systems, because the systems have mutual dynamical interactions on their dynamic responses as may be readily seen in Eqs. (47) and (48), but the presented method can be extended to the case of a composite plate with several massspringdamper attachments. However, we analyzed the case of a composite plate with only twoattached massspringdamper systems, for the sake of simplicity and brevity.
4.1. Verification of the Results
Since the present problem has not been investigated before, a special example of a plate with a moving load but without attached mass already treated by Malekzadeh et al. [25] is reexamined, for verification purposes. The considered plate is a symmetrically laminated [0/90/90/0] composite plate with the following specifications:
Thickness of each layer is 1 mm.
The present results are compared in Fig. 2 with results of the formulation proposed by Malekzadeh et al. [25], for the following parameters of the moving load:
Fig. 2 reveals that the results are coincident.
4.2. The Effect of M, C and K on the Forced Vibration
In Fig. 3 the effects of the magnitude of the suspended mass on the vibration of the plate treated in the previous section are depicted. As may be deduced from Fig. 3, the amplitude of the vibration for both the free (after ) and forced (before ) vibrations of the plate increases when the magnitude of the suspended mass is increased. This happens mainly because the exerted force by the attached system on the plate increases as the weight of the suspended mass increases. The amplitude of the free vibration increases with the suspended mass, almost linearly. The effects of the magnitude of suspended mass on the vibration frequencies of the middle point of the plate are not quite clear in this figure. In section 43, the effect of the suspended mass on the frequencies of the vibrations is explained in more details.
If a damper is included in parallel with the spring, the amplitude of vibration of the middle point of the plate decays as time elapses (see Fig. 4). As may be deduced from Fig. 4 when coefficient of the damper increases, the decay in the amplitude of vibration happens with a higher rate, as it has been already expected.
It is known that the vibration response of a system may be considered as a superposition of the effects of the various mode shapes (mode superposition). For a continuous system, infinite natural frequencies and consequently, infinite vibration modes may be defined. Fig. 4 shows that the effects of the higher vibration modes have been superimposed on the main response. Moreover, as Fig. 4 implies, the amplitudes of the higher vibration modes are considerably smaller, due to their smaller reversal times. Therefore, the effects of these higher modes are secondary ones. One of the superiorities of the present Laplace transform approach is that it retains all the frequency contents, as Fig. 4 confirms.
The effect of the stiffness of the spring is shown in Fig. 5. It is obvious that increasing the spring stiffness reduces the amplitude of vibration slightly. This fact stems from the evidence that the average stiffness of the system of the plate with attached massspring system increases when the stiffness of the spring increases. The effect of the spring stiffness on the vibration frequency is investigated in section 4.3 in more details.
4.3. The Effects of M and K on the Frequencies of the Vibration
In this section the effects of the parameters of the massspring system (M and K) on the frequencies of the vibrations of the system are investigated. For this purpose imaginary parts of roots of characteristic equation of the system are analyzed.
Figure 2. A comparison between the present results and the results of Malekzadeh et al. [25] for a laminated plate [0/90/90/0] with moving load but without the suspended mass ( )
Figure 3. The effects of the suspended mass on the time history of the transverse vibration of the middle point of the laminated [0/90/90/0] plate with attached massspring system ( )
Figure 4. The effects of the damper on time history of the transverse vibration of the middle point of the laminated [0/90/90/0] plate with attached massspringdamper system ( )
The characteristic equation of the system can be simply deduced from denominator of the right side of Eq. (27) as follows
(49) 
The imaginary components of roots of Eq. (49) are frequencies of vibration of the system. Table 1 reports the effects of the magnitude of the suspended mass on the frequencies of the transverse vibration of the plate of the previous section for a relatively small spring stiffness ( ) and in Table 2, depicts the effects for a larger stiffness ( ). Comparing the results of Table 1 and Table 2 shows that in the larger stiffness, the changes in the magnitude of the suspended mass results in relatively greater changes in the frequencies of the vibration of the plate. This is mainly because the stiffness of the spring is the coupling parameter between the plate and the massspring system and when its magnitude experiences an increase, the coupling parameter becomes stronger and the same changes in mass results in a relatively larger change in the frequencies of the vibration. On the other hand, Table 2 shows that attaching a massspring system to a plate results in a considerable increase in the first two frequencies of vibration and an increase in the magnitude of the attached mass, results in a relatively smaller decrease in the frequencies. This is because when we attach a massspring system to a plate, the total energy of system jumps to a higher level and its increase appears in higher frequencies of vibration and when magnitude of the suspended mass increases, the average inertia of the system increases and leads to a small reduction in the first two frequencies. On the other hand, comparing results of Tables 1 and 2 shows that in larger spring stiffness, the frequencies of vibration are higher.
4.4. The Effects of P_{0}, V_{x} and b_{0}
In Fig. 6 the effects of velocity and magnitude of the moving load on the time history of vibration of the composite plate and the suspended mass are illustrated. As it can be seen from Fig. 6, due to the presence of a damper parallel to the spring with damping coefficient of , the amplitude of vibration reduces for all magnitudes of the velocity and the moving load.
Both the velocity and magnitude of the moving load affect the maximum deflection of the middle point of the plate, as it is expected. When the velocity or magnitude of the moving load increases, the maximum deflection increases. The effect of the path line of the moving load on the time history of transverse deflection of middle point of the plate is shown in Fig. (7). It can be deduced from Fig. (7) that the maximum deflection is higher when the path line passes through the center of the plate, mainly because when the load path is closer to the edges of the plate, the supports compensate for the effects of the external load to some extent.
Figure 5. The effect of the spring stiffness on the time histories of the transverse vibration of the middle point of the [0/90/90/0] laminated plate and the suspended massspringdamper system (
Table 1. The effect of the suspended mass on the first four frequencies of the transverse vibration of the laminated [0/90/90/0] plate with attached massspring system ( )
Natural frequencies of the plate 
Frequencies of the massspring system 
Mass 

attached 
separated 

548.442 
456.143 
289.881 
137.116 
 
 
(Plate without attached massspring) 
548.444 
456.143 
290.005 
137.247 
7.061 
7.071 

548.444 
456.143 
290.005 
137.247 
4.465 
4.472 

548.444 
456.143 
290.005 
137.247 
3.157 
3.162 
Table 2. The effect of the suspended mass on the first four frequencies of the transverse vibration of the laminated [0/90/90/0] plate ( )
Natural frequencies of the plate 
Frequencies of the massspring system 
Mass 

attached 
separated 

548.442 
456.143 
289.881 
137.116 
 
 
(Plate without attached massspring) 

548.445 
456.144 
303.411 
151.019 
61.333 
70.710 

548.445 
456.144 
302.950 
149.349 
39.286 
44.722 

548.445 
456.144 
302.803 
148.843 
27.884 
31.623 

(a) 
(b) 
Figure 6. The effect of velocity and magnitude of the moving load on (a) the time history of transverse vibration of the middle point of the laminated [0/90/90/0] plate and (b) the suspended mass (
Figure 7. The Effect of the path of the moving load on the time history of transverse vibration of the middle point of the laminated [0/90/90/0] plate and the suspended massspringdamper system (
)
4.5. The Effects of the Positions of the Attached MassSpring Systems on the Frequencies of the Vibrations
The Effect of the positions of the attached massspring systems on the frequencies of vibrations of a composite plate with twoattached massspring system can be inferred from Table 3. As can be deduced from this table, the general conclusion of section 4.4 regarding distance of the moving load from the plate edge is valid for the frequencies as well, i.e., when the attached systems are closer to the center of plate, the frequencies are higher (in comparison with the first row of Tables 1 and 2 which shows the frequencies of a plate without any attachments).
4.6. The Vibration of the Suspended Mass
In this section, the vibration of the suspended mass, attached to the laminated [0/90/0/90/0/90/0/90] composite plate by means of a spring and a damper, is investigated. The thickness of the plate is ( for each layer) and the density of the material of the plate is which imply that the total mass of the plate is . Other properties related to the composite plate are calculated as follows
In Fig. 8, the effect of magnitude of the suspended mass on the time history of its vibration when attached to such a plate is plotted.
It may be deduced from Fig. 8 that an increase in the magnitude of the suspended mass alters its vibration pattern considerably, i.e., increases the amplitude and decreases the frequency of its vibration. The effect of damper’s coefficient on the time history of vibration of the middle point of the plate (where the massspringdamper system is attached) and the suspended mass is illustrated in Fig. 9.
In the present study, the dynamic response and vibration analysis of a composite plate with attached massspringdamper system under moving load is accomplished analytically for the first time. Employing the CPT description to explain the behavior of the composite laminated plate, exact results are extracted for the time histories, the frequencies of vibrations of the plate and the suspended mass, and a comprehensive analysis of the effects of various parameters is presented and illustrated graphically. The results are verified by the results available in literature for a special case and an excellent concordance is noticed.
Table 3. The effect of the position of the attached massspring systems on the first four frequencies of transverse vibration of the laminated [0/90/90/0] plate with twoattached massspring systems ( )
Natural frequencies of the plate 
Frequencies of the massspring system 
Locations on the attached system 

attached 
separated 

548.443 
456.143 
291.142 
140.970 
9.684 
10 

549.771 
456.925 
291.113 
138.403 
9.824 
10 

549.105 
456.538 
290.506 
140.342 
9.728 
10 

549.110 
457.323 
290.499 
139.054 
9.802 
10 

548.444 
457.725 
289.881 
139.709 
9.780 
10 
Figure 8. The effect of the suspended mass on the time history vibration of the suspended mass for the laminated [0/90/0/90/0/90/0/90] plate with attached massspringdamper system ( )
Figure 9. The effect of the damper coefficient on the time histories of vibrations of (a) the attachment point of the laminated [0/90/0/90/0/90/0/90] plate with attached massspringdamper system and (b) the suspended mass (
)
Apart from the novelties presented in the modeling and solution stages, some of the drawn practical conclusions can be summarized as:
References
[1] Wong JY. Theory of Ground Vehicles. 3^{rd} edition, John Wiley & Sons Inc.; 2001.
[2] Jazar RN. Vehicle Dynamics: Theory and Applications. Springer; 2008.
[3] Ellis BR, Ji T. HumanStructure Interaction in Vertical Vibrations. Proc Inst Civ Eng, Struct Build, 1997; 122: 1–9.
[4] Snowdon JC. Vibration of Cantilever Beams to Which Dynamic Absorbers Are Attached. J Acoust Soc Am 1966; 39: 878–886.
[5] Ranjan V, Ghosh MK. Forced Vibration Response of Thin Plate with Attached Discrete Dynamic Absorbers. ThinWall Struct 2005; 43: 1513–1533.
[6] Turhan O. On the Fundamental Frequency of Beams Carrying a Point Mass: Rayleigh Approximation Versus Exact Solutions. J Sound Vib 2000; 230(2): 449–459.
[7] Kompaz O, Telli S. Free Vibration of a Rectangular Plate Carrying Distributed Mass. J Sound Vib 2002; 251: 39–57.
[8] Wong WO. The Effect of Distributed Mass Loading on Plate Vibration Behavior. J Sound Vib 2002; 252: 577–583.
[9] Chiba M, Sugimoto T. Vibration Characteristics of a Cantilever Plate with Attached SpringMass System. J Sound Vib 2003; 260: 237–263.
[10] Li QS. An Exact Approach for Free Vibration Analysis of Rectangular Plates with Lineconcentrated Mass and Elastic LineSupport. Int J Mech Sci 2003; 45: 669–685.
[11] Zhou D, Ji T. Free Vibration of Rectangular Plates with Continuously Distributed SpringMass. Int J Solids Struct 2006; 43: 6502–6520.
[12] Éshmatov BK, Khodzhaev DA. Dynamic Stability of a Viscoelastic Plate with Concentrated Masses. Int Appl Mech 2008; 44(2): 208–216.
[13] Khodzhaev DA, Éshmatov BK. Nonlinear Vibrations of a Viscoelastic Plate with Concentrated Masses. J Appl Mech Tech Phys 2007; 48(6): 905–914.
[14] Ciancio PM, Rossit CA, Laura PAA. Approximate Study of the free Vibrations of a Cantilever Anisotropic Plate Carrying a Concentrated Mass. J Sound Vib 2997; 302: 621–628.
[15] Alibeigloo A, Shakeri M, Kari MR. Free Vibration Analysis of Antisymmetric Laminated Rectangular Plates with Distributed Patch Mass Using ThirdOrder Shear Deformation Theory. Ocean Eng 2007; 35: 183–190.
[16] Watkins RJ, Santillan S, Radice J, Barton Jr O. Vibration Response of an Elastically PointSupported Plate with Attached Masses. ThinWall Struct 2010; 48: 519–527.
[17] Malekzadeh K, Tafazoli S, Khalili SMR. Free Vibrations of Thick Rectangular Composite Plate with Uniformly Distributed Attached Mass Including Stiffness Effect. J Compos Mater 2010; 44: 2897–2918.
[18] Amabili M, Carra S. Experiments and Simulations for LargeAmplitude Vibrations of Rectangular Plates Carrying Concentrated Masses. J Sound Vib 2012; 331: 155–166.
[19] Agrawal OP, Stanisic MM, Saigal S. Dynamic Responses of Orthotropic Plates Under Moving Masses. Eng Arch 1988; 58: 9–14.
[20] Taheri MR, Ting EC. Dynamic Response of Plates to Moving Loads Finite Element Method. Comput Struct 1990; 34: 509–521.
[21] Zaman MM, Taheri R., Alavappillaix A. Dynamic Response of a Thick Plate on Viscoelastic Foundation to Moving Loads. Int J Numer Ana Meth Geomech 1991; 15: 627–647.
[22] de Faria AR, Oguamanam DCD. Finite Element Analysis of the Dynamic Response of Plates Under Traversing Loads Using Adaptive Meshes. ThinWall Struct 2004; 42: 1481–1493.
[23] Wu JJ. Dynamic analysis of a rectangular plate under a moving line load using scale beams and scaling laws. Comput Struct 2005; 83: 1646–1658.
[24] Sun L. Dynamics of plate generated by moving harmonic loads. J Appl Mech 2005; 72: 772–777.
[25] Malekzadeh P, Fiouz AR, Razi H. ThreeDimensional Dynamic Analysis of Laminated Composite Plates Subjected to Moving Load. Compos Struct 2009; 90: 105–114.
[26] Reddy JN. Mechanics of Laminated Composite Plate: Theory and Analysis. CRC Press; 1997.
[27] Spiegel MR. Advanced Mathematics for Engineers and Scientists. McGrawHill; 1971.