Table of Contents Author Guidelines Submit a Manuscript
Shock and Vibration
Volume 2014, Article ID 894213, 20 pages
Research Article

Ground Vibration Isolation of Multiple Scattering by Using Rows of Tubular Piles as Barriers

1Zhejiang University City College, Hangzhou 310015, China
2Zhejiang Tongji Vocational College of Science and Technology, Hangzhou 311231, China
3Institute of Geotechnical Engineering, Zhejiang University, Hangzhou 310058, China
4MOE Key Laboratory of Soft Soils and Geo-Environmental Engineering, Zhejiang University, Hangzhou 310058, China
5Zhejiang Institute of Hydraulics and Estuary, Hangzhou 310020, China
6Zhejiang Provincial Key Laboratory of Hydraulics Disaster Prevention and Mitigation, Hangzhou 310020, China

Received 19 June 2013; Accepted 9 March 2014; Published 7 May 2014

Academic Editor: Vadim V. Silberschmidt

Copyright © 2014 Miao-miao Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


A new formal solution for the multiple scattering of plane harmonic waves by a group of arbitrary configuration tubular piles in an elastic total space is derived. Each order of scattering satisfies prescribed boundary conditions at the interface of tubular piles, which is delivered as the sum of incident and scattering waves. The first order performs the scattering wave by each scattered pile and the subsequent orders resulted from the excitation of each pile of first order of scattering from the remaining tubular piles. Advanced scattering orders can be regarded as the same manners. Several series of scattering coefficients are figured out with the aids of addition theorem so that the exact steady-state solution for the scattered displacement and stress is obtained. Particularly, when internal diameter of tubular piles tends to be infinitely small, it degenerates to a solid pile problem. By imposing the normalized displacement amplitudes and transmissibility indices, the influences of specific parameters such as scattering orders, internal and external diameter ratio of piles, pile material rigidity, position and distances between tubular pile and pile rows, and pile numbers are discussed. Certain recommended conclusions have been drawn as the guidelines of practical engineering design for discontinuous barrier of tubular piles.

1. Introduction

The progress of modern society and development brings about a great deal of environmental problems. Traffic, machine operation, and infrastructure construction induced vibration is one of those which is becoming a most talked topic in recent years. Generally, vibration pollutions give rise to not only the uncomfortable residential experiences around neighborhood, but also the overruns of high-tech plants that have strict limit states. In numerous methods of vibration isolation, installing a barrier between vibration sources and the protected equipment or construction appears to be the ideal means. For the past few decades, trenches for screening as research hotspots have already been discussed. Trenches regardless of whether they are in-filled or not serve as a normal perfect way of practical approaches due to their convenience of construction and effective results emerged from the continuous barriers.

However, a series of disadvantages such as high water table levels and field soil instability set limits to the screening in certain circumstances. Thereafter, the alternative method of vibration isolation which is usually considered as the discontinuous barriers by rows of piles is widely adopted in soft soil area. Piles can be driven deeply into soil and arranged in any desirable formation creating wave barriers. Numerous relevant literature researches about using solid piles have been presented during the last several years.

Traditionally from experimental methods, Barkan [1] proposed sheet piles field tests in order to examine the performances creatively. Richart et al. [2] employed piles that served as the isolation barrier soon afterwards. Woods et al. [3] investigated the isolation effectiveness of a row of piles as screening by holography technique. They raised a criterion of designing that the piles’ diameters should be no less than 0.6 times Rayleigh wavelengths and the net separations between piles should be no more than 0.4 times. Haupt [4] carried out a series of model scale tests on the vibration isolation of various measures including concrete walls in a laboratory ground. Afterwards, Liao and Sangrey [5] inspected the isolating effectiveness by model piles when sound waves propagated in fluid media. They concluded a point of view that flexible piles played better roles than stiff ones. Boroomand and Kaynia [6] presented an approximate analytical foundation model to examine the efficiency of rigid and elastic piles as isolator in half space. Kani and Hayakawa [7] put use to field tests on PC-wall piles to evaluate the reduction of ground vibration.

Unlike the classic empirical approaches, numerical methods were widely utilized recently. Especially the boundary element method seemed to show a good strategy for dynamics response. More detailed descriptions on the barrier systems may be found. Kattis et al. [8, 9] put forward a 3D frequency domain BEM analyzing isolated traits of a row of piles. Several patterns of piles which included tubular, solid with circular and square cross-section have been taken as the research objects. The tubular piles were regarded as long cylindrical cavities.

