Journal of Sound and Vibration (1986) 106(2), 333-345

DYNAMIC RESPONSE OF CIRCULAR PLATES IN CONTACT

WITH A FLUID SUBJECTED TO GENERAL DYNAMIC
PRESSURES ON A FLUID SURFACE

K. NAGAYA AND K. NAGAI

Department of Mechanical Engineering, Gunma University, Kiryu, Gunma 376, Japan

(Received 16 March 1985, and in revised form 5 June 1985)

This paper is concerned with a method for solving dynamic response problems of a
circular plate in contact with a fluid whose surface is excited by general dynamic pressures.
By utilizing the Fourier expansion and the Laplace transform methods, the expression for
the dynamic response of displacement is obtained in a general form which is applicable
to general dynamic pressures. As applications, numerical calculations have been carried
out for three types of sinusoidal, trapezoidal and explosive pressures. The results obtained
in a certain type of impact pressure are compared with the exact ones.

1. INTRODUCTION
Plates are important elements of machines and structures, and their vibration problems
have been investigated in detail by many investigators [l-lo]. The vibration of structures
including plates in contact with a fluid also is of importance in many fields, and its study
has been reported [ll]. Kubota and Suzuki [12] analyzed plate vibration problems by
making use of the theory of potential flow, neglecting the viscosity of the fluid. However,
in this study, they neglected the displacement in the radial direction of the circular plate,
and applied elementary ring theory without consideration of the circular plate boundary
conditions. They also performed experimental tests, but in their experiment the boundary
conditions were inside edge clamped and outside edge free, different from the theoretical
boundary conditions, so that comparisons between theoretical and experimental results
could be made only for the effects of fluid inertias. Nagaya and Takeuchi [ 131 considered
the bottom plate of a tank, and analyzed the forced vibrations of plates of arbitrary shape
in contract with a fluid by utilizing the equation of motion of the plate and the Navier-
Stokes equation, with inclusion of the viscosity of the fluid. In this study harmonic
vibration problems due to the base excitation were investigated, but studies of transient
response problems have not been investigated thoroughly. Situations in which the fluid
is vaporizing in the tank, and the excitation is due to various dynamic pressures are often
observed in practice. In such problems the dynamic pressure is a complicated function
of time, and hence the analysis cannot be developed by using the general integral transform
procedure.
For this situation, this paper presents an analysis of transient response problems of
circular plates in contact with a fluid, whose surface is excited by general dynamic
pressures. For treating complicated dynamic loads, Nagaya suggested a method in which
the general dynamic load was expanded into Fourier series and the differential equations
were solved by use of the Laplace transform methods. By using this procedure he analyzed
the transient response of towers and piles excited by waves [14] or earthquakes [ 151. In
the present paper this method is extended to problems of a circular plate in contact with
333
0022-460X/86/080333 + 13 %03.00/O @ 1986 Academic Press Inc. (London) Limited
334 K. NAGAYA AND K. NAGAI

a fluid, and the theoretical solutions for general transient loads are derived. The results
obtained are applicable to any such problem for which the numerical load versus time
data are available for inclusion in the numerical calculation process. As applications,
numerical calculations have been carried out for sinusoidal, trapezoidal and explosive
impact pressures, and the dynamic behaviours of the plate have been investigated in
detail. In order to prove the correctness of this analysis, the results obtained in the case
of sinusoidal pressures are compared with the exact ones.

2. METHOD OF ANALYSIS

In this analysis a bottom plate of a tank is considered to provide an analytical model.

As shown in Figure l(a) the bottom plate, of circular shape, is welded into the rigid
cylinder and is on an elastic foundation. Since the welded edge of the plate is elastic,
there is a bending slope with the welded edge indicated by a spiral spring, the analytical
model is shown as Figure l(b). If one considers a viscous fluid, the problem must be

(0) (b)

Figure 1. Definition sketch for the title problem.

solved by using the Navier-Stokes partial differential equation. In this case, as the fluid
vibrates along the z direction and the side wall,of the tank is rigid, the velocity of the
fluid in the radial direction is negligible. Under this assumption, with the origin of
cylindrical co-ordinates r, 0 and z at the center of the circular plate, and upon applying
the equation of continuity, the Navier-Stokes equation of a non-compressible and viscous
fluid is
av a% 1 au 1 a2v
-=_Lap+v,
at p'az ( -g+--+‘T,
rar r ae >
(1)

where v is the velocity of the fluid, p’ is the density, p is the pressure, v’is the kinematic
viscosity, and t is the time. Under the initial condition (v) r=O= 0, the Laplace transforma-
tion f(r, 8, z, s) = (Tf(r, 0, z, t) ems’dt of equation (1) yields
a~/az=~{v*-(i/~')~}~, (2)
where V2 is the two-dimensional Laplacian operator, and y = p’~’ is the coefficient of
viscosity. When the surface of the fluid is excited by a general dynamic uniformly
 PLATE SUBJECT TO TRANSIENT FLUID PRESSURE 335

distributed pressure, equation (2) can be written as

~=/1{v*-(1/z+}ii(Z+h’)+ Q, (3)

where h’ is the depth of the fluid and Q is the pressure applied to the fluid surface. The
equation of motion for the bottom plate on a Winkler foundation with a uniform initial
tension S due to the pressure of the tluid is

DV2V2w - SV2w +pha2w/at2+ kw = (~)~zz~,, 14)

where D = Eh3/ 12( 1 - v’) is the flexural rigidity, w is the bending displacement, p is the
density, h is the thickness, E is the Young’s modulus, Y is Poisson’s ratio of the plate,
respectively, and k is the spring constant per unit area of the elastic foundation. Applying
the Laplace transformation to equation (4) under the initial conditions ( w),=~ = 0 and
(aw/dt),,, = 0 gives the transformed equation. The velocity of the plate must be equal to
the velocity of the fluid over the plate surface. Hence, substituting this condition, i.e.,
.sM1= ( u),=,~, into equation (3) yields

(P)~_,, = ph’s{V2 - (l/ v’)s}w + Q. 15)

Since the pressure of the fluid over the plate surface acts on the plate, by substituting
equation (5) into the transform of equation (4), one obtains the following partial differen-
tial equation for the displacement w of the plate:

(6)

For a solid plate, the homogeneous solution of equation (6) can be written by

*I= i ? En[AjnJn(Qr)+CjnIn(P~)I@jm (7)

j=1 n=O

cr2=(q*+A4y-~, p2=(7y+A4y2+~,

T =t{(slD)+(~h'slD)l, A4= -{(ph.s2/D)+(p’h’s2/D)+(k/D)},

F,, = l/2 (n = O), E, = 1 (n ?= l), @,” = cos (no), Q2” = sin (n6). (8)
Ajn and C,, are constants of integration, J,(ar) is the Bessel function of the first kind of
order n, and I,(@-) is the modified Bessel function of the first kind, both with complex
arguments. The particular solution of equation (6) can be found by making use of
Lagrange’s method of variation of parameters. The general solution i+ is given by the
addition of the homogeneous and the particular solutions,

* = i ? EnlAjnJn(ar)+Cj,In(Pr)l@jn
j=1 n=O

+tL(Pr)L(PS) - K&)L(PS)l (9)

Here Y,(ar) is the Bessel function of the second kind of order n, and K,(/?r) is the
modified Bessel function of the second kind.
336 K. NAGAYA AND K. NAGAI

Since the pressure Q may be a complicated function of time, it may be difficult to

apply the Laplace transformation to the pressure Q directly. Hence, in this paper, the
pressure Q is expanded into the following Fourier series [14,15]:

Qn= 00 $, Em[qjnmcos (2r~ny) + qi*nmsin (2rmr)], (11)

4jnm = 2
J0
’(Qjn/Qo) COS(2rmy) dyp *
qjnm = 2 l(Qjn/ Qo) sinOrmy)
J0
dy,

E, = l/2 (m =O), .a,=1 (m>l), Y = tl To, (12)

Q. is a reference pressure which is utilized to obtain dimensionless expressions, To is a
reference time and y is the dimensionless time. Applying the Laplace transformation to
equation (11) yields

qjnmS + $Tnn(2~ml TO)

ojn=Qo f urn s2+c2Tm,TJ2 (13)
m=O 0

and substituting equation (13) into equation (9) gives

* = ? ? a,[Aj,J,(ar) + Cj,I,(Br)]@j”
j=ln=O

- Y,(cu)J,(af)l+ (2I~)[L(Pr)L(PT) - L(Pr)L(PS)lIS d5 @‘jn. (14)

When the edge of the plate is restrained against rotation and the translation displacement
is fixed, the boundary conditions are

(G)r=R=O, (fir>,,, = k’(iG/dr),=,, (1%

where R is the radius of the circular plate, and k’ is the spring constant against rotation
per unit length along the plate edge. The boundary condition that the velocity of the fluid
is zero [(v),,, = 0] at the side wall mugt be satisfied. In this analysis the solution (equation
(14)) satisfies the condition of continuity between the plate and the fluid, and hence the
solution satisfies the boundary condition of the fluid at the side wall automatically. If
the fluid surface is excited by a dynamic uniformly distributed pressure, the fluid pressure
(P)*,~ has no dependence on the direction 0, and the particular solution in equation (14)
becomes [ - Qo/(Da2P2)CZxo c,[q,,,s+qZ(27rm/ To)]/[s2-P(2rrm/ To)‘]& Then one may
consider the case of n = 0 in equation (14), and substituting equation (14) into equation
(15) yields

q,,,s + qZ(2nml To)

Jo(aR)A,o+lo(/3R)CIo=2Q0 f
Da2p2 ,,,EOEm s2+ (2nm/ To)2 ’

{a2Jo(cyR)-[a(l-vy)/R-ak’/D]J,(aR)}A,,,

-{P’Io(PR)-Ml- v)lR-pk’lDlI,(PR)}C,o=O. (16)

In equations (16), since qi,,,, and q$,, have no terms with respect to j, n, these are replaced
by qm and qz respectively. The constants of integration are determined from equations
 PLATE SUBJECT TO TRANSIENT FLUID PRESSURE 337

(16), and hence the displacement is given by

*
ti = ; GlQo[4mS+ qnl(27w To)1 [n, _ nl,
(17)
m=O Dcx*p*[S*+ (27rm/ To)Q

A = Io(PR)& + JoCaR)&, 4 = Io( + Jo('Y~)&,

B,=a2Jo(crR)-[a(l-v)/R-ak’/D]J,(aR),

B2=P2MPRHP(1- v)IR-W/W,W). (18)

The co-ordinates of the singular points are necessary in order to perform the Laplace
inversion transformation. In this case, they are

s = *2iTm/ To, s = f i( Dt4/phR4)li2, (19)

where i = &i and ,$ is the eigenvalue of the circular plate. The residues for s = + 2irrm/ T,,
are

RI = f ia,FG’[(-iq, + qz)f(i@) e-*jirmy

rn=O

-(iq, + qz)f(+) e2inmy ]/{2(1+H)[(27rm)2-KG2/(1+H)]} (20)

and the residues for s = *i( Dt4/phR4)l’* are written by making use of the Ith eigenvalue
6,

R,, = f 2ia,FG&(iG~:q, +2rmq*,)

m=O

x (AT),=,, eif:Gy [5:(1+H)-K1[(2?rm)2-G21fl(~)~_~,}, (21)

where the dimensionless expressions are

H = p’h’lph, K = R4kl D, S* = R’S/ D,

E = v’ph’/ D, G= (D/phR4)1’2To,

F=R’Q,ID, K’= Rk’l D, Q = (2nm/ G)l/‘, f= r/R,

(Y2=(aR)2+,*2+A*4)1/2-,,*, ~=(~R)2+,*2+~*2)“2+V*,

f(p)=A:‘lA*-1, A:= I,(P;-)Bf+J,(ar)B;, A*= 10(~)B~+Jo(~)B~,

BT=c?~J,(~)-[~~(~-v)-(YK’]J,((U), B,*=pI,(&-[&l- v)+K’]I,@),

n*=[S*+i(p*(eH)“*]/2, A*4=q4(l+H)-K. (22)

The Ith eigenvalue obtained from A = 0 has a conjugate value so that 6 = 5; + i&’ and
and & = 5; -i&‘, where 5; is the real part and &’ is the imaginary part. Hence the
dimensionless displacement w/R as a function of time is the sum of the residues,

wlR=R,+ f (R,,+&,). (23)

I=1

where R2, is given by replacing 6 by .$ in equation (21).

338 K. NAGAYA AND K. NAGAI

3. NUMERICAL EXAMPLES

As examples numerical calculations have been carried out for the three cases of
sinusoidal, trapezoidal (shown in Figure 2(b)) and explosive (shown in Figure 2(a))
impact pressures. In the numerical calculation, the equations of the curve of the dynamic
pressure versus time are required in order to obtain the Fourier coefficients (see equations

Figure 2. Geometry of the dynamic impact pressures versus time. (a) Explosive type; (b) trapezoidal type.
0, Fourier, N = 20; -, exact.

(12)). However it is difficult to obtain those equations for complex dynamic loads,
especially from experimental results for the pressures. Hence, in this paper, the Fourier
coefficients have been obtained by using the dynamic pressures at many sample points,
and assuming the equation between an arbitrary sample point and the next one to be a
straight line. When the Fourier coefficients q,,, and q*, have been determined, equation
(23) gives the numerical values of the transient response of the circular plate in contact
with the fluid, but the number of residues and the truncation number of the Fourier series,
N, in equation (23) must be decided according to the convergence of the series solution.
In the present numerical example, the residues up to the third mode have been considered,
with N up to 15 in the case of a sinusoidal pressure and N up to 20 in the other cases,
to give sufficient accuracy. Figures 2(a) and (b) show the comparisons between the curves
considered and the results of the Fourier approximation with N = 20. The Fourier series
solutions show good agreement with the curves although in Figure 2(b) there are slight
discrepancies the sharp corners. To verify the present approximate results, the exact
dimensionless displacement of the circular plate for a typical case of a sinusoidal pressure
 PLATE SUBJECT TO TRANSIENT FLUID PRESSURE 339

0 02 04 06 08 IO
Y

Figure 3. Dynamic displacement response at the center of the circular plate subjected to a sinusoidal pressure
PO= sin (2~7). K = 0.25, S* = 0, H = 0, E = 0, F = 0.1, K' = 104,Y = 0.3. -, Exact; 0, present results.

[Q/ Q,, = sin (27r7)] is shown in Figures 3 and 4 (see the Appendix). The present analytical
results are in very good agreement with the exact ones.
The displacement response curves for the dynamic pressures of Figure 2(a) type are
shown in Figures 5-10. In the figures, the effects of the time T,, which is the load duration
(see Figures 2(a) and (b)) is shownin Figure 5. It can be seen that the dynamic response
curves are similar to those for an impact pressure when G[ = (D/phR4)“* T,] is large,
but when G is small, the dynamic response curves are different from the impact pressure
curves because of the influence of the vibration of both the fluid and the plate. Figure 6
shows the response of the displacement versus the fluid inertia H[ = p’h’/ph). In this
case the natural frequencies and the displacements decrease with increasing values of the

12 -1

:
0
5 08
e
1

04

0 02 04 06 00 10
r/R

Figure 4. Modes for a sinusoidal pressure PO = sin (277~). K = 0.25, S* = 0, H = 0,E = 0,F = 0.1,K' = IO",
v = 0.3, y = 0.28. -, Exact; 0, present results.
340 K. NAGAYA AND K. NAGAI

OL 1 I / i i
0 02 04 06 08 10

Figure 5. Dynamic displacement response at the center of the circular plate due to an explosive impact
pressure as shown in Figure 2(a). K = 0.25, S* = 0, H = 0,E = 0,F = 0.1, K’= 104, v = 0.3.

fluid inertia H. It is interesting to see that in Figure 6 the first mode vibration of the plate
appears in the curve for H = 0 and G = 20, but does not appear in Figure 5( H = 0, G = 50).
Figure 7 depicts the dependence of the displacement on the viscosity of the fluid. As the
viscosity E[ = v’@‘/D] increases, the damping increases and the natural frequencies
decrease slightly, but the effect of the viscosity can be neglected in case of fluids with
small viscosity like water and oil (about E = lo-l6 through lO_“). Figure 8 shows the
displacement response of the plate restrained by the spring against rotation around its
edge. One can see that when the fluid becomes thick, the response becomes steady for
large spring constants, Figure 9 shows the effects of the initial tension S*[ = R’S/D] on

Figure 6. Effects of fluid inertia on the dynamic displacement response at the center of the circular plate to
the pressure as shown in Figure 2(a). K = 0.25, S* = 0, E = 0, G = 20, F = 0.1, K'= 104, v = 0.3.
 PLATE SUBJECT TO TRANSIENT FLUID PRESSURE 341

Figure 7. Effects of viscosity of fluid on the dynamic displacement response at the center of the circular plate
to the pressure as shown in Figure 2(a). K = 0.25, S* = 0, H = 200, G = 20, F = 0.1, K’ = 104, Y = 0.3.

the displacement. In this case the natural frequencies increase with the tension, of course.
Figure 10 shows the relation between the displacement and the rigidity of the Winkler
foundation. In the figure the natural frequencies increase with the spring constants.
The numerical results for a dynamic pressures of Figure 2(b) type are shown in Figures
11 and 12. The eigen period of the first mode of the circular plate in contact with the
fluid is r, and hence the appropriate dimensionless expression is r/ T, = 2rr/(&G). By
use of this expression it can be seen that in Figures 5-12 the periods of the response
curves agree with the eigen periods of the first mode for small values of G. This implies
that the circular plates vibrate in the principal mode of vibration. Since the inertia force
is small for large values of G, the response curve has the same shape as that for an impact
pressure, as shown in Figure 5 (G = 50). However the curve for G = 50 in Figure 11 has
some small oscillations and the period is different from that of the first vibration mode.
In this case, there are some errors due to lack of convergence and the inevitable numerical
calculation errors. When G has a comparatively small value, the error becomes small

I I
-6-O
0 02 0.4 06 08 10

Figure 8. Effects of springs against rotation around the edge on the dynamic displacement response at the
center of the circular plate to the pressure as shown in Figure 2(a). K = 0.25, S* = 0, H = 200, E = 0, G = 20,
F = 0.1, V = 0.3.
342 K. NAGAYA AND K. NAGAI

30

2-o

; 1.0
b
X

cc
2 0

-1.0

-2.0 L I / I I J
0 0.2 0.4 O-6 0.0 10
Y

Figure 9. Effects of initial tensions on the dynamic displacement response at the center of the circular plate
to the pressure as shown in Figure 2(a). K = 0.25, H = 0, E = 0, G = 5.0, F = 0.1, K’= 104,Y = 0.3.

and curves with sufficient accuracy, which include the effects of the first mode vibration,
can be obtained, as shown in the case of H = 0 in Figure 12.

4. CONCLUSION

An analysis for determining the dynamic response of a circular plate in contact with
a fluid whose surface is excited by general dynamic pressures has been presented. The
results may be summarized as follows.

2.0

F I.0
sz
x
e
B 0

-1.0

Figure10. Effects of Winkler’s foundation on the dynamic displacement response at the center of the circular
plate to the pressure as shown in Figure 2(a). S* = 0, H = 0, E = 0, G = 5.0, F = 0.1, K’= 104, Y = 0.3. -,
K = 5.0; - - -, K = 50; - - -, K = 100.
 PLATE SUBJECT TO TRANSIENT FLUID PRESSURE 343
30 ._

20

rz- 10
b

Y
9 0
1

-10

-2 0 il w v v
I i /
0 02 04 06 08 i-0
Y

Figure 11. Dynamic displacement response at the center of the circular plate subjected to a trapezoidal
impact pressure as shown in Figure 2(b). K = 0.25, S* = 0, H = 0,E = 0,F = 0.1,K' = 104, Y = 0.3.

(1) By utilizing the equation of motion of the plate and the Navier-Stokes equation
with inclusion of the fluid viscosity, an expression for the transient response has been
obtained in a general form which is applicable to general transient pressures. The analysis
has been developed by a combination of Fourier expansion and Laplace transform
methods.
(2) As applications, numerical calculations have been carried out for three cases, of
sinusoidal, trapezoidal and explosive impact pressures, and the characteristics of the
dynamic response for those impact pressures have been investigated in detail.
(3) In order to prove the correctness of the analysis, the results obtained in the case
of the sinusoidal impact pressure have been compared with those obtained by the exact
method. It is clear that the present analytical results agree well with the exact ones, and
that the present results possess sufficient accuracy.

Figure 12. Effects of fluid inertia on the dynamic displacement response at the center of the circular plate
to the pressure as shown in Figure 2(b). K = 0.25, S* = 0, e=O, G=20, F=O.l, K’=104, v=O.3.
344 K. NAGAYA AND K. NAGAI

(4) Since this method involves no restrictions on the shape of the dynamic pressures,
the equation obtained-is applicable to any dynamic response problems for circular plates,
when the numerical data of the particular pressure curve of interest is used in the numerical
calculation process.

REFERENCES
1. A. W. LEISSA 1969 NASA SP-160. Vibration of plates, 41-159.
2. H. D. CONWAY 1960 Transactions of the American Society of Mechanical Engineers, Journal
of Applied Mechanics 28, 288-291. The bending and Aexural vibrations of simply supported
polygonal plates by point-matching.
3. T. IRIE, G. YAMADA and Y. NARITA 1977 Bulletin of the Japan Society of Mechanical Engineers
20, 10851092. Free vibration of a rectangular plate supported on the sides and some segments.
4. T. IRIE, G. YAMADA and Y. NARITA 1978 Bulletin of the Japan Society of Mechanical Engineers
21, 1696-1702. Free vibration of clamped polygonal plates.
5. J. C. M. YU 1971 Journal of the Acoustical Society of America 49, 781-785. Application of
conformal mapping and variational method to the study of natural frequencies of polygonal
plates.
6. P. A. A. LAURA and R. H. GUTIERREZ 1976 Journal of Sound and Vibration 48, 327-332.
Fundamental frequency of vibration of clamped plates of arbitrary shape subjected to a
hydrostatic state of in-plane stresses.
7. P. A. A. LAURA, R. H. GUTIERREZand D. S. STEINBERG1977 Journal of Sound and Vibration
55,49-53. Vibration of simply-supported plates of arbitrary shape carrying concentrated masses
and subjected to a hydrostatic state of in-plane stresses.
8. J. MAZUMDAR and J. R. COLEBY 1976 Journal of Sound and Vibration 45,503-512. Simplified
approach to the vibration analysis of elastic plates due to sonic boom.
9. D. BUCCO, J. MAZUMDAR and G. SVED 1979 Journal of Sound and Vibration 67, 253-262.
Vibration analysisof plates of arbitraryshape-a new approach.
10. K. NAGAYA 1980 Transactions of the American Society of Mechanical Engineers, Journal of
Apphed Mechanics 47, 620-626. Dynamic response of a plate with arbitrary shape.
11. R. F. BLEVINS 1984 Structural Vibrations in a Fluid. New York: Van Nostrand.
12. Y. KUBOTA and T. SUZUKI 1984 Transactions of the Japan Society of Mechanical Engineers C
50, 243-247. Added mass effect on disc vibrating in fluid, (In Japanese.)
13. K. NAGAYA and J. TAKEUCHI 1984 Journalof the AcousticalSocietyofAmerica 75,1511-1518.
Vibration of a plate with arbitrary shape in contact with a fluid.
14. K. NAGAYA 1985 American Society of Civil Engineers Journal of Waterway Part Coastal and
Ocean Engineering 111, 337-353. Response of variable cross-sectionalmembersto waves.
15. K. NAGAYA and Y. HAI 1985 Journal of Sound and Vibration 103, 119-138. Seismic response
of underwater members of variable cross-section.

APPENDIX

The exact solution for a circular plate whose surface is excited by a sinusoidal dynamic
uniformly distributed pressure PO sin (at) is obtained as follows.
The equation of motion of the circular plate on a Winkler foundation is

DV2V2w + pha2w/at2+ kw = POsin (at). (Al)

Applying the Laplace transformation to equation (Al) under the initial conditions (w),=,, =
0 and (aw/at),,, = 0 gives

V2V2~-J4P = Poa/[D(s2+a2)], where 14= -(phs’/D+k/D). (A2, A3)

The general solution of equation (A2) is obtained as

&a
‘=jil“go En[AjnJn(O'r)+ CjnIn(lr)]Qjn - D14(s2 + a’)’ (A4)
 PLATE SUBJECT TO TRANSIENT FLUID PRESSURE 345

When the edge of the plate is fixed, the boundary conditions are

(W),=R =o, (?M/?Ir),=, = 0. (AS)

Since the circular plate is excited by a uniformly distributed pressure, the terms with
n = 0 must be chosen in equation (A4). Then equations (A5) yield

A,,JdW) + GoI, = 2E,~l[~54(~2+ a*)13 A,,J,@R) - C,,I,(lR) =O. (A6)

From equations (A6), one has

w={~,a/D~4(s2+a2>}{(n:/n’)-1>, (A7)

A: = II(5R)Jo(5r) + J1(5RM5r), A’= 1,(5R)Jo(5R)+J,(5R)Io(5R). (Ag)

Performing the Laplace inversion transformation and substituting a = 27r/ T,, gives the
displacement of the circular plate as

w/R=Rft f (R;,+l?,*,), (A9)

/=I
where I?$ is given by replacing & ( = &+i&‘) by & ( = 51 -i.$‘) in equation (All ),
R* = iF*G’[f*(i(P*) e-*i”y-f*(cp*) e2imY] 4irr&GF*(A:*),,s, eiG’fv
I
2[ (274’ - KG21
, R%=
(S:- WN27d2- ~‘S:1(d~‘*/‘&4.=,,’
(AlO, Al 1)

K = R4k/ D, G = ( D/phR4)“*To, F*h R3Po/D,

Cp*= (2~r/G)“~, c4=((gR)4=~4-K, J= r/R,

f*(cp)=A:*/A’*-1, A:* = I,(<IJo<&->+

J,<~IL,(&L
A’* = IAi?Jo(~I+ J,k%(h. (AL?)