HOW TO PREDICT THE VIBRATION EFFECT DUE TO MACHINE OPERATION ON ITS SURROUNDINGS

Similar documents
NON-LINEAR DYNAMIC BEHAVIOR OF THIN RECTANGULAR PLATES PARAMETRICALLY EXCITED USING THE ASYMPTOTIC METHOD, PART 2: COMPUTATION OF THE PHASE ANGLE

Keywords: pile, liquefaction, lateral spreading, analysis ABSTRACT

TORSIONAL VIBRATION SUPPRESSION IN AUTOMATIC TRANSMISSION POWERTRAIN USING CENTRIFUGAL PEN- DULUM VIBRATION ABSORBER

POWER DENSITY OPTIMIZATION OF AN ARRAY OF PIEZOELECTRIC HARVESTERS USING A GENETIC ALGORITHM

Methods for detecting fatigue cracks in gears

Comparative study on different walking load models

NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS

Prediction of the Excitation Force Based on the Dynamic Analysis for Flexible Model of a Powertrain

Flexible Pipes in Trenches with Stiff Clay Walls

u y

Participation Factors. However, it does not give the influence of each state on the mode.

REFINED STRAIN ENERGY OF THE SHELL

I Poles & zeros. I First-order systems. I Second-order systems. I E ect of additional poles. I E ect of zeros. I E ect of nonlinearities

Estimating Laterally Loaded Pile Response

KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS

Soil Dynamics Prof. Deepankar Choudhury Department of Civil Engineering Indian Institute of Technology, Bombay

Dispersion relation of surface plasmon wave propagating along a curved metal-dielectric interface

Chapter a. Spring constant, k : The change in the force per unit length change of the spring. b. Coefficient of subgrade reaction, k:

Modelling a Partly Filled Road Tanker during an Emergency Braking

Journal of System Design and Dynamics

Chapter 2 Introductory Concepts of Wave Propagation Analysis in Structures

A PIEZOELECTRIC BERNOULLI-EULER BEAM THEORY CONSIDERING MODERATELY CONDUCTIVE AND INDUCTIVE ELECTRODES

Visualizing Free and Forced Harmonic Oscillations using Maple V R5*

Robust Predictive Control of Input Constraints and Interference Suppression for Semi-Trailer System

Dynamic Analysis on Vibration Isolation of Hypersonic Vehicle Internal Systems

FINITE ELEMENT SIMULATION OF BACKWARD EXTRUSION DIES USING SHRINK FITTED RINGS

The directivity of the forced radiation of sound from panels and openings including the shadow zone

Higher order theory for analytic saddle point approximations to the Ρ Ρ and Ρ Ś reflected arrivals at a solid/solid interface

Uniformly best wavenumber approximations by spatial central difference operators: An initial investigation

Estimation of dynamic behavior and energy efficiency of thrust hybrid bearings with active control

Design of Isolated Bridges from the Viewpoint of Collapse under Extreme Earthquakes

Session 5: Review of Classical Astrodynamics

5.5 The concepts of effective lengths

BENDING INDUCED VERTICAL OSCILLATIONS DURING SEISMIC RESPONSE OF RC BRIDGE PIERS

PRELIMINARY SELECTION OF BASIC PARAMETERS OF DIFFERENT TORSIONAL VIBRATION DAMPERS INTENDED FOR USE IN MEDIUM-SPEED DIESEL ENGINES

Wave Drift Force in a Two-Layer Fluid of Finite Depth

Implementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis

Modified Quasi-Static, Elastic-Plastic Analysis for Blast Walls with Partially Fixed Support

DETERMINING SHEAR ELASTICITY MODULUS BY TORSIONAL PENDULUM

Dynamic Soil Structure Interaction

VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES

Lower bound solutions for bearing capacity of jointed rock

EE 508 Lecture 13. Statistical Characterization of Filter Characteristics

Evaluating Process Capability Indices for some Quality Characteristics of a Manufacturing Process

1 University of Edinburgh, 2 British Geological Survey, 3 China University of Petroleum

Time Domain Calculation of Vortex Induced Vibration of Long-Span Bridges by Using a Reduced-order Modeling Technique

Radar Dish. Armature controlled dc motor. Inside. θ r input. Outside. θ D output. θ m. Gearbox. Control Transmitter. Control. θ D.

A High Precision Double Tubed Hydrostatic Leveling System for Accelerator Alignment Applications

Optical Fibres - Dispersion Part 1

