Parallel Volume Integral Equation Method for Two-Dimensional Elastodynamic Analysis

2019 ◽  
Vol 16 (06) ◽  
pp. 1840025
Author(s):  
Jungki Lee ◽  
Hogwan Jeong

The parallel volume integral equation method (PVIEM) is applied for the analysis of two-dimensional elastic wave scattering problems in an unbounded isotropic solid containing various types of multiple multilayered anisotropic inclusions. It should be noted that the volume integral equation method (VIEM) does not require the use of the Green’s function for the anisotropic inclusion — only the Green’s function for the unbounded isotropic matrix is needed. A detailed analysis of the SH wave scattering problem is presented for various types of multiple multilayered orthotropic inclusions. Numerical results are presented for the elastic fields at the interfaces for square and hexagonal packing arrays of various types of multilayered orthotropic inclusions in a broad frequency range of practical interest. Standard parallel programming was used to speed up computation in the VIEM. The PVIEM enables us to investigate the effects of single/multiple scattering, fiber packing type, fiber volume fraction, single/multiple layer(s), multilayer’s shapes and geometry, isotropy/anisotropy, and softness/hardness of various types of multiple multilayered anisotropic inclusions on displacements and stresses at the interfaces of the inclusions and far-field scattering patterns. Also, powerful capabilities of the PVIEM for the analysis of general two-dimensional multiple scattering problems are investigated.

2015 ◽  
Vol 2015 ◽  
pp. 1-48 ◽  
Author(s):  
Jungki Lee ◽  
Hyechang Lee ◽  
Hogwan Jeong

The parallel volume integral equation method (PVIEM) is applied for the analysis of elastic wave scattering problems in an unbounded isotropic solid containing multiple multilayered anisotropic elliptical inclusions. This recently developed numerical method does not require the use of Green’s function for the multilayered anisotropic inclusions; only Green’s function for the unbounded isotropic matrix is needed. This method can also be applied to solve general two- and three-dimensional elastodynamic problems involving inhomogeneous and/or multilayered anisotropic inclusions whose shape and number are arbitrary. A detailed analysis of the SH wave scattering is presented for multiple triple-layered orthotropic elliptical inclusions. Numerical results are presented for the displacement fields at the interfaces for square and hexagonal packing arrays of triple-layered elliptical inclusions in a broad frequency range of practical interest. It is necessary to use standard parallel programming, such as MPI (message passing interface), to speed up computation in the volume integral equation method (VIEM). Parallel volume integral equation method as a pioneer of numerical analysis enables us to investigate the effects of single/multiple scattering, fiber packing type, fiber volume fraction, single/multiple layer(s), multilayer’s shape and geometry, isotropy/anisotropy, and softness/hardness of the multiple multilayered anisotropic elliptical inclusions on displacements at the interfaces of the inclusions.


2018 ◽  
Vol 35 (1) ◽  
pp. 432-476 ◽  
Author(s):  
Jungki Lee ◽  
Hogwan Jeong

Purpose The purpose of this paper is to calculate near field and far field scattering of SH waves by multiple multilayered anisotropic circular inclusions using parallel volume integral equation method (PVIEM) quantitatively. Design/methodology/approach The PVIEM is applied for the analysis of elastic wave scattering problems in an unbounded solid containing multiple multilayered anisotropic circular inclusions. It should be noted that this numerical method does not require the use of the Green’s function for the inclusion – only the Green’s function for the unbounded isotropic matrix is needed. This method can also be applied to solve general elastodynamic problems involving inhomogeneous and/or anisotropic inclusions whose shape and number are arbitrary. Findings A detailed analysis of the SH wave scattering problem is presented for multiple multilayered orthotropic circular inclusions. Numerical results are presented for the displacement fields at the interfaces and the far field scattering patterns for square and hexagonal packing arrays of multilayered circular inclusions in a broad frequency range of practical interest. Originality/value To the best of the authors’ knowledge, the solution for scattering of SH waves by multiple multilayered anisotropic circular inclusions in an unbounded isotropic matrix is not currently available in the literature. However, in this paper, calculation of displacements on interfaces and far field scattering patterns of multiple multilayered anisotropic circular inclusions using PVIEM as a pioneer of numerical modeling enables us to investigate the effects of single/multiple scattering, fiber packing type, fiber volume fraction, single/multiple layer(s), the multilayer’s geometry, isotropy/anisotropy and softness/hardness.


Author(s):  
Jungki Lee

A number of analytical techniques are available for the stress analysis of inclusion problems when the geometries of inclusions are simple (e.g., cylindrical, spherical or ellipsoidal) and when they are well separated [9, 41, 52]. However, these approaches cannot be applied to more general problems where the inclusions are anisotropic and arbitrary in shape, particularly when their concentration is high. Thus, stress analysis of heterogeneous solids or analysis of elastic wave scattering problems in heterogeneous solids often requires the use of numerical techniques based on either the finite element method (FEM) or the boundary integral equation method (BIEM). However, these methods become problematic when dealing with elastostatic problems or elastic wave scattering problems in unbounded media containing anisotropic and/or heterogeneous inclusions of arbitrary shapes. It has been demonstrated that the volume integral equation method (VIEM) can overcome such difficulties in solving a large class of inclusion problems [6,10,20,21,28–30]. One advantage of the VIEM over the BIEM is that it does not require the use of Green’s functions for anisotropic inclusions. Since the elastodynamic Green’s functions for anisotropic media are extremely difficult to calculate, the VIEM offers a clear advantage over the BIEM. In addition, the VIEM is not sensitive to the geometry or concentration of the inclusions. Moreover, in contrast to the finite element method, where the full domain needs to be discretized, the VIEM requires discretization of the inclusions only.


Sign in / Sign up

Export Citation Format

Share Document