FG Statik und Dynamik
Refine
Year of publication
Document Type
- Scientific journal article peer-reviewed (26)
- Conference Proceeding (18)
- Conference publication not peer-reviewed (15)
- 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 (3)
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
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.
Cyclic loading may cause elastic-plastic strains to accumulate if a ratcheting mechanism is present. After a number of cycles, strain accumulation may cease so that a state of either elastic or plastic shakedown is reached. Determination of strains and other quantities in the state of shakedown (post-shakedown quantities) by means of incremental analyses is costly. Direct methods such as the Simplified Theory of Plastic Zones (STPZ) aim at providing estimates of the post-shakedown quantities, bypassing cycle-by-cycle analyses. If geometric effects such as stress stiffening play a role, determination of accumulated strains even becomes more complicated. The STPZ has been further developed in order to account for the combination of plasticity and stress stiffening with respect to elastic shakedown. The theory is described and illustrated using examples such as a pipe bend subjected to cyclic in-plane bending. The implications of cyclic as opposed to constant stress stiffness are discussed. The effect of stress stiffening on ratcheting interaction diagrams (RID), separating regions of elastic and plastic shakedown in the space of loading parameters, is discussed.
The behavior of elastic-plastic structures under cyclic loading can be determined by incremental elastic-plastic analyses where a given load histogram is analysed cycle-by-cycle until shakedown is achieved. Many cycles may be required for this, if a ratchet mechanism is causing plastic strains to accumulate in each cycle until either elastic or plastic shakedown is achieved. The complexity of the elastic-plastic response of a structure is further increased if geometric effects are present. It may be very costly to get the accumulated strains and strain ranges in the state of shakedown by incremental analyses so that simplified or direct methods have been developed as an alternative. The Simplified Theory of Plastic Zones (STPZ) has proven itself for estimating the required quantities in the state of elastic and plastic shakedown within the framework of the 1st order theory, i.e., if the equilibrium conditions are satisfied for the undeformed structure. It is shown in this paper, how the STPZ can be expanded to capture 2nd order effects introduced by the equilibrium of the deformed structure. Some examples are used to demonstrate its applicability and the quality of the results, but also its limitations with respect to determining the accumulated strains and elastic-plastic strain ranges.
The main feature of isogeometric analysis is the usage of a coherent geometry description for design and analysis. In most cases, Non-Uniform Rational B-splines (NURBS) are used within the frame of the finite element method, which combines the most common geometrical modeling concept with the most common analysis method in structural mechanics. NURBS are a versatile tool for geometric modeling, and in order to define complex geometric structures, a multitude of tensor-product NURBS patches is required. Typical software in Computer-Aided design can manage smoothness requirements across common interfaces of neighboring patches. For the finite element method, a different type of coupling across interfaces must be achieved. Due to the nature of tensor product NURBS, the simple concept of mutual refinement and subsequent coupling by shared degrees of freedom is prohibitively costly and, in some cases, even not possible. A multitude of coupling methods has been proposed over the last years. The most common concepts are known as mortar methods. In particular, the dual mortar method has been shown to yield very efficient computations. A recent paper by the authors has provided an isogeometric mortar method with mathematically proven optimal convergence of the stress errors over the entire domain. We use dual basis functions, which have support only on one interface and avoid interrelations between different interfaces. Models with a large number of intersecting interfaces can be handled. However, the basis functions have full support on the interfaces. In our current contribution, we propose the use of approximate dual basis functions with the advantage of having local support on the interfaces. These functions fulfill the duality only in an approximate way, but still guarantee the optimal degree for the convergence of the mortar method. Since the duality is not fulfilled, an additional lumping of the mortar matrix is introduced. The error of this lumping can be analyzed mathematically and is not significant in comparison to the global approximation error of the finite element method. The use of the approximate dual basis functions restores the local support of basis functions along the interface while the convergence properties remain intact. Numerical examples show the convergence behavior for simple and complex models.
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.
Many pressure vessel and piping components have to withstand high internal pressures and are therefore thick-walled so that geometric effects such as stress stiffening need not be accounted for. However, thin-walled, or moderately thick structures may be sensitive to these effects. Design Codes such as the ASME Boiler and Pressure Vessel Code usually provide little guidance on when they are to be accounted for. In general, the effects of stress stiffening are difficult to estimate even for experienced engineers and can only be estimated by detailed finite element analyses. In the opinion of the authors, this effect deserves more attention. This is particularly true for simplified elastic-plastic methods for fatigue and ratcheting assessment of structures subjected to cyclic loading.
Cyclic loading may cause elastic-plastic strains to accumulate if a ratcheting mechanism is present. After a number of cycles, strain accumulation may cease so that a state of either elastic or plastic shakedown is reached. Determination of accumulated strains, strain ranges and other quantities in the state of shakedown (post-shakedown quantities) by means of incremental analyses is costly, in particular if geometric effects such as stress stiffening play a role. Direct methods aim at providing estimates of the post-shakedown quantities, bypassing cycle-by-cycle analyses. These methods claim to deliver the post-shakedown quantities with high accuracy and low computational effort.
The Simplified Theory of Plastic Zones can account for the combination of plasticity and stress stiffening. The theory is described and illustrated by examples. The Simplified Theory of Plastic Zones (STPZ) has proven itself for estimating the post-shakedown quantities in the state of elastic and plastic shakedown within the framework of the 1st order theory, i.e. if the equilibrium conditions are satisfied for the undeformed structure.
In this paper, the results of elbows subjected to various loading parameters are compared, considering and neglecting stress stiffening. Thus, the results show the influence 2nd order effects can generate. It is further shown that the STPZ can capture 2nd order effects introduced by the equilibrium of the deformed structure. Some examples are used to demonstrate its applicability and the quality of the results, e.g. for a pipe bend subjected to cyclic in-plane bending.