Most of the existing theoretical researches concerning only one row of piles subjected to the incident waves were still not consummate. Avilés and Sánchez-Sesma [10, 11] analyzed the properties of one row solid piles acted as the isolator of SH, SV, and P waves in an elastic medium. The numerical investigation showed that rigid piles were superior to the flexible ones. The most valid diameter of piles was between 0.25 and 1 wave shear wavelengths. Later on, they developed a theoretical model to study the amplitude reduction behind one row of cylindrical solid piles when Rayleigh waves are incident. They concluded that the optimum efficiency would be gotten when the diameter of piles equaled 0.25 wavelengths and the pile length equaled 2 wavelengths. Gao et al.  [12, 13] employed the half-space Green’s function to inspect 3D multirows of square piles as barrier in a far field. Xu et al. [14, 15] initiated single-row tubular piles as screening barrier utilizing wave function expansion method under the circumstance of P and SV waves’ incidence. Cai et al. [1618] studied the wave scattering by single-row piles in poroelastic soil with the aids of Fourier-Bessel series expansion. Tsai et al. [19] applied a 3D BEM in the frequency domain to investigate the screening effectiveness of circular piles in a row for a massless square foundation subjected to harmonic vertical loading. Xia et al. [20] delivered the analytical solutions of an arbitrary configuration of rigid solid piles as wave scatterers based on the solutions of acoustic and electromagnetic waves scattering by immovable rigid cylinders.

However, in the circumstance that stiffness of piles is much smaller than soil, namely, elastic piles case, the screening effectiveness is often more deteriorated than rigid ones; thus tubular piles are substitutes for solid piles as discontinuous barrier. Not only are they a cost-saving and convenient constructed solution of materials compared to rigid ones, but also they cut off the medium of wave propagation derived from their hollowed formations, even rivaling the continuous barrier.

Some theories of acoustic and electromagnetic waves scattered by shells in elastic medium as the basis of tubular piles issue have been accumulated. More recently, Huang and Gaunaurd  [21, 22] denoted a solution of scattering of coherent acoustic waves by two spherical shells. Robert et al. [23] proposed a way of solving the scattering problems by finite cylindrical cavities. Le Bas et al. [24]studied the multiple scattering and fission phenomena by multishell in water.

Nevertheless, to date, the theory of multiple scattering by tubular piles as discontinuous barrier is still immature and deficient. Assumption of scattering orders has remained to be discussed as single from the literatures mentioned above, in which the scatterers are far apart so that their excitation can be taken as the incident waves [25]. However, separations between piles and pile-rows are small enough to guarantee the overall effect actually. In addition, little attention has been paid to analyze tubular piles in elastic medium in which the soil is treated as an isotropic homogeneous solid. Other methodologies have been most adopted by frequency domain boundary element, let alone considered of more general arrangement of piles, most of whose cross-section configurations trends to be identical and the arrangement is linear.

This paper presents a creative method of calculating multiple scattering by an arbitrary configuration and diameter of tubular piles as discontinuous screen-scatters. The response of tubular piles to this excitation is the so-called first-order scattering. In calculating the second-order scattering from one of the scatterers, the first-order radiation of all the others is taken to be the secondary excitation, and the fields of all these remaining tubular piles are translated to the origin of that particular pile from which the second-order scattering is to be calculated. Following this procedure for all the objects present, a recursion relation is obtained for higher order scattered fields in terms of the scattered field of the single, component object, and thus successively higher order corrections are made to the single scattering solution for the configuration. With the aids of addition theorems, boundary conditions both in and out tubular piles have been implemented for each pile and scattering fields of the remaining piles can be expressed with respect to the origin of piles for which the boundary conditions applied. Numerical analyses have been taken to confirm the multiple scattering calculating valid and reasonable for engineering applications.

2. The Geometry and Theoretical Model

The geometric position of the scatterers chosen for the theoretical analysis of the problem is shown in Figure 1; the cylindrical tubular piles are located in an isotropic, homogeneous, elastic infinite medium, of which and are Lamé constants and is matrix mass density. The cylinders extend to infinity with zero displacement along the pile shaft so that the problem can be simplified as a two-dimensional, plane-strain problem.

Figure 1: Geometric model and analysis coordinate systems.

The reference coordinate system begins with the origin 0 and is set up as a rectangular Cartesian coordinate or a polar coordinate . The first tubular pile above the-axis is defined as number 0; the next ones are identified by the numbers , ,  , and so forth. The piles are with the external radii , , and and internal radii , , and , respectively. The point stands for an arbitrary position with respect to the origin. The distances between point and pile-centers are , , and so forth; the angles independent of each cylinder are , , and so forth; their center-to-center spaces are , , and so forth.

The multiple scattering procedure can be described as follows. When a series of steady-state harmonic waves propagate from infinite distance with incident angle , the first order of scattering is regarded as the most common “single scattering” hypothesis in which all the scatterings are calculated by translating the primary incident plane waves to the origin of each pile and then summing all the responses of the first excitations. On the basis of the first scattering, second order of scattering is excited by the first scattering waves which perform as the “incident” waves, and the fields of all these remaining piles are translated to the origin of that specific pile. It gives rise to a subsequent recursion relationship for successive orders so that the multiple scattering physical processes made sense.

Equation of motion for the displacement vector at the point is expressed as [26] where is the gradient operator with , of polar coordinate or of rectangular Cartesian coordinate system.

The displacement vector could be attained by introducing a scalar and a vector potential , , separately, where and are functions of or with time dependence exp().