Magnetostrictive Dynamic Vibration Absorber (DVA) for Passive and Active Damping

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)

MODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL

Development of self-adaptively loading for planetary roller traction-drive transmission

Meshless Methods for Scientific Computing Final Project

645. Active control of the process and results of treatment

Optimal Design of Truss Structures Using a Neutrosophic Number Optimization Model under an Indeterminate Environment

Evaluation of straightening capacity of plate roll straightener

351. Massarsch, K.R. and Fellenius, B.H., Engineering assessment of ground vibrations caused by impact pile driving. Geotechnical Engineering

Damage Identification from Power Spectrum Density Transmissibility

CONDUCTION SHAPE FACTOR MODELS FOR 3-D ENCLOSURES

AN EVALUATION OF A SIMPLE DYNAMICAL MODEL FOR IMPACTS BETWEEN RIGID OBJECTS

ON THE PREDICTION OF EXPERIMENTAL RESULTS FROM TWO PILE TESTS UNDER FORCED VIBRATIONS

Pulse Propagation in Optical Fibers using the Moment Method

Optimization of Gear Design and Manufacture. Vilmos SIMON *

Lecture 1.2 Units, Dimensions, Estimations 1. Units To measure a quantity in physics means to compare it with a standard. Since there are many

The effect of dynamic bending moments on the ratchetting behavior of stainless steel pressurized piping elbows

COMPARISON OF FREQUENCY DEPENDENT EQUIVALENT LINEAR ANALYSIS METHODS

Chapter 1 Fundamentals

P043 Anisotropic 2.5D - 3C Finite-difference Modeling

The student will experimentally determine the parameters to represent the behavior of a damped oscillatory system of one degree of freedom.

GROUND VIBRATION FROM RIGID FOUNDATION BY BEM-TLM

Bulletin of the Seismological Society of America

Chapter 6: Sound Wave Equation

Study on a Ship with 6 Degrees of Freedom Inertial Measurement System and Related Technologies

Computational Acoustic Attenuation Performance of Helicoidal Resonators Comparable to Experiment

ABSTRACT. Ossama Salem Ali Ahmed, Ph.D., interaction of the pile with the surrounding soils under vibratory loading will modify

Vibrations of Machine Foundations

Evaluation of Seismic Earth Pressures at the Passive Side

On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm

Response of tall Steel Moment Resisting Frames (SMRF) to near-fault pulse-like ground motion

Study on Characteristics of Sound Absorption of Underwater Visco-elastic Coated Compound Structures

Position Control of Induction Motors by Exact Feedback Linearization *

Prob. 1 SDOF Structure subjected to Ground Shaking

Effect of accumulators on cavitation surge in hydraulic systems

Active Impact Sound Isolation with Floating Floors. Gonçalo Fernandes Lopes

Study on the Electromagnetic Force Affected by Short-Circuit Current in Vertical and Horizontal Arrangement of Busbar System

Use of Transformations and the Repeated Statement in PROC GLM in SAS Ed Stanek

Continuous Steel Casting System Components

Rocking behaviour of a rigid foundation with an arbitrary embedment

Determination of Pile Bearing Capacity By In Situ Tests

VIBRATIONS OF SHALLOW SPHERICAL SHELLS AND GONGS: A COMPARATIVE STUDY

Numerical Modelling of Dynamic Earth Force Transmission to Underground Structures

On Gravity Waves on the Surface of Tangential Discontinuity

MODELING OF THE DYNAMIC RESPONSE OF A FRANCIS TURBINE

Sloshing Effects in Half-Full Horizontal Cylindrical Vessels Under Longitudinal Excitation

MTH234 Chapter 15 - Multiple Integrals Michigan State University

Temperature, current and doping dependence of non-ideality factor for pnp and npn punch-through structures

Design Procedures For Dynamically Loaded Foundations

y p 2 p 1 y p Flexible System p 2 y p2 p1 y p u, y 1 -u, y 2 Component Breakdown z 1 w 1 1 y p 1 P 1 P 2 w 2 z 2

The time and space characteristics of magnetomotive force in the cascaded linear induction motor

Transcription:

