S. K. Datta and T. Chakraborty Department of Mechanical Engineering and CERES University of Colorado Boulder, Colorado 80309

Similar documents
Detection of Protective Coating Disbonds in Pipe Using Circumferential Guided Waves

Quasi-Rayleigh Waves in Butt-Welded Thick Steel Plate

Quantitative Crack Depth Study in Homogeneous Plates Using Simulated Lamb Waves.

EWGAE 2010 Vienna, 8th to 10th September

Long Range Ultrasonic Testing - Case Studies

Rayleigh Wave Interaction and Mode Conversion in a Delamination

Electronic Noise Effects on Fundamental Lamb-Mode Acoustic Emission Signal Arrival Times Determined Using Wavelet Transform Results

GUIDED WAVES FOR DAMAGE MONITORING IN PLATES FOR NOTCH DEFECTS

Corrosion detection under pipe supports using EMAT Medium Range Guided Waves

Time Reversal FEM Modelling in Thin Aluminium Plates for Defects Detection

Ultrasonic Guided Wave Testing of Cylindrical Bars

Study on Propagation Characteristics of Ultrasonic Guided Wave for EMAT Sensor

Guided Wave Travel Time Tomography for Bends

In-Situ Damage Detection of Composites Structures using Lamb Wave Methods

Effect of Attenuation on Inspection Range and Sensitivity in Long- Range Guided Wave NDT of Coated and Buried Pipes

High contrast air-coupled acoustic imaging with zero group velocity Lamb modes

Excitation and reception of pure shear horizontal waves by

A Simulation Study of Attenuation Factors in a Gas Pipeline Guided Wave Testing

DETECTION OF CORROSION IN BOTTOM PLATES OF GAS AND OIL TANKS USING GUIDED ULTRASONIC WAVES AND ELECTROMAGNETIC ULTRASONIC (EMAT) TRANSDUCERS

Guided Wave in Engineering Structures Using Non-Contact Electromagnetic Acoustic Transducers A Numerical Approach for the Technique Optimisation.

Mode mixing in shear horizontal ultrasonic guided waves

Selective Excitation of Lamb Wave Modes in Thin Aluminium Plates using Bonded Piezoceramics: Fem Modelling and Measurements

Penn State University ESM Ultrasonics R&D Laboratory Joseph L. Rose Research Activities

arxiv:physics/ v1 [physics.optics] 28 Sep 2005

THE LONG RANGE DETECTION OF CORROSION IN PIPES USING LAMB WAVES

Non-Destructive Method Based on Rayleigh-Like Waves to Detect Corrosion Thinning on Non- Accessible Areas

G. Hughes Department of Mechanical Engineering University College London Torrington Place London, WClE 7JE, United Kingdom

ON LAMB MODES AS A FUNCTION OF ACOUSTIC EMISSION SOURCE RISE TIME #

Professor Emeritus, University of Tokyo, Tokyo, Japan Phone: ;

A NEW APPROACH FOR THE ANALYSIS OF IMPACT-ECHO DATA

Use of Lamb Waves High Modes in Weld Testing

Passive Polymer. Figure 1 (a) and (b). Diagram of a 1-3 composite (left) and a 2-2 composite (right).

A Numerical study on proper mode and frequency selection for riveted lap joints inspection using Lamb waves.

VERSATILE USAGE OF ELECTROMAGNETIC ACOUSTIC TECHNOLOGIES FOR IN-LINE INSPECTION OF AGEING PIPELINES

Fig. 1 Feeder pipes in the pressurized heavy water reactor.

Application of Ultrasonic Guided Waves for Characterization of Defects in Pipeline of Nuclear Power Plants. Younho Cho

PRACTICAL ASPECTS OF ACOUSTIC EMISSION SOURCE LOCATION BY A WAVELET TRANSFORM

Testing of Buried Pipelines Using Guided Waves

Determination of the width of an axisymmetric deposit on a metallic pipe by means of Lamb type guided modes

EFFECTS OF LATERAL PLATE DIMENSIONS ON ACOUSTIC EMISSION SIGNALS FROM DIPOLE SOURCES. M. A. HAMSTAD*, A. O'GALLAGHER and J. GARY

Measurement of phase velocity dispersion curves and group velocities in a plate using leaky Lamb waves

Properties of Interdigital Transducers for Lamb-Wave Based SHM Systems

