24SIMPLE EXPRESSIONS OF THE DYNAMIC STIFFNESS OF GROUPED PILES⎛α⎞ ⎛⎞2R0R0k ≅ ⎜ +⎟ =⎜ +⎟s2πξkα1µLa0.4 µ La2πξkα20. 8πξkα1(2.19)⎝ α1La⎠ ⎝ La⎠For <strong>the</strong> <strong>dynamic</strong> case (ω →∞), non-dimensional frequency a 0converges on:ω R0a0= i = ia(2.20)vT pEquation (2.16) is thus approximated by:k ≅ L i ⋅ 2πξα ⋅ a + 0. πξ α(2.21)s( )µa k 18k1From Equations (2.19) and (2.21), soil <strong>stiffness</strong> will presumably be approximated as:⎧⎛R⎫0⎞k ⎨⎜⎟s≅ µ La2πξkα2+ 0.8πξkα1+ i ⋅ 2πξkα1⋅ a⎬(2.22)⎩⎝La⎠⎭Even without <strong>the</strong> soil above <strong>the</strong> active pile length, <strong>the</strong> pile group exhibits its own<strong>stiffness</strong>, k (Figure 2.7b), which is described as:kg4EIp µ L0α3gα3= α3 3= µ L3 3 0LaLaα0122≅ (2.23)⎛ d ψ ( ζ )where, ∫ ⎟ ⎞α3=⎜ dζ20 ⎝ dζ⎠Both ksand k gsustain <strong>the</strong> massamong <strong>the</strong> <strong>piles</strong>. This massmgLa02 20ψ z)mg<strong>of</strong> <strong>the</strong> embedded pile group with soil caughtmgis approximated by:20≅ ∫ ρ πR( dz = ρ πRL α(2.24)ssaTherefore <strong>the</strong> overall <strong>stiffness</strong> k xx<strong>of</strong> <strong>the</strong> pile cap for sway motion is given as:kxx = ks + kg −mgω 2 (2.25)From Equations (2.22), (2.23) and (2.24), Equation (2.25) is rewritten as:⎡⎪⎧R ⎛⎞⎪⎫⎤0α⎢⎨⎜32k⎟xx≅ µ La2πξkα2+⎬ + ⋅ ⋅ − ⋅ ⎥⎢0.8πξkα1+3i 2πξkα1a πα1a (2.26)⎣⎪⎩La⎝ α0 ⎠⎪⎭⎥⎦Substituting Equation (2.3) into Equation (2.26), one obtains:⎡⎪⎧R ⎛⎞⎪⎫⎤0α⎢⎨⎜32k⎟xx≅ µ L0 2πξkα2+⎬ + ⋅ ⋅ − ⋅ ⎥⎢0.8πξkα0α1+2i 2πξkα0α1a πα0α1a (2.27)⎣⎪⎩L0⎝α0 ⎠⎪⎭⎥⎦It is now obvious that k xxin Equation (2.27) has <strong>the</strong> following <strong>simple</strong> form withfrequency-independent <strong>stiffness</strong> k 0, and damping and mass parameters c 0and m 0respectively:2≅ k + i ⋅ c ⋅ a − m ⋅(2.28)where,k xx0 00a1
SIMPLE EXPRESSIONS OF THE DYNAMIC STIFFNESS OF GROUPED PILES 25k0R0c0m0= c1+ c2, = c3andµ L0L0µ L0µ L0= c4(2.29a)-(2.29c)with c1= 2πξ kα2,α3c2= 0. 8πξ kα0α1+2α0, c3= 2πξ kα0α1and c4= πα0α1.(2.30a)-(2.30d)The above equations show some important features <strong>of</strong> <strong>the</strong> pile cap <strong>stiffness</strong>. Among <strong>the</strong>parameters specified in <strong>the</strong> above equations, c1, c2, c3and c 4are dependent on <strong>the</strong>shape function ψ ( ζ ), which may not differ drastically in different soil-pile systems aslong as <strong>piles</strong> exhibit a flexible nature, and k / L 0µ0alone includes a term proportionalto 0 0. Equation (2.28) was derived with <strong>the</strong> intention <strong>of</strong> showing what could be <strong>the</strong>most important key parameters that determine . The assumption taken to derive <strong>the</strong>equation is good enough for this purpose, but certainly is an oversimplification <strong>of</strong> reality.Since soils below active pile lengths are not allowed to deform at all, <strong>the</strong> assumption isliable to lead to overestimation <strong>of</strong> <strong>the</strong> <strong>stiffness</strong> parameter k 0and underestimation <strong>of</strong><strong>the</strong> damping parameter c0. Therefore, parameters c 1, c2, c3and c 4were obtainednot directly from Equations (2.30a)-(2.30d), but in such a way that <strong>the</strong> overall errorwould be minimized for <strong>the</strong> variety <strong>of</strong> soils and pile parameters examined. Theparameters that have been considered are: 1) pile parameters such as group-pile <strong>stiffness</strong>,EI ( = npEpIp), and active pile length ratio, L0 / L; and 2) soil parameters includingshear modulus µ and material damping D . In this discussion, only a homogeneoussoil pr<strong>of</strong>ile with a square arrangement <strong>of</strong> <strong>piles</strong> is considered. The best fit <strong>of</strong> <strong>the</strong> valuesfrom Equation (2.28) to rigorous solutions <strong>of</strong> k xxwas obtained by setting c 1, c 2, c 3and c 4at 2π , π /2, 2π and π /4, respectively. Some representative cases areshown in Figures 2.8a-2.8f.The present <strong>simple</strong> <strong>expression</strong> <strong>of</strong> k xx(Equation (2.28)) allows <strong>the</strong> effects <strong>of</strong> overallsite non-linearity to be reflected by simply replacing <strong>the</strong> shear modulus <strong>of</strong> <strong>the</strong> intact soil,µ , with <strong>the</strong> complex modulus, µ'( 1 + iD ); this describes equivalent-linear features <strong>of</strong><strong>the</strong> soil experiencing <strong>dynamic</strong> loading, and is obtained from shear-modulus-reductionand damping ratio curves <strong>of</strong> <strong>the</strong> soil. This manipulation, however, causes <strong>the</strong> <strong>stiffness</strong>and damping parameters k 0and c 0in Equations (2.29a) and (2.29b) to be slightlydependent on frequency as:k0⎛ R0π ⎞= ⎜2π+ ⎟ −2πD⋅ a(2.31a)µ ' L ⎝ L 2 ⎠00c0⎛ R0π π 2⎞= 2π+ ⎜2π+ − a ⎟ D(2.31b)µ ' L0⎝ L02 4 ⎠When <strong>the</strong> effect <strong>of</strong> D cannot be ignored in Equations (2.31a) and (2.31b), appropriatekxx