ICSV14 Cairns Australia 9-1 July, 007 HOW TO PREDICT THE VIBRATION EFFECT DUE TO MACHINE OPERATION ON ITS SURROUNDINGS Alexandre Augusto Simões 1, Márcio Tadeu de Almeida and Fabiano Ribeiro do Vale Almeida Abstract 1 Deartment of Mechanical Engineering, UNIS Av. Cel. José Alves, 56, Vila Pinto, Brazil Deartment of Mechanical Engineering, Federal University of Itajubá Av. BPS, 1303, Brazil mtaev@uol.com.br The satisfactory erformance of a machine is concerned with its installation, i.e., with a welldesigned foundation and soil investigation. The vibration amlitude of the machinefoundation system must be ket under ermissible values generally secified by the machine s manufacturers and natural frequencies must be determined to avoid resonance. There may be situations when amlitudes of the machine may be within the accetable limits for satisfactory machine oeration but the vibrations may affect other machines around, be harmful to adjacent structures and recision equiment and annoying to ersons. Here in will be shown how to redict in a simle way the effect of the vibration due to the machine oeration on its surroundings considering the foundation as a vibration source on the surface of an elastic medium. 1. INTRODUCTION All machines generated unbalanced dynamic loads (forces and moments) that induce vibrations. These dynamic loads are caused by machine oeration condition such as wear, imbalance, misalignment, etc, and must be ket under a safety limit by the maintenance ersonal for the satisfactory oeration of the entire lant. Another machine s vibration roblem is caused by resonance. The resonance condition is reached when one or more frequencies of the exciting loads are near or equal to any of the natural frequencies of the machine-foundation-soil system. Unfortunately the resonance roblem can t be readily solved because a structural modification will be necessary to bring the vibration level to acceted values.the vibration generated by the machine oeration condition is then transmits through the soil by the block foundation what may affect the vicinity in a bad way. A vibrating footing is then a source of wave generation in the soil. Herein the resonse of a block foundation embedded in the soil and excited by a vertical harmonic load will be outlined and used to redict vibration resonse of the foundation. Also the roblem of the vibration roagation through the soil due to concentrated dynamic load acting on the soil surface will be treated to

ICSV14 9-1 July 007 Cairns Australia estimated the decrease in amlitude of the vibration with distance from the source and to evaluated the effects of the vibration on the erformance of other equiments in the vicinity and structures around. VIBRATION PROPAGATION THROUGH THE SOIL DUE TO PERIODIC CONCENTRATED LOAD The roblem of the surface wave roagation in an elastic medium (elastic half-sace) was first studied by Lord Rayleigh [7] and is known as the Rayleigh wave. Later Lamb [6] studies the resonse of the medium as it was excited by oscillating vertical force or ulse loading acting at a oint or along a line on the surface and within the body. See Figure (1). Lamb s solutions for calculation of vertical surface dislacement at a distance r from the concentrated dynamic load on the surface of the medium may be written as following, Pe iωt w Gr 1 + ( F if ) = (1) in which, P = amlitude of the concentrated oscillating force alied to surface area, ω = circular frequency of force alication, G = shear modulus of the ground, r = radial distance from the load, i = imaginary number, F 1, F = dimensionless Dynamic Boussinesq Dislacement Functions and given by the following equations, Dynamic Concentrated Load alied to soil surface soil Point located at a distance r from the source r Wave roagation through the soil Wave attenuation Figure 1. Dynamic concentrated load alied to soil surface and wave roagation through the soil. s 0.5 a x(s x ) YO (ax) F1 = πykdyo (ay) + π 0.5 0 (x 1) + 4x (s x ) (1 x ) 0.5 x(s x ) K (ax) O + π dx 0.5 0.5 0 (x 1) 4x (s x ) (1 x ) 0.5 dx + 1 0 3 0.5 4x (x s )(1 x ) YO (ax) dx 4 4 (x 1) 16x (x s )(x 1) () s 0.5 1 3 0.5 a x(s x ) J O(ax) 4x (x s )(1 x ) JO(ax) = πykdjo (ay) + + π dx 0.5 0.5 dx 4 4 0 (x 1) + 4x (s x ) (1 x ) 0 (x 1) 16x (x s )(x 1) F (3) in which, a is the dimensionless frequency and is defined by the following equation where ρ is the soil mass density, ρ a = ωr (4) G

