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)
- Interpolation of rotations (2)
- Locking (2)
- Post-shakedown quantities (2)
- Reissner–Mindlin shell (2)
- Reissner–Mindlin shells (2)
- STPZ (2)
- Scaled boundary finite element method (2)
- ratcheting (2)
- Adjusted approximation spaces (1)
- Airy equation (1)
- Approximate dual basis functions for NURBS (1)
- Approximate duals (1)
- Bemessung (1)
- Bézier extraction (1)
- Bézier projection (1)
- Cloth Simulation (1)
- Collocation method (1)
- Computerunterstütztes Verfahren (1)
- Continuity (1)
- Cooling Towers (1)
- Coupling of non-conforming meshes (1)
- Deformation Gradient (1)
- Domain Decomposition (1)
- Dual basis functions for NURBS (1)
- Efficient Integration (1)
- Elbow (1)
- Explicit dynamics (1)
- Ferroelectric material (1)
- Ferroelectrics (1)
- Ferromagnetic material (1)
- Finite element method (1)
- Finland (1)
- Fließfigur (1)
- Geometric effects (1)
- Geometrical nonlinear Reissner–Mindlin shell (1)
- Geometrically nonlinear Reissner–Mindlin shell (1)
- Ginzburg–Landau equation (1)
- Hybrid collocation-Galerkin method (1)
- Hybrid collocation–Galerkin method (1)
- Integration schemes (1)
- Interpolation of Rotations (1)
- Interpolation of the director (1)
- Ke (1)
- Lahti (1)
- Mass lumping (1)
- Mass lumping scheme (1)
- Mixed Formulations (1)
- Mortar method (1)
- Mortar methods (1)
- Multiferroic heterostructure (1)
- N-body Problem (1)
- Numerical integration (1)
- Numerical simulations (1)
- Optimal convergence (1)
- Phase field method (1)
- Phase-field modeling (1)
- Plasticity (1)
- Plastizieren (1)
- Platte (1)
- Proceedings of the 8th World Conference on Timber Engineering (1)
- Progressive buckling (1)
- Reissner-Mindlin Shell Analysis (1)
- Rotation interpolation (1)
- Scaled Boundary FEM (1)
- Second order effects (1)
- Shape Matching (1)
- Shell not requiring drilling rotation stabilization (1)
- Simplified theory of plastic zones (1)
- Solidifcation (1)
- Spectral elements (1)
- Spline Basis Functions (1)
- Strain range (1)
- Structural vibrations (1)
- Surface represented solid (1)
- Three-dimensional displays (1)
- Traglast (1)
- Transverse isotropy (1)
- Transverse shear locking (1)
- Treatment of kinks (1)
- Twice Yield (1)
- Vereinfachte Fließzonentheorie (1)
- Wave propagation (1)
- Weak Substitution Method (1)
- Wrinkle Vector Field (1)
- Zarka's method (1)
- adapted approximation orders (1)
- caled boundary finite element method (1)
- continuous director vector interpolation (1)
- convective basis systems (1)
- direct methods (1)
- direct shakedown analysis (1)
- domain coupling (1)
- domain decomposition (1)
- dual basis functions (1)
- element formulation (1)
- explicit dynamics (1)
- fatigue (1)
- finite element method (1)
- implicit dynamics (1)
- isogeometric analysis (1)
- limit load (1)
- mortar methods (1)
- multi-parameter loading (1)
- multi-patch connection (1)
- multi-scale (1)
- nonlinear isogeometric analysis (1)
- nonlinear material behavior (1)
- phase-field (1)
- progressive Deformation (1)
- simplified theory of plastic zones (1)
- stress stiffening (1)
- theory of porous media (1)
- two-dimensional elasticity (1)
- wcte 2004 (1)
- weak substitution method (1)
- zyklische Belastung (1)
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.
The focus of this study is on providing an efficient and highly accurate mass lumping scheme for explicit dynamic calculations in isogeometric analysis (IGA). Computer-Aided design models make use of higher-order polynomials, usually Non-Uniform Rational Splines (NURBS), to build up the geometry model. In IGA they are chosen as shape functions to take advantage of their high continuity and to keep the description of geometry exact.
Thus, compared to standard Finite Element Method (FEM), a smaller number of elements is required to gain results of the same quality level. In general, explicit time integration methods require a huge number of time steps to obtain numerical stable results of reasonable quality. Hence, diagonal mass matrices are preferred to reduce the computational costs within each time step. Row-sum lumping and diagonal scaling method are two well-known techniques, which have been developed as mass lumping schemes for standard FEM.
Unfortunately, these standard FEM mass lumping schemes deteriorate convergence rates in IGA dynamics significantly. With rising order of the basis functions, accuracy of the results decreases. The reason for this is that shape functions of higher order lead to more dense mass matrices and thus the error caused by lumping the mass matrices increases. Therefore, new mass lumping schemes have to be developed, which are more suitable for the use of higher-order shape functions.
Several types of dual basis functions (duals), which have already been used in the isogeometric mortar method, are considered for a new mass lumping scheme. Keeping the initial NURBS as shape functions, the duals are applied as test functions for the variational form. Depending on the chosen type of duals, diagonal or banded mass matrices are obtained without any additional lumping. Compared to lumping the original mass matrices based on NURBS as test functions, applying additional row-sum lumping to the banded matrices causes significantly smaller errors.
The presented new mass lumping scheme is a very promising approach, in which duals are used as test functions in IGA dynamics with explicit time integration, e.g. the Central Differences Method. The implementation into existing codes can be easily done. The assembled system matrices, if based on NURBS as shape and test functions, only have to be multiplied with a transformation matrix. Numerical examples show that if using this dual lumping scheme within IGA methods, explicit dynamic analysis can be performed efficiently with higher-order shape functions.
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.
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.
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 aim of this study is to provide a new mass lumping scheme for explicit dynamic calculations in isogeometric analysis (IGA).
Shape functions for IGA methods are taken from the Computer-Aided design (CAD) model. The use of higher-order polynomials, usually Non-Uniform Rational B-Splines (NURBS), ensures an exact geometry description [1]. Thus, in comparison to standard finite element method (FEM) the number of elements and therefore the computational costs can be lowered. In addition to that the convergence rate raises equally with the polynomial order. Increasing the degree of shape functions,
less elements are necessary for results of the same quality.
In explicit dynamics, mass lumping schemes are commonly applied to reduce computational costs by using diagonal mass matrices. Well-known techniques like row-sum lumping or diagonal scaling method were developed for dynamic analysis with standard FEM. Unfortunately, they are neither suitable for higher-order shape functions in IGA, nor in the spectral element method [2]. It is not possible to take advantage of increasing the polynomial order to lower the number of elements
in explicit dynamic IGA, because the error caused by lumping the mass matrix also increases. Thus, other mass lumping schemes, suitable for higher-order shape functions, have to be developed. As they are already used for isogeometric mortar method [3], several types of dual basis functions (duals) are studied for a new mass lumping scheme. Using duals as test functions in IGA leads to diagonal consistent mass matrices [4] or banded matrices, depending on the chosen type of duals.
In case of the banded matrices, additional row-sum lumping causes smaller errors than lumping the original mass matrix with standard NURBS as test functions. For IGA dynamics with explicit time integration, e.g. the Central Difference Method, using duals as test functions is a very promising approach. They can be easily implemented in existing methods through multiplying the already assembled system matrices with a transformation matrix based on the underlying NURBS curve. Applying the dual lumping scheme, an explicit dynamic analysis
could also be performed efficiently using IGA methods with higher-order shape functions.
In isogeometric analysis (IGA), which was founded by Hughes et al. [1], the geometry representation is used for the analysis as well. Hence, due to the exact description of the geometry, analysis results can be improved [1, 2]. Therefore, different kinds of splines, like non-uniform rational B-splines (NURBS) [3], are used as shape functions for the discretizations. In linear elasticity problems, for standard formulations, shear locking phenomena can occur due to the different orders of the derivatives in the unbalanced strain-displacement relation. This effect can be reduced using shape functions of higher order, causing an increased computational effort. For low-order formulations, the degrees of shape functions can be adapted accordingly in order to counteract this effect. For an isogeometric displacement-stress mixed Reissner-Mindlin shell formulation such adapted approximation spaces were investigated in [4], for instance. In [5], additionally to a proper choice of shape function spaces, convective coordinates are employed in the derivation of isogeometric shell formulations.
In this contribution, a convective displacement-based isogeometric formulation is introduced, wherein the displacements in the different surface directions are approximated independently using appropriate approximation orders. Therefore, two different meshes are generated from the NURBSdescribed geometry representation, employing order elevation solely for one of the two surface directions. Thus, the order is elevated in opposite directions for these two meshes. The two different possibilities of order elevation are investigated. Furthermore, the use of different convective basis systems is studied. This comprises convective basis systems in each control point, computed according to [6], as well as convective basis systems determined from the local geometry direction in each integration point. The achieved results are compared to whose of a two-dimensional displacementstress mixed formulation presented in [7].
Isogeometric analysis and mixed finite element methods offer promising opportunities to enhance analysis results for complex problems like incompressible elasticity and are able to cope with different locking phenomena. In this contribution, a mixed two-field isogeometric formulation with independent approximations for displacements and stresses is derived, and its ability to counteract different types of locking is investigated using two examples. Furthermore, the influence of the continuity of the stress shape functions on the accurancy of results and convergence behaviour is shown.
An isogeometric element formulation for linear two-dimensional elasticity based on the Airy equation
(2022)
The aim of this work is to derive a formulation for linear two-dimensional elasticity using just one degree of freedom. This degree of freedom is used to directly discretize the Airy bipotential equation, which requires higher order basis functions. Isogeometric structural analysis is based on shape functions of the geometry description in Computer-Aided design software.
These shape functions can easily fulfill the continuity requirement of the bipotential equation.
Thus, an Airy element formulation can be obtained through isogeometric methods. In this contribution Non-Uniform Rational B-splines are used to discretize the domain and to solve the occurring differential equations. Numerical examples demonstrate the accuracy of the evolved formulation for a quadratic plate under different load situations.
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.
We investigate the mortar finite element method for second order elliptic boundary value problems on domains which are decomposed into patchesk with tensor-product NURBS parameterizations. We follow the methodology of IsoGeometric Analysis (IGA) and choose discrete spaces Xh,k on each patch k as tensor-product NURBS spaces of the same or higher degree as given by the parameterization. Our work is an
extension of Brivadis et al. (Comput Methods Appl Mech Eng 284:292–319, 2015) and highlights several aspects which did not receive full attention before. In particular, by choosing appropriate spaces of polynomial splines as Lagrange multipliers, we obtain
a uniform infsup-inequality. Moreover, we provide a new additional condition on the discrete spaces Xh,k which is required for obtaining optimal convergence rates of the mortar method. Our numerical examples demonstrate that the optimal rate is lost if this condition is neglected.
Isogeometric Analysis (IGA), which was founded by Hughes et al. [1], tries to unify Computer-aided design (CAD) and finite element analysis (FEA) by using the same model for geometry representation and analysis. Therefore, in contrast to common finite element analysis, non-uniform rational B-splines (NURBS) and other kinds of splines are used as shape functions instead of the usual polynomials. Due to the exact representation of the geometry, analysis results can be improved [1, 2]. For mixed formulations, stresses and/or strains or pressures are approximated independently and in addition to the usual displacement approximation. Using, this kind of method is more robust and offers more accurate results. Hence, such formulations can be employed to solve incompressible elasticity problems, for instance [3]. Recent investigations, e.g. [4], have already combined isogeometric analysis and mixed formulations in order to benefit from the advantages of both methods.
In this contribution, a mixed isogeometric method is proposed in order to improve the analysis results and counteract different locking phenomena. Therefore, spline basis functions are used and the displacement shape functions of a two-dimensional isogeometric plane stress and plane strain element are supplemented by independent stress shape functions. These additional stress shape functions are chosen to be of one order lower compared to the displacement shape functions, but with adapted continuity. Evaluating the errors of the calculated displacements and stresses for a plane stress and an incompressible elasticity plane strain problem, it is shown that the proposed method can lead to an improved accuracy of results compared to a standard isogeometric formulation. Furthermore, the influence of different continuities on the convergence behavior and 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 the opinion of the authors, this effect deserves more attention. Therefore, the purpose of this paper is to illuminate the effect of stress stiffening by investigating some examples, with particular attention to elastic-plastic strain ranges and the plastic strain range enhancement factor Ke used for fatigue analyses.
The simplified theory of plastic zones (STPZ) was mainly developed to determine strain ranges and accumulated strains in the state of shakedown at cyclic loading between prescribed levels of loading. Kinematic hardening is an indispensable feature of the STPZ. The plastic limit load, however, is defined for monotonic loading and elastic–plastic material behavior without hardening. Simply assigning a zero value or a numerically very low value of the tangent modulus when applying the STPZ is generally not possible due to arising numerical instabilities. It is, therefore, not immediately obvious how the STPZ can be used to determine the maximum load level that can be applied to a structure without developing a kinematic mechanism. This paper describes the theory and the analysis steps required and provides some illustrative examples. Typically, between one and three linear elastic analyses and some local calculations are required to provide either the exact value or at least a reasonable estimate of a range of the plastic limit load, as well as of the associated stress and strain fields and displacements that are not provided by classical limit analysis.
Cyclic and over-elastic loading can lead to an accumulation of plastic strains. If there is a cyclic load, which is driven by a single parameter, the lifecycle design can be very costly in terms of computational effort. If more than one cyclic load parameter is to be taken into account, which is then a multi-parameter loading, this task can become even more complex and costly. To solve this problem efficiently, different techniques are proposed. One of these techniques is based on step-by-step calculations of the strain ranges for a reduced set of loadings. Once these strain ranges are known, the accumulated state for each individual load case can be estimated using the Simplified Theory of Plastic Zones (STPZ), which requires just a few linear elastic analyses. It is shown that cyclic loads, which occur in intervals, can be replaced by interval-free calculations, which reduce the computational effort enormously. All these techniques lead to a procedure, which delivers good estimations in terms of post-shakedown quantities with very low computational effort compared to incremental step-by-step calculations. The results of the STPZ are presented by an example. A thick-walled cylinder is loaded with a constant axial force and subjected to cyclic shear and cyclic internal pressure. In general, for structures exhibiting ratcheting, hundreds or more load cycles must be analysed via step-by-step calculations until the shakedown state is reached. Using the STPZ, post-shakedown quantities, including strain ranges and accumulated strains can be estimated efficiently and the structure can be designed according to the rules of the ASME Codes. The computational effort and the quality of the results of the STPZ are compared with a step-by-step calculation.
In this paper, we propose a new template fitting method that can capture fine details of garments in target 3D scans of dressed human bodies. Matching the high fidelity details of such loose/tight-fit garments is a challenging task as they express intricate folds, creases, wrinkle patterns, and other high fidelity surface details. Our proposed method of non-rigid shape fitting – FoldMatch – uses physics-based particle dynamics to explicitly model the deformation of loose-fit garments and wrinkle vector fields for capturing clothing details. The 3D scan point cloud behaves as a collection of astrophysical particles, which attracts the points in template mesh and defines the template motion model. We use this point-based motion model to derive regularized deformation gradients for the template mesh. We show the parameterization of the wrinkle vector fields helps in the accurate shape fitting. Our method shows better performance than the stateof-the-art methods. We define several deformation and shape matching quality measurement metrics to evaluate FoldMatch on synthetic and real data sets.
We develop a mixed geometrically nonlinear isogeometric Reissner–Mindlin shell element for the analysis of thin-walled structures that leverages Bézier dual basis functions to address both shear and membrane locking and to improve the quality of computed stresses. The accuracy of computed solutions over coarse meshes, that have highly non-interpolatory control meshes, is achieved through the application of a continuous rotational approach. The starting point of the formulation is the modified Hellinger–Reissner variational principle with independent displacement, membrane, and shear strains as the unknown fields. To overcome locking, the strain variables are interpolated with lower-order spline bases while the variations of the strain variables are interpolated with the corresponding Bézier dual bases. Leveraging the orthogonality property of the Bézier dual basis, the strain variables are condensed out of the system with only a slight increase in the bandwidth of the resulting linear system. The condensed approach preserves the accuracy of the non-condensed mixed approach but with fewer degrees of freedom. From a practical point of view, since the Bézier dual basis is completely specified through Bézier extraction, any spline space that admits Bézier extraction can utilize the proposed approach directly.
Patch coupling in isogeometric analysis of solids inboundary representation using a mortar approach
(2020)
This contribution is concerned with a coupling approach for nonconforming NURBS patches in the framework of an isogeometric formulation for solids inboundary representation. The boundary representation modeling technique inCAD is the starting point of this approach. We parameterize the solid according to the scaled boundary finite element method and employ NURBS basis functions for the approximation of the solution. Therefore, solid surfaces consist of several sections, which can be regarded as patches and discretized independently. The main objective of this study is to derive an approach for the connection of independent sections in order to allow for local refinement and thus an accurate and efficient discretization of the computational domain. Non-conforming sections are coupled with a mortar approach within a master-slave framework. The coupling of adjacent sections ensures the equality of mutual deformations along the interface in a weak sense and is enforced by constraining the NURBS basis functions on the interface. We apply this approach to non-linear problems in two dimensions and compare the results with conforming discretizations.
Die isogeometrische Methode definiert sich durch den Einsatz einer einheitlichen Geometriebeschreibung für Entwurf und Berechnung.
Insbesondere bei der Berechnung dünnwandiger Flächentragwerke ist durch die Verwendung der exakten Geometrie ein großer Gewinn an Genauigkeit und Zuverlässigkeit möglich. Um neben hoher Genauigkeit auch hohe Effizienz zu erreichen, wird in diesem Beitrag die Verwendung effizienter numerischer Integrationsmethoden für die Berechnung der Steifigkeitsmatrix untersucht.
Modeling of electric field-induced magnetization switching in multi- ferroic heterostructures
(2019)
Multiferroic heterostructures consist of materials with either pronounced ferroelectric or ferromagnetic effect. The combination of both types of material, be it in layers, columns or inclusions, potentially yields a significant magneto‐electric coupling effect even at room temperature. The magnetization in the ferromagnetic material can be controlled by the application of electric fields to the ferroelectric material. In this contribution a linear elastic continuum formulation is coupled with a phase field formulation for the polarization and magnetization in the ferroelectric and the ferromagnetic layer, respectively. The strain transfer at the interface of the layers yields a magneto‐electric coupling effect within the heterostructures. The finite element method is used to discretize the arising differential equations. A numerical example provides a proof of concept for the simulation of the magneto‐electric coupling effect in multiferroic heterostructures.
Metastable austenitic steels can undergo phase transformation. As an allotrope two crystal configurations are of interest: the softer austenitic parent phase and the martensitic phases. Here, the bain orientation relationship leads to distinct orientations for the martensitic variants with a different transformation strain [7]. A phase field approach is used to model the transformation, where a multi‐valued order parameter urn:x-wiley:16177061:media:PAMM201900465:pamm201900465-math-0001 identifies the austenitic parent phase and the martensitic variants. This allows to define bulk and surface energies as regularized functions in terms of the order parameter and its gradient. The kinetics of the martensitic transformation are temperature dependent. Temperatures below an equilibrium temperature favour the growth of the martensitic phase, whereas temperatures above the equilibrium favour the austenitic phase. Approaching the equilibrium temperature slows down the transformation [5]. In this work we consider a static crack under mode I and mode II loading.
Isogeometric analysis fosters the integration of design and analysis by using the geometry description of the CAD system also for the numerical analysis. Hereby, the use of NURBS surfaces is common but entails the need for a coupling of non-conforming patches. The use of mortar methods allows a coupling which requires neither additional variables nor empirical parameters. In this contribution dual basis functions are used in order to obtain an accurate and efficient mortar method.
Shell elements for slender structures based on a Reissner-Mindlin approach struggle in pure bending problems. The stiffness of such structures is overestimated due to the transversal shear locking effect.
Here, an isogeometric Reissner-Mindlin shell element is presented, which uses adjusted control meshes for the displacements and rotations in order to create a conforming interpolation of the pure bending compatibility requirement. The method is tested for standard numerical examples.
In this contribution, an isogeometric Reissner-Mindlin shell element is presented, which uses adjusted approximation spaces for the displacements and the rotations in order to avoid transversal shear locking effects. These locking effects arise especially in pure bending problems, with decreasing thickness of the shell. Their origin lies in the not conforming interpolations of the displacements and rotations in the formulation of the compatibility requirements.
One possibility to overcome this difficulty would be to increase the polynomial degree of the used NURBS shape functions, as it is proposed in [1]. However, the locking effects are not completely eliminated and the computational time for the formation of the stiffness matrix increases significantly with rising polynomial degrees. For lower polynomial degrees, there exist only a few effective concepts for the prevention of locking. Beir˜ao da Veiga and his group proposed one of these methods for the elimination of transversal shear locking effects for plate problems [2]. They suggested the implementation of different control meshes, with adjusted polynomial degrees for the interpolation of the displacements and the rotations. In this way transversal shear locking effects due to the coupling of shear strains and curvature are avoided.
The isogeometric concept still holds because the reference geometry for the displacements and the rotations is the same and only the refinement is different.
This method is now extended to an isogeometric Reissner-Mindlin shell formulation. The used shell element is derived from continuum theory. Nodal basis systems are computed with a globalL2 fit and thus the interpolated director vectors of the reference configuration coincide best possible with the normal vectors. The interpolation of the current director vector is performed using a full SO(3) update, so the nodal rotations are interpolated. This deviates from standard rotation interpolation procedures, where the current director vector is interpolated. The more accurate interpolation used in this contribution leads to more accurate results, see [3]. The accuracy and efficiency of the shell element are examined for some standard benchmark examples.
Heterostructures of ferroelectric and ferromagnetic layers are commonly used to obtain electromagnetic effects. The elastic coupling between the layers is widely acknowledged as the main mechanism responsible for the electro-magneto interaction. Within this contribution we study the coupling of ferroelectric and ferromagnetic layers with well-defined interfaces. The intention is to simulate the switching of the magnetization with the help of electric fields, which has been studied experimentally in [1]. Each layer is simulated by using mechanically coupled phase field modeling, whereby the approaches presented in [2] and [3] will be used. The strains in each layer depend on the direction of the polarization/magnetization. A mismatch of these strains will be compensated by local deformations at the interface as the coupling results from the coherent deformation at the interface. This leads to the possibility to alter the magnetization direction by changing the electric polarization and vice verse. Numerical simulations will illustrate the evolution of the ferroic microstructures with a focus on the strain coupling and the resulting interactions between layers.
In this contribution we present a coupling approach for non-conforming NURBS patches in the framework of an isogeometric formulation in boundary representation [1]. In order to fit the boundary representation modeling technique in CAD, we follow the idea of the scaled boundary finite element method [2] for the parameterization of the solid. Thus, two-dimensional solid surfaces are partitioned into sections in relation to a central point, the scaling center. Each section is parameterized with a circumferential parameter along the boundary and a radial scaling parameter in the interior of the domain. We employ NURBS basis functions for the approximation of the solution in both parametric directions. The approximation in scaling direction is flexible and allows for local refinement within the computational domain.
For the coupling of non-conforming sections we employ a mortar approach as presented in [3]. We establish a master-slave relation for the interface control points between adjacent sections based on the equality of mutual deformations along the interface. The coupling is realized in a weak manner by constraining the basis functions, which means that the NURBS basis functions of the slave side are related to those of the master side so that the solution of both sections along the interface is equal in a weak sense. This approach is derived numerically and its application to nonlinear problems is straightforward.
We study several numerical examples of linear and nonlinear problems in solid mechanics and compare the results with conforming computations.
The efficiency of the NURBS-based isogeometric Galerkin method highly depends on the chosen numerical integration scheme for the evaluation of the stiffness matrix and the residual vector. In an efficient implementation, the computational costs for the formation of the stiffness matrix scale almost linearly with the total number of integration points. Furthermore, also the accuracy of the solution
depends on the chosen integration scheme, since spurious locking effects increase with the number of integration points. Gauss integration is probably the most common approach for numerical integration in NURBS-based IGA, whereby the higher continuity in IGA can be exploited for a reduction of the number of integration points, see [2]. More sophisticated integration schemes have been proposed in [3,4], where a nonlinear equation is solved in order to compute integration weights and quadrature points. Using the minimum possible number of quadrature points, the nonlinear equation aims at integrating the univariate B-spline spaces exactly. Thus, the methods can be considered as optimal.
In this contribution, an implementation of the optimal integration method [3] for general NURBS surfaces with non-uniform knot vectors is presented. The integration method is applied to the Poisson problem, two-dimensional plane stress problems and isogeometric shell elements as proposed in [5].
The numerical examples of this study focus on the accuracy of the integration for NURBS surfaces with non-uniform weights. A further very important aspect is the comparison of computational costs between the different methods. The alleviation of locking effects, which comes along with a reduction of the number of integration points, is studied for the shell elements proposed in [5]. Furthermore, a comparison to the locking-free formulation of [6] is given.
Numerical methods for the modeling of the magnetization vector in multiferroic heterostructures
(2017)
Multiferroic heterostructures are commonly used to obtain electro‐magnetic coupling effects. Thereby, the ferroelectric layer is used to control the magnetization in the ferromagnetic layer. The coupling between the layers is obtained by the mechanical coupling between the layers, which have well‐defined interfaces. Within this contribution we use phase field models to define the polarization and magnetization in the ferroelectric and ferromagnetic layers, respectively. A coupling between polarization/magnetization and strains in each layer in combination with coherent deformations at the interface yields an electromagnetic coupling within the entire heterostructure. Numerical formulations for the interpolation of the polarization vector are well‐defined in the literature. However, the establishment of a consistent numerical formulation for the ferromagnetic layer, where the length of the magnetization vector has to be constant, remains a difficult task. We propose a new numerical approach for the consistent treatment of the ferromagnetic layer and provide numerical simulations which illustrate the electromagnetic coupling effect.
Magneto-electric coupling along interfaces between ferroelectric and ferromagneticmaterials allows the magnetization in ferromagnetic layers to be controlled by electrical fields.The coupling effect is mainly due to the deformation coupling between the layers. Thus, furthercoupling effects at the interfaces are neglected in this contribution. Phase field formulationsare used to model the polarization and magnetization in the ferroelectric and ferromagneticlayers, respectively. A coupling between the phase field and the strains is introduced in eachlayer in combination with a mechanical coupling at the interface. The numerical formulationfor the ferroelectric layer is taken from literature. A special focus is set on the discretization ofthe length-constrained order parameter in the ferromagnetic layer. We show a new approachto enforce this constraint and provide a numerical simulation which illustrates the magneto-electric coupling effect.
Interpolation of Rotations and Coupling of Patches in Isogeometric Reissner– Mindlin Shell Analysis
(2015)
This work is concerned with the development of an efficient and robust isogeometric Reissner–Mindlin shell formulation. The basic assumption of shell theories is a dimensional reduction of the three-dimensional continuum to a two-dimensional surface embedded in the three-dimensional space. Consequently, the geometry is described by a reference surface in combination with a director vector field, which defines the expansion in the thickness direction. The main objective of isogeometric analysis is to use the same model description for design and analysis. Thin-walled structures are usually defined by a reference surface and an associated thickness in industrial design software. Thus, the usage of isogeometric shell elements can avoid costly conversions to volumetric geometry descriptions.The usage of NURBS surfaces (Non-Uniform Rational B-splines) possibly yields high continuity between elements. This requires a rethinking of all concepts used in conventional shell elements, which base on linear Lagrange basis functions. The shell formulation presented in this work is derived from the continuum theory and uses an orthogonal rotation described by Rodrigues' tensor to compute the current director vector. Large deformations and finite rotations can be described accurately. The discretization requires nodal director vectors which interpolate the normal vector as exact as possible. A new method for the definition of nodal basis systems and nodal director vectors is derived. Basing on this, a criterion for the automatic assignment of the correct number of rotational degrees of freedom for each node is proposed. This allows stable computations of geometries with kinks while requiring neither the usage of drilling rotation stabilization nor manual user interaction. The main part of this work is the derivation of various concepts for the interpolation of the current director vector, which is a function of the rotational state. The respective concepts differ in the quantity which is actually interpolated and in the chosen update formulation for the rotations. The influence of each concept on the global deformation convergence behavior is assessed with the help of numerical examples. The results suggest that proper convergence behavior for all orders of NURBS basis functions can only be attained if interpolated director vectors are rotated. Concepts of this type are more accurate and expensive than concepts which rotate nodal director vectors. But the higher computational effort pays off for geometries with arbitrary curvature and for basis functions of higher order. Geometries with kinks require a multiplicative rotational update formulation for concepts that rotate interpolated director vectors.Three different integration rules are considered within the numerical examples. Besides full and reduced Gauss integration also a new non-uniform Gauss integration concept following Adam et al. (2015) is assessed. A special focus is put on the interaction between the chosen rotational concept and the integration scheme. The reduction of the number of integration points from full to reduced integration slightly alleviates locking effects. The further reduction entailed by non-uniform integration significantly reduces locking effects in some examples. But this only yields higher accuracy if a concept which rotates interpolated director vectors is chosen. The reduction of locking deteriorates the accuracy of deformation results in other cases. The efficiency of the presented shell formulation is compared to standard shell formulations in terms of computational costs to attain a pre-defined error level. The most effective combination of integration scheme and rotational concept is shown to be competitive to standard shell formulations.A further main concern of this work is the derivation of a mortar-type method for the coupling of non-conforming NURBS surface patches. Methods to handle non-conforming patches without mutual refinement are essential for an efficient application of NURBS-based isogeometric analysis. The proposed method bases on a substitution relation, which is derived from the weak fulfillment of the equality of mutual displacements along the interface. A static condensation can be performed with the help of this substitution relation in order to attain a coupled global system of equations. The variational formulation is not altered and the global stiffness matrix remains positive definite. Numerical examples show the applicability of the method. A comparison to reference results and to computations with the Lagrange multiplier method is given. The applicability of the coupling method for the presented Reissner-Mindlin shell formulation is shown with the help of two nonlinear examples.
This work presents a framework for the computation of complex geometries containing intersections of multiple patches with Reissner–Mindlin shell elements. The main objective is to provide an isogeometric finite element implementation which neither requires drilling rotation stabilization, nor user interaction to quantify the number of rotational degrees of freedom for every node. For this purpose, the following set of methods is presented. Control points with corresponding physical location are assigned to one common node for the finite element solution. A nodal basis system in every control point is defined, which ensures an exact interpolation of the director vector throughout the whole domain. A distinction criterion for the automatic quantification of rotational degrees of freedom for every node is presented. An isogeometric Reissner–Mindlin shell formulation is enhanced to handle geometries with kinks and allowing for arbitrary intersections of patches. The parametrization of adjacent patches along the interface has to be conforming. The shell formulation is derived from the continuum theory and uses a rotational update scheme for the current director vector. The nonlinear kinematic allows the computation of large deformations and large rotations. Two concepts for the description of rotations are presented. The first one uses an interpolation which is commonly used in standard Lagrange-based shell element formulations. The second scheme uses a more elaborate concept proposed by the authors in prior work, which increases the accuracy for arbitrary curved geometries. Numerical examples show the high accuracy and robustness of both concepts. The applicability of the proposed framework is demonstrated.
A NURBS based hybrid collocation-Galerkin method for the analysis of boundary represented solids
(2015)
The paper is concerned with a new numerical method, NURBS based hybrid collocation–Galerkin method (NURBS-HCGM), to solve the in-plane motion problem of elastic solids. It combines the merits of the so-called scaled boundary finite-element method (SB-FEM) and the isogeometric collocation method. For the analysis, the boundary scaling technique of SB-FEM is adopted. It leads to a formulation, where only the boundary of a structure is discretized. Here, the NURBS basis functions are employed for the description of the geometry of the boundary as well as for the approximation of the displacements at the boundary. This is in accordance with the boundary representation modeling technique, which is commonly employed in computer aided design software. The inner domain is described by a radial scaling parameter. Applying the weak form only in circumferential direction the governing partial differential equations of elasticity are transformed to an ordinary differential equation (ODE) of Euler type, where the unknown displacements are a function of the radial scaling parameter. In the present work a NURBS based collocation scheme is introduced to solve this equation. NURBS basis functions are suggested for the approximation of the displacements in scaling direction. The higher continuity provided by NURBS allows to use collocation to solve the ODE directly instead of using the weak form in scaling direction. The proposed approach is validated by comparison with the eigenvalue solution of the ODE. It is remarked that the eigenvalue solution is restricted to linear problems, whereas the proposed method could be extended to nonlinear problems. In general, the presented formulation will allow to model patches bounded by an arbitrary number of contour boundaries. The accuracy of the proposed approach is analyzed and estimated with respect to analytical solutions. The computational cost is investigated with the help of numerical examples and is compared to isogeometric Galerkin approach.
In this contribution, a mortar‐type method for the coupling of non‐conforming NURBS (Non‐Uniform Rational B‐spline) surface patches is proposed. The connection of non‐conforming patches with shared degrees of freedom requires mutual refinement, which propagates throughout the whole patch due to the tensor‐product structure of NURBS surfaces. Thus, methods to handle non‐conforming meshes are essential in NURBS‐based isogeometric analysis. The main objective of this work is to provide a simple and efficient way to couple the individual patches of complex geometrical models without altering the variational formulation. The deformations of the interface control points of adjacent patches are interrelated with a master‐slave relation. This relation is established numerically using the weak form of the equality of mutual deformations along the interface. With the help of this relation, the interface degrees of freedom of the slave patch can be condensated out of the system. A natural connection of the patches is attained without additional terms in the weak form. The proposed method is also applicable for nonlinear computations without further measures. Linear and geometrical nonlinear examples show the high accuracy and robustness of the new method. A comparison to reference results and to computations with the Lagrange multiplier method is given.
This paper presents a numerical method to solve the three-dimensional elasticity problem of surface represented solids. A surface oriented formulation is derived, in which the three-dimensional solid is described by its boundary surfaces and a radial scaling center. Scaling the boundary surfaces with respect to the scaling center yields a parameterization of the complete solid. The definition of the boundary surfaces is sufficient for the description of the solid. Thus, the formulation conforms ideally to the boundary representation modeling technique used in CAD. In the present approach, no tensor-product structure of three-dimensional objects is exploited to parameterize the physical domain. The weak form of the equation of motion and the weak form of the Neumann boundary conditions are enforced only at the boundary surfaces. It leads to a transformation of the governing partial differential equations of elasticity to an ordinary differential equation (ODE) of Euler type. Solving the ODE leads to the displacement in the radial scaling direction. In the present approach, the isogeometric Galerkin finite element method is employed to describe the geometry of the boundary surfaces and also to approximate the displacement response of the boundary surfaces. It exploits two-dimensional NURBS objects to parameterize the boundary surfaces. The final Euler type ODE is solved by the NURBS based collocation method. The displacement response in the radial scaling direction is approximated by the NURBS basis functions. Hence, the method can be extended to nonlinear problems. The accuracy of the method is validated against analytical solutions. In general, the proposed formulation is able to model solids bounded by an arbitrary number of surfaces.
An efficient and robust rotational formulation for isogeometric Reissner–Mindlin shell elements
(2016)
This work is concerned with the development of an efficient and robust isogeometric Reissner–Mindlin shell formulation for the mechanical simulation of thin-walled structures. Such structures are usually defined by non-uniform rational B-splines (NURBS) surfaces in industrial design software. The usage of isogeometric shell elements can avoid costly conversions from NURBS surfaces to other surface or volume geometry descriptions. The shell formulation presented in this contribution uses a continuous orthogonal rotation described by Rodrigues’ tensor in every integration point to compute the current director vector. The rotational state is updated in a multiplicative manner. Large deformations and finite rotations can be described accurately. The proposed formulation is robust in terms of stable convergence behavior in the nonlinear equilibrium iteration for large load steps and geometries with large and arbitrary curvature, and in terms of insensitivity to shell intersections with kinks under small angles. Three different integration schemes and their influence on accuracy and computational costs are assessed. The efficiency and robustness of the proposed isogeometric shell formulation is shown with the help of several examples. Accuracy and efficiency is compared to an isogeometric shell formulation with the more common discrete rotational concept and to Lagrange-based finite element shell formulations. The competitiveness of the proposed isogeometric shell formulation in terms of computational costs to attain a pre-defined error level is shown.
In isogeometric analysis, identical basis functions are used for geometrical representation and analysis. In this work, non‐uniform rational basis splines basis functions are applied in an isoparametric approach. An isogeometric Reissner–Mindlin shell formulation for implicit dynamic calculations using the Galerkin method is presented. A consistent as well as a lumped matrix formulation is implemented. The suitability of the developed shell formulation for natural frequency analysis is demonstrated by a numerical example. In a second set of examples, transient problems of plane and curved geometries undergoing large deformations in combination with nonlinear material behavior are investigated. Via a zero‐thickness stress algorithm for arbitrary material models, a J2‐plasticity constitutive law is implemented. In the numerical examples, the effectiveness, robustness, and superior accuracy of a continuous interpolation method of the shell director vector is compared with experimental results and alternative numerical approaches.
We discuss the application of non-uniform rational B-splines (NURBS) in the scaled boundary finite element method (SBFEM) for the solution of wave propagation problems at rather high frequencies. We focus on the propagation of guided waves along prismatic structures of constant cross-section. Comparisons are made between NURBS-based discretizations and high-order spectral elements in terms of the achievable convergence rates. We find that for the same order of shape functions, NURBS can lead to significantly smaller errors compared with Lagrange polynomials. The difference becomes particularly important at very high frequencies, where spectral elements are prone to instabilities. Furthermore, we analyze the behavior of NURBS for the discretization of curved boundaries, where the benefit of exact geometry representation becomes crucial even in the low-frequency range.
This contribution defines and compares different methods for the computation of dual basis functions for B-splines and Non-Uniform Rational B-splines (NURBS). They are intended to be used as test functions for the isogeometric mortar method, but other fields of application are possible, too. Three different concepts are presented and compared. The first concept is the explicit formula for the computation of dual basis functions for NURBS proposed in the work of Carl de Boor. These dual basis functions entail minimal support, i.e., the support of the dual basis functions is equal to the support of the corresponding B-spline basis functions. In the second concept dual basis functions are derived from the inversion of the Gram matrix. These dual basis functions have global support along the interface. The third concept is the use of approximate dual basis functions, which were initially proposed for the use in harmonic analysis. The support of these functions is local but larger than the support of the associated B-spline basis functions. We propose an extension of the approximate dual basis functions for NURBS basis functions. After providing the general formulas, we elaborate explicit expressions for several degrees of spline basis functions. All three approaches are applied in the frame of the mortar method for the coupling of non-conforming NURBS patches. A method which allows complex discretizations with multiple intersecting interfaces is presented. Numerical examples show that the explicitly defined dual basis functions with minimal support severely deteriorate the global stress convergence behavior of the mechanical analysis. This fact is in accordance with mathematical findings in literature, which state that the optimal reproduction degree of arbitrary functions is not possible without extending the support of the dual basis functions. The dual basis functions computed from the inverse of the Gram matrix yield accurate numerical results but the global support yields significantly higher computational costs in comparison to computations of conforming meshes. Only the approximate dual basis functions yield accurate and efficient computations, where neither accuracy nor efficiency is significantly deteriorated in comparison to computations of conforming meshes. All basic cases of T-intersections and star-intersections are studied. Furthermore, an example which combines all basic cases in a complex discretization is given. The applicability of the presented method for the nonlinear case and for shell formulations is shown with the help of one numerical example.