# Wave Dispersion and Propagation in Linear Peridynamic Media

Authors: Xiaolong Zhang; Zhenying Xu; Qing Yang
Journal: Shock and Vibration
 (2019)

---

## Abstract

We detail the linear peridynamic wave equation with a nonlocal integral form based on the linear peridynamic and dynamic theory. Wave dispersion in an infinite maraging steel material is obtained by analyzing the linear peridynamic wave equation. The dispersion curves, group velocity, phase velocity, and other wave parameters of the shear and longitudinal waves in an infinite media are obtained using numerical methods. We obtained the optimal calculation parameters by analyzing the weight function, horizon, mesh size, and other numerical calculation parameters on the dispersion curve. We simulated the propagation of waves in an infinite media by applying these parameters into the peridynamic wave equation. We conclude that the wave model can generate waves that propagate in all directions with initial loads. The wavefront is an ellipsoid.


---

## Body

## 1. Introduction

Ultrasonic waves with unidirectionality and penetration properties are widely used in nondestructive testing [1–4]. In order to study the propagation of ultrasonic waves, the community usually uses finite element analysis methods to solve elastic dynamic equations in continuous media [5–8]. However, the derivative function does not exist at the defect when there is a defect (discontinuity) in the media [3, 4]. It is necessary to remesh the grid to make the defect outside of the boundary [9]. This method can simulate ultrasonic wave propagation in a media with the defect; however, the stresses at the crack tip are mathematically singular because of their undefined spatial derivatives. To study ultrasonic wave propagation in the media, we used the integral equation to describe the wave displacement change with location and time—named “the peridynamic method” [10, 11]. This method falls into the category of nonlocal models [12, 13] because particles separated by a finite distance can interact with each other [14, 15]. The peridynamic theory of solid mechanics has been proposed as a means of treating discontinuous media via a mathematical model that does not require a smooth distribution of mass or differentiability of the deformation [10]. The starting point of the theory is that the internal forces acting on a material point are determined via interactions between the point and all others within a finite distance from it. The peridynamic model uses integrals rather than differential equations so that the mathematical structure does not break down when a discontinuity occurs. Rather, the fracture is treated as a natural outcome of deformation that emerges according to the equations of motion and the constitutive model [16].The theory of peridynamic waves was first proposed by Silling [10] in 2000, with a study of the propagation of linear stress waves and a discussion of wave dispersion. In 2003, Silling et al. [17] discussed the propagation of waves in a discontinuities bar. Later, Zimmermann [18] explored many features of the peridynamic theory including certain aspects of wave motion, material stability, and numerical solution techniques. In 2007, Silling et al. [19] gave the relationship between stress, displacement, and location under small deformation conditions making it easier to solve the microelastic peridynamic wave equation. Silling [20] gave the linearization theory of peridynamic in 2010. This theory paves the way for the study of linear elastic waves, making the peridynamic wave equation easier to solve. Bažant et al. [21] systematically studied the peridynamic stress wave in 2016 and gave the bond-based peridynamic stress wave theory and dispersion analysis. In the same year, Silling [22] studied solitary waves in elastic media and found that the solitary wave velocity is greater than that of the elastic wave. Butt et al. [23] conducted a detailed study of the peridynamic stress wave under linear elastic conditions in 2017. They then analyzed the state-based peridynamic wave dispersion and studied the mesh size, horizon, and weight function influence on dispersion.In this paper, we gave the wave propagation equation a nonlocal integral form based on the research of linear peridynamic theory. The dispersion matrix [10] of the wave is given from the linear peridynamic theory, and the dispersion of the longitudinal and shear waves in the infinite media is analyzed. The phase velocity and group velocity of the waves in the infinite peridynamic media are given. The effects of weight function, horizon, and mesh size on the dispersion are analyzed. The peridynamic simulation is used in the infinite media to analyze the difference between the shear wave and the longitudinal wave. By analyzing the wave model, we conclude that the wavefront of the single-point excited wave is an ellipsoid.


## 2. Basic Equations of Linear Peridynamic Media