KAERI Feeder Tube Inspection Using EMAT Generated Circumferential Guided Waves

EMAT Application on Incoloy furnace Tubing Ramamohan Reddy M (ASNT Level III UT, PCN Level III UT,PAUT&TOFD)

ULTRASONIC GUIDED WAVE ANNULAR ARRAY TRANSDUCERS FOR STRUCTURAL HEALTH MONITORING

Simulation of ultrasonic guided wave inspection in CIVA software platform

Guided Wave Inspection of Supported Pipe Locations Using Electromagnetic Acoustic Transducers

Ultrasonic Testing using a unipolar pulse

Ultrasonic Guided Wave Applications

FATIGUE CRACK DETECTION IN METALLIC MEMBERS USING SPECTRAL

18th World Conference on Non-destructive Testing, April 2012, Durban, South Africa

Acoustic emission signal attenuation in the waveguides used in underwater AE testing.

Keywords: Ultrasonic Testing (UT), Air-coupled, Contact-free, Bond, Weld, Composites

describe sound as the transmission of energy via longitudinal pressure waves;

CHARACTERISTICS AND APPLICATIONS OF ELECTROMAGNETIC SURFACE WAVE TRANSDUCERS

ASSESSMENT OF WALL-THINNING IN CARBON STEEL PIPE BY USING LASER-GENERATED GUIDED WAVE

Piezoelectric transducer excitation for guided waves propagation on pipeline with flexural wave modes

Multiple crack detection of pipes using PZT-based guided waves

Laser Vibrometer Measurement of Guided Wave Modes in Rail Track. Sensor Science and Technology, CSIR Material Science and Manufacturing,

Assessment of lamination defect near the inner surface based on quasi-symmetric circumferential Lamb waves

DEVELOPMENT OF STABILIZED AND HIGH SENSITIVE OPTICAL FI- BER ACOUSTIC EMISSION SYSTEM AND ITS APPLICATION

In service application of EMAT in Boiler Water Wall Tubes and High Temperature Components

ABSTRACT 1. INTRODUCTION

Finite Element Modeling and Simulation of Ultrasonic Guided Wave Propagation using Frequency Response Analysis

Ultrasonic Guided Waves for NDT and SHM

Co-Located Triangulation for Damage Position

ULTRASONIC GUIDED WAVE FOCUSING BEYOND WELDS IN A PIPELINE

ACOUSTO-ULTRASONIC EVALUATION OF HYBRID COMPOSITES USING

DEVELOPMENT OF ULTRASONIC WAVE NONDESTRUCTIVE INSPECTION ROBOT WITHOUT COUPLING MEDIUM USING EMAT

SOME OBSERVATIONS ON RAYLEIGH WAVES AND ACOUSTIC EMISSION IN THICK STEEL PLATES #

Quick Assessment of the Anomalies in Concrete Structure Using Dispersive Characteristic of Surface wave

CRACK SIZING USING A NEURAL NETWORK CLASSIFIER TRAINED WITH DATA OBTAINED FROM FINI1E ELEMENT MODELS

DISBOND DETECTION AND CHARACTERIZATION USING HORIZONT ALL Y

SPARSE ARRAY TOMOGRAPHY SYSTEM FOR CORROSION EXTENT MONITORING H. Bian, H. Gao, J. Rose Pennsylvania State University, University Park, PA, USA

RECENT ADVANCEMENTS IN THE APPLICATION OF EMATS TO NDE

Experimental and theoretical investigation of edge waves propagation and scattering in a thick plate with surface-breaking crack-like defect

REFLECTION AND TRANSMISSION OF LAMB WAVES AT DISCONTINUITY IN PLATE Z. Liu NDT Systems & Services AG, Stutensee, Germany

Low wavenumber reflectors

FATIGUE CRACK CHARACTERIZATION IN CONDUCTING SHEETS BY NON

MODELLING AND EXPERIMENTS FOR THE DEVELOPMENT OF A GUIDED WAVE LIQUID LEVEL SENSOR

ULTRASONIC GUIDED WAVES FOR AGING WIRE INSULATION ASSESSMENT

USE OF GUIDED WAVES FOR DETECTION OF INTERIOR FLAWS IN LAYERED

THE EXTRACTION METHOD FOR DISPERSION CURVES FROM SPECTROGRAMS USING HOUGH TRANSFORM

