Document Type : Research Article
Authors
^{1} Department of Aerospace Eng.
^{2} department of SGC
^{3} S.G. center of research
Abstract
Keywords

Mechanics of Advanced Composite Structures 5 (2018) 173 – 185


Semnan University 
Mechanics of Advanced Composite Structures journal homepage: http://MACS.journals.semnan.ac.ir 
An Analytical Approach to Thermoelastic Bending of Simply Supported Advanced Ribbed Composite Plates
M. Shahravi^{a}, S. Fallahzade^{b}, M. Mokhtari^{a,c}*,
^{a} Department of Aerospace Engineering, Maleke Ashtar University of Technology, Tehran, Iran
^{b} Department of Mechanical Engineering, Tehran, Iran
^{c} Telecommunication department, Iranian space research center, Space research center, Tehran, Iran
Paper INFO 

ABSTRACT 
Paper history: Received 20170830 Received in revised form 20180112 Accepted 20180129 
In the present paper, an analytical approach is used to study the thermal deflections of a simply supported composite plate with a beamlike stiffener. The results for a plate–beam system exposed to a sinusoidal thermal load is used to study the effects of the low Earth orbit (LEO) thermal conditions on the composite plates, which have been used in the structure of satellites and spacecraft. To solve the governing equations of the system, the Laplace transform method for the time domain is used with the Navier series expansions. As the employed method is completely analytical, the results are exact.




Keywords: Ribbed composite plate Thermoelastic bending Laplace transform Advanced composite 