The continuum peridynamic wave theory is in contrast to the classical continuum theory. Any material pointx in the reference configuration is acted upon by forces due to the deformation of all the material points x′ within the neighborhood of finite radius δ centered at x. The radius δ is called the horizon, and the material points within this neighborhood of x in the reference configuration are called the family of x.A material pointx at time t has a force density bx,t acting on it, and the motion is given by the following equation [24]:(1)ρxu¨x,t=Lux,t+bx,t,∀x∈β,t≥0,(2)Lux,t=∫βT¯x,tx′−x−T¯x′,tx−x′dVx′,where β is the reference configuration of the body, ρ is the density in the reference configuration, u is the displacement, and b is the body force density. Equation (2) is a function of displacement that represents the internal force density (per unit volume) exerting on x by other body points. Since force vector state T¯x,t depends only on the deformation of the family of x, and defining the vector x′−x=ξ to be a bond, we assume the following:(3)ξ>δ⇒Tξ=0.Under the assumption of a small deformation that does not need a small deformation gradient, the peridynamic equation (1) is accurately approximated by the following linear integrodifferential equation:(4)ρxu¨x,t=∫NxCx,quq,t−ux,tdVq+bx,t,where Nx is a neighborhood of x. As discussed in [20], the radius of Nx is in general 2δ, where δ is the horizon. Cx,q is a tensor-valued function which is named “the micromodulus function.”Similar to the linear classical theory, the linear peridynamic theory concerns small deformations. However, the damage and fracture that may be included in a linearized model, as noted above, does not require small displacement gradients [20]. Under this assumption of smallness, the peridynamic equation of equilibrium reduces to a linear integral equation. The linearized theory is more compatible with implicit solvers than the full theory. It is also much easier to obtain theoretical results about the properties of the linear equations such as well posedness than those of the nonlinear model. The properties of waves are also much easier to analyze within the linearized theory.


## 3. Wave Dispersion in Linear Peridynamic Media

The one-dimensional linear peridynamic wave model is popular because of its simplicity. Butt et al. [23] focused on the dispersion properties of a state-based linear peridynamic model. They derived the dispersion relation for one-, two-, and three-dimensional cases and specifically investigated the effect of the horizon, mesh size, and weighting function on the dispersion. Here, we focus on the dispersion characteristics of waves based on linear peridynamics. To study the peridynamic wave dispersion, we provide the general form of a plane wave:(5)ux,t=Aeikx−ωt.Of these,A is the vibration amplitude and u is the displacement, with k as the wavenumber and w as the circular frequency. When the external force density in (4) is zero, we can obtain the following equation by substituting (5) into (4):(6)−ω2ρxux,t=∫NxCξuq,t−ux,tdVq.Here,ξ=q−x is a bond between two material points. Changing the form of equation (6) leads to(7)−ω2ρxux,t=∫NxCξeikξ−1dVξ⋅ux,t.The Euler transform in (7) is as follows:(8)ω2ρxux,t=∫NxCξ1−coskξ−isinkξdVξ⋅ux,t.The micromodulus functionCξ is an even function, and sinkξ is an odd function and integrated into the asymmetric domain Nx. Thus, we can get the following:(9)∫NxCξisinkξdVξ=0.Thus, we conclude from (8) and (9) that(10)ω2ρxux,t=∫NxCξ1−coskξdVξ⋅ux,t.We conclude from equation (10) that the dispersion of the peridynamic wave is only related to the micromodulus function, density, and horizon. These parameters are inherent properties of the material.The micromodulus functionCξ=Cx,q can be written as(11)Cx,q=15μm2−9κm2∫Hpϖp−xϖp−qp−x⊗p−qdVp+2rq−x,(12)rξ=15μmϖξM⊗M.Of these,M=ξ/ξ, m=∫Hωξξ2dVξκ is the bulk modulus, μ is the shear modulus, and ϖ is the weight function. The weight function and horizon selection are under intense investigation (equation (11)). If the selected horizon and weight function are consistent with the properties of the material, then the simulation results are close to the experimental results. Otherwise, a large error will occur. Equations (10) and (11) can lead to the following:(13)ω2ρxux,t=Mk,ξux,t,(14)Mk,ξ=∫NxCξ1−coskξdVξ.From (11) and (12), we found that Mk,ξ is a matrix (dispersion matrix). This matrix first appeared in [10].This work assumes that a wave propagates along the wave vectorsn=1,0,0T in a Cartesian coordinate system. At this time, the stress wave propagates along x1, which can be seen as a one-dimensional rod model. Thus, we can obtain ξ=ξ1. The expansion of (13) is(15)ω2ρxu1x,tu2x,tu3x,t=M11M12M13M21M22M23M31M32M33u1x,tu2x,tu3x,t.When we study a longitudinal wave propagating alongx1, we can obtain only the displacement of u1x,t≠0. Equation (15) leads to the following:(16)ω112ρxu1x,t=M11u1x,t=∫NxC11ξ1−coskξ1dVξu1x,t.Equation (16) shows that the longitudinal wave along the direction x1 is related to the material parameter C11ξ; thus, we can obtain the longitudinal wave dispersion equation by simplifying (16) to be(17)ω11=1ρx∫NxC11ξ1−coskξ1dVξ.Similarly, the displacement of the horizontal shear wave propagating in the directionx1 only has u2x,t≠0. A combination of (14) and (15) shows that the dispersion is related to the material parameters C22ξ. The dispersion equation is as follows:(18)ω12=1ρx∫NxC22ξ1−coskξ1dVξ.Similarly, as the displacement of the vertical shear wave propagates alongx1, only u3x,t≠0. Combination of (14) and (15) shows that the dispersion is related to the material parameters C33ξ. The dispersion equation is as follows:(19)ω13=1ρx∫NxC33ξ1−coskξ1dVξ.Equations (11) and (12), and the symmetry of the spatial integral indicate that the micromodulus function matrix in an infinite media is a symmetric matrix; the values on the diagonal are equal. Therefore, (18) and (19) indicate the same shear wave dispersion equation with only different vibration directions.


