Construction of Dynamic Fundamental Solutions for Piezoelectric Solids

1995 ◽  
Vol 48 (11S) ◽  
pp. S222-S229 ◽  
Author(s):  
Naum Khutoryansky ◽  
Horacio Sosa

Fundamental solutions are derived within the framework of transient dynamic, three-dimensional piezoelectricity. The purpose of the article is to show alternate integral representations for such solutions. Thus, a representation over the unit sphere in accordance to a methodology based on the plane wave decomposition is provided. It is shown, however, that more efficient representations from a computational point of view can be achieved through appropriate coordinate transformations. Hence, representations of the fundamental solutions over surfaces of slowness are provided as novel alternatives to more classical approaches. The computational benefits of these new representations are displayed through a numerical example involving a transversely isotropic piezoelectric solid.

Author(s):  
Naum Khutoryansky ◽  
Horacio Sosa

Abstract Fundamental solutions for transient dynamic piezo-electricity are derived through the plane wave decomposition and represented in three alternative manners, namely over the unit sphere, over the material’s slowness surface and over a line of the latter. The computational virtues of the uni-dimensional integral representation are exposed through a numerical example concerning a transversely isotropic piezoelectric ceramic subjected to unitary impulsive applied forces.


Geophysics ◽  
1984 ◽  
Vol 49 (6) ◽  
pp. 767-770 ◽  
Author(s):  
R. F. Stöckli

The ray‐tracing problem is considered the solution to a minimum travel time problem for media where each layer may be regarded as a transversely isotropic homogeneous solid. The wave surface‐wavefront at t = 1 s, corresponding to a wave generated at the point source, associated with each layer’s anisotropy is approximated by surfaces which are not more difficult to handle, from a computational point of view, than ellipsoidal surfaces. These approximating surfaces are those used in ray‐tracing computation; a ray being a true ray approximation is thus obtained.


2020 ◽  
Vol 10 (3) ◽  
pp. 1033 ◽  
Author(s):  
Pierre Massé ◽  
Thibaut Carpentier ◽  
Olivier Warusfel ◽  
Markus Noisternig

Directional room impulse responses (DRIR) measured with spherical microphone arrays (SMA) enable the reproduction of room reverberation effects on three-dimensional surround-sound systems (e.g., Higher-Order Ambisonics) through multichannel convolution. However, such measurements inevitably contain a nondecaying noise floor that may produce an audible “infinite reverberation effect” upon convolution. If the late reverberation tail can be considered a diffuse field before reaching the noise floor, the latter may be removed and replaced with an extension of the exponentially-decaying tail synthesized as a zero-mean Gaussian noise. This has previously been shown to preserve the diffuse-field properties of the late reverberation tail when performed in the spherical harmonic domain (SHD). In this paper, we show that in the case of highly anisotropic yet incoherent late fields, the spatial symmetry of the spherical harmonics is not conducive to preserving the energy distribution of the reverberation tail. To remedy this, we propose denoising in an optimized spatial domain obtained by plane-wave decomposition (PWD), and demonstrate that this method equally preserves the incoherence of the late reverberation field.


2014 ◽  
Vol 709 ◽  
pp. 113-116 ◽  
Author(s):  
Leonid Igumnov ◽  
I.P. Маrkov ◽  
A.A. Belov

Direct boundary element method formulation for transient dynamic linear piezoelectricity is presented. Integral representations of Laplace transformed dynamic piezoelectric fundamental solutions are used. Laplace domain BEM solutions inverted in real time by the stepping method. Numerical example of transient piezoelectric analysis is presented.


2017 ◽  
Vol 743 ◽  
pp. 153-157 ◽  
Author(s):  
Leonid A. Igumnov ◽  
Ivan Markov

In this paper, the direct boundary element method in the Laplace domain is applied for the solution of three-dimensional transient dynamic problems of anisotropic elasticity in multi-connected domains. The formulation is based upon the integral representations of anisotropic dynamic fundamental solutions. As numerical example the problem of an anisotropic elastic prismatic solid with cubic cavity is investigated.


2011 ◽  
Vol 2011 ◽  
pp. 1-22 ◽  
Author(s):  
L. Godinho ◽  
P. Amado-Mendes ◽  
A. Pereira

