FG Statik und Dynamik
Refine
Year of publication
Document Type
- Scientific journal article peer-reviewed (28)
- Conference Proceeding (18)
- Conference publication not peer-reviewed (16)
- Conference publication peer-reviewed (11)
- Article (6)
- Part of a book (chapter) (3)
- Scientific journal article not peer-reviewed (2)
- Book (2)
- Book (publisher) (2)
- Doctoral thesis (2)
Way of publication
- Open Access (4)
Keywords
- Isogeometric analysis (13)
- NURBS (9)
- Isogeometric Analysis (3)
- Mortar Method (3)
- NURBS basis functions (3)
- Ratcheting (3)
- Shakedown (3)
- Stress stiffening (3)
- Cyclic strain accumulation (2)
- Dual basis functions (2)
Institute
Der Beitrag vergleicht die Genauigkeit von Reissner-Mindlin-Schalenelementen zwischen der isogeometrischen Methode und der spektralen Elemente-Methode. Während die erste durch die hohe Kontinuität zwischen den Elementen eine sehr hohe Genauigkeit in Bezug auf Anzahl der Freiheitsgrade aufweist, besticht die letztere durch eine einfachere Formulierung und eine bessere Konditionierung auch für sehr hohe Ansatzordnungen.
A spectral finite element Reissner–Mindlin shell formulation with NURBS-based geometry definition
(2024)
A curved non-isoparametric Reissner–Mindlin shell element is developed for analyzing thin-walled structures. The standard kinematic description of the element requires the calculation of the director vector. To address this demand accurately, similar to isogeometric analysis (IGA), the geometry is defined by utilization of the non-uniform rational B-splines (NURBS) imported directly from computer-aided design (CAD) files. Then, shape functions of the Legendre spectral element method (SEM) are used to interpolate the displacements. Consequently, the shell director vector and Jacobian of the transformation are calculated properly according to the presented formulation. On the other hand, in Legendre SEM combined with Gauss–Lobatto–Legendre quadrature, the integration points and the element nodes coincide. Thus, the easily computable local coordinate systems at the integration points can be used directly as nodal basis systems. A separate calculation of nodal basis systems at control points, which is the source of either complexity or error in IGA shells, is not required. Given the condition number of the stiffness matrix in the developed method, super high-order elements can also be used. Very high order p-refined elements are used in addition to h-refinement of the mesh to show the capability of higher order elements to analyze problems without mesh refinement. The validity and convergence rate of the method are investigated and verified through various cases of h- and p-refinement in challenging obstacle course problems.
An efficient mass lumping scheme for isogeometric analysis based on approximate dual basis functions
(2024)
In this contribution, we propose a new mass lumping scheme for explicit dynamics in isogeometric analysis (IGA). To this end, an element formulation based on the idea of dual functionals is developed. Non-Uniform Rational B-splines (NURBS) are applied as shape functions and their corresponding dual basis functions are applied as test functions in the variational form, where two kinds of dual basis functions are compared. The first type are approximate dual basis functions (AD) with varying degree of reproduction, resulting in banded and diagonally-dominant mass matrices. Dual basis functions derived from the inversion of the Gram matrix (IG) are the second type and already yield diagonal mass matrices. We will show that it is possible to apply the dual scheme as a transformation of the resulting system of equations based on NURBS for both – shape and test functions. Hence, it can be easily implemented into existing IGA routines and it is also promising to retain the accuracy known from similar formulations without mass lumping. Applying additional row-sum lumping to the mass matrices is either not necessary for IG or the caused loss of accuracy is lowered to a reasonable magnitude in the case of AD. Numerical examples show a significantly better approximation of the dynamic behavior for the dual lumping scheme compared to standard NURBS approaches making use of conventional row-sum lumping. In a nutshell, applying IG yields accurate numerical results but fully populated stiffness matrices occur, which are entirely unsuitable for explicit dynamic simulations, while combining AD and row-sum lumping leads to efficient dynamical computations, with respect to effort and accuracy.
In order to overcome locking effects that especially occur for lower order finite element formulations, different methods can be employed. This can be conducted using mixed formulations or adapted approximation orders, for instance. Hence, in order to tackle shear locking that is caused by non-matching interpolation degrees in the shear strain equation, an irreducible and a mixed Reissner-Mindlin plate formulation with accordingly adapted conforming discretizations are derived within the scope of this contribution. In addition, non-uniform rational B-splines (NURBS) are employed therefore, in order to benefit from the properties and refinement strategies offered by isogeometric analysis (IGA) and to achieve more accurate results. The effect of various combinations of interpolation orders on the convergence behavior and the ability to alleviate locking is investigated for both the irreducible and the mixed isogeometric plate formulation and examined for a benchmark example. This is also supplemented by investigations on the stability of the considered variants, tested by the existence of the correct number of zero-energy modes.
Für eine Lebensdauervorhersage veränderlich belasteter Tragwerke, etwa des Anlagen- und Maschinenbaus sowie des Bauingenieurwesens, werden die zyklisch akkumulierten Verzerrungen und ggf. auch die elastisch-plastischen Dehnschwingbreiten benötigt. Die Vereinfachte Fließzonentheorie (VFZT) ist eine direkte Methode, die Abschätzungen dieser und aller anderen mechanischen Größen im elastischen und im plastischen Einspielzustand liefert. Das vorliegende Buch stellt die VFZT ausführlich dar und legt Wert darauf, dass sich nicht nur Wissenschaftler, sondern auch in der Praxis tätige Ingenieure sowie Studierende höherer Semester ein Bild von den Möglichkeiten und Grenzen machen können. Zahlreiche Abbildungen und Anwendungsbeispiele unterstützen das Verständnis.
On the use of mixed basis function degrees within a convective isogeometric element formulation
(2023)
Solving linear elasticity problems using standard finite element methods, different locking phenomena can occur. In order to counteract these effects, mixed methods or formulations using higher approximation orders can be employed, for instance. As a consequence, this leads to an increased computational effort. Hence, a selective elevation of orders in decisive directions within a purely displacement-based element formulation is proposed in this contribution. Within isogeometric analysis (IGA), the geometry is discretized using non-uniform rational B-splines (NURBS), which simultaneously represent the basis functions for the analysis. Due to the fact that the geometry can be preserved exactly during analysis, this can increase the accuracy of results. In this contribution, convective basis systems that are aligned with the local geometry are employed combined with selective order elevation. The required convective basis systems are interpolated from those determined in each control point.
Employing isogeometric analysis (IGA), the geometry of structures is discretized by non-uniform rational B-splines (NURBS), which simultaneously provide the basis functions for the corresponding analysis as well. Due to the high continuity within patches, a unique local convective basis system, which corresponds to the local geometry directions of the mesh, can be defined in every point.
Locking Phenomena, which strongly affect purely displacement-based low order elements, can be counteracted by employing particular methods or higher polynomial degrees; both options increase the resulting computational effort significantly. Using order elevation only in specific directions has the potential to optimize the ratio between locking counteraction and computational costs.
In this contribution, the use of directed deformations based on convective basis systems in each control point is proposed for a displacement-based isogeometric formulation with specifically adapted orders. Therefore, distinct meshes for the interpolation of the displacements in each direction are generated based on the initial geometry. Subsequently, the order of every mesh is elevated in only one direction. Consequently, different possibilities for the combination of order elevations have to be examined.
This procedure is conducted for a two-dimensional linear elasticity problem. The benefit of a convective formulation with direction dependent degrees is shown by a comparison to a standard isogeometric formulation. Furthermore, the impact of the directions selected for order elevation on the accuracy of the results is investigated.