### 3.1. Longitudinal Wave Dispersion

This peridynamic method was used to simulate the wave dispersion in maraging steel. The following terms were used: material densityρx=8000kg/m3, bulk modulus κ=1.59×1011pa, shear modulus μ=8.2×1010pa, mesh size h=1mm, and horizon δ=3.015h. If there was a vibration at the center of the infinite media, then the longitudinal wave propagates along x1within the infinite media, and we can obtain the discrete equation from dispersion equation (17) as(20)ω11=1ρ∑i=1NxC11ξi1−coskξi1ΔVξi,where(21)C11=15μm2−9κm2∑j=1p,j≠iϖpj−xϖpj−qipj1−x1pj1−qi1ΔVj+30μmϖξqi1−x1ξ2,(22)m=∑k=1Hxωξkξk2ΔVξk.The literature [5] shows that the relationship between the wavenumber, phase velocity, and circular frequency is as follows:(23)Cp=ωk.The group velocity can be found from the phase velocity via formula (24):(24)Cg=dωdk.Figures1–3 plots the dispersion curve, phase velocity, and group velocity, respectively, from (20)–(24).Figure 1
Longitudinal wave dispersion.Figure 2
Phase velocity of longitudinal wave.Figure 3
Group velocity of longitudinal wave.Figure1 shows that the circular frequency of the longitudinal wave increases with the wavenumber, but the rate of circular frequency gradually decreases. Figure 2 shows that the phase velocity of the longitudinal wave reaches a maximum of 5539.8 m/s when the wavenumber is near zero, which is consistent with the longitudinal wave velocity of 5500 m/s in maraging steel. The phase velocity gradually decreases with increasing wavenumber. The group velocity is equal to the phase velocity when the wavenumber is near zero (Figure 3). At this point, there is no dispersion, but the longitudinal wave begins to disperse as the wavenumber increases, and the group velocity decays faster than the phase velocity.The analysis above shows that the longitudinal wave dispersion is related to the material parameters. Equation (21) indicates that the material parameters are related to the weight function, horizon, and mesh size. Next, we separately study the effects on the dispersion curve. Given a material point x at the origin, the radius of the horizon domain is δ=3.015h and the mesh size is h=1mm. The dispersion curves at different weight functions are shown in Figure 4.Figure 4
Longitudinal wave dispersion curve affected by a weight function.In Figure4, we see that the dispersion curve of the wave function is basically the same under the weight functions ω1 and ω2, but the rate of circular frequency of the dispersion curve is negative when the wavenumber is large. This indicates that the group velocity of the wave is negative, which does not match the actual situation. Therefore, the weight function ω3 is better for wave simulation when the wavenumber is large. The dispersion under the three weight functions is similar when the wavenumber is small, and the calculation is simple under ω2.When the mesh size ish=1mm, the weight function is exp−ξ2/δ2, and the radius of the horizon domain is δ=2mm,3mm,5mm,and7mm. The dispersion image of the longitudinal wave is shown in Figure 5.Figure 5
Longitudinal wave dispersion curve affected by the horizon.The horizon has less influence on dispersion when the wavenumber is small (Figure5). When the wavenumber is large, the dispersion becomes more obvious as the radius of the horizon increases. When the radius of the horizon is 5 mm and 7 mm, the group velocity is negative, which is inconsistent with the actual situation. The calculation is complicated in both cases. A horizon radius of 2 mm and 3 mm can better reflect the dispersion phenomenon of longitudinal waves.When the horizon isδ=3mm, the weight function is exp−ξ2/δ2 and the mesh size is h=1.5mm, 1.0mm, 0.6mm, and 0.5mm. The dispersion curve of the longitudinal wave is shown in Figure 6.Figure 6
Longitudinal wave dispersion curve affected by mesh size.Figure6 shows that the mesh size has little effect on the dispersion when the wavenumber is small, but the mesh size is more pronounced when the wavenumber is large. When the mesh size is h1=1.5mm, the group velocity is negative, which is not consistent with the actual data. Although the dispersion is not obvious when h3=0.6mm and h=0.5mm, the computation is complex.