Substituting to the wave equation, (2) turns into

Those are known as the Helmholtz equation [27], in which , , , and are used to denote wave velocities of P and S waves of the matrix. is the so-called circular frequency of harmonic waves, with the corresponding wave numbers and , respectively.

2.1. SH Wave Incidence

In the propagation of SH wave (the direction vertically polarized), the displacement potential accords with the Helmholtz equations mentioned above: Subscript of means the wave numbers of soil.

The incident SH wave displacement potential of unit magnitude is taken as where inc indicates incident wave, is a scalar coefficient, , and = time, with being suppressed. The vector gives the direction of incidence and is the position vector of far field point with respect to origin 0.

SH wave’s potentials due to the th order scattering from th pile are defined as : Thus the total scattering wave field from all piles with number is constituted as By applying the separation variables method, the scattered field under the polar reference system could be delivered as where the superscript sc denotes scattering waves, is the first kind of Hankel function with order , and is the scattering undetermined coefficient of th pile varying with the multiple order . The total waves become

The tubular elastic piles problem differs from the solid or rigid piles. Therefore the solutions could be readily extended to more general boundary conditions, and expressions could also be obtained for the field inside any pile when it exists.

In consideration of interior transmitted waves which consist of two parts of waves, external and internal propagations, the “transmitted” wave field in the interior of the th pile is expressed as [28] where and are undetermined transmitted coefficients for the th pile subject to continuity boundary conditions and is defined as the second kind of Hankel function with order . The subscripts of refer to pile and the rest of the variables refer to the same meaning as mentioned above.

As known to all, incident SH wave function series expansion equation is In order to transform generalized wave functions to a specific pile’s coordinate , , through employing the boundary conditions, incident SH wave is expanded to a Fourier-Bessel series in cylindrical wave functions [28]: where , in which the time factor and displacement amplitudes in common are both omitted.

The expansion equation of scattering wave function is presented in the same way as listed in (8).

For the sake of solving the undetermined complex coefficients, firstly, It is assumed that the interfaces of soil and piles complied with elastic boundary conditions, which is considered the displacement and stress of interface are continuous and the radial stress of inner boundary is free:where , and defines the shear modulus of soil and pile with the subscripts and , separately.

Substitute (8), (10), and (12) into (13a), (13b), and (13c); the first order of scattering and transmitted coefficients is obtained:where we set The subscripts and of mass density refer to soil and pile.

With regard to the expression of , it is set that when internal radius , so that , ; meanwhile ; it is inferred that . As a result, degenerates to , which is the same solution form of fixed rigid solid pile [20].

It is demonstrated that the single scattering coefficients are more general. On the other hand, when , thereby Hankel functions and ; (15) turns into

Obviously, it degrades a cylindrical solid elastic pile problem [29]; the solutions of tubular pile problem are verified to be correct.

Since the th order of scattering from th pile is also the excitation for the th order scattering from the th, given above with respect to origin , the th order have to be translated to the origin . This translation of the origins for the functions can be accomplished by Graf’s addition theorem [30] for the cylindrical wave functions:

It is taken no account of transmitted waves’ superposition as secondary wave source for the result that inner transmitted waves are influenced a little as the excitation of other piles’ transmitted waves. Therefore the multiple scattering () boundary conditions become where ,  and .

After a substitution course, the th order undetermined complex coefficients are given aswhere we set

It is noted that the operator is expressed in the same form as . When shear modulus , the operator of th order scattering coefficient ; meanwhile internal radius and Hankel functions and . The same degeneration process in which it transforms to fixed rigid solid and elastic solid piles’ scattering complex coefficients is received. The theoretical solution of total scattering and transmitted waves could be determined finally.

2.2. P and SV Waves Incidence

A coupling scattering phenomenon under P or SV waves incidence occurs in front of the barrier, and it contains both P and SV waves. There are of no difference scattering waves functions form between solid pile and tubular pile. However, transmitted waves include both dispersed outward and converged inward stationary waves, which are expressed as expansions ,, , and , where , ,, and are transmitted coefficients of th order and , are defined as wave numbers of P and S waves.

Thus the wave function expansions of incident, scattering, and transmitted waves are

The radial and tangential displacements are decomposed:

The stress tensors of P and SV waves are expressed as where is unit tensor.

The radial and tangential stress of stress components involved arewhere the subscripts of all variables refer to the th pile.

Assuming as the first scattering, the boundary of soil-pile interface at satisfies the elastic conditions:where   and.

Substitute (22)–(24b) and (27a) to (28f); the first scattering coefficients and and transmitted coefficients , , , and meet the relations below: where stands for matrix, is the scattering or transmitted coefficients matrix, and presents a vector like thiswhere elements and   indicate ; ; and . All the elements in , , and are listed in the Appendix.

Matrices can be solved by Cramer rule only with substituting vector into every column of and renewing them so that the first scattering and transmitted coefficients are , , , , , and .

The same degeneration process will be prone to be received when and ; transmitted coefficients  , , , and are equal to 0, which represents that there is no transmitted field when screen consists of rigid solid piles [20]; on the other side, when , it is equivalent to a cavity problem.