© 2018 Published by Semnan University Press. All rights reserved. 
Fibrous composites are widely used in many engineering structures that are subjected to severe thermal environments. The composites have many attractive properties, including temperature resistance and a low thermal coefficient of expansion. In fibrous composite laminates, the coefficients of thermal expansion in the direction of the fibers are usually much smaller than those in the transverse direction [1]. Therefore, several plate theories have been developed by various researchers to predict the correct bending behavior of composite laminates under mechanical/thermal loads. These plate theories have been reviewed in various studies [26]. Reddy [7] presented a thermal analysis of laminated composite plates using Kirchhoff’s classical plate theory (CPT), Mindlin’s firstorder shear deformation theory (FSDT), and Reddy’s higherorder shear deformation theory (HSDT) [8]. Tauchert [9] investigated the stationary twodimensional temperature, stress, and displacement distributions for a simply supported slab consisting of bonded orthotropic layers. Khdeir and Reddy [10] developed refined plate theories to study the thermal stresses and deformations of crossply rectangular laminates using the statestress approach. Savoia and Reddy [11] solved the transient heat conduction equation for a given temperature distribution across the thickness of laminates for a threedimensional (3D) stress analysis of a square laminate subjected to a sudden uniform temperature change.
Rohwer et al. [12] obtained the thermal stresses in laminated composite plates using the higherorder theories. Carrera and Ciuffreda [13] obtained a closedform solution for the thermal analysis of laminated composite plates using a unified formulation. A refined FSDT was presented and used for the thermal analysis of laminated structures without using a shear correction factor [14]. Zenkour [15] used the parabolic and trigonometric shear deformation theories for the thermomechanical analysis of crossply laminated composite plates. Zenkour et al. [16] extended the thermomechanical analysis of laminated composite plates resting on an elastic foundation. Zhen et al. [17] used the global–local higherorder theory to carry out a thermal analysis of a fourlayered symmetric crossply laminated composite plate subjected to the actual temperature field. Gao and Zhao [18] developed the refined plate theory, which was based on the thermoelasticity theory, for the thermoelastic bending analysis of rectangular plates. Khdeir [19] obtained an exact solution for the thermoelastic bending analysis of crossply laminated arches with arbitrary boundary conditions. Kant and Shiyekar [20] developed an HSDT with a Taylor seriestype expansion in the thickness direction of the displacements to analyze composite and sandwich plates under thermal loading. Noda et al. [21] provided various solution methods to solve 3D heat conduction problems; the temperature profiles for 1D, 2D, and 3D linear elastic bodies, such as beams, plates, and shells, were obtained by solving the respective Fourier heat conduction equation. However, a review of the literature showed that such solution techniques are rarely employed for the thermoelastic analysis of multilayered plates.
In recent years, the problems associated with thermomechanical deflections, bending, buckling, and vibrations of composite and functionally graded (FG) sandwich plates have been analyzed and solved by some authors using fourvariable trigonometric shear deformation theories. Tounsi et al. [22] used a refined trigonometric shear deformation theory to analyze the thermoelastic bending of FG sandwich plates. Zidi et al. [23] used a fourvariable refined plate theory for the bending analysis of functionally graded material (FGM) plates under hygrothermomechanical loading. Beldjelili et al. [24] used a fourvariable trigonometric plate theory to analyze the hygrothermomechanical bending behavior of sigmoid (S)FGM plates resting on various elastic foundations. Bouderba et al. [25, 26] developed a simple shear deformation theory to analyze the thermal stability of FG sandwich plates. Chikh et al. [27] performed a thermal buckling analysis of crossply laminated plates using a simplified HSDT.
Structures consisting of composite plates stiffened by a set of beams form a class of structural elements that have practical importance in various engineering applications, such as aircrafts and ships. Because aerospace and marine vehicles are subjected to thermal and dynamic loads, confident predictions of the natural frequencies and amplitudes of the vibrations of the structural components are essential for preventing excessive vibration levels, which may result in fatigue failure [2830].
The addition of stiffeners to composite plates complicates the dynamic analysis and, thus, simplifying assumptions have to be made to facilitate a solution to the problem. Many analytical and numerical methods have been proposed to study the vibrations of ribstiffened plates. The proposed approaches include the orthotropic [31] and grillage [32] models, the Lagrange multiplier formalism [33], the Rayleigh–Ritz method [34, 35], the finite difference method [36], the finite element method [37, 38], the differential quadrature method [39], and the meshless method [40]. Hygrothermomechanical bendin analysis of SFGM plates on variable elastic foundation studied by Beldjelili et al [41].
According to the literature, the thermomechanical bending analysis of ribbed composite plates has never been performed. Thus, the present study is the first attempt in using an analytical approach to solve the problem of thermal bending of a composite plate with a beamlike stiffener. To study the results, an applied case that occurs frequently in satellite structure design problems and has been never solved was considered. Assumptions of the CPT were considered to describe the motion of the composite plate [7], while the motion of the stiffener was assumed to follow the Euler–Bernoulli beam equation. To determine the relationship between the governing equations of the plate and stiffener, a compatibility equation was used [29]. In addition, the quasistatic equation of heat conduction was used to describe the temperature field of the plate. The plate and the beam were simply supported along the mechanical boundaries to allow a Fouriertype expansion of the Levy solution to be applied to trace the deflection of the composite plate in terms of the spatial coordinates. The thermal boundary condition on one side of the plate was considered to be in a sinusoidal form, which is common in LEO satellite applications. The temperature equation with the assumption of a small strain rate was uncoupled from the equation of motion of the plate by omitting the coupling term of displacement from the temperature equation. The quasistatic equation of heat conduction was solved separately using the Laplace transform method. Furthermore, the Laplace transform method was utilized to solve the system of governing equations of motion of the plate. Finally, the method of residues (residue theorem [41]) was applied to transform the algebraic results of the complex domain (sdomain) into the time domain. The solution method that is used in the present study is fully analytical and the results are exact within the framework of the CPT. Thus, this paper comprehensively investigates and graphically represents the effects of various parameters, such as the dimension of the crosssection area of the stiffener and the orbital temperature parameters, on the time history bending of a composite plate and stiffener.
2.1. Equations of Motion of a Composite Plate with a Stiffener
According to the CPT, the governing equations of motion of the composite plate may be written as follows in the absence of inplane forces [7]:
(1) 
where and
(2) 

(3) 
where
(4) 

(5) 

(6) 
where and are measured from the midplane of the plate. The stiffener can be modeled as an Euler–Bernoulli beam parallel to the yaxis as shown in Fig. 1 assuming that the inplane displacements and have no direct effect on the motion of the stiffener; the deflections of the plate are small; and the torsional interaction moment between the beam and the plate is negligible.
Therefore, the governing equation of the beamlike stiffener can be written as follows [34]:
(7) 
where is the deflection of the stiffener and is the interaction force between the beam and the plate. Considering Equation (7) and neglecting the effect of the torsional interaction moment, the distributed transverse load can be expressed as [22]:
(8) 
Substituting Equations (2), (3), and (8) into Equation (1) leads to the governing equations of the plate in terms of the displacement components:
(9) 
Figure 1 Schematic of a composite with parallel stiffeners
(10) 