ICSV14 9-1 July 007 Cairns Australia ν is the Poisson s ratio of the soil and the s arameter is given by the equation below, 1 ν s =, (5) ( 1 ν) J O, is the Bessel function of the first kind and order zero; Y O, is the Bessel function of second kind and order zero;.5 x 10-6 1.5 vertical dislacement (mm) 1 0.5 0-0.5 0 5 10 15 0 5 30 35 40 45 50 distance from source (m) Figure. Surface dislacement of an elastic medium due to an oscillating concentrated load with amlitude P=1(N) and a circular frequency ω = 6.83 (rad/s) or f=10(hz). Soil data: G= 3.x10 7 (N/m ); mass density ρ=1800 (kg/m 3 ); Poisson s ratio ν=0.3. Rayleigh velocity V R =13.7 (m/s) and wavelength λ=1.4 (m). K O, is the Modified Bessel function of order zero; y is the ositive real root of the equation below and is a function of Poisson s ratio ν only, (y 0.5 0.5 1) 4y (y s ) (y 1) = 0 (6) K d, is a function of the Poisson s ration ν and the value for y of Equation (6) and is exressed by Equation (7). K d 0.5 (y s ) = 4 4y(4y 3y (s + 1) + s ) 8y(y 1) 0.5 0.5 (y s ) (y 1) (7) A lot of the surface dislacement caused by a unit-oscillating load with a circular frequency ω = 6,8 (rad/s) is shown in Figure (). In this figure is ossible to see how the amlitude of the surface wave decreases with the distance from source. The dimensionless functions F 1 and F were numerically evaluated using the traezoidal method. The surface wave, also known as Rayleigh wave, travels at a seed given by the equation below. R V S y V = (8) in which, y is the value obtained using the Equation (6); V S is the shear velocity given by the Equation (9).

ICSV14 9-1 July 007 Cairns Australia F(t)=Pe iωt H Rigid Block Foundation Side Layer G S, ρ S h r O G, ρ Figure 3. Schematic of a block foundation embedded in a soil and subjected to a eriodic force F(t). V = G S (9) ρ 3. VIBRATION OF A RIGID BLOCK FOUNDATION EMBEDDED ON THE SOIL The rigid block that suorts machines may be of various kinds and is usually embedded on the soil as shown in Figure (3). The roblem of a vibrating rigid circular footing on the surface of an elastic medium was examined first by Reissner [10] and after by Sung [1] and Quinlan [9]. The solutions were modified for the urose to obtain a solution for vertical vibrations similar to the equation for damed vibrations of a single degree of freedom system. See Figure (4). The vertical stiffness,, and daming,, for an embedded foundation are given by the following equations below (Prakash and Puri [8], Novak and Beredugo [], Richart et al [11]). G S = GrO C1 + S1 (10) G ro GrO G S = C + S ω (11) G ro Here, C 1 and C are dimensionless arameters derived from the soil reaction at the base of the foundation, and S 1 and S are arameters derived from the reactions acting on the sides of the block foundation. G and ρ are resectively the shear modulus and the mass density of the soil layer under the base of the foundation; G S and ρ S are the shear modulus and the mass density of the layer surrounding the sides of the block foundation as shown in Figure (3). The arameter is the embedment deth. If the footing is circular, r O is the radius of the footing. In the case of a rectangular, square, or other footing shae of area A b the equivalent radius r O is given by Equation (1). The arameters C 1 and C deend on the deth H of the soil, Poisson s ratio ν, and the dimensionless frequency a O given by Equation (13). The arameters S 1 and S deend only on the dimensionless frequency a O. A b r O = (1) π a ωr O O = (13) VS

ICSV14 9-1 July 007 Cairns Australia F(t)=Psin(ωt) Mass M = Mass of the Block foundation lus the mass of suorted equiment. Figure 4. Schematic of a one degree-of-freedom model of a machine-foundation system. The vertical amlitude w O of the steady-state resonse of a single-degree-of-freedom system with resect to the harmonic forcing function F(t) is given by Equation (14). w O P kzs = (14) (1 Ω ) + ( ξω) in which, P is the amlitude of the harmonic forcing function F(t) as shown in Figure (3) and Ω is the frequency ratio defined by the following exression, ω ω Ω = = ω n (15) M where ω n is the natural frequency and M is the total mass, i. e., the mass of the inertial block and the equiment. ξ is the daming ratio defined by the exression below, c ξ = ZS (16) ω n M 4. DEPENDENCE OF THE AMPLITUDE OF SOIL VIBRATION ON DISTANCE FROM THE SOURCE OF WAVES A vibrating machine foundation generates waves in the soil. The energy of the oscillating footing is carried away by a combination of waves: the comression wave, shear wave and the Rayleigh wave. The energy density of each one of these waves decreases with distance from the footing. Figure () shows the decrease in amlitude of the wave with increased distance when a concentrated vertical load is alied at the surface of an elastic medium. The decrease in amlitude, as shown in Figure (), can be exressed by the suggested exression below (Richart et al [11], Barkan [1]). d w = w1 (17) χ where, w is the unknown amlitude at a distance χ from the source; w 1 is the known amlitude at distance d from source;