Similarly, the elastic boundary conditions of th order are

where ,,   and.

Thus the formula of th order () and ()th coefficients is displayed: where the implications of and are no different from former ones, defines a realized th coefficients matrix, and is vector. Expressions yieldwhere the subscripts of elements are ; and ,  . The rest of the elements are attached in the Appendix. The degeneration process is omitted for abbreviation.

For incidence of SV waves, the only distinction is to substitute incident P wave for SV wave functions (12) at the right-hand side.

The first coefficients expressions of incident SV waves meetwhere detailed expressions also can be got afterwards and the th order solutions are no longer different from P waves.

3. Numerical Analysis on Screening Effectiveness by Rows of Tubular Piles

In this study, the surrounding soil is considered as a homogeneous, isotropic medium, the piles’ lengths are assumed infinite, and herein this issue goes to a total-space problem. Assuming a series of time-harmonic plane waves with incident angle which propagates perpendicularly to this pile group, a numerical calculation model of one row parallel tubular piles with circular cross-section is depicted in Figure 2(a), a rectangular coordinate system is placed from the left, the length of barrier system is , and every tubular pile has identical inner diameter and external diameter , while the separations from center to center of piles are equally spaced with . Some other notations are given as follows:shear modulus ratio ( and are shear modulus of pile and soil, separately),mass density ratio .

Figure 2: The reference system and analysis model of incident plane waves scattering by rows of piles .

Frequencies of incident P and S wave are normalized aswhere the first letter in subscripts of frequency or wave length   represents incident P or S waves. The second letter in subscripts and the ones in wavenumbers refer to soil or pile. Wavenumbers in different medium meet the relations below:where and are Poisson ratios of soil and tubular pile.

Other two-row piles with hexagon appearance are also shown in Figure 2(b), with the pile-row distance . To testify the correctness of multirow applications, one can take the single-row as a special situation of Figure 2(b) while the space between rows of piles    degenerates to 0. Computations over a wide range of complex arguments require sufficient accuracy only if the Bessel function expansion order truncated at no less than 8 [29]. Besides, normalized displacement amplitudes , , and are employed in numerical calculation, so as all other geometrical parameters.

3.1. Influences of Scattering Orders

Figure 3 illustrates several curves of displacement amplitudes versus distance under different scattering orders. varies from 1 to 4. From these figures, it is observed that the same amplitude amplified phenomenon emerges at the front of isolation barrier. The numerical results indicate that a much smaller increment comes up when orders of scattering increase. However, the disparity between single and multiple orders cannot be neglected for the reason that there is coherent phase interference in the successive orders. will be used in the numerical analysis owing to these two curves of and are almost superposed with the incidence of both S and P. It might be attributed to the presence of the dissipation effect. This result, which is consistent with the result of Sun [29] using the same multiple scattering hypothesis in rigid solid pile, is a convincing fact of this research method. It is interpreted that when the elastic wave propagates through the pile barrier, there exists an interference effect due to the wave scattering. The scattering waves are stronger in the center of the pile barrier than at the edge, where more wave energy is transmitted to the screened zone because of a finite size of the barrier. These infinite series’ convergence resulted from the scattered waves’ diminishing along with the higher order of scattering. A comparison between Figures 3(a) and 3(b) indicates that gapping away each tubular pile from front to back properly causes a prominent improvement of displacement amplitudes reduction along the midline of barrier. In other words, the double-row system decreases no less than 2 times incident waves’ amplitude at the range of under the same number of piles as single-row. A further detailed parametric study has been taken as follows.

Figure 3: Normalized displacement amplitudes curves along midline behind the isolation system under different scattering orders by rows of tubular piles ,  ,   ,   .
3.2. Influences of Different Position behind the Barrier

In order to investigate the displacement representation of different spot at the back of pile row, the variations of and with varied positions versus -axis by single or double-rows of tubular piles are presented in Figure 4. It is demonstrated that the closer distance it is away from the midline of barrier, the superior responses of amplitudes reduction it will get. In the examination of single elastic row tubular piles in Figures 4(a) and 4(c) under the incident SH and P waves, similar characters versus are received. In the center, it reacts 1/3 displacement amplitudes of the barrier edge within the bounds of . However, delicate distinction can be differentiated in various positions for due to oscillation verging to steady for distance. Either the rise of shear modulus ratio or pile row distances (a process of tubular piles degrading to rigid solid piles) achieves the remarkable decrease of amplitude in Figures 4(b) and 4(d). The best isolated zone moves back to . These appearances agree with subsequent parametric studies.

Figure 4: Normalized displacement amplitudes at different positions behind the barrier , , .
3.3. Influences of Tubular Pile Wall Thickness Ratio

