TY - GEN A1 - Kikis, Georgia A1 - Dornisch, Wolfgang A1 - Klinkel, Sven T1 - A method for the elimination of shear locking effects in an isogeometric Reissner-Mindlin shell formulation T2 - Proceedings of the 7th GACM Colloquium on Computational Mechanics for Young Scientists from Academia and Industry, October 11-13, 2017 in Stuttgart, Germany N2 - 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. Y1 - 2017 U6 - https://doi.org/10.18419/opus-9334 SP - 500 EP - 503 PB - Institute for Structural Mechanics University of Stuttgart CY - Stuttgart ER - TY - GEN A1 - Kikis, Georgia A1 - Dornisch, Wolfgang A1 - Klinkel, Sven T1 - Isogeometric Reissner-Mindlin shell analysis - adjusted approximation spaces for the reduction of shear locking effects T2 - VI International Conference on Isogeometric Analysis N2 - 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. Y1 - 2017 UR - http://congress.cimne.com/iga2017/admin/files/fileabstract/a74.pdf PB - CIMNE CY - Barcelona ER - TY - GEN A1 - Chasapi, Margarita A1 - Dornisch, Wolfgang A1 - Klinkel, Sven T1 - Patch coupling in isogeometric analysis of solids inboundary representation using a mortar approach T2 - International Journal for Numerical Methods in Engineering N2 - 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. KW - Isogeometric analysis KW - Scaled Boundary FEM KW - Mortar Method Y1 - 2020 U6 - https://doi.org/10.1002/nme.6354 SN - 1097-0207 VL - 121 IS - 14 SP - 3206 EP - 3226 ER - TY - GEN A1 - Dornisch, Wolfgang A1 - Klinkel, Sven A1 - Simeon, Bernd T1 - Isogeometric Reissner–Mindlin shell analysis with exactly calculated director vectors T2 - Computer Methods in Applied Mechanics and Engineering N2 - An isogeometric Reissner–Mindlin shell derived from the continuum theory is presented. The geometry is described by NURBS surfaces. The kinematic description of the employed shell theory requires the interpolation of the director vector and of a local basis system. Hence, the definition of nodal basis systems at the control points is necessary for the proposed formulation. The control points are in general not located on the shell reference surface and thus, several choices for the nodal values are possible. The proposed new method uses the higher continuity of the geometrical description to calculate nodal basis system and director vectors which lead to geometrical exact interpolated values thereof. Thus, the initial director vector coincides with the normal vector even for the coarsest mesh. In addition to that a more accurate interpolation of the current director and its variation is proposed. Instead of the interpolation of nodal director vectors the new approach interpolates nodal rotations. Account is taken for the discrepancy between interpolated basis systems and the individual nodal basis systems with an additional transformation. The exact evaluation of the initial director vector along with the interpolation of the nodal rotations lead to a shell formulation which yields precise results even for coarse meshes. The convergence behavior is shown to be correct for k-refinement allowing the use of coarse meshes with high orders of NURBS basis functions. This is potentially advantageous for applications with high numerical effort per integration point. The geometrically nonlinear formulation accounts for large rotations. The consistent tangent matrix is derived. Various standard benchmark examples show the superior accuracy of the presented shell formulation. A new benchmark designed to test the convergence behavior for free form surfaces is presented. Despite the higher numerical effort per integration point the improved accuracy yields considerable savings in computation cost for a predefined error bound. KW - Isogeometric analysis KW - Reissner–Mindlin shell KW - NURBS KW - Interpolation of the director Y1 - 2013 U6 - https://doi.org/10.1016/j.cma.2012.09.010 SN - 0045-7825 VL - 253 SP - 491 EP - 504 ER - TY - GEN A1 - Klinkel, Sven A1 - Chen, Lin A1 - Dornisch, Wolfgang T1 - A NURBS based hybrid collocation-Galerkin method for the analysis of boundary represented solids T2 - Computer Methods in Applied Mechanics and Engineering N2 - 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. KW - Collocation method KW - NURBS basis functions KW - Scaled boundary finite element method KW - Hybrid collocation–Galerkin method Y1 - 2015 U6 - https://doi.org/10.1016/j.cma.2014.10.029 SN - 0045-7825 VL - 284 SP - 689 EP - 711 ER - TY - GEN A1 - Kikis, Georgia A1 - Dornisch, Wolfgang A1 - Klinkel, Sven T1 - Adjusted approximation spaces for the treatment of transverse shear locking in isogeometric Reissner–Mindlin shell analysis T2 - Computer Methods in Applied Mechanics and Engineering N2 - Transverse shear locking is an issue that occurs in Reissner–Mindlin plate and shell elements. It leads to an artificial stiffening of the system and to oscillations in the stress resultants for thin structures. The thinner the structure is, the more pronounced are the effects. Since transverse shear locking is caused by a mismatch in the approximation spaces of the displacements and the rotations, a field-consistent approach is proposed for an isogeometric degenerated Reissner–Mindlin shell formulation. The efficiency and accuracy of the method is investigated for benchmark plate and shell problems. A comparison to element formulations with locking alleviation methods from the literature is provided. KW - Isogeometric analysis KW - Reissner–Mindlin shells KW - Transverse shear locking KW - Adjusted approximation spaces Y1 - 2019 U6 - https://doi.org/10.1016/j.cma.2019.05.037 SN - 0045-7825 VL - 354 SP - 850 EP - 870 ER - TY - GEN A1 - Sobota, Paul M. A1 - Dornisch, Wolfgang A1 - Müller, Ralf A1 - Klinkel, Sven T1 - Implicit dynamic analysis using an iso- geometric Reissner–Mindlin shell formulation T2 - International Journal for Numerical Methods in Engineering N2 - 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. KW - isogeometric analysis KW - Reissner–Mindlin shell KW - NURBS KW - nonlinear material behavior KW - Structural vibrations KW - implicit dynamics KW - continuous director vector interpolation Y1 - 2017 U6 - https://doi.org/10.1002/nme.5429 SN - 1097-0207 VL - 110 IS - 9 SP - 803 EP - 825 ER - TY - GEN A1 - Chen, Lin A1 - Klinkel, Sven A1 - Dornisch, Wolfgang T1 - Hybrid collocation-Galerkin approach for the analysis of surface represented 3D-solids employing SB-FEM T2 - Computer Methods in Applied Mechanics and Engineering N2 - 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. KW - Surface represented solid KW - NURBS basis functions KW - Scaled boundary finite element method KW - Hybrid collocation-Galerkin method Y1 - 2015 U6 - https://doi.org/10.1016/j.cma.2015.07.004 SN - 0045-7825 VL - 295 SP - 268 EP - 289 ER - TY - GEN A1 - Dornisch, Wolfgang A1 - Müller, Ralf A1 - Klinkel, Sven T1 - An efficient and robust rotational formulation for isogeometric Reissner–Mindlin shell elements T2 - Computer Methods in Applied Mechanics and Engineering N2 - 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. KW - Isogeometric analysis KW - Geometrically nonlinear Reissner–Mindlin shell KW - NURBS KW - Interpolation of rotations KW - Numerical integration Y1 - 2016 U6 - https://doi.org/10.1016/j.cma.2016.01.018 SN - 0045-7825 VL - 303 SP - 1 EP - 34 ER - TY - GEN A1 - Dornisch, Wolfgang A1 - Vitucci, Gennaro A1 - Klinkel, Sven T1 - The weak substitution method – An application of the mortar method for patch coupling in NURBS-based isogeometric analysis T2 - International Journal for Numerical Methods in Engineering N2 - 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. KW - nonlinear isogeometric analysis KW - NURBS KW - domain decomposition KW - weak substitution method KW - multi-patch connection KW - mortar methods Y1 - 2015 U6 - https://doi.org/10.1002/nme.4918 SN - 1097-0207 VL - 103 IS - 3 SP - 205 EP - 234 ER -