### 3.2. Shear Wave Dispersion

The material parameters are consistent with the longitudinal wave, and a shear wave propagates along the directionx1 with a load at the center of the infinite media. The dispersion equation is indicated by (18) or (19). The two dispersion equations are consistent. By discretizing the dispersion equations, we can obtain(25)ω12=1ρ∑i=1NxC22ξi1−coskξi1ΔVξi,C22=15μm2−9κm2∑j=1p,j≠iϖpj−xϖpj−qipj2−x2pj2−qi2ΔVj+30μmϖξqi2−x2ξ2.The termm is given in (22).From (22)–(24), the dispersion curve, phase velocity, and group velocity are given in Figures 7–9, respectively.Figure 7
Shear wave dispersion.Figure 8
Phase velocity of the shear wave.Figure 9
Group velocity of the shear wave.Figures7–9 indicate that the properties of the shear wave such as dispersion and wave velocity are similar to those of the longitudinal wave. The circular frequency increases with the wavenumber, but the rate of circular frequency gradually decreases. The phase velocity and group velocity of the shear wave reach 3216.4 m/s when the wavenumber nears zero. This agrees with the shear wave velocity of 3200 m/s in maraging steel. As the wavenumber increases, the phase velocity gradually slows down and starts to disperse and the group velocity decays faster than the phase velocity.Through the above analysis, it can be observed that the dispersion of the shear wave and the longitudinal wave is affected by similar material parameters, such as the weight function, horizon, and mesh size. Given a material pointx at the center of the origin, the radius of the horizon domain is δ=3.015h and the mesh size is h=1mm, the dispersion curves at different weight functions are shown in Figure 10.Figure 10
Shear wave dispersion curve affected by the weight function.When the mesh size ish=1mm, the weight function is exp−ξ2/δ2 and the radius of the horizon domain is δ=2mm, 3mm, 5mm, and 7mm; the dispersion of the shear wave is shown in Figure 11.Figure 11
Shear wave dispersion curve affected by the horizon.When the horizon isδ=3mm, the weight function is exp−ξ2/δ2 and the mesh sizes are h=1.5mm, 1.0mm, 0.6mm, and 0.5mm; the dispersion curve of the shear wave is shown in Figure 12.Figure 12
Shear wave dispersion curve affected by the mesh size.Figures10–12 indicate that the dispersion of the shear wave is influenced by the weight function, horizon, and mesh size. It is the same as the longitudinal wave. The dispersion curve is roughly the same under the three weight functions, except that the dispersion of the wave function is more pronounced under the weight functions ω1 and ω2. When the wavenumber is low, the horizon has little effect on the dispersion. When the wavenumber is large, the dispersion becomes more obvious with horizon increases. A small mesh size leads to a more obvious dispersion. However, the mesh size is relatively small when the horizon is large, and there are more material points involved in the horizon—this increases the complexity of the calculation.Figures4 and 10 show that the dispersion with different weight functions has little difference when the wavenumber is small. However, when the wavenumber is around 1000 rad/m, the group velocity becomes negative with weight functions ω1 and ω2. This indicates that the wave propagates in the inverse direction; this phenomenon does not occur in the real world. Thus, the weight function ω3 is better for simulation. Figures 5 and 11 show that the dispersion is noticeable as the horizon domain increases. The wave group velocity becomes negative when the horizon is greater than three times the mesh sizes—especially in the longitudinal wave. The peridynamic wave equation will converge to the classical Navier equation when the horizon is small. The horizon should be selected between two and three times the mesh sizes. In Figures 6 and 12, we see that the dispersion is weakened with decreasing mesh size. When the mesh size is small, there will be more material points on the horizon to complicate the calculations. The analysis above shows that the horizon should be three times the mesh sizes [24]. Therefore, in this paper, the weight function should be exp−ξ2/δ2, the horizon should be δ=3mm, and the mesh size should be h=1mm. This makes the wave simulation agree with natural reality.