The tubular pile wall thickness ratio is a unique characteristic to investigate screening effect of the barrier, which has not been taken into account in previous studies. Figure 5 depicts the curves of the midline amplitudes reduction along with the variation of tubular piles’ wall thickness ratio . It can be observed that the maximum reduction occurs in the range behind the isolation system. The displacement amplitudes decrease to 0.3 in Figure 5(a) due to the wall thickness attenuating to 1/10 () of external diameter, which amounts to 3/5 displacement amplitudes of solid piles barrier . The parallel advantages can be perceived in Figures 5(c) and 5(d) at the range backing to when it is of incident P waves. For double-row piles, evident decrease of displacement is obtained while the thickness wall reduces. Besides, comparing Figure 5(a) with Figure 5(b) or Figure 5(c) with Figure 5(d), it appears that double-row reflects better than single-row. With the same material utilization amount, tubular piles present to have more efficiently isolated achievements than traditional solid piles do. In this case, tubular piles save much cost and have high application value in practical engineering. Thus ratio is taken for analyses below.

Figure 5: Normalized displacement amplitudes curves along midline behind the isolation system under different ratios of internal and external diameters of tubular piles , , , .
3.4. Influences of Shear Modulus Ratio

Different from the rigid solid piles, the tubular piles have another significant factor called shear modulus ratio which evaluates the flexibility of pile body materials and soil properties. Figure 6 presents a series of contours following the changed from 10 to infinite when it is under both SH and P waves incidence. It is noted that all optimum screening regions are at the center of the isolated scope. Those figures reveal that the minimum values of the normalized displacement amplitudes decrease with the increasing of shear modulus ratio . As varied to 100 when it is of SH waves incidence, the best performance area draws back to other than when equals 10. The isolation effectiveness performs to reduce 80% displacement amplitudes meanwhile. When grows from 1000 to +∞, it can be observed that further increasing of shear modulus ratio appears to have little impact on the screening effectiveness. The shapes of the contour containing are approximately the same. That is to say, it almost equates with the solid piles under the circumstance that is greater than 1000. The degeneration characteristics are proved that the numerical algorithm is valid. In other words, elastic tubular piles are far better formation of discontinuous barrier than rigid solid piles. The same tendency can be received in Figures 6(e)6(h) when it is of P waves incidence.

Figure 6: Contours of normalized displacement amplitudes by double-row tubular piles , , , , .
3.5. Influences of Separations between Each Tubular Pile

Figure 7 depicts the variation of spacing between the tubular piles behind the isolation system. Generally, for all the curves at the distances from the pile unity, they reveal the consistent variation trends in which the closer it gets near the core behind the barrier, the superior reduction of displacement amplitudes it manifests. The normalized value of separation varies from 3.0 to 5.5, respectively. One can observe that considerable decrease appears with the separation rises. With 3 times of pile diameter, Figure 7 suggests that the value of or acquires no less than 50% incident wave is isolated compared to the pile-row with the normalized separation of 5.0. The rational explanation is that the vibration isolation barrier behaves as an integrated obstacle when the separation is within a certain limit. Along with the separations from pile to pile setting farther away, the pile group displays nearly a series of independent piles instead of a whole unit. However, the effectiveness verges to convergence at the fringe of barrier. As expected, it is noticed that a far field point at -axis except varying puts up better performance at nearer limits than a remote distance , whether SH or P waves incidence that is also a conformance resulted from Figure 3.   = 3.0 is recommended as the designing criterion in practical engineering to ensure distinguished screening effects.

Figure 7: Normalized displacement amplitudes curves with different separations of tubular piles .

Another significant factor to measure the effectivity of barrier as isolation system is the so-called transmissibility index (TI) [10]: where and are defined as the total wave field, and those with subscripts 0 are demonstrated incident wave as well. Therefore, TI varies between 0 and 1. The transmissibility indices with different positions distributed in the rear of barrier are illustrated by dash, dot dash, and solid line. The -axis indicates that separations of tubular pile vary in the extent of . As expected initially, despite carrying on rigid tubular piles , which are even equivalent with rigid solid piles, or elastic piles , TI keeps a relatively lower range when pile spacing locates in certain limits . It is a proof that discontinuous barrier plays a role of integrity other than those expressed as disperse individuals. However, on the other side of view, Figures 8(c) and 8(d) display a more preferable result in rigid tubular pile than under the circumstances of elastic ones in Figures 8(e) and 8(f), which is an affirmation with Section 3.4. Neither nearer nor farther confines could bring about an ideal screening effectiveness only when their positions of far field points are within . Double-row tubular piles (Figures 8(c) and 8(d), ) evidently reveal outstripping isolating characteristics other than single-row (Figures 8(a) and 8(b), ). Consequently, it is another vigorous identification that can be utilized as the most optimum spacing to ensure the desired effect of pile rows.

Figure 8: Transmissibility index curves with different separations of tubular piles ,   .
3.6. Influences of Distances between Tubular Pile Rows