(11) 
Note that Equations (7) and (9)–(11) are the governing equations of motion of a composite plate with a beamlike stiffener. These equations must be solved in conjunction with the continuity condition at the interface between the plate and the beam, which can be expressed as follows [22]:
(12) 
2.2 Heat Conduction and Temperature Distribution
In the present study, the quasistatic equation of heat conduction will be considered only for the plate. It has been assumed that the presence of the beam does not affect the heat conduction process. Thus, the quasistatic equation of heat conduction for a laminated composite plate can be written as follows [35]:
(13) 
Assuming that the thermal properties of the plate are constant in any direction, Equation (13) can be expressed as:
(14) 
where is the coefficient of thermal diffusivity and . The boundary and initial conditions of Equation (14) are as follows:
(15) 
To solve Equation (14) with the prescribed boundary conditions of Equation (15), the following definition of the finite Fourier transform may be used in the direction [41]:
(16) 
The inverse transform may be obtained from
(17) 
Employing the integration by parts rule, the following relation is derived based on the definition presented in Equation (16):
(18) 
Applying the finite Fourier transform in the direction of by making use of Equation (18) on Eqs. (14) and (15) leads to the following:
(19) 

(20) 
To satisfy the boundary condition, the double sinusoidal expansion of will be used, which can be presented as follows:
(21) 
Substituting from Equation (21) into Equation (19) results in the following:
(22) 
Applying the Laplace transform to Equation (22) gives the Laplace transform version of as follows:
(23) 
where and . Making use of the inverse Fourier transform of Equation (17) with Equation (21) results in the following:
(24) 
Therefore,
(25) 
where
(26) 
In the following sections, Equations (4) and (5) in the form of which has been presented in Equation (25), will be used to calculate the thermal force vector and the thermal moment vector .
The simply supported boundary conditions for the classical linear plate theory are as follows:
(27) 

(28) 
The boundary conditions in Equation (27) may be satisfied using the following Naviertype solutions:
(29) 
where and . may also be expressed in the following double Fourier series:
(30) 
where
(31) 
Substituting from Equation (8) into Equation (31) leads to the following:
(32) 
where
(33) 
After the displacements of Equation (29) are first inserted into Equations (2) and (3), and then inserted into Equation (28), the result reveals that the Navier solution of Equation (29) exists only if .
It follows that the Navier solution for simply supported boundary conditions applies to plates with the following characteristics: a single generally orthotropic layer; symmetrically laminated plates with multiple orthotropic layers; and antisymmetric crossply laminated plates, which include the former cases as special cases. Substituting the Navier solution of Equation (29) into the governing Equations (9–11) leads to the following equation in terms of the unknown coefficients of the Navier solution :
(34) 

(35) 

(36) 
where .
Expanding the thermal force vectors and the thermal moment vectors in terms of the double Fourier sries leads to:
(38) 

(39) 
where
(40) 

(41) 
Note that has been defined by Equation (26). The temperaturerelated terms in the governing Equations (9)–(11) can then be written as follows:
(42) 

(43) 

(44) 
The particular form of the solution of the problem, which is based on the double Fourier series, requires that and be defined to achieve an analytical solution. Thus, the configuration of the composite plate must be as follows:
(45) 

(46) 
The conditions in Equations (45) and (46) are automatically satisfied for the following conditions: singlelayer plates with a generally orthotropic layer; symmetrically laminated plates with multiple orthotropic layers; and antisymmetric crossply laminated plates. In order to include and , the temperature distribution should be expanded in a double cosine series. Then , , ,and must be equal to zero [7]. Substituting , , and from Equations (42)–(44) into Equations (34)–(36) results in the following governing equations of motion of the plate:
(47) 
in which
(48) 
Equations (7) and (48) are the governing equations of the system shown in Fig. (1), and must be solved simultaneously. The lateral deflection of the beam may be presented as follows:
(49) 
Substituting from Equation (49) into Equation (7) results in the following equation for the motion of the stiffener:
(50) 
where
(51) 
Equations (47) and (50) are the coupled equtations governing the motion of the system of a composite plate with a beamlike stiffener, which must be solved in conjunction with the continuity condition that has been expressed in Equation (7). In order to solve Equations (47) and (50), the Laplace transform may be employed. Assuming that all of the initial conditions are zero, after the Laplace transform is applied, the form of the governing equations will be as follows:
(52) 

