Benjamin Marussig

Benjamin Marussig, Jürgen Zechner, Gernot Beer et al.

An isogeometric boundary element method for problems in elasticity is presented, which is based on an independent approximation for the geometry, traction and displacement field. This enables a flexible choice of refinement strategies, permits an efficient evaluation of geometry related information, a mixed collocation scheme which deals with discontinuous tractions along non-smooth boundaries and a significant reduction of the right hand side of the system of equations for common boundary conditions. All these benefits are achieved without any loss of accuracy compared to conventional isogeometric formulations. The system matrices are approximated by means of hierarchical matrices to reduce the computational complexity for large scale analysis. For the required geometrical bisection of the domain, a strategy for the evaluation of bounding boxes containing the supports of NURBS basis functions is presented. The versatility and accuracy of the proposed methodology is demonstrated by convergence studies showing optimal rates and real world examples in two and three dimensions.

Gernot Beer, Benjamin Marussig, Jürgen Zechner et al.

In this work a novel approach is presented for the isogeometric Boundary Element analysis of domains that contain inclusions with different elastic properties than the ones used for computing the fundamental solutions. In addition the inclusion may exhibit inelastic material behavior. In this paper only plane stress/strain problems are considered. In our approach the geometry of the inclusion is described using NURBS basis functions. The advantage over currently used methods is that no discretization into cells is required in order to evaluate the arising volume integrals. The other difference to current approaches is that Kernels of lower singularity are used in the domain term. The implementation is verified on simple finite and infinite domain examples with various boundary conditions. Finally a practical application in geomechanics is presented.

Gernot Beer, Benjamin Marussig, Jürgen Zechner

In this work a novel method for the analysis with trimmed CAD surfaces is presented. The method involves an additional mapping step and the attraction stems from its sim- plicity and ease of implementation into existing Finite Element (FEM) or Boundary Element (BEM) software. The method is first verified with classical test examples in structural mechanics. Then two practical applications are presented one using the FEM, the other the BEM, that show the applicability of the method.

Guilherme Henrique Teixeira, Nepomuk Krenn, Peter Gangl et al.

Topology optimization is a valuable tool in engineering, facilitating the design of optimized structures. However, topological changes often require a remeshing step, which can become challenging. In this work, we propose an isogeometric approach to topology optimization driven by topological derivatives. The combination of a level-set method together with an immersed isogeometric framework allows seamless geometry updates without the necessity of remeshing. At the same time, topological derivatives provide topological modifications without the need to define initial holes [7]. We investigate the influence of higher-degree basis functions in both the level-set representation and the approximation of the solution. Two numerical examples demonstrate the proposed approach, showing that employing higher-degree basis functions for approximating the solution improves accuracy, while linear basis functions remain sufficient for the level-set function representation.

Benjamin Marussig

Trimming is a core technique in geometric modeling. Unfortunately, the resulting objects do not take the requirements of numerical simulations into account and yield various problems. This paper outlines principal issues of trimmed models and highlights different analysis-suitable strategies to address them. It is discussed that these concepts not only provide important computational tools for isogeometric analysis, but can also improve the treatment of trimmed models in a design context.

Teoman Toprak, Michael Loibl, Guilherme H. Teixeira et al.

In the field of scientific computing, one often finds several alternative software packages (with open or closed source code) for solving a specific problem. These packages sometimes even use alternative methodological approaches, e.g., different numerical discretizations. If one decides to use one of these packages, it is often not clear which one is the best choice. To make an informed decision, it is necessary to measure the performance of the alternative software packages for a suitable set of test problems, i.e. to set up a benchmark. However, setting up benchmarks ad-hoc can become overwhelming as the parameter space expands rapidly. Very often, the design of the benchmark is also not fully set at the start of some project. For instance, adding new libraries, adapting metrics, or introducing new benchmark cases during the project can significantly increase complexity and necessitate laborious re-evaluation of previous results. This paper presents a proven approach that utilizes established Continuous Integration tools and practices to achieve high automation of benchmark execution and reporting. Our use case is the numerical integration (quadrature) on arbitrary domains, which are bounded by implicitly or parametrically defined curves or surfaces in 2D or 3D.

Benjamin Marussig, Jürgen Zechner, Gernot Beer et al.

We explore extended B-splines as a stable basis for isogeometric analysis with trimmed parameter spaces. The stabilization is accomplished by an appropriate substitution of B-splines that may lead to ill-conditioned system matrices. The construction for non-uniform knot vectors is presented. The properties of extended B-splines are examined in the context of interpolation, potential, and linear elasticity problems and excellent results are attained. The analysis is performed by an isogeometric boundary element formulation using collocation. It is argued that extended B-splines provide a flexible and simple stabilization scheme which ideally suits the isogeometric paradigm.

Jürgen Zechner, Benjamin Marussig, Gernot Beer et al.

In this paper the isogeometric Nyström method is presented. It's outstanding features are: it allows the analysis of domains described by many different geometrical mapping methods in computer aided geometric design and it requires only pointwise function evaluations just like isogeometric collocation methods. The analysis of the computational domain is carried out by means of boundary integral equations, therefor only the boundary representation is required. The method is thoroughly integrated into the isogeometric framework. For example, the regularization of the arising singular integrals performed with local correction as well as the interpolation of the pointwise existing results are carried out by means of Bezier elements. The presented isogeometric Nyström method is applied to practical problems solved by the Laplace and the Lame-Navier equation. Numerical tests show higher order convergence in two and three dimensions. It is concluded that the presented approach provides a simple and flexible alternative to currently used methods for solving boundary integral equations, but has some limitations.