The influence of multirow separation distance between tubular pile rows is investigated in Figure 9. Those figures of meshes represent intuitionistic isolation characters where the back of the discontinuous barrier is located. Aside from ubiquitous amplification of displacement amplitudes near the piles, one can perceive a sharp decrease presented in a certain range, almost shaping like low-lying hollows concentrated at the front. In the case of in Figure 9(a), the screening effectiveness performs as a single-row which keeps at an average of 0.7~0.8. With the larger separations between pile rows in Figure 9(c) it takes, the reduces to 0.4~0.6. It is a valid fact that suggests retaining particular spacing distances of tubular pile rows, a much more remarkable effect it will get without increasing the number of piles. Similar meshes appearance is shown in Figures 9(b) and 9(d) with under the incidence of P waves. grows from 0.2~0.4 to 0.6~0.7 accompanied with rises. The increasing gap results in interaction phenomenon weakened so that the pile group no longer serves as a whole entirety. That is to say, not only pile separation but also row distance is the significant factor to form thickness of the obstacle to recede both incident SH and P waves. More prospective effectiveness on condition that multirow separation distance is will be guaranteed.

Figure 9: Normalized displacement amplitudes meshes behind the isolation system under different row-separations of tubular piles .
3.7. Influences of Tubular Pile Numbers

The numbers of piles varied from 6 to 12 along the midline behind pile rows are illustrated in Figure 10. In general, displacement amplitudes along the midline go to smaller when the numbers of piles increase. and damp at a more distant bound with more tubular piles. To some extent, a 12-pile barrier screens off 70–80% incident P waves while 6-pile barrier screens no more than 50–60% of it. Aside from the phenomena observed, an approximate result has been obtained in Figure 10(b), in which tending to a certain large value make the tubular pile coincide with rigid solid pile case. This degradation suggests a supportive fact of analytic solution of tubular pile from another point of view.

Figure 10: Normalized displacement amplitudes curves along midline behind the isolation system under different numbers of tubular piles , ,   ,   .

Numbers of rows, namely, the thickness of barrier, are an important factor affecting isolation of tubular piles. A treble-row shaped hexagon is illustrated in Figure 11, which displayed 12 piles versus 8 ones before. Comparing Figure 6(a) with Figure 12(a), it is observed that with growing of 50% pile numbers, which is double the thickness of screen, the minimum value of treble-row piles reduces to 0.1. Incident SH waves could been isolated no less than 40% at most area behind barrier. The minimum value of double-row is 0.5 at certain distance comparatively, and most of the region behind the screen incident waves is insulated only 30% at most region behind the screen. Analogous appearances can be received from Figure 6(e) and Figure 12(b) when it is of incident P waves.

Figure 11: The reference system and analysis model of incident plane waves scattering by multirows of piles (hexagon arrangement, ).
Figure 12: Contours of normalized displacement amplitudes by multirows of tubular piles ,   ,   ,   .
3.8. Influences of Configuration of Tubular Piles

This paper presents a methodology of calculating scattering by arbitrarily arranged tubular piles. Different settings lead to diverse performances of isolation effect. Figure 13 shows the rectangular configuration reference system and Figure 14 depicts contours of normalized displacement amplitudes by incident SH and P waves. With the same numbers of tubular piles , it can be perceived that no less than 40% of incident SH waves are screened when the barrier is hexagon shaped versus 20% of the incident waves isolated when the barrier is rectangular arranged from Figures 6(a) and 14(a), which increases double efficiency of isolation while it is hexagon. An obvious contrast can be discovered in Figures 14(b) and 6(e) under the incidence of P waves. Therefore, rectangular shaped system actually broadens the separations between the tubular piles with the same width of barrier so that less interaction and interference will occur. Reasonable arrangement is acceptable to meet the requirement of vibration attenuation in practical engineering.

Figure 13: The reference system and analysis model of incident plane waves scattering by multirows of tubular piles (rectangular arrangement, ).
Figure 14: Contours of normalized displacement amplitudes by double-rows of tubular piles (rectangular arrangement, , , , , , ).

4. Conclusions