(53) 
where
(54) 
Substituting from Equation (53) into Equation (52) leads to
(55) 
The Laplace transform version of Equation (12) is
(56) 
Based on the Fourier series expansions, the following relations are valid:
(57) 

(58) 
From Equations (55) and (56) we have
(59) 
where
(60) 
The third row of Equation (59) can be written as
(61) 
where
(62) 

(63) 
Substituting Equation (61) into Equation (58) results in the following:
(64) 
Replacing with on both sides of Equation (64) leads to the following expression for :
(65) 
where
(66) 

(67) 
Substituting from Equation (65) into Equation (59), and then substituting the results into the Laplace transform version of Equation (29) leads to the Laplace transform version of the displacement components , and . By making use of the inverse Laplace transform, the Navier solution (29) for , and may be obtained in the time domain. The inverse Laplace transform of , and is determined using Maple software and the residue theorem. In order to successfully accomplish the process of applying the inverse Laplace transform to any rational function using the residue theorem [37], the roots of the denominator must be found. In the case of Equation (30), the most challenging part of the denominator is . However, the series of converges rapidly by the growth in the number of terms. Thus, the roots of can be found for any desired degree of accuracy and can be used in the implementation of the residue theorem.
In this the section, different aspects of the current problem and the effects of various parameters, including the periodic thermal load and the stiffener, on the results will be discussed in detail using some numerical examples. In this regard, consider a symmetrically laminated [0/90/90/0] composite plate with the following specifications:
The thickness of each layer is 1 mm and the dimensions of this plate are the mass of the plate would be . Thus, according to Equation 6, the nonzero coefficients of this plate are:
The plate has a stiffener with a crosssection area of (Fig. 1), which is attached to the plate in It has the following mechanical and thermal properties:
where is the thermal diffusivity coefficient of the composite plate with the assumption of homogeneity and is the thermal expansion coefficient of the composite plate.
Furthermore, suppose that a sinusoidal temperature has been applied to the upper surface of the plate; in this case, the sinusoidal temperature is applied to the stiffened composite plates of the outer surface of the structure of a LEO satellite. In the equation, is the amplitude of the temperature and is the frequency of the temperature function. In the case of an LEO satellite, it takes about 90 minutes (an orbital period) to complete a revolution in its orbit around Earth. The effects of and on the temperature of the midpoint of the stiffened plate are shown in Fig. 2, which result from Equations (24) and (25). Tracing the temperature variation through the thickness suggests that using a linear approximation of temperature for many applied purposes would result in an acceptable precision. However, in this study, the exact solution of the thermal conduction equation has been used to achieve higher levels of accuracy.
Figure 2. Effects of the orbital temperature parameters of and on the temperature at the midpoint of the composite plate ( ) with a stiffener ( )
Fig. 3 shows the effect of the height of the stiffener, which is an important parameter, on the time history deflection of the midpoint of the plate. As can be deduced from Fig. 3, using a stiffener with a relatively small height will not diminish the amplitude of the deflection of the plate to a considerable extent. Therefore, to decrease the amplitude considerably, the height of the stiffener (beam) must be comparable with the plate’s thickness (i.e., more than 2 times larger). Fig. 4 shows the effect of the stiffener’s height on the maximum deflection of the plate. This illustration also shows that a stiffener with a relatively small height does not significantly affect the maximum amplitude of the plate and that the sensitivity of the amplitude to the height of the stiffener increases as the height of the stiffener increases. However, the case for the width ( ) of the crosssection area of the plate is quite different. For example, Fig. 4 shows the effect of the stiffener’s width on the maximum deflection of the plate. As can be seen in the figure, the decrease in the maximum deflection of the plate with respect to has almost a linear pattern. This is because of the moment of inertia of the crosssection area of the beam with respect to its neutral axis ( ), which is a linear function of and a cubic function of This means that for a constant crosssection area, the beam should have a high and a low to reduce the maximum deflection of the beam. Of course we have to keep in mind that a beam with very low may cause high degrees of stress concentration, which must be avoided for design purposes (see Fig. 5).
Figure 3. Effect of the height of the stiffener ( ) on the timehistory deflection of the midpoint of the plate for
Figure 4. Effect of of the stiffener on the maximum amplitude of the deflection of the midpoint of the plate for b_{1}=0.01 m
In Fig. 6, the effect of and on the deflection of the midpoint of the composite plate ( ) with is shown. Note that all the figures have been produced under the assumption that the inner space of the satellite is at a constant temperature. However, this is an optimistic assumption for the electronic devices that are inside the satellite and a pessimistic assumption for the plate deflection because it results in greater deflections. In reallife situations, such an assumption cannot be fulfilled because the inner temperature of the satellite will probably change with the environment as the satellite moves in front of the sun in its orbit. Thus, the actual deflections will be less than the deflections presented in this paper.
In the present paper, an analysis of the thermal deflections of a simply supported ribbed composite plate with application in LEO satellite structures is accomplished for the first time. By employing the CPT description to explain the behavior of the composite laminated plate, exact results are extracted for the time histories and an analysis of the effects of various parameters is presented and illustrated graphically.
In addition to the novelties presented in the modeling and solution stages, some of the practical conclusions that have been drawn may be summarized as follows:
Figure 5. Effect of of the stiffener on the maximum amplitude of deflection of the midpoint of the plate for
Figure 6. Effects of the orbital temperature parameters of and on the deflection of the midpoint of the composite plate ( ) with a stiffener ( ).
References
[1] Sayyad AS, Ghugal YM, Mhaske BA. A FourVariable Plate Theory for Thermoelastic Bending Analysis of Laminated Composite Plates, Journal of Thermal Stresses 2015; 38: 904–925.
[2] Tauchert TR, Thermally Induced Flexure, Buckling and Vibration of Plates. Appl. Mech. Rev. 1991; 44(8):347–360.
[3] Argyris J, Tenek L. Recent Advances in Computational Thermostructural Analysis of Composite Plates and Shells with Strong Nonlinearities. Appl. Mech. Rev. 1997; 50(5): 285–306.
[4] Kant T, Swaminathan K. Estimation of Transverse/Interlaminar Stresses in Laminated Composites—A Selective Review and Survey of Current Developments. Compos. Struct. 2000; 49: 65–75.
[5] Reddy JN, Arciniega RA. Shear Deformation Plate and Shell Theories: From Stavsky to Present. Mech. Adv. Mater. Struct. 2004; 11(6): 535–582.
[6] Wanji C, Zhen W. A Selective Review on Recent Development of DisplacementBased Laminated Plate Theories. Recent Pat. Mech. Eng. 2008; 1:29–44.
[7] Reddy JN, Mechanics of Laminated Composite Plates, CRC Press, Boca Raton, 1997.
[8] Reddy JN, A Simple Higher Order Theory for Laminated Composite Plates, ASME J. Appl. Mech. 1984; 51: 745–752.
[9] Tauchert TR. Thermoelastic Analysis of Laminated Orthotropic Slabs. J. Therm. Stresses 1980; 3: 117–132.
[10] Khdeir AA, Reddy JN, Thermal Stresses and Deflections of Cross Ply Laminated Plates Using Refined Plate Theories. J. Therm. Stresses 1991; 14: 419–439.
[11] Savoia M, Reddy JN. ThreeDimensional Thermal Analysis of Laminated Composite Plates, Int. J. Solids Struct. 1995; 32: 593–608.
[12] Rohwer K, Rolfes R, Sparr H. HigherOrder Theories for Thermal Stresses in Layered Plates, Int. J. Solids Struct. 2001; 38: 3673–3687.
[13] Carrera E, A Ciuffred. ClosedForm Solutions to Assess MultilayeredPlate Theories for Various Thermal Stress Problems. J. Therm. Stresses 2004; 27: 1001–1031.
[14] Fares ME, Zenkour AM. Mixed Variational Formula for the Thermal Bending of Laminated Plates. J. Therm. Stresses. 1999; 22: 347–365.
[15] Zenkour AM. Analytical Solution for Bending of CrossPly Laminated Plates under ThermoMechanical Loading. Compos. Struct. 2004; 65: 367–379.
[16] Zenkour AM, Allam MNM, Radwan AF. Bending of CrossPly Laminated Plates Resting on Elastic Foundations under ThermoMechanical Loading. Int. J. Mech. Mater. Des. 2013;9: 239–251.
[17] Zhen W, Cheng YK, Lo SH, Chen W. Thermal Stress Analysis For Laminated Plates Using Actual Temperature Field. Int. J. Mech. Sci. 2007; 49: 1276–1288.
[18] Gao Y, Zhao B. The Refined Theory of Thermoelastic Rectangular Plates. J. Therm. Stresses 2007; 30: 505–520.
[19] Khdeir AA. An Exact Solution for the Thermoelastic Deformations of CrossPly Laminated Arches with Arbitrary Boundary Conditions. J. Therm. Stresses. 2011; 34: 1227–1240.
[20] Kant T, Shiyekar SM. An Assessment of a Higher Order Theory for Composite Laminates Subjected to Thermal Gradient. Compos. Struct. 2013; 96: 698–707.
[21] Noda N, Hetnarski RB, Tanigawa Y, Thermal Stresses, 2^{nd} Edition, Taylor & Francis: New York; 2003.
[22] Tounsi A, Houari MSA, Benyoucef S. A refined trigonometric shear deformation theory for thermoelastic bending of functionally graded sandwich plates. Aerospace Sci. Tech. 2013; 24: 209220.
[23] Zidi M, Tounsi A, Houari MSA, Bég OA. Bending analysis of FGM plates under hygrothermomechanical loading using a four variable refined plate theory. Aerospace Sci. Tech. 2014; 34: 2434.
[24] Beldjelili Y, Tounsi A, Mahmoud SR. Hygrothermomechanical bending of SFGM plates resting on variable elastic foundations using a fourvariable trigonometric plate theory. Smart Struct. Syst., Int. J. 2016; 18(4): 755786.
[25] Bouderba B, Houari MS, Tounsi A, Mahmoud SR. Thermal stability of functionally graded sandwich plates using a simple shear deformation theory. Structural Engineering, Mechanics 2016; 58(3): 397422.
[26] Bousahla AA, Benyoucef S, Tounsi A, Mahmoud SR. On thermal stability of plates with functionally graded coefficient of thermal expansion. Struct. Eng. Mech. 2016; 60(2): 313335.
[27] Chikh A, Tounsi A, Hebali H, Mahmoud SR. Thermal buckling analysis of crossply laminated plates using a simplified HSDT, Smart Structures Systems 2017; 19(3): 289297.
[28] Dozio L, Ricciardi M. Free vibration analysis of ribbed plates by a combined analytical–numerical method. Journal of Sound and Vibration 2009; 319: 681–697.
[29] Ney SF, Kulkarni GG. On the transverse free vibration of beamslab type highway bridges. Journal of Sound and Vibration 1972; 21: 249–261.
[30] Attia A, Tounsi A, Bedia EA, Mahmoud SR. Free vibration analysis of functionally graded plates with temperaturedependent properties using various four variable refined plate theories. Steel and composite structures 2015; 18(1):187212.
[31] Balendra T, Shanmugam NE. Free vibration of plate structures by grillage method. Journal of Sound and Vibration 1985; 99: 333–350.
[32] Dowell EH. Free vibrations of an arbitrary structure in terms of component modes. Journal of Applied Mechanics 1972; 39: 727–732.
[33] Liew KM, Xiang Y, Kitipornchai S, Meek JL. Formulation of Mindlin–Engesser model for stiffened plate vibration. Computer Methods in Applied Mechanics and Engineering 1995; 120(3–4): 339–353.
[34] Berry A, Locqueteau C. Vibration and sound radiation of fluidloaded stiffened plates with consideration of inplane deformation. Journal of the Acoustical Society of America 1996; 100 (1): 312–319.
[35] Asku G, Ali R. Free vibration analysis of stiffened plates using finite difference method. Journal of Sound and Vibration 1976; 48: 15–25.
[36] Bhimaraddi A, Carr AJ, Moss PJ. Finite element analysis of laminated shells of revolution with laminated stiffeners. Computers and Structures 1989; 33: 295–305.
[37] Harik IE, Guo M. Finite element analysis of eccentrically stiffened plates in free vibration. Computers and Structures 1993; 49: 1007–1015.
[38] Zeng H, Bert CW. A differential quadrature analysis of vibration for rectangular stiffened plates. Journal of Sound and Vibration 2001; 241: 247–252.
[39] Rao SS. Vibration of Continuous Systems. Wiley; Hoboken: 2007.
[40] Hetnarski RB, Eslami MR, Gladwell GML, Thermal stresses: advanced theory and applications. Springer: 2009.
[41] Spiegel MR. Advanced Mathematics for Engineers and Scientists. New York: McGrawHill; 1971.
[42] Sneddon IN, The use of integral transforms, McGrawHill; 1972.