Y. Li and R. B. Thompson Ames Laboratory Iowa State University Ames, Iowa 50011

Implementation of electromagnetic acoustic resonance in pipe inspection

Pipe Testing Using Guided Waves

Damage Detection in Stiffened Composite Panels Using Lamb Wave

THE USE OF MAGNETOSTRICTIVE EMAT TRANSDUCERS ON OXIDE SCALED BOILER TUBES

Research on An Inspection Method for De-bond Defects in Aluminum. Skin-Honeycomb Core Sandwich Structure with Guided Waves

APPLICATIONS OF GUIDED WAVE PROPAGATION ON WAVEGUIDES WITH IRREGULAR CROSS-SECTION. Zheng Fan

Air Coupled Ultrasonic Inspection of Steel Rubber Interface

APPLICATION OF ULTRASONIC GUIDED WAVES FOR INVESTIGATION OF COMPOSITE CONSTRUCTIONAL COMPONENTS OF TIDAL POWER PLANTS

PIEZOELECTRIC WAFER ACTIVE SENSORS FOR STRUCTURAL HEALTH MONITORING STATE OF THE ART AND FUTURE DIRECTIONS

Sensitivity analysis of guided wave characters for transducer array optimisation on pipeline inspections

NONDESTRUCTIVE EVALUATION OF ADHESIVE BONDS USING LEAKY LAMB WAVES* Cecil M. Teller and K. Jerome Diercks. Yoseph Bar-Cohen and Nick N.

Detection of Cracks at Rivet Holes in Thin Plates Using Lamb-Wave Scanning

NUMERICAL MODELING OF AIR-COUPLED ULTRASOUND WITH EFIT. D. E. Chimenti Center of Nondestructive Evaluation Iowa State University Ames, Iowa, USA

ULTRASONIC SIGNAL CHARACTERIZATIONS OF FLAT-BOTTOM HOLES IN

FEASIBILITY OF THICKNESS MAPPING USING ULTRASONIC GUIDED WAVES. Pierre Belanger

A Numerical Study of Depth of Penetration of Eddy Currents

Transcription:

LAMB WAVE MODES IN COAL-TAR-COATED STEEL PLATES INTRODUCTION N. C. Banik, M. Land, and G. L. Puckett Research and Engineering T. D. Williamson, lnc. P.O. Box 2299 Tulsa, Oklahoma 74101 S. K. Datta and T. Chakraborty Department of Mechanical Engineering and CERES University of Colorado Boulder, Colorado 80309 In order to study the feasibility of using ultrasonic Lamb wave modes for detection and sizing of corrosion-related flaws in buried steel pipelines, we have calculated Lamb wave modes and performed numerous experiments on steel plates coated on one side with coal-tar enamel. The purpose of this paper is to discuss these theoretical and experimental results. To protect a pipeline from corrosion, it is coated on the external surface with a thick layer of coating. The most frequent1y used coatin~ is coal-tar enamel. It is usually applied in the factory, but over-the-ditch applicatlon is not uncommon, especially in older lines. The nominal thickness of coating may vary from one-quarter to one-half of the steel-wall thickness, but little attention is paid to assure a constant thickness or bonding strength of coating. In the lines which have been in the ground for years, the state of coating can vary considerably pver the length of the line. Lamb waves excited and propagating in the steellayer will be modified in the varying coating conditions as well as in the presence of corrosion-related flaws in steel. The question is whether or not the effects of coating variation are discernible from those due to corrosion. As a first step, we have studied the effects of coal-tar enamel coating on flat steel plates. We have done an exact calculation of Lamb wave modes in coated steelelates for various thicknesses of coating. The method used is similar to that descnbed by Jones [1]. We have also done numerous experiments on steel plates and pipes coated with coal-tar enamel as a function of the thickness of coating. THEORY The geometry of the two-iayered medium is shown in Figure 1. The steel and coal-tar layers are bonded together rigidly at an interface. The thickness of each layer is finite and uniform; otherwise, each layer is elastically isotropic and homogeneous. Lamb waves propagate along the horizontal X-axis, and elastic displacements are in the X-Z plane. The calculation of Lamb wave modes in this medium is straightforward. Seeking solutions of the elastic wave equation in each 1345