A multiple scattering method different from the single scattering of calculating elastic waves scattering by an arbitrary configuration and radius of parallel cylindrical tubular piles is proposed. By means of constructing transmitted waves field inside the tubular piles, applying elastic pile-soil interface boundary conditions, and employing wave function expansion and coordinates conversion theory, it acquires scattering and transmitted complex coefficients of each order. Total scattering and transmitted wave fields are determined consequently. Example analysis of single and multirow tubular piles, which are focused on scattering orders, types of incident waves, materials of screen, and geometrical parameters have been discussed.(1)An amplification phenomenon of displacement amplitudes behind the tubular pile screen stimulated by the incident plane waves is conspicuous. Displacement amplitudes increment appears to be gradually decreasing with the orders of scattering increasing and is nearly superposed far apart the barrier especially. It is concerned that displacement amplitudes are converged when truncates at 3 to 4. No matter what kind of incident waves, the computation of multiple scattering method is much smaller than the way based on single scattering hypothesis. It is presented to be the basis of pile barrier optimal design in practical engineering.(2)The most prominent amplification of displacement is discovered centered behind the pile barrier. Meanwhile, the minimum values of that appear exactly behind the amplification zone and the amplitudes of two sides are apparently varying small. It indicates that a similar manifestation like continuous barrier characteristics is achieved by tubular piles.(3)Internal and external radius ratio or called thickness of tube wall is a significant influencing factor of screening effectiveness which is distinct from solid piles. Whether single or multiple rows of them, a reduction of thickness makes displacement amplitudes weakened remarkably. When internal radius approaches to 0, that is, degenerating to elastic solid piles, the decrement of displacement is less than that of tubular piles of the same external radius. As a result, tubular pile is an ideal shape of much beneficial and practical substitute of discontinuous barrier for solid ones. is proved to be the optimum ratio to acquire better isolation result.(4)The shear modulus ratio of materials of piles and soil, or say the stiffness of barrier, influences the function of screening directly. An obvious magnification emerges with the pile-soil shear modulus ratio increasing; the optimum isolation region falls back and the minimum displacement amplitudes decrease so that the isolation effectiveness improves evidently. In addition, tubular pile groups with are regarded as fixed rigid piles for the normalized displacement amplitudes contour is tantamount to that of .(5)Variation of separation between piles has much impact on screening of incident waves. It is more sensitive to displacement increment at central area than at two sides behind the barrier. TI which is called transmitted index is introduced to evaluate the isolation efficiency. The transmitted waves become less when the separations of single-row piles and double-row reduce. The best separation of tubular piles tends to be .(6)Row-to-row distance needs to keep certain range because pile rows are isolated as discrete barrier gradually when the distance increases. With the same numbers of tubular piles, double-row performs more superiorly than single-row tubular piles. It works sufficiently when distance is about .(7)Increasing either of width or thickness of barrier gives rise to improvement of the effectiveness of screening. The rectangle configuration widens the separation of tubular piles and reduces the chances of multiple scattering coherent phases’ intervention. Compared with rectangle arrangement, hexagon configuration with the same numbers of tubular piles is preferred.