ICSV14 9-1 July 007 Cairns Australia F(t)=Pe iωt [N] Machine mass Machine mass Z 1 (t) Precision machine or sensitive equiment Z 4 (t) Z (t) k a c a k c Z 3 (t) Figure 5. Schematic of active (left side) and assive (right side) isolation roblem. 5. EFFECT ON SURROUNDING OF WAVES FROM INDUSTRIAL SOURCES If a machine is directly attached to its block foundation, the vibration generated due to the normal oerating condition of the machine is transmitted through the soil and may affect the erformance of other machines in the vicinity or be harmful to adjacent structures even at large distance. On the other hand, it is sometimes ossible to isolate the source of vibration from its surroundings (active isolation), or to isolate a device in such a way that the motion of its base caused by a disturbing source in the vicinity is reduced to accetable limits (assive isolation). See Figure (5).Herein, the roblem of how to redict the effect due to the machine oeration on its vicinity was analyzed considering three cases. First, the roblem was treated without the resence of isolation mounts, i.e., without active and assive isolators. In this case the machine was considered directly attached to the foundation and a single degree of freedom model was used to determine the amlitude of vibration of the machine and its foundation. The roblem was solved alying Equation (14) with Z 1 =Z = w O. See Figure (6) for more details. Equation (17) was modified to take in account the vibration attenuation with distance due to disersion of the wave energy. The calculated vibration amlitude of the machine block foundation was considered as the known vibration amlitude at a distance d = r 0. In this way, the Equation (17) becomes, ro Z3 = w (18) 1 ro + b In this case w 1 =Z 1 =Z and Z 3 =Z 4 due to the absence of assive and active isolators. A distance of b=1.5(m) was considered. See Figure (5) or Figure (6).In the second case only active isolation was considered and a two-degree of freedom system was required to solve this roblem. See Figures (5) and (7). In this second case is necessary to consider w 1 =Z in Equation (18) and as in the first case Z 3 =Z 4 due to the absence of assive isolators. F(t)=Pe iωt [N] b [m] Precision machine or sensitive equiment Z 1 (t)=z (t) Machine mass M hi Z 3 (t)=z 4 (t) b [m] Figure 6. Schematic of the first case treated, i.e., without the resence of active and assive isolation mounts.

ICSV14 9-1 July 007 Cairns Australia F(t)=Pe iωt Z 1 (t) k a c a Figure 7. Schematic of the two degree of freedom of the system showing the machine with its active isolators and the block foundation with its equivalent soil stiffness and daming. Finally both the active and assive isolators were alied. In this last case the transmissibility dislacement ratio must be calculated using Equation (19) below (Inman [5]). Z 4 Z 3 1+ (ξ Ω ) = (1 ) Ω + (ξ Ω In Equation (19), ξ is the daming ratio of the assive isolator and Ω is the frequency ratio of the frequency of the oscillatory base motion ω b to the natural frequency of the assive isolation system ω n. The frequency of the base motion is equal to the frequency of the excitation, ω b =ω. Table 1 resents the data required to the comutation. The results of the comutation of each one of the cases is shown in Table. The force amlitude used to excited the system was P=1500(N) at ω=188.5 (rad/s). The values to the natural frequency and daming ratio of the active and assive isolators were chosen to be equals to 5(Hz) and 0.05(/), resectively. 6. CONCLUSION In the industry the roblem of vibration is always of first concern and must be ket below ermissible values for satisfactory oeration of the entire lant. The goal of this aer was to redict the influence of the vibration in the vicinity due to oeration condition of the machine in a simle way and aid the engineer to get some revious design requirements or remedial measures, in a low cost manner, by roer selection of the block size and contact area, active and assive isolators. The roblem of the vibration roagation through the soil due to concentrated dynamic load acting on the soil surface was treated and the results showed that the vibration attenuation with distance due to disersion of the wave energy may be aroximated by Equation (17) and used to estimated the vibration level that affect structures and delicate machines around as well cause discomfort to workers. ) 1/ Table 1. Soil, Block Size and Mass, Machine data, dynamic soil arameters Shear modulus G=3.14x10 7 (N/m ) Soil Data Mass density ρ=1650 (kg/m 3 ) Poisson s ratio ν=0.5 (/) Lengthxwidthxthickness 1.50mx0.80mx0.60m Block Size Embedded deth =0.5m and its mass Block foundation Mass M B = 178 (kg)* Machine Data Equiment Mass M E = 300 (kg) Seed 1800 rm (30Hz) r O 0.6m Dynamic Soil 1.34x10 8 (N/m) Parameters 7.8x10 5 (Ns/m) * The block foundation mass was calculated using the concrete mass density equal to 400 kg/m 3 Z (t) (19)