## 3.1. Longitudinal Wave Dispersion

This peridynamic method was used to simulate the wave dispersion in maraging steel. The following terms were used: material densityρx=8000kg/m3, bulk modulus κ=1.59×1011pa, shear modulus μ=8.2×1010pa, mesh size h=1mm, and horizon δ=3.015h. If there was a vibration at the center of the infinite media, then the longitudinal wave propagates along x1within the infinite media, and we can obtain the discrete equation from dispersion equation (17) as(20)ω11=1ρ∑i=1NxC11ξi1−coskξi1ΔVξi,where(21)C11=15μm2−9κm2∑j=1p,j≠iϖpj−xϖpj−qipj1−x1pj1−qi1ΔVj+30μmϖξqi1−x1ξ2,(22)m=∑k=1Hxωξkξk2ΔVξk.The literature [5] shows that the relationship between the wavenumber, phase velocity, and circular frequency is as follows:(23)Cp=ωk.The group velocity can be found from the phase velocity via formula (24):(24)Cg=dωdk.Figures1–3 plots the dispersion curve, phase velocity, and group velocity, respectively, from (20)–(24).Figure 1
Longitudinal wave dispersion.Figure 2
Phase velocity of longitudinal wave.Figure 3
Group velocity of longitudinal wave.Figure1 shows that the circular frequency of the longitudinal wave increases with the wavenumber, but the rate of circular frequency gradually decreases. Figure 2 shows that the phase velocity of the longitudinal wave reaches a maximum of 5539.8 m/s when the wavenumber is near zero, which is consistent with the longitudinal wave velocity of 5500 m/s in maraging steel. The phase velocity gradually decreases with increasing wavenumber. The group velocity is equal to the phase velocity when the wavenumber is near zero (Figure 3). At this point, there is no dispersion, but the longitudinal wave begins to disperse as the wavenumber increases, and the group velocity decays faster than the phase velocity.The analysis above shows that the longitudinal wave dispersion is related to the material parameters. Equation (21) indicates that the material parameters are related to the weight function, horizon, and mesh size. Next, we separately study the effects on the dispersion curve. Given a material point x at the origin, the radius of the horizon domain is δ=3.015h and the mesh size is h=1mm. The dispersion curves at different weight functions are shown in Figure 4.Figure 4
Longitudinal wave dispersion curve affected by a weight function.In Figure4, we see that the dispersion curve of the wave function is basically the same under the weight functions ω1 and ω2, but the rate of circular frequency of the dispersion curve is negative when the wavenumber is large. This indicates that the group velocity of the wave is negative, which does not match the actual situation. Therefore, the weight function ω3 is better for wave simulation when the wavenumber is large. The dispersion under the three weight functions is similar when the wavenumber is small, and the calculation is simple under ω2.When the mesh size ish=1mm, the weight function is exp−ξ2/δ2, and the radius of the horizon domain is δ=2mm,3mm,5mm,and7mm. The dispersion image of the longitudinal wave is shown in Figure 5.Figure 5
Longitudinal wave dispersion curve affected by the horizon.The horizon has less influence on dispersion when the wavenumber is small (Figure5). When the wavenumber is large, the dispersion becomes more obvious as the radius of the horizon increases. When the radius of the horizon is 5 mm and 7 mm, the group velocity is negative, which is inconsistent with the actual situation. The calculation is complicated in both cases. A horizon radius of 2 mm and 3 mm can better reflect the dispersion phenomenon of longitudinal waves.When the horizon isδ=3mm, the weight function is exp−ξ2/δ2 and the mesh size is h=1.5mm, 1.0mm, 0.6mm, and 0.5mm. The dispersion curve of the longitudinal wave is shown in Figure 6.Figure 6
Longitudinal wave dispersion curve affected by mesh size.Figure6 shows that the mesh size has little effect on the dispersion when the wavenumber is small, but the mesh size is more pronounced when the wavenumber is large. When the mesh size is h1=1.5mm, the group velocity is negative, which is not consistent with the actual data. Although the dispersion is not obvious when h3=0.6mm and h=0.5mm, the computation is complex.