Elements of matrices are where in which the subscripts of and can be or .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The authors are grateful to the anonymous reviewers for their helpful advice and valuable comments. This work was supported by the National Natural Science Foundation of China under Grant (no. 51309201), Public Welfare Scientific Research Foundation of the Ministry of Water Resource of China (no. 201201043), and Post-doctoral Merited Research Foundation in Zhejiang Province of China (no. Bsh1202040).


  1. D. D. Barkan, Dynamics of Bases and Foundation, McGraw-Hill, New York, NY, USA, 1962.
  2. F. E. Richart, J. R. Hall, and R. D. Woods, Vibrations of Soils and Foundations, Prentice-Hall, Englewood Cliffs, NJ, USA, 1970.
  3. R. D. Woods, N. E. Barnett, and R. Sagesser, “Holography, a new tool for soil dynamics,” Journal of Geotechnical Engineering, vol. 100, no. 11, pp. 1231–1247, 1974. View at Google Scholar
  4. W. A. Haupt, “Model test on screening of surface waves,” in Proceedings of the 10th International Conference on Soil Mechanics and Foundation Engineering, vol. 3, pp. 215–222, Stockholm, Sweden, 1981.
  5. S. Liao and D. A. Sangrey, “Use of piles as isolation barriers,” Journal of the Geotechnical Engineering Division, vol. 104, no. 9, pp. 1139–1152, 1978. View at Google Scholar · View at Scopus
  6. B. Boroomand and A. M. Kaynia, “Stiffness and damping of closely spaced pile groups,” in Proceedings of the 1st International Conference on Soil Dynamics and Earthquake Engineering, pp. 490–501, Karlsruhe, Germany, September 1991. View at Scopus
  7. Y. Kani and K. Hayakawa, “Simulation analysis about effects of a PC wall-pile barrier on reducing ground vibration,” in Proceedings of the 13th International Offshore and Polar Engineering Conference, pp. 677–682, Honolulu, Hawaii, USA, 2003.
  8. S. E. Kattis, D. Polyzos, and D. E. Beskos, “Vibration isolation by a row of piles using a 3-D frequency domain BEM,” International Journal for Numerical Methods in Engineering, vol. 46, no. 5, pp. 713–728, 1999. View at Google Scholar · View at Scopus
  9. S. E. Kattis, D. Polyzos, and D. E. Beskos, “Modelling of pile wave barriers by effective trenches and their screening effectiveness,” Soil Dynamics and Earthquake Engineering, vol. 18, no. 1, pp. 1–10, 1999. View at Publisher · View at Google Scholar · View at Scopus
  10. J. Aviles and F. J. Sanchez-Sesma, “Piles as barriers for elastic waves,” Journal of Geotechnical Engineering, vol. 109, no. 9, pp. 1133–1146, 1983. View at Google Scholar · View at Scopus
  11. J. Avilés and F. J. Sánchez-Sesma, “Foundation isolation from vibration using piles as barriers,” Journal of Engineering Mechanics, vol. 114, no. 11, pp. 1854–1870, 1988. View at Google Scholar
  12. G. Y. Gao, Z. Y. Li, C. Qiu, and Z. Q. Yue, “Three-dimensional analysis of rows of piles as passive barriers for ground vibration isolation,” Soil Dynamics and Earthquake Engineering, vol. 26, no. 11, pp. 1015–1027, 2006. View at Publisher · View at Google Scholar · View at Scopus
  13. Z.-Y. Li, G.-Y. Gao, C. Qiu, and X.-J. Yang, “Analysis of multi-row of piles as barriers for isolating vibration in far field,” Chinese Journal of Rock Mechanics and Engineering, vol. 24, no. 21, pp. 3990–3995, 2005. View at Google Scholar · View at Scopus
  14. P. Xu, T.-D. Xia, and X.-M. Zhou, “Study on effect of barrier of a row of hollow pipe piles on isolation of incident plane SV waves,” Chinese Journal of Geotechnical Engineering, vol. 29, no. 1, pp. 131–136, 2007. View at Google Scholar · View at Scopus
  15. P. Xu, X.-M. Zhou, and T.-D. Xia, “Discontinuous barrier used a row of elastic piles for incident elastic waves,” Journal of Vibration Engineering, vol. 20, no. 4, pp. 388–395, 2007. View at Google Scholar · View at Scopus
  16. Y.-Q. Cai, G.-Y. Ding, and C.-J. Xu, “Screening of plane S waves by an array of rigid piles in poroelastic soil,” Journal of Zhejiang University: Science A, vol. 9, no. 5, pp. 589–599, 2008. View at Publisher · View at Google Scholar · View at Scopus
  17. Y.-Q. Cai, G.-Y. Ding, and C.-J. Xu, “Amplitude reduction of elastic waves by a row of piles in poroelastic soil,” Computers and Geotechnics, vol. 36, no. 3, pp. 463–473, 2009. View at Publisher · View at Google Scholar · View at Scopus
  18. Y.-Q. Cai, G.-Y. Ding, C.-J. Xu, and J. Wang, “Vertical amplitude reduction of Rayleigh waves by a row of piles in a poroelastic half-space,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 33, no. 16, pp. 1799–1821, 2009. View at Publisher · View at Google Scholar · View at Scopus
  19. P.-H. Tsai, Z.-Y. Feng, and T.-L. Jen, “Three-dimensional analysis of the screening effectiveness of hollow pile barriers for foundation-induced vertical vibration,” Computers and Geotechnics, vol. 35, no. 3, pp. 489–499, 2008. View at Publisher · View at Google Scholar · View at Scopus
  20. T.-D. Xia, M.-M. Sun, C. Chen, W.-Y. Chen, and X. Ping, “Analysis on multiple scattering by an arbitrary configuration of piles as barriers for vibration isolation,” Soil Dynamics and Earthquake Engineering, vol. 31, no. 3, pp. 535–545, 2011. View at Publisher · View at Google Scholar · View at Scopus
  21. H. Huang and G. C. Gaunaurd, “Acoustic scattering of a plane wave by two spherical elastic shells,” The Journal of the Acoustical Society of America, vol. 98, no. 4, pp. 2149–2156, 1995. View at Publisher · View at Google Scholar · View at Scopus
  22. H. Huang and G. C. Gaunaurd, “Acoustic scattering of a plane wave by two spherical elastic shells above the coincidence frequency,” The Journal of the Acoustical Society of America, vol. 101, no. 5, pp. 2659–2668, 1997. View at Publisher · View at Google Scholar · View at Scopus
  23. S. Robert, J.-M. Conoir, H. Franklin, and F. Luppé, “Resonant elastic scattering by a finite number of cylindrical cavities in an elastic matrix,” Wave Motion, vol. 40, no. 3, pp. 225–239, 2004. View at Publisher · View at Google Scholar · View at Scopus
  24. P.-Y. le Bas, F. Luppé, J.-M. Conoir, and H. Franklin, “N-shell cluster in water: multiple scattering and splitting of resonances,” The Journal of the Acoustical Society of America, vol. 115, no. 4, pp. 1460–1467, 2004. View at Publisher · View at Google Scholar · View at Scopus
  25. L. Rayleigh, “On the passage of waves through apertures in plane screens and allied problems,” Philosophical Magazine, vol. 43, no. 263, pp. 259–272, 1897. View at Publisher · View at Google Scholar
  26. A. C. Eringen and E. S. Suhubi, Elastodynamics, Academic Press, 1974.
  27. J. D. Achenbach, Wave Propagation in Elastic Solids, North-Holland, Amsterdam, The Netherlands, 1973.
  28. Y. H. Pao and C. C. Mow, Diffraction of Elastic Waves and Dynamic Stress Concentrations, Crane and Russak, New York, NY, USA, 1973.
  29. M. M. Sun, Research on multiple scattering of elastic waves by rows of discontinuous barrier [Ph.D. thesis], Zhejiang University, 2012.
  30. G. N. Watson, A Treatise on the Theory of Bessel Functions, Cambridge University Press, Oxford, UK, 1966.