Search

Eigenvalues of parameter-dependent quadratic eigenvalue problems form eigencurves. The critical points on these curves, where the derivative vanishes, are of practical interest. A particular example is found in the dispersion curves of elastic waveguides, where such points are called zero-group-velocity (ZGV) points. Recently, it was revealed that the problem of computing ZGV points can be modeled as a multiparameter eigenvalue problem (MEP), and several numerical methods were devised. Due to their complexity, these methods are feasible only for problems involving small matrices. In this paper, we improve the efficiency of these methods by exploiting the link to the Sylvester equation. This approach enables the computation of ZGV points for problems with much larger matrices, such as multi-layered plates and three-dimensional structures of complex cross-sections., Comment: 18 pages

We present a semi-analytical approach to compute quasi-guided elastic wave modes in horizontally layered structures radiating into unbounded fluid or solid media. This problem is of relevance, e.g., for the simulation of guided ultrasound in embedded plate structures or seismic waves in soil layers over an elastic half-space. We employ a semi-analytical formulation to describe the layers, thus discretizing the thickness direction by means of finite elements. For a free layer, this technique leads to a well-known quadratic eigenvalue problem for the mode shapes and corresponding horizontal wavenumbers. Rigorously incorporating the coupling conditions to account for the adjacent half-spaces gives rise to additional terms that are nonlinear in the wavenumber. We show that the resulting nonlinear eigenvalue problem can be cast in the form of a multiparameter eigenvalue problem whose solutions represent the wave numbers in the plate and in the half-spaces. The multiparameter eigenvalue problem is solved numerically using recently developed algorithms.

The scaled boundary finite element method (SBFEM) has recently been employed as an efficient means to model three-dimensional structures, in particular when the geometry is provided as a voxel-based image. To this end, an octree decomposition of the computational domain is deployed and each cubic cell is treated as an SBFEM subdomain. The surfaces of each subdomain are discretized in the finite element sense. We improve on this idea by combining the semi-analytical concept of the SBFEM with certain transition elements on the subdomains' surfaces. Thus, we avoid the triangulation of surfaces employed in previous works and consequently reduce the number of surface elements and degrees of freedom. In addition, these discretizations allow coupling elements of arbitrary order such that local p-refinement can be achieved straightforwardly.

The dispersion curves of (elastic) waveguides frequently exhibit crossings and osculations (also known as veering, repulsion, or avoided crossing). Osculations are regions in the dispersion diagram where curves approach each other arbitrarily closely without ever crossing before veering apart. In semi-analytical (undamped) waveguide models, dispersion curves are obtained as solutions to discretized parameterized Hermitian eigenvalue problems. In the mathematical literature, it is known that such eigencurves can exhibit crossing points only if the corresponding matrix flow (parameter-dependent matrix) is uniformly decomposable. We discuss the implications for the solution of the waveguide problem. In particular, we make use of a simple algorithm recently suggested in the literature for decomposing matrix flows. We also employ a method for mode tracing based on approximating the eigenvalue problem for individual modes by an ordinary differential equation that can be solved by standard procedures., Comment: accepted version, published in Ultrasonics

The Reynolds equation, combined with the Elrod algorithm for including the effect of cavitation, resembles a nonlinear convection-diffusion-reaction (CDR) equation. Its solution by finite elements is prone to oscillations in convection-dominated regions, which are present whenever cavitation occurs. We propose a stabilized finite-element method that is based on the variational multiscale method and exploits the concept of orthogonal subgrid scales. We demonstrate that this approach only requires one additional term in the weak form to obtain a stable method that converges optimally when performing mesh refinement., Comment: accepted version published in CMAME

Dispersion curves of elastic waveguides exhibit points where the group velocity vanishes while the wavenumber remains finite. These are the so-called zero-group-velocity (ZGV) points. As the elastodynamic energy at these points remains confined close to the source, they are of practical interest for nondestructive testing and quantitative characterization of structures. These applications rely on the correct prediction of the ZGV points. In this contribution, we first model the ZGV resonances in anisotropic plates based on the appearance of an additional modal solution. The resulting governing equation is interpreted as a two-parameter eigenvalue problem. Subsequently, we present three complementary numerical procedures capable of computing ZGV points in arbitrary nondissipative elastic waveguides in the conventional sense that their axial power flux vanishes. The first method is globally convergent and guarantees to find all ZGV points but can only be used for small problems. The second procedure is a very fast, generally-applicable, Newton-type iteration that is locally convergent and requires initial guesses. The third method combines both kinds of approaches and yields a procedure that is applicable to large problems, does not require initial guesses and is likely to find all ZGV points. The algorithms are implemented in "GEW ZGV computation" (doi: 10.5281/zenodo.7537442)., Comment: Added Appendices A, B and D. Some clarifications in the text

Typical areas of application of explicit dynamics are impact, crash test, and most importantly, wave propagation simulations. Due to the numerically highly demanding nature of these problems, efficient automatic mesh generators and transient solvers are required. To this end, a parallel explicit solver exploiting the advantages of balanced octree meshes is introduced. To avoid the hanging nodes problem encountered in standard finite element analysis (FEA), the scaled boundary finite element method (SBFEM) is deployed as a spatial discretization scheme. Consequently, arbitrarily shaped star-convex polyhedral elements are straightforwardly generated. Considering the scaling and transformation of octree cells, the stiffness and mass matrices of a limited number of unique cell patterns are pre-computed. A recently proposed mass lumping technique is extended to 3D yielding a well-conditioned diagonal mass matrix. This enables us to leverage the advantages of explicit time integrator, i.e., it is possible to efficiently compute the nodal displacements without the need for solving a system of linear equations. We implement the proposed scheme together with a central difference method (CDM) in a distributed computing environment. The performance of our parallel explicit solver is evaluated by means of several numerical benchmark examples, including complex geometries and various practical applications. A significant speedup is observed for these examples with up to one billion of degrees of freedom and running on up to 16,384 computing cores.

Advanced transition elements are of utmost importance in many applications of the finite element method (FEM) where a local mesh refinement is required. Considering problems that exhibit singularities in the solution, an adaptive hp-refinement procedure must be applied. Even today, this is a very demanding task especially if only quadrilateral/hexahedral elements are deployed and consequently the hanging nodes problem is encountered. These element types, are, however, favored in computational mechanics due to the improved accuracy compared to triangular/tetrahedral elements. Therefore, we propose a compatible transition element - xNy-element - which provides the capability of coupling different element types. The adjacent elements can exhibit different element sizes, shape function types, and polynomial orders. Thus, it is possible to combine independently refined h- and p-meshes. The approach is based on the transfinite mapping concept and constitutes an extension/generalization of the pNh-element concept. By means of several numerical examples, the convergence behavior is investigated in detail, and the asymptotic rates of convergence are determined numerically. Overall, it is found that the proposed approach provides very promising results for local mesh refinement procedures., Comment: 51 pages, 44 figures, 4 tables