## 3.2. Shear Wave Dispersion

The material parameters are consistent with the longitudinal wave, and a shear wave propagates along the directionx1 with a load at the center of the infinite media. The dispersion equation is indicated by (18) or (19). The two dispersion equations are consistent. By discretizing the dispersion equations, we can obtain(25)ω12=1ρ∑i=1NxC22ξi1−coskξi1ΔVξi,C22=15μm2−9κm2∑j=1p,j≠iϖpj−xϖpj−qipj2−x2pj2−qi2ΔVj+30μmϖξqi2−x2ξ2.The termm is given in (22).From (22)–(24), the dispersion curve, phase velocity, and group velocity are given in Figures 7–9, respectively.Figure 7
Shear wave dispersion.Figure 8
Phase velocity of the shear wave.Figure 9
Group velocity of the shear wave.Figures7–9 indicate that the properties of the shear wave such as dispersion and wave velocity are similar to those of the longitudinal wave. The circular frequency increases with the wavenumber, but the rate of circular frequency gradually decreases. The phase velocity and group velocity of the shear wave reach 3216.4 m/s when the wavenumber nears zero. This agrees with the shear wave velocity of 3200 m/s in maraging steel. As the wavenumber increases, the phase velocity gradually slows down and starts to disperse and the group velocity decays faster than the phase velocity.Through the above analysis, it can be observed that the dispersion of the shear wave and the longitudinal wave is affected by similar material parameters, such as the weight function, horizon, and mesh size. Given a material pointx at the center of the origin, the radius of the horizon domain is δ=3.015h and the mesh size is h=1mm, the dispersion curves at different weight functions are shown in Figure 10.Figure 10
Shear wave dispersion curve affected by the weight function.When the mesh size ish=1mm, the weight function is exp−ξ2/δ2 and the radius of the horizon domain is δ=2mm, 3mm, 5mm, and 7mm; the dispersion of the shear wave is shown in Figure 11.Figure 11
Shear wave dispersion curve affected by the horizon.When the horizon isδ=3mm, the weight function is exp−ξ2/δ2 and the mesh sizes are h=1.5mm, 1.0mm, 0.6mm, and 0.5mm; the dispersion curve of the shear wave is shown in Figure 12.Figure 12
Shear wave dispersion curve affected by the mesh size.Figures10–12 indicate that the dispersion of the shear wave is influenced by the weight function, horizon, and mesh size. It is the same as the longitudinal wave. The dispersion curve is roughly the same under the three weight functions, except that the dispersion of the wave function is more pronounced under the weight functions ω1 and ω2. When the wavenumber is low, the horizon has little effect on the dispersion. When the wavenumber is large, the dispersion becomes more obvious with horizon increases. A small mesh size leads to a more obvious dispersion. However, the mesh size is relatively small when the horizon is large, and there are more material points involved in the horizon—this increases the complexity of the calculation.Figures4 and 10 show that the dispersion with different weight functions has little difference when the wavenumber is small. However, when the wavenumber is around 1000 rad/m, the group velocity becomes negative with weight functions ω1 and ω2. This indicates that the wave propagates in the inverse direction; this phenomenon does not occur in the real world. Thus, the weight function ω3 is better for simulation. Figures 5 and 11 show that the dispersion is noticeable as the horizon domain increases. The wave group velocity becomes negative when the horizon is greater than three times the mesh sizes—especially in the longitudinal wave. The peridynamic wave equation will converge to the classical Navier equation when the horizon is small. The horizon should be selected between two and three times the mesh sizes. In Figures 6 and 12, we see that the dispersion is weakened with decreasing mesh size. When the mesh size is small, there will be more material points on the horizon to complicate the calculations. The analysis above shows that the horizon should be three times the mesh sizes [24]. Therefore, in this paper, the weight function should be exp−ξ2/δ2, the horizon should be δ=3mm, and the mesh size should be h=1mm. This makes the wave simulation agree with natural reality.


