scholarly journals MATHEMATICAL MODELING OF ELASTIC WAVE PROPAGATION IN A BODY WITH A DEFECT

2021 ◽  
Vol 3 (1) ◽  
pp. 67-71
Author(s):  
Olena Stankevych ◽  

The solution of the dynamic problem of calculation the wave field of displacements on the surface of an elastic half-space caused by the opening of an internal crack under the action of torsional forces is presented. Based on the solutions of the boundary integral equations, the nature of the change in the amplitude-frequency characteristics of elastic oscillations on the surface of a rigid body depending on the size of the defect is shown.

2021 ◽  
Vol 3 (1) ◽  
pp. 45-50
Author(s):  
Olena Stankevych ◽  
◽  
Nazar Stankevych ◽  

The dynamic problem of the displacement field in an elastic half-space caused by the time-steady displacement of the surfaces of the system of disc-shaped coplanar cracks is solved. The solutions are obtained by the method of boundary integral equations. The dependences of elastic displacements on the surface of the half-space on the wave number, the number of defects and the depths of their occurrence are constructed.


Geophysics ◽  
1987 ◽  
Vol 52 (6) ◽  
pp. 765-771 ◽  
Author(s):  
B. Kummer ◽  
A. Behle ◽  
F. Dorau

We have constructed a hybrid scheme for elastic‐wave propagation in two‐dimensional laterally inhomogeneous media. The algorithm is based on a combination of finite‐difference techniques and the boundary integral equation method. It involves a dedicated application of each of the two methods to specific domains of the model structure; finite‐difference techniques are applied to calculate a set of boundary values (wave field and stress field) in the vicinity of the heterogeneous domain. The continuation of the near‐field response is then calculated by means of the boundary integral equation method. In a numerical example, the hybrid method has been applied to calculate a plane‐wave response for an elastic lens embedded in a homogeneous environment. The example shows that the hybrid scheme enables more efficient modeling, with the same accuracy, than with pure finite‐difference calculations.


2020 ◽  
Vol 36 (6) ◽  
pp. 749-761
Author(s):  
Y. -Y. Ko

ABSTRACTWhen the Symmetric Galerkin boundary element method (SGBEM) based on full-space elastostatic fundamental solutions is used to solve Neumann problems, the displacement solution cannot be uniquely determined because of the inevitable rigid-body-motion terms involved. Several methods that have been used to remove the non-uniqueness, including additional point support, eigen decomposition, regularization of a singular system and modified boundary integral equations, were introduced to amend SGBEM, and were verified to eliminate the rigid body motions in the solutions of full-space exterior Neumann problems. Because half-space problems are common in geotechnical engineering practice and they are usually Neumann problems, typical half-space problems were also analyzed using the amended SGBEM with a truncated free surface mesh. However, various levels of errors showed for all the methods of removing non-uniqueness investigated. Among them, the modified boundary integral equations based on the Fredholm’s theory is relatively preferable for its accurate results inside and near the loaded area, especially where the deformation varies significantly.


2017 ◽  
Vol 743 ◽  
pp. 158-161
Author(s):  
Andrey Petrov ◽  
Sergey Aizikovich ◽  
Leonid A. Igumnov

Problems of wave propagation in poroelastic bodies and media are considered. The behavior of the poroelastic medium is described by Biot theory for partially saturated material. Mathematical model is written in term of five basic functions – elastic skeleton displacements, pore water pressure and pore air pressure. Boundary element method (BEM) is used with step method of numerical inversion of Laplace transform to obtain the solution. Research is based on direct boundary integral equation of three-dimensional isotropic linear theory of poroelasticity. Green’s matrices and, based on it, boundary integral equations are written for basic differential equations in partial derivatives. Discrete analogue are obtained by applying the collocation method to a regularized boundary integral equation. To approximate the boundary consider its decomposition to a set of quadrangular and triangular 8-node biquadratic elements, where triangular elements are treated as singular quadrangular. Every element is mapped to a reference one. Interpolation nodes for boundary unknowns are a subset of geometrical boundary-element grid nodes. Local approximation follows the Goldshteyn’s generalized displacement-stress matched model: generalized boundary displacements are approximated by bilinear elements whereas generalized tractions are approximated by constant. Integrals in discretized boundary integral equations are calculated using Gaussian quadrature in combination with singularity decreasing and eliminating algorithms.


2019 ◽  
Vol 10 (1) ◽  
pp. 182-198
Author(s):  
A. Aimi ◽  
L. Desiderio ◽  
M. Diligenti ◽  
C. Guardasoni

Abstract Starting from a recently developed energetic space-time weak formulation of the Boundary Integral Equations related to scalar wave propagation problems, in this paper we focus for the first time on the 2D elastodynamic extension of the above wave propagation analysis. In particular, we consider elastodynamic scattering problems by open arcs, with vanishing initial and Dirichlet boundary conditions and we assess the efficiency and accuracy of the proposed method, on the basis of numerical results obtained for benchmark problems having available analytical solution.


Sign in / Sign up

Export Citation Format

Share Document