Jtam.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 47, 1, pp. 127-141, Warsaw 2009 INFLUENCE OF ELASTIC MATERIAL COMPRESSIBILITY ON PARAMETERS OF THE EXPANDING SPHERICAL STRESS WAVE. I. ANALYTICAL SOLUTION TO THE PROBLEM Edward Włodarczyk Military University of Technology, Faculty of Mechatronics, Warsaw, Poland e-mail: Edward.Włodarczyk@wat.edu.pl Mariusz Zielenkiewicz Military Institute of Armament Technology, Zielonka, Poland e-mail: m.zielenkiewicz@chello.pl We investigated the influence of elastic material compressibility on parameters of the expanding spherical stress wave. Thematerial compressibility is represented by Poisson’s ratio ν. The stress wave is generated by pressure created inside the sphe- rical cavity. The isotropic elastic material surrounds this cavity. Analytical closed form formulae determining thedynamical state ofmechanical parameters (displace- ment, particle velocity, strains, stresses, andmaterial density) in thematerial have been derived. These formulae were obtained for surge pressure p(t) = p0 = const inside the cavity. From analysis of these formulae it results that Poisson’s ratio ν substantially influences the course of material parameters in space and time. All parameters intensively decrease in space together with increase of the Langrangian coordinate r. On the contrary, these parameters oscillate versus time around the- ir static values. These oscillations decay with a lapse of time. We can mark out two ranges of the parameter ν values in which vibrations of the parameters are damped with a different degree. Thus, a decrease in Poisson’s ratio in the range ν ¬ 0.4 causes an intense decay of oscillation of parameters. On the other hand, in the range 0.4 < ν < 0.5, i.e. in quasi-compressible materials the damping of parameters vibrations is very low. In the limiting case when ν = 0.5, i.e. in the incompressible material damping vanishes, and the parameters harmonically oscil- late around their static values. The abnormal behaviour of the material occurs in the range 0.4< ν ¬ 0.5. In this case an insignificant increment of Poisson’s ratio causes considerable an increase of the parameters vibration amplitude. The speci- fic influence of Poisson’s ratio on the parameters of the expanding spherical stress wave in elastic media is the main result of this paper. As we see it, this fact may be the contribution supplementing the description of properties of the expanding spherical stress wave in elastic media. Key words: expansion of spherical stress wave, isotropic elastic material, dynamic load 128 E. Włodarczyk, M. Zielenkiewicz 1. Introduction In some theoretical analyses of dynamics of liners driven by explosives, the compressibility of their materials is neglected. On one hand, this assumption makes analytical solution tomany boundary value problems possible (Walters and Zukas, 1989; Cole, 1948; Kaliski et al., 1992; Gurney, 1943, 1947; Taylor, 1961; Trębiński et al., 1988a,b, 1989a; Włodarczyk and Zielenkiewicz, 2008). On the other hand, this simplification in a real physical systemneutralises the wavy course of the process in this system. On the contrary, results of theore- tical analyses (Lambour andHarley, 1965; Vidart et al., 1965; Knoepfel, 1970; Trębiński et al., 1989b) and experimental studies (Gimenez et al., 1985;Deren- towicz et al., 1984) show that wave phenomena have substantial quantitative and qualitative influence on the driving process of solids. Bearing inmind the results of these publications and theneeds of explosion mechanics, the influence of the elastic material compressibility on parameters of the expanding spherical stress wave has been theoretically investigated in this paper. The stress wave has been generated by pressure dynamically cre- ated inside the spherical cavity. The isotropic linear-elasticmaterial surrounds this cavity. Thematerial compressibility is represented by Poisson’s ratio ν. The paper consists of two parts. An analytical solution to the considered problem and its introductory analysis are placed in the first part. In turn, the vast quantitative and qualitative analysis of mechanical parameters (displace- ment, particle velocity, strains, stress, andmaterial density) of the expanding spherical stresswave in the compressible isotropic linear-elasticmaterial versus the Lagrangian coordinate r and time t are presented in the second part. The results of analysis presented in this paper can be used, among other things, to research spherical ballistic casings (Włodarczyk and Zielenkiewicz, 2008). From our point of view, the results of this analysis are a modest con- tribution of knowledge to the theory of propagation of stress waves in solids. 2. Formulation of the problem Propagation of the expanding spherical stress wave in an unbounded elastic material is considered. Thematerial is isotropic and compressible. The stress wave has been generated by a time–dependent pressure created inside the spherical cavity (Fig.1). The symbol r0 placed in Fig.1 denotes the initial radius of the cavity. Influence of elastic material compressibility... 129 Fig. 1. Physical scheme of the boundary value problem Taking into account spherical symmetry of the problem, it can be consi- dered as a one-dimensional boundary value problem. Independent variables of the problem are the Langrangian coordinate r and time t. The states of stress and strain in thematerial are represented by the following components: σr – radial stresses, σϕ =σθ – circumferential stresses, εr – radial strain and εϕ = εθ – circumferential (tangential) strains. The rest of components of the stress and the strain tensors are equal to zero in the considered coordinate system. According to the linear elasticity theory (Nowacki, 1970), we have εr = ∂u ∂r εϕ = εθ = u r υ= ∂u ∂t (2.1) and σr =2µεr +λ(εr +2εϕ)= (2µ+λ)εr +2λεϕ =(2µ+λ) ∂u ∂r +2λ u r (2.2) σϕ =2µεϕ+λ(εr +2εϕ)= 2(µ+λ)εϕ+λεr =2(µ+λ) u r +λ ∂u ∂r where u is the radial displacement, υ is the radial particle velocity, λ and µ are Lame’s constants, namely λ= νE (1+ν)(1−2ν) µ= E 2(1+ν) (2.3) In turn, the symbols E and ν denote Young’s modulus and Poisson’s ratio, respectively. 130 E. Włodarczyk, M. Zielenkiewicz For an element of the linear-elastic material, the equation of motion can be written in the form ρ0 ∂2u ∂t2 = ∂σr ∂r + 2(σr −σϕ) r (2.4) where ρ0 is the initial density of the material. Through substitution of expressions (2.2) into Eq. (2.4) and simple trans- formations, we obtain ∂2u ∂t2 = a2 ∂2u ∂r2 +2a2 (1 r ∂u ∂r − u r2 ) (2.5) where a2 = 1−ν (1+ν)(1−2ν) a20 a 2 0 = E ρ0 (2.6) The quantity a denotes the velocity of stress wave propagation. Equation (2.5) has been solved for the following boundary conditions u(r,t)≡ 0 for r= r0+at (2.7) and σr(r,t)=    (2µ+λ) ∂u(r,t) ∂r +2λ u(r,t) r =−p(t) for r= r0 p(t)> 0 σr(r,t)≡ 0 for r=∞ (2.8) According to themass conservation law, the equation ofmediumcontinuity can be written as follows ( 1+ u r )2( 1+ ∂u ∂r ) = ρ0 ρ(r,t) (2.9) where ρ(r,t) is current density of the material. For small strains, Eq. (2.9) can be transformed into the form 1+ ∂u ∂r +2 u r = ρ0 ρ (2.10) The structure of the analytical solution to the problem formulated above has been presented below. Influence of elastic material compressibility... 131 3. General analytical solution to the problem The general solution of Eq. (2.5) has the following form (Achenbach, 1975; Broberg, 1956; Graff, 1975; Hopkins, 1960; Kaliski et al., 1992) u(r,t) = ϕ′(r−r0−at) r − ϕ(r−r0−at) r2 (3.1) where the symbol ϕ′ denotes derivative of the function ϕ with respect to its argument. For ϕ(∞) 6=∞ and ϕ′(∞) 6=∞ solution (3.1) fulfils also boundary con- dition (2.8)2. Boundary condition (2.7) and solution (3.1) yield ϕ′(0)=ϕ(0)= 0 (3.2) where r−r0−at=0 is the equation of the wave front propagating from the cavity (Fig.1). The numerical values of the independent variables r and t are contained within the intervals r0 ¬ r¬∞ t­ r−r0 a (3.3) Expression (3.1) and condition (2.7)1 yield ϕ′′(x0)−2hϕ′(x0)+ 2h r0 ϕ(x0)=− 1 E (1+ν)(1−2ν) 1−ν r0p ( −x0 a ) (3.4) where h= 1−2ν 1−ν 1 r0 ­ 0 x0 =−at (3.5) Straightforward integration of Eq. (3.4) together with initial conditions (3.2) results in ϕ(x0)=− (1+ν) √ 1−2ν E r20 x0 ∫ 0 p (y−x0 a ) ehy sin(ωy) dy (3.6) where ω= √ 1−2ν (1−ν)r0 The function ϕ and its derivatives ϕ′ and ϕ′′ uniquely determine all pa- rameters of the considered problem. 132 E. Włodarczyk, M. Zielenkiewicz If the spherical cavity surface is loaded by the pressure p = p0 = const created in a statical way, then Eq. (2.5) can be written in the form d2us dr2 +2 (1 r dus dr − us r2 ) =0 (3.7) with boundary conditions σr(r0)= (2µ+λ) dus dr ∣ ∣ ∣ r=r0 +2λ us r0 =−p0 p0 > 0 (3.8) σr(∞)= 0 The general solution to Eq. (3.7) has the form us(r)=Cr+ D r2 (3.9) From relationships (2.3), boundary conditions (3.8) and solution (3.9) it follows that C =0, and D=(1+ν)r30p0/(2E). Finally, the parameters determining the statical state of the material are defined by the formulae us(r)= 1+ν 2 p0 E r0 (r0 r )2 εrs(r)=−(1+ν) p0 E (r0 r )3 εϕs(r)= 1+ν 2 p0 E (r0 r )3 σrs(r)=−p0 (r0 r )3 σϕs(r)= p0 2 (r0 r )3 (3.10) 4. Particular solution to the problem for surge pressure inside the cavity Integration of Eq.(3.6) for p(t) = p0 = const and differentiation of the func- tion ϕ(x) yield ϕ(x) =−1+ν 2 r30 p0 E [1+ √ 1−2νehx sin(ωx)− ehxcos(ωx)] ϕ′(x)=−(1+ν) √ 1−2νr20 p0 E ehx sin(ωx) (4.1) ϕ′′(x)=−(1+ν)(1−2ν) 1−ν r0 p0 E ehx[ √ 1−2ν sin(ωx)+cos(ωx)] Influence of elastic material compressibility... 133 where ωx= √ 1−2ν 1−ν ( r r0 −1− √ 1−ν (1+ν)(1−2ν) a0t r0 ) (4.2) hx= 1−2ν 1−ν ( r r0 −1− √ 1−ν (1+ν)(1−2ν) a0t r0 ) In order to simplify the quantitative analysis of the stresswave parameters, the following dimensionless quantities have been introduced ξ= r r0 η= a0t r0 U = u r0 Us = us r0 V = υ a0 R= ρ ρ0 Sr = σr p0 Srs = σrs p0 Sϕ = σϕ p0 Sϕs = σϕs p0 P = p0 E (4.3) The dimensionless variables ξ and η, according to relationships (3.3), (2.6) and (4.3) are contained within the intervals 1¬ ξ¬∞ η­ √ (1+ν)(1−2ν) 1−ν (ξ−1) (4.4) The above-mentioned formulae determining the stress wave parameters, expressed by dimensionless quantities (4.3), can be written as follows U(ξ,η) = 1+ν 2 P ξ2 {1− [ √ 1−2ν(2ξ−1)sin(ωx)+cos(ωx)]ehx} V (ξ,η) = P ξ {[ (1−2ν) √ 1+ν 1−ν − √ 1−ν2 1 ξ ] sin(ωx)+ + √ (1+ν)(1−2ν) 1−ν cos(ωx) } ehx εr(ξ,η) =−(1+ν)P { 1 ξ3 + [√ 1−2ν (1−2ν 1−ν 1 ξ −2 1 ξ2 + 1 ξ3 ) sin(ωx)+ − ( 1 ξ3 − 1−2ν 1−ν 1 ξ ) cos(ωx) ] ehx } εϕ(ξ,η)= 1+ν 2 P ξ3 {1− [ √ 1−2ν(2ξ−1)sin(ωx)+cos(ωx)]ehx } (4.5) R(ξ,η) = 1 1+εr(ξ,η)+2εϕ(ξ,η) 134 E. Włodarczyk, M. Zielenkiewicz Sr(ξ,η) =− 1 ξ3 {1+(ξ−1)[ √ 1−2ν(ξ−1)sin(ωx)+(ξ+1)cos(ωx)]ehx} Sϕ(ξ,η) = 1 2ξ3 { 1+ [√ 1−2ν ( − 2ν 1−ν ξ2−2ξ+1 ) sin(ωx)+ − ( 2ν 1−ν ξ2+1 ) cos(ωx) ] ehx } where ωx= √ 1−2ν 1−ν (ξ−1)− 1√ 1−ν2 η hx= 1−2ν 1−ν (ξ−1)− √ 1−2ν 1−ν2 η (4.6) On the cavity surface, i.e. for ξ=1, from these formulae we obtain U(1,η) = 1+ν 2 P [ 1− √ 2(1−ν)exp ( − √ 1−2ν 1−ν2 η ) sin ( − 1√ 1−ν2 η+α ) ] V (1,η) = √ 1−ν2P exp ( − √ 1−2ν 1−ν2 η ) sin ( 1√ 1−ν2 η+αV ) εr(1,η) =−(1+ν)P [ 1− √ 2ν2 1−ν exp ( − √ 1−2ν 1−ν2 η ) sin ( − 1√ 1−ν2 η+α ) ] εϕ(1,η) =U(1,η) (4.7) R(1,η) = [ 1− (1+ν)(1−2ν) √ 2 1−ν P exp ( − √ 1−2ν 1−ν2 η ) · ·sin ( − 1√ 1−ν2 η+α ) ]−1 Sr(1,η) =−1 − boundary condition Sϕ(1,η) = 1 2 − ν+1 √ 2(1−ν) exp ( − √ 1−2ν 1−ν2 η ) sin ( − 1√ 1−ν2 η+α ) where tanα= 1√ 1−2ν tanαV = √ 1−2ν ν (4.8) In the limiting case when ν = 0.5, i.e. for an incompressible material, formulae (4.6) and (4.8) yield ωx=− 2√ 3 η hx=0 α= π 2 αV =0 (4.9) Influence of elastic material compressibility... 135 In this case from formulae (4.5) we obtain U(ξ,η) = 3 4 P 1 ξ2 ( 1− cos 2√ 3 η ) V (ξ,η) = √ 3 2 P 1 ξ2 sin 2√ 3 η εr(ξ,η) =− 3 2 P 1 ξ3 ( 1− cos 2√ 3 η ) εϕ(ξ,η)= 3 4 P 1 ξ3 ( 1− cos 2√ 3 η ) (4.10) R(ξ,η) =1 Sr(ξ,η) =− 1 ξ3 [ 1+(ξ2−1)cos 2√ 3 η ] Sϕ(ξ,η) = 1 2ξ3 [ 1− (2ξ2+1)cos 2√ 3 η ] The expressions analogous to formulae (4.10) have been obtained for a spherical casing inWłodarczyk and Zielenkiewicz (2008), namely U(ξ,η)= 3 4 P β3 β3−1 1 ξ2 [ 1− cos ( 2√ 3 √ 1+ β+1 β2 η )] V (ξ,η) = √ 3 2 P β3 β3−1 √ 1+ β+1 β2 1 ξ2 sin ( 2√ 3 √ 1+ β+1 β2 η ) εr(ξ,η) =− 3 2 P β3 β3−1 1 ξ3 [ 1− cos ( 2√ 3 √ 1+ β+1 β2 η )] εϕ(ξ,η) = 3 4 P β3 β3−1 1 ξ3 [ 1− cos ( 2√ 3 √ 1+ β+1 β2 η )] R(ξ,η)= 1 (4.11) Sr(ξ,η) =− 1 ξ3−1 { (β ξ )3 −1+ − [(β ξ )3− (β2+β+1) (β ξ −1 ) −1 ] cos ( 2√ 3 √ 1+ β+1 β2 η )} Sϕ(ξ,η)= 1 2(β3−1) { 2+ (β ξ )3+ − [(β ξ )3+2(β2+β+1) (β ξ −1 ) +2 ] cos ( 2√ 3 √ 1+ β+1 β2 η )} 136 E. Włodarczyk, M. Zielenkiewicz where β = r1/r0. The symbols r1 and r0 denote the outer and inner radii of the casing, respectively. From the analysis of formulae (4.10) and (4.11) it follows that by their means, for β­ 5, we obtain comparable numerical results. The differences are smaller than 1%. 5. Introductory analysis of the problem We consider parameters of the stress wave which has been created by surge pressure p(t)= p0 = const inside the spherical cavity. From introductory ana- lysis of the formulae derived above it results that all parameters of the stress wave intensively decrease in space together with an increase of the distance from the centre of the system.On the contrary, these parameters oscillate ver- sus time around their static values in the respective spherical sections of the medium. Exemplary courses of the relative displacement of the cavity surface, U(1,η)/P, versus the dimensionless time, η = a0t/r0, for a few values of Poisson’s ratio ν are depicted in Fig.2. It is well known that the quantity ν represents compressibility of the medium. As it can be seen, the parameter ν substantially influences the course of function U(1,η)/P versus η. Fig. 2. Calculated relative displacement of the cavity surface U(1,η)/P versus dimensionless time η for selected values of Poisson’s ratio ν Influence of elastic material compressibility... 137 Wecanmarkout two ranges of thevalues ν inwhichvibrationof the cavity surface is dampedwith a different degree. Thus, a decrease of the parameter ν in the range ν ¬ 0.4 causes an intense decay of the cavity surface vibration. For these numerical values ν, the displacement of the cavity surface attains its statical value, i.e. Us/P =(1+ν)/2, during one cycle of the vibration (Fig.2). On the other hand, in the range 0.4 < ν < 0.5, i.e. in quasi-compressible media the vibration damping is very low. In the limiting case, when ν = 0.5, i.e. in the incompressible medium, the damping vanishes and the cavity surface harmonically vibrates around its static position Us/P = 0.75 with the constant amplitude AU =0.75 (Fig.2). It is necessary to take into account abnormal behaviour of the materials in the range 0.4 < ν ¬ 0.5. In this case, an insignificant increment of the parameter ν causes a considerable increase of the vibration amplitude of the cavity surface (Fig.3). For example, an increment of the Poisson’s ratio by ∆ν = 0.5− 0.4 = 0.1 causes a variation of the relative displacement of the cavity surface by ∆U/P =1.5−0.7=0.8. Fig. 3. Calculated relative displacement of the cavity surface U(1,η)/P versus dimensionless time η and Poisson’s ratio ν Theanalysis presented above applies tomotion of the cavity surface ξ=1. Respective spherical sections of thematerial for ξ > 1displace in ananalogous way. The displacements of these sections are adequately reduced (Fig.4). The material density changes insignificantly in both ranges of Poisson’s ratio (Fig.5). The maximal increment of the relative density ∆R is smaller than 0.2%. From presented analysis it follows that the considered material 138 E. Włodarczyk, M. Zielenkiewicz Fig. 4. Calculated relative displacement U(ξ,η)/P versus dimensionless variables ξ and η for Poisson’s ratio ν=0.3 Fig. 5. Calculated relative density R versus dimensionless time η for ξ=1 and selected values of Poisson’s ratio ν is compressible even though its density changes very little. As it is known, the compressibility measure of a linear-elastic material is Poisson’s ratio. The density change is not a measure of the material compressibility. The phenomenon of specific influence of Poisson’s ratio ν on the space r and time t variability of parameters of the spherical stress wave expanding in elastic media presented above is the main result of this paper. As we see it, this paper may be a contribution supplementing the description of properties Influence of elastic material compressibility... 139 of the spherical stress wave expanding in elastic media and has important significance for technical problems. According to the author’s knowledge, this phenomenon has not been described in the available literature yet. The vast quantitative and qualitative analysis of the mechanical parame- ters (displacement, particle velocity, strains, stresses and material density) of the expanding spherical stress wave in the compressible isotropic elasticmate- rial versus the Langrangian coordinate r and time is presented in the second part. References 1. Achenbach J.D., 1975, Wave Propagation in Elastic Solids, North-Holland Publ. Co., American Elsevier, Amsterdam-NewYork 2. Broberg K.B., 1956, Shock Waves in Elastic and Elastic-Plastic Media, Sto- kholm 3. Cole R.H., 1948,Underwater Explosions, PrincetonUniversity Press, Prince- ton, New Jersey 4. Derentowicz H., Dybek W., et al., 1984, Generation strong shock waves into deuterium by the Voitenki’s explosive generator, Biul. WAT, 33, 5 [in Polish] 5. Gimenez P., et al., 1985, EOS of detonation products obtained from sym- metrical deflection of liners investigated by laser interferometers techniques, In: The 8-th International Symposium on Detonation, Albuquerque 6. Graff K.F., 1975,Wave Motion in Elastic Solids, University Press, Oxford 7. Gurney R.W., 1943, The initial velocities of fragments from bombs, shells and grenades,BRL Report 405, September 8. Gurney R.W., 1947, Fragmentation of bombs, shells and grenades,BRL Re- port 635, March 9. Hopkins H.G., 1960, Dynamic expansion of spherical cavities in metals, In: Progress in Solid Mechanics, Sneddon J.N., Hill R. (Edit.), vol II, North- Holland Publishing Company, Amsterdam 10. Kaliski S., Rymarz C., Sobczyk K., Włodarczyk E., 1992, Waves, El- sevier, Amsterdam 11. Knoepfel H., 1970, Pulsed High Magnetic Fields, North-Holland Publishing Company, Amsterdam-London 140 E. Włodarczyk, M. Zielenkiewicz 12. LambourB.O.,Harley J.E., 1965,Proc the 4-th Symposium onDetonation, Washington 13. Nowacki W., 1970,Theory of Elasticity, PWN,Warszawa [in Polish] 14. Taylor G.I., 1961, Analysis of the explosion of a long cylindrical bomb deto- nated at one end, In:The Scientific Papers of Sir Geoffrey IngramTaylor, Vol. III, Batchelor G.K. (Edit.), University Press, Cambridge 15. Trębiński R., Trzciński W.A., Włodarczyk E., 1988a, Theoretical ana- lysis of the process of driving a cylindrical liner by the products of grazing detonation, J. Tech. Phys., 29, 3-4 16. TrębińskiR.,TrzcińskiW.A.,WłodarczykE., 1988b,Theoretical analy- sis of the influence of an external envelopeon the process of driving a cylindrical liner by the products of grazing detonation, J. Tech. Phys., 29, 3-4 17. Trębiński R., Trzciński W.A., Włodarczyk E., 1989a, Theoretical ana- lysis of the influence of the external medium on the launching process of a cylindrical liner by the products of grazing detonation, J. Tech. Phys., 30, 2 18. Trębiński R., Trzciński W.A., Włodarczyk E., 1989b, Theoretical ana- lysis of the influence of material compressibility of a cylindrical liner on the process of launching it by products of grazing detonation, J. Tech. Phys., 30, 3-4 19. Vidart A., Beatrix P., et al., 1965,Proc the 4-th Symposium on Detona- tion, Washington 20. Walters W.P., Zukas J.A., 1989, Fundamentals of Shaped Charges, Wiley and Sons, NewYork 21. Włodarczyk E., Zielenkiewicz M., 2008, The dynamics of a thick-walled spherical casing loaded with a time depending internal pressure, Journal of Theoretical and Applied Mechanics, 46, 1 Teoretyczna analiza wpływu ściśliwości ośrodka sprężystego na parametry ekspandującej kulistej fali naprężenia. I. Analityczne rozwiązanie problemu Streszczenie Zbadano wpływ ściśliwości materiału, reprezentowanej przez współczynnik Po- issona ν, na parametry ekspandującej kulistej fali naprężeń, wywołanej ciśnieniem wewnątrz kulistej kawerny w nieskończonym izotropowym ośrodku sprężystym.Wy- prowadzono zamknięte wzory określające dynamiczny stan mechanicznych parame- trów (przemieszczenia, prędkości przemieszczenia, odkształcenia, naprężenia oraz gę- stości materiału) w ośrodku dla stałego ciśnienia p(t) = p0 = const przyłożonego Influence of elastic material compressibility... 141 w sposób nagły. Z analizy tych wzorówwynika, że współczynnik Poissona ν wpływa zasadniczo na przebiegi tych parametróww czasie i przestrzeni.Wszystkie parametry maleją intensywnie wraz ze wzrostem współrzędnej Lagrange’a r. Z drugiej strony, parametry te oscylują w czasie wokół ich wartości statycznych. Oscylacje te z cza- sem zanikają. Możemy wyróżnić dwa przedziały wartości parametru ν, w których tłumienie oscylacji zachodzi z różną intensywnością. Tak więc dla ν ¬ 0.4 tłumienie jest bardzo intensywne. Natomiast w przedziale 0.4 < ν < 0.5, tzn. w materiałach quasi-ściśliwych, poziom tłumienia jest bardzo niski. W przypadku granicznym dla ν=0.5, tzn. dlamateriałównieściśliwych, tłumienie zanika i parametrymechaniczne oscylują wokół ich wartości statycznych. Ponadto w przedziale 0.4<ν ¬ 0.5 można zaobserwować anormalne zachowanie materiału. Niewielki przyrost wartości współ- czynnika Poissona skutkuje znaczącymwzrostem amplitudy drgań parametrów.Opis szczególnego wpływu współczynnika Poissona na parametry ekspandującej kulistej fali naprężenia jest głównym rezultatem niniejszej pracy. Według autorówmoże ona być uzupełnieniem opisu właściwości ekspandującej kulistej fali naprężeń w ośrodku sprężystym. Manuscript received January 18, 2008; accepted for print July 16, 2008