## 4. Numerical Simulation of Peridynamic Wave

In order to study the ultrasonic wave propagation in an infinite media, we studied a cubic piece of maraging steel with the following properties:0.1m×0.1m×0.1m, material density ρx=8000kg/m3, bulk modulus κ=1.59×1011pa, shear modulus μ=8.2×1010pa, mesh size h=1mm, horizon radius δ=3.015h, and time step Δt=10−8s. A displacement of uy=10−10cos−760000πt was placed along y at the center of the cube. Based on the peridynamic simulation method [25], the displacement in direction y on the xoy plain can be obtained in Figure 13.When the displacement load is applied in the directiony at the center of the cube maraging steel material on the surface xoy, the load generated a longitudinal wave that propagates in the direction y. The shear wave propagates in the direction x. The other directions are superpositions of the two waveforms. By analyzing the waveform at t=5μs and t=10μs in Figure 13, we can obtain the wave propagation in the direction y. It is faster than that in the direction x, which is consistent with the observation that the velocity of the longitudinal wave is greater than that of the shear wave. Thus, the wavefront on the xoy surface is an elliptical surface. When t=15μs in Figure 13, the waveform is reflected with the free boundary leading to a waveform transform; the incident wave and the reflected wave are superimposed resulting in the wavelength shortening and an increase in the wavenumber. This phenomenon is even more pronounced when t=20μs in Figure13.Figure 13
Displacement of a material point in directiony on the surface xoy. (a) t = 5 μs. (b) t = 10 μs. (c) t = 15 μs. (d) t = 20 μs.
(a)
(b)
(c)
(d)Ultrasonic waves propagate in all directions in the cube. The two different shear waves can be obtained at the same time on the surfacexoz, and they are both displaced in the direction y. These two shear wave vibrations are shown in Figure 14.Figure 14
Vibration on thexoz plane at different times. (a) t = 5 μs. (b) t = 10 μs.
(a)
(b)The ultrasonic waves propagate in the directionsx and z with the same wave velocity when given a load in the direction y. Figure 14 shows that the vibrations in each direction on xoz are basically the same; the wavefront is a circular surface. The wave propagates in an infinite media, and the longitudinal wave propagates along the y-axis. Here, the horizontal shear wave propagates along the x-axis, and the vertical shear wave propagates along the z-axis. The wavefront in the infinite media is an ellipsoid because the longitudinal wave and shear wave have different velocities.


## 5. Conclusions

In this paper, the integral equations for the wave are given by the peridynamic theory. By studying the linearization of the wave peridynamic equation, the dispersions of the longitudinal and shear waves are given in the infinite media. The phase velocity and the group velocity of the wave are obtained, which are consistent with the wave velocity in classical mechanics. The dispersion curve is influenced by the weight function, horizon, and mesh size. When the weight function isexp−ξ2/δ2, the horizon is three times the mesh size and the dispersion curve can respond to the wave dispersion in the maraging steel. It is better than the dispersion curve with other parameters. The peridynamic theory evaluated the numerical simulation of the wave. Under initial loads at direction y, we found that the longitudinal wave was excited in the direction of y, which on the other main axis was a shear wave. The longitudinal and shear wave have different wave velocities, and thus the wavefront in the cube is an ellipsoid.



---

*Source: 9528978-2019-06-09.xml*