This paper presents a mode-tracing approach for elastic guided waves based on analytically computed derivatives and includes a study of interesting phenomena in the dispersion curve representation. Numerical simulation is done by means of the Scaled Boundary Finite Element Method (SBFEM). Two approaches are used to identify the characteristics of the resulting wave modes: Taylor approximation and Padé approximation. Higher order differentials of the underlying eigenvalue problem are the basis for these approaches. Remarkable phenomena in potentially critical frequency regions are identified and the tracing approach is adapted to these regions. Additionally, a stabilization of the solution process is suggested.
When modeling the propagation of elastic guided waves in plates or cylinders, Finite Element based numerical methods such as the Scaled Boundary Finite Element Method (SBFEM) or the Semi-Analytical Finite Element (SAFE) Method lead to an eigenvalue problem to be solved at each frequency. For the particular case of shear horizontal modes in a homogeneous plate or torsional modes in a homogeneous cylinder, the problem can be drastically simplified. The eigenvalues become simple functions of the frequency, while the eigenvectors are constant. The current contribution discusses how this behavior is represented in the numerical formulation and derives the expressions for the eigenvalues and eigenvectors as well as the dynamic stiffness matrix of infinite elastic waveguides.
In this contribution, we present an efficient approach for the transient and time-causal modeling of guided waves in viscoelastic cylindrical waveguides in the context of ultrasonic material characterization. We use the scaled boundary finite element method (SBFEM) for efficient computation of the phase velocity dispersion. Regarding the viscoelastic behavior of the materials under consideration, we propose a decomposition approach that considers the real-valued frequency dependence of the (visco-)elastic moduli and, separately, of their attenuation. The modal expansion approach is utilized to take the transmitting and receiving transducers into account and to propagate the excited waveguide modes through a waveguide of finite length. The effectiveness of the proposed simulation model is shown by comparison with a standard transient FEM simulation as well as simulation results based on the exact solution of the complex-valued viscoelastic guided wave problem. Two material models are discussed, namely the fractional Zener model and the anti-Zener model; we re-interpret the latter in terms of the Rayleigh damping model. Measurements are taken on a polypropylene sample and the proposed transient simulation model is used for inverse material characterization. The extracted material properties may then be used in computer-aided design of ultrasonic systems.
In this paper, an approach is presented to model the propagation of elastic waves and their interaction with defects in plate structures. The formulation is based on the Scaled Boundary Finite Element Method (SBFEM), a general semi-analytical method requiring the discretization of boundaries only. For a homogeneous finite or infinite plate section, only the through-thickness direction of the plate is discretized. To describe a defect, the full boundary of a short plate section of irregular shape is discretized. High-order spectral elements are employed for the discretization. The formulation for infinite plates can model the transmission into an unbounded domain exactly. Results are compared with conventional Finite Element Analyses in both time domain and frequency domain. The presented approach allows for the simulation of complex reflection and scattering phenomena using a very small number of degrees of freedom while the mesh consists of one-dimensional elements only.
An approach is presented to model elastic waveguides of arbitrary cross-section coupled to infinite solid media. The formulation is based on the scaled boundary-finite element method. The surrounding medium is approximately accounted for by a dashpot boundary condition derived from the acoustic impedances of the infinite medium. It is discussed under which circumstances this approximation leads to sufficiently accurate results. Computational costs are very low, since the surrounding medium does not require discretization and the number of degrees of freedom on the cross-section is significantly reduced by utilizing higher-order spectral elements.
In this paper, a method to determine the complex dispersion relations of axially symmetric guided waves in cylindrical structures is presented as an alternative to the currently established numerical procedures. The method is based on a spectral decomposition into eigenfunctions of the Laplace operator on the cross-section of the waveguide. This translates the calculation of real or complex wave numbers at a given frequency into solving an eigenvalue problem. Cylindrical rods and plates are treated as the asymptotic cases of cylindrical structures and used to generalize the method to the case of hollow cylinders. The presented method is superior to direct root-finding algorithms in the sense that no initial guess values are needed to determine the complex wave numbers and that neither starting at low frequencies nor subsequent mode tracking is required. The results obtained with this method are shown to be reasonably close to those calculated by other means and an estimate for the achievable accuracy is given.
The simulation of guided waves in plate structures and cylinders coupled to infinite fluids is addressed. The approach is based on the Scaled Boundary Finite Element Method. Only a straight line is discretized that represents the through-thickness direction or the radial direction. The surrounding fluid is accounted for by employing a damping boundary condition that is based on the analytical description of the radiation impedance. Since the radiation impedance is a function of the wavenumber in the waveguide, an iterative solution procedure is applied. The algorithm is highly efficient while the results are in agreement with the Global Matrix Method.
This paper addresses the computation of dispersion curves and mode shapes of elastic guided waves in axisymmetric waveguides. The approach is based on a Scaled Boundary Finite Element formulation, that has previously been presented for plate structures and general three-dimensional waveguides with complex cross-section. The formulation leads to a Hamiltonian eigenvalue problem for the computation of wavenumbers and displacement amplitudes, that can be solved very efficiently. In the axisymmetric representation, only the radial direction in a cylindrical coordinate system has to be discretized, while the circumferential direction as well as the direction of propagation are described analytically. It is demonstrated, how the computational costs can drastically be reduced by employing spectral elements of extremely high order. Additionally, an alternative formulation is presented, that leads to real coefficient matrices. It is discussed, how these two approaches affect the computational efficiency, depending on the elasticity matrix. In the case of solid cylinders, the singularity of the governing equations that occurs in the center of the cross-section is avoided by changing the quadrature scheme. Numerical examples show the applicability of the approach to homogeneous as well as layered structures with isotropic or anisotropic material behavior.
In this paper we propose an algorithm to compute specific parts of the dispersion curves for elastic waveguides. The formulation is based on an axisymmetric representation of the Scaled Boundary Finite Element Method, where the wavenumbers of propagating modes are obtained as solutions of a Hamiltonian eigenvalue problem. The novel solution procedure involves tracing selected modes over a given frequency range and computing the corresponding solutions by means of inverse iteration. The resulting algorithm is applied in the context of material characterization, where the efficiency of the computation is crucial.
Numerical simulation of ultrasonic guided waves using the scaled boundary finite element method
(2012)
The formulation of the Scaled Boundary Finite Element Method is applied for the computation of dispersion properties of ultrasonic guided waves. The cross-section of the waveguide is discretized in the Finite Element sense, while the direction of propagation is described analytically. A standard eigenvalue problem is derived to compute the wave numbers of propagating modes. This paper focuses on cylindrical waveguides, where only a straight line has to be discretized. Higher-order elements are utilized for the discretization. As examples, dispersion curves are computed for a homogeneous pipe and a layered cylinder.
In this paper, a numerical approach for the computation of dispersion relations for three-dimensional waveguides with arbitrary cross-section is proposed. The formulation is based on the Scaled Boundary Finite Element Method (SBFEM). It is an extension of the approach previously derived for plate structures. It is shown that the wavenumbers of guided waves in a waveguide can be obtained as the eigenvalues of the Z matrix, which is well known in the SBFEM. The Hamiltonian properties of this matrix are utilized to derive an efficient way to compute the group velocities of propagating waves as eigenvalue derivatives. The cross-section of the waveguide is discretized using higher-order spectral elements. It is discussed in detail how symmetry axes can be utilized to reduce computational costs. In order to sort the solutions at different frequencies, a mode-tracking algorithm is proposed, based on the Padé expansion.
In this paper a numerical approach, based on the Scaled Boundary Finite Element Method (SBFEM), is described to obtain dispersion relations for propagating modes in wave guides. While the formulation is developed for plate structures, it can easily be extended to wave guides with arbitrary cross-section. The cross-section is discretized in the Finite Element sense while all equations remain analytical in the direction of propagation. The wave numbers of all propagating modes are obtained as the solutions of a standard eigenvalue problem. The group velocities can be calculated accurately as the eigenvalue derivatives. The use of higher-order elements drastically increases the efficiency and accuracy of the computation. This approach can be used for wave guides with arbitrary distribution of material parameters.
Der Einsatz geführter Wellen für die
zerstörungsfreie Prüfung mit Ultraschall eröffnet neue Möglichkeiten,
räumlich ausgedehnte Bauteile mit begrenzter Zugänglichkeit
auf ihre Integrität zu prüfen und gewinnt daher
zunehmend an Bedeutung. Dieser Artikel behandelt die
physikalischen Grundlagen der Schallausbreitung. Deren Verständnis
bildet die Grundlage für die Entwicklung geeigneter
Prüfsysteme. An Beispielen werden verschiedene Möglichkeiten
zur Simulation der Schallausbereitung vorgestellt. Aktuelle
Lösungsansätze zur Prüfung von plattenförmigen Strukturen
und von Rohrleitungen werden beschrieben, wobei besonders
auf die Sensortechnik und die speziellen Anforderungen
an die Prüfhardware eingegangen wird. -----------------------------------------------------------------------------------------------------------------------------------
Guided waves are widely used for non-destructive
testing using ultrasound. Recently, the method has become increasingly
important for integrity tests of spatially extended
components with limited accessibility. This article discusses the
basic physics of the sound propagation of guided waves. Their
understanding forms the basis for the successful development
of adapted inspection systems. Examples for simulating the
wave propagation using different methods are presented. Current
approaches for the inspection of plate-like structures and
pipelines are described with focus on sensor technology and the
specific requirement on the test hardware.