ICSV14 9-1 July 007 Cairns Australia Table. Amlitudes Z 1, Z, Z 3, Z 4 of the vibratory movements of the system shown in Fig. 5 to three secific cases. First Case: without the resence of isolation mounts Z 1 =Z =w O 10(µm) Z 3 =Z 4 5.4(µm) Second Case: active isolation Z 1 140(µm) Z 0.3(µm) Z 3 =Z 4 0.16 (µm) Third Case: active and assive isolation Z 1 140(µm) Z 0.3(µm) Z 3 0.16 (µm) Z 4 0.0053 (µm) The vibration roblem of footings resting on or embedded on the soil was treated here considering the solutions given in literature to a circular footing with radius r O. For the case of a square and rectangular contact area is accetable to use solutions of a rigid circular contact area as was done herein for aroximate resonse calculations (Prakash and Puri [8], Novak and Beredugo [], Richart et al [11]). Solutions for irregular shaes of the foundation are available and must be found in the literature (Gazetas [4]). Also Models with several degrees of freedom may be used for redicting the interaction of the machine with its foundation. The interaction changes the rotor-bearing critical seeds redicted by the manufacturer what may causes excessive vibration of the system (Chen [3]). The results obtained and resented in Table shows that when a machine is rigidly bolted to floor or its block foundation, the vibration of the machine itself may be reduced comared to the case when active isolators are used. The use of active isolators reduces the vibration transmitted through the soil what will reduce the bad effects of the vibration on its vicinity. So, when active mounts are rovided the vibration transmitted to the soil will be reduced but this may cause significant vibration to the machine during its oeration or even during its starting and stoing stages. The use of assive isolation hels to reduce the effect of the vibration and is also ossible to rovide isolation using wave barriers or trench barriers (Prakash and Puri [8], Richart et al [11]). REFERENCES [1] D.D. Barkan, Dynamic of Bases and Foundations, McGraw Hill Book Comany, Inc, 196. [] Y.O. Beredugo, M. Novak, Vertical Vibration of Embedded Footings, Journal of the Soil Mechanics and Foundations Division, 197. [3] H.M. Chen, S.B. Malanoski, Fan/Foundation Interaction A Simlified Calculation Procedure, Journal of Engineering for Power, Vol. 103,. 805-810, 1981. [4] G. Gazetas, Foundation Vibrations, Foundation Engineering Handbook, nd Edition, H. Y. Fang Ed., Van Nostrand Reinholds, Chater 15,. 553-593, 1991. [5] D.J. Inman, Engineering Vibration, Prentice-Hall International, Inc, 1994. [6] H. Lamb, On the Proagation of Tremors Over the Surface of An Elastic Solid, Philosohical Transactions of the Royal Soc., London, 1904. [7] Lord Rayleigh, On Waves Proagated Along the Plane Surface of An Elastic Solid, London Mathematical Soc. Proceedings, 1885. [8] S. Prakash, V.K. Puri, Foundation for Machines Analysis and Design, John Wiley & Sons, Inc, 1988. [9] P.M. Quinlan, P. M., The Elastic Theory of Soil Dynamics, Symosium on Dynamic Testing of Soils, ASTM Secial Technical Publication No 156,. 3-34, 1953. [10] E. Reissner, Stationare, Axialsymmetriche durchein Schuttelude Masse erregte Schwingungen eines Homogenen Elastischen Halbraumes, Igenieur-Archiv, Vol. 7, 1936. [11] F.E. Richart, R.D. Woods, J. Hall, Vibrations of Soil and Foundations, Prentice Hall, Inc, 1970. [1] T.Y. Sung, Vibrations in Semi-Infinite Solids due to Periodic Surface Loading, Symosium on Dynamic Testing of Soils, ASTM Secial Technical Publication No 156,. 35-64, 1953.