z,~ ~~,~ Steel h1 Coal Tar :: x Fig. 1. Geometry of layered steel and coating medium. layer and invoking boundary conditions at the top and bottom surfaces and at the interface, one obtains an 8 x 8 determinantal equation as the dispersion relation. We have reduced the size of the determinantal equation to 4 X 4 using four of the eight boundary conditions. The equation has been solved for real and positive values of frequency as a function of the wave number. An efficient numerical iterative scheme has been used; otherwise, the solutions obtained are exact. The elastic parameters used in the calculation are given in Table 1. We have kept the thickness of steel plate fixed at 1.0 cm and varied the thickness of coating from O to.3 cm. TABLE 1 Chosen Elastic Properties of the Layers Longitudinal Shear Thickness Densit) Wave Velocity Wave Velocity (cm) (gm/cc (m/s) (m/s) Steel 1.0 7.86 5960 3235 Coal-tar 0-0.3 1.20 2500 1000 Figures 2a-2d show the calculated dispersion curves for the first severalloworder Lamb wave modes in steel in the presence of four different thicknesses of the coal-tar layer. The ordinate F is the angular frequency multiplied by the thickness of steel and divided by the shear wave velocity in steel. The abscissa K is the wavenumber multiplied by the thickness of steel. AlI propagating modes lying in the range 0-15 for F and K are shown. In Figure 2a, the coating thickness is zero and so the modes seen are the standard symmetric and antisymmetric modes in bare steel plates. There are eight such modes in the chosen range. In Figure 2b, the modes correspond to a steel plate with a thin, O.5-mm layer of coal-tar coating. The presence of even such a thin coating layer causes a substantial modification in the modal structure. First of ali, there is an extra mode in the same range of frequency and wavenumber. The mode marked 7 in Figure 2b does not correspond to any mode in Figure 2a. Second, at large wave numbers, many of the modes are pulled downward in frequency, tending towards Rayleigh and shear modes in coating rather than in steel. For example, mode 1 in Figure 2b behaves similarly as mode 1 in Figure 2a for K up to about 6.0, but afterwards deviates substantialiy downward, becoming the Rayleigh wave mode in coating. Third, at some points, degenerate symmetric and antisymmetric modes in the bare plates are split up in the presence of coating, and nondegenerate modes are perturbed due to nearby modes. In the presence of coating, the symmetry is broken--a mode may have both symmetric and antisymmetric components. Fourth, the critical frequency of higher order modes are lowered to some extent. In the limit of long wavelengths, the 1346

FRAC _ 0.00 FRAC... 0.05 Frequency vs. Wavenumber 15r-~~~,...,.;17-,""':"~"""7--r~7'"~-;T--r-:A 1. IL 1. 15 K K FRAC _ 0.20 FRAC... 0.30.... (c) K (d) K Fig. 2. Lamb wave modes in bare and coal-tar coated (b-d) steel plates. The vertical and horizontal axes are reduced frequency and wavenumber as defined in text. layered medium behaves like a single medium with elastic properties which are weighted averages of the elastic properties of the two media, causing a lowering of the critical frequency of steel plate modes in the presence of elastically soft coal-tar. As we increase the thickness of coating, the magnitudes of the above four characteristics increase. More and more modes show up in the viewing window of frequency and wavenumber; modes tend to become shear and Rayleigh wave modes in coating at lower wavenumbers; critical frequencies of higher order modes are further lowered; and more perturbations of modes take place. These effects are clearly seen in Figure 2d which are Lamb wave modes when the thickness of the coating layer is 0.3 cm. This thickness is quite common in pipelines. An intermediate case is shown in Figure 2c in which the coating thickness is 0.2 cm. The effect of lowering the frequency as a function of wavenumber will result in the lowering of the phase and the group velocity. The phase velocity is the ratio of the frequency and the wavenumber. The group velocity is obtained from the tangents of the dispersion curves. We have calculated the group velocity curves for modes shown in Figures 2a-2d. These are given in Figures 3a-3d where we have plotted the ratio of the group velocity to the shear wave velocity in steel as a function of the reduced frequency F. Clearly, there is a general trend of lowering the group velocity, especially at large values of the wavenumber. Comparing Figures 3a and 3b, we find that the first antisymmetric mode in the bare plate, asymptotical- 1347

FRAC_o.oo 2.o.---~~.,...:G:;:.r.:..:ou=rp~v..;:.e..;:.loc.:;:Ity~V':":S'...:,.F.:...:rec.:;q.:..:ue:..:,n.:..:cy'--.-~~---, llmc -o.. 2.or-~~~G_r_ou-Tp_v_e_IOC_I..:::..ty_v_s._F_req-;-Ue_n,...;cy'-,-~~---, FRAC,., 0.20 FRAC:: 0.30 ~or-~~~g~r...:,. Tup~v~e_IOC~I~~v_s_.F_re,q~u_en_c~y.---~~ ur-~~~g_ro_u,p_v_el_oc_i~~v_s._f_re--,qu_e...,n~cy~~~---, (c) (d) Fig. 3 Group velocity curves corresponding to modes shown in Figures 2a-2d respectively. ly, goes to the Rayleigh wave mode in steel, but in the presence of 0.05 cm of coating, it is substantially pulled downward at K values larger than 10, perhaps going to the Rayleigh wave mode in coating. That this mode asymptotically becomes the Rayleigh wave mode in coating has been verified by increasing the range of K. On increasing the thickness of coating (see Figures 3c and 3d), we find that more and more modes have group velocities moving towards the shear wave and Rayleigh wave speeds in coating and do so at lower values of wavenumber. In fact, as we increase the wavenumber, we see that the group velocity of a mode initially tends to the shear or the Rayleigh wave speed in steel and then swings towards the shear and the Rayleigh wave speeds in coating. This phenomenon, called "terracing" by Mindlin [2] is, according to him, "due to the coupling of the equivoluminal mode with the lower dilatational mode." 1348

To summarize the effects of coating as observed in the theoretical calculation, we note the following: (1) More and more new modes which become Rayleigh and shear modes in coating at short wavelengths appear with increasing coating thickness. (2) Phase and group velocities of bare-plate modes generally decrease with increasing coating thickness, but at an arbitrary operating point they may also increase because of new modes. (3) The critical frequency of a bare-plate mode may de crease with increasing coating thickness. EXPERIMENT In the inspection of pipelines, transmitters and receivers are placed on the inner steel surf ace of the pipe. The new modes which have group velocities much lower than the Rayleigh wave speed in steel cannot propagate in steel. They will propagate in coating, however. As a Lamb mode is excited in steel, during its propagation, a part of its energy willleak into the coating, creating one or more propagating modes in the coating. The number of such modes wiii depend partly on the coating thickness. As a result, there will be an apparent attenuation of the received Lamb wave signal. This apparent attenuation will increase with coating thickness untii the thickness is several times the wavelength. Because of the slowing of the group velocity, the arrival time of the received signal wiii be delayed with coating thickness, causing a phase delay in the signal at a given time. If a new mode appears at the operating frequency and wavenumber, there is also a possibility of a phase advancement. There may also be a detectable shift in the amplitude spectrum towards lower frequencies, but this will depend on frequency response of the transducer. Thus, large changes in both amplitude and phase of the received Lamb wave signal are expected with variation in the thickness of coating. We have performed a large number of experiments with meander-coil, electromagnetic-acoustic transducers (EMAT's) at various combinations of operating frequency and wavelength and on coated plate samples with different thicknesses of coal-tar coating. Experiments were done in the pitch-catch transmission mode. A detailed comparison of theory with experiment is not possible at this time. We present some sample data. We prepared aflat-plate sample, half of which was coated with a layer of 0.3cm thick coal-tar enamel and the other half remained uncoated. The thickness of steel was 0.954 cm. We performed several experiments on this plate sample, movin~ transducers from the uncoated to the coated area, keeping other operating conditions fixed during each experiment. The results shown in Figures 4 and 5 were taken with EMAT transducers at two separate operating frequencies and wavelengths but with the fixed receiver-transmitter distance of 29 cm. Figure 4a shows the Lamb wave EMAT signal transmitted and received on the uncoated area of the plate. The transducers had peak frequencies at 93 KHz and coii wavelengths at 2 cm. The Lamb wave mode seen is the AO mode. In Figure 4b, we see the same signa! on the coated area. The changes in amplitude (note the difference in the vertical scale in the Figures) and phase of the signal are de ar. The amplitude is reduced by a factor of about 2.3 and the phase is delayed by about 80 degrees in the presence of coating. 1349

10~ ~ ~ to... 5V r- (\ f\ fi f\ fi. f'1 fi h \ Lfl /'\ 1\ 1\ f I I 1 1 ~ \ II ~ ~ I I 1\1 IV - \. r-' V V - - 70.0ua 70.0 ua Fig. 4. Lamb wave signal at a frequency of 93 KHz on the uncoated area of the plate and on the coated area. In Figure 5a, we show the EMAT signal on the uncoated area of the plate when the transmitter frequency is 270 KHz and the coii wavelength is 1 cm. The first burst of signals arrivmg at about 90 microsecond is predominantly the AO mode. The second burst is SO mixed with some higher order modes. In Figure 5b, the signal seen was obtained when we moved the transducers to the coated area. We first notice that the slow SO and higher order modes seen in the uncoated area are almost gone in the coated area. That this is not due to simple attenuation has been verified by decreasing the distance between the transmitter and the receiver. We believe that the generation efficiency ofthese higher order modes at the particular operating frequency and wavelength is much less in the coated plate than in the uncoated plate. Second, the amplitude of the fast AO mode has been greatly reduced in the coating, by a factor of about 8. Third, the phase of the AO mode has been delayed by a large amount. We also measured the amplitude and the group velocity as a function of the thickness of coating. For this purpose, a plate sample was made having eight sections of differing coating thicknesses. The transducers were moved from one section to the next. In each section we monitored the peak aml?litude and the group velocity of the received signal. In Figure 6, we show the vanation in the amplitude L1~V ~ ~ l~v A~ ~ IA h nl~ ~A~ An -I\.t.. '" L,\.A InA~J d... '" V IV 'v IVV I~ IJV ' vv yv~ VII v 'V' IV'V' J" 75.0 us.. -- 75.0 us Fig. 5. Lamb wave signal at a frequency of 270 KHz on the uncoated and coated areas of the plate sample. 1350

of the received signal as a function of coating thickness. The amplitude was normalized with respect to the amplitude in the uncoated area. Figure 6a shows the results for the mode at 100 KHz and 2 cm coii wavelength. In Figure 6b, we see the results for the AO and SO modes at 270 KHz and 1 cm coii wavelength. We find that the amplitude decreases exponentially with increased coating thickness. As a function of the thickness/wavelength ratio, the amplitude of the AO mode drops essentially at the same rate at both the low and the high frequency. The effect on the SO mode is more severe than on the AO mode. In Flgure 7, we show the ratio of the measured group velocity to the shear wave velocity in steel for the AO mode at 270 KHz. We find that the group velocity gradually decreases with increased coating thickness. We have done numerous other experiments at different combinations of operating frequency and wavelength of the transducer. We have found that higher order modes are generally more severely modified in coating than the lowest order AOmode. Amplltude VI. Coatlng Thlcknell 1.2.,.-...:":..:-c.:.:.: ' -:...::'.. ::.:-=AO"-- ----, Amplltude VI. Coatlng Thlckness 1.2,.- " -_.. ;...'-_2_70_-.;...A_OO_50 --. i 1.. s ) u I DAO... sa --- --- 0.2 o."... 0.2 DA O Fig. 6. The variation in the peak amplitude ofthe AO mode at 100 KHz. Group Veloclty VI. Coatlng Thlckness '.02 r--...;".:..-.. ::..;'..-"--=O'-'Ao;..::-=- ----,.... ~......... ~-~~~-----.-~~~--r--,--,...--l o 0.02 0.G4 OM 0.01 0.1 0.12 0,14 CootIno-_1 The same for the AO mode (squares) and SO mode (pluses) at 270 KHz. Fig. 7. The ratio of the measured group velocity to the shear wave speed in steel as a function of coating thickness for the AO mode at 270 KHz. 1351

CONCLUSION The theoretical and experimental results presented above clearly show that effects of coating on Lamb wave amplitude and phase are drastic. Therefore, it may be difficult to use Lamb wave modes for detection and sizing of corrosion-related flaws in pipelines where the variation in the coating condition is large. Additional theoretical and experimental work, especially with simulated corrosion flaws, is needed to make a definitive judgment about the prospect of using Lamb waves in the NDE of buried pipelines. REFERENCES 1. J. P. Jones, in Journal of Appljed Mecbanjcs 31, 213 (1964). 2. R. D. Mindlin, in An Introductjon to tbe Matbematjcal TbeoI}' of Vjbratjous of El astic Plates, prepared for the U. S. Army Signal Corps of Engineering Laboratories, Project 142B (1955). 1352