Several numerical and analytical models have been used to study underwater acoustics problems. The most accurate and realistic models are usually based on the solution of the wave equation using a variety of methods. Here, a hybrid numerical-analytical model is proposed to address the problem of underwater sound scattering by an elastic shell structure, which is assumed to be circular and that is buried in a fluid seabed bellow a water waveguide. The interior of the shell is filled with a fluid that may have different properties from the host medium. The analysis is performed by coupling analytical solutions developed both for sound propagation in the waveguide and in the vicinity of the circular hollow pipeline. The coupling between solutions is performed using the method of fundamental solutions. This strategy allows a compact description of the propagation medium while being very accurate and highly efficient from the computational point of view.


This paper is concerned with the study of transient response of a transversely isotropic elastic half-space under internal loadings and displacement discontinuities. Governing equations corresponding to two-dimensional and three-dimensional transient wave propagation problems are solved by using Laplace–Fourier integral transforms and Laplace−Hankel integral transforms, respectively. Explicit general solutions for displacements and stresses are presented. Thereafter boundary-value problems corresponding to internal transient loadings and transient displacement discontinuities are solved for both two-dimensional and three-dimensional problems. Explicit analytical solutions for displacements and stresses corresponding to internal loadings and displacement discontinuities are presented. Solutions corresponding to arbitrary loadings and displacement discontinuities can be obtained through the application of standard analytical procedures such as integration and Fourier expansion to the fundamental solutions presented in this article. It is shown that the transient response of a medium can be accurately computed by using a combination of numerical quadrature and a numerical Laplace inversion technique for the evaluation of integrals appearing in the analytical solutions. Comparisons with existing transient solutions for isotropic materials are presented to confirm the accuracy of the present solutions. Selected numerical results for displacements and stresses due to a buried circular patch load are presented to portray some features of the response of a transversely isotropic elastic half-space. The fundamental solutions presented in this paper can be used in the analysis of a variety of transient problems encountered in disciplines such as seismology, earthquake engineering, etc. In addition these fundamental solutions appear as the kernel functions in the boundary integral equation method and in the displacement discontinuity method.


2010 ◽  
Vol 26 (4) ◽  
pp. 483-491
Author(s):  
Y. C. Shiah ◽  
W. X. Sun

ABSTRACTDue to a lack of closed-form solutions for three dimensional anisotropic bodies, the computational burden of evaluating the fundamental solutions in the boundary element method (BEM) has been a research focus over the years. In engineering practice, transversely isotropic material has gained popularity in the use of composites. As a degenerate case of the generally anisotropic material, transverse isotropy still needs to be treated separately to ease the computations. This paper aims to investigate the computational efficiency of the BEM implementations for 3D anisotropic, transversely isotropic, and isotropic bodies. For evaluating the fundamental solutions of 3D anisotropy, the explicit formulations reported in [1,2] are implemented. For treating transversely isotropic materials, numerous closed form solutions have been reported in the literature. For the present study, the formulations presented by Pan and Chou [3] are particularly employed. At the end, a numerical example is presented to compare the computational efficiency of the three cases and to demonstrate how the CPU time varies with the number of meshes.


Geophysics ◽  
1982 ◽  
Vol 47 (10) ◽  
pp. 1375-1401 ◽  
Author(s):  
Sven Treitel ◽  
P. R. Gutowski ◽  
D. E. Wagner

A point‐source seismic recording can be decomposed into a set of plane‐wave seismograms for arbitrary angles of incidence. Such plane‐wave seismograms possess an inherently simple structure that make them amenable to existing inversion methods such as predictive deconvolution. Implementation of plane‐wave decomposition (PWD) takes place in the frequency‐wavenumber domain under the assumption of radial symmetry. This version of PWD is equivalent to slant stacking if allowance is made for the customary use of linear recording arrays on the surface of a three‐dimensional medium. An imaging principle embodying both kinematic as well as dynamic characteristics allows us to perform time migration of the plane‐wave seismograms. The imaging procedure is implementable as a two‐dimensional filter whose independent variables are traveltime and angle of incidence.


Sign in / Sign up

Export Citation Format

Share Document