Kirchhoff-Love shell representation and analysis using triangle configuration B-splines

This paper presents the application of triangle configuration B-splines (TCB-splines) for representing and analyzing the Kirchhoff-Love shell in the context of isogeometric analysis (IGA). The Kirchhoff-Love shell formulation requires global $C^1$-continuous basis functions. The nonuniform rational B-spline (NURBS)-based IGA has been extensively used for developing Kirchhoff-Love shell elements. However, shells with complex geometries inevitably need multiple patches and trimming techniques, where stitching patches with high continuity is a challenge. On the other hand, due to their unstructured nature, TCB-splines can accommodate general polygonal domains, have local refinement, and are flexible to model complex geometries with $C^1$ continuity, which naturally fit into the Kirchhoff-Love shell formulation with complex geometries. Therefore, we propose to use TCB-splines as basis functions for geometric representation and solution approximation. We apply our method to both linear and nonlinear benchmark shell problems, where the accuracy and robustness are validated. The applicability of the proposed approach to shell analysis is further exemplified by performing geometrically nonlinear Kirchhoff-Love shell simulations of a pipe junction and a front bumper represented by a single patch of TCB-splines

Isogeometric Kirchhoff–Love shell formulations for general hyperelastic materials

We present formulations for compressible and incompressible hyperelastic thin shells which can use general 3D constitutive models. The necessary plane stress condition is enforced analytically for incompressible materials and iteratively for compressible materials. The thickness stretch is statically condensed and the shell kinematics are completely described by the first and second fundamental forms of the midsurface. We use C1-continuous isogeometric discretizations to build the numerical models. Numerical tests, including structural dynamics simulations of a bioprosthetic heart valve, show the good performance and applicability of the presented methods

A comparison of smooth basis constructions for isogeometric analysis

In order to perform isogeometric analysis with increased smoothness on complex domains, trimming, variational coupling or unstructured spline methods can be used. The latter two classes of methods require a multi-patch segmentation of the domain, and provide continuous bases along patch interfaces. In the context of shell modeling, variational methods are widely used, whereas the application of unstructured spline methods on shell problems is rather scarce. In this paper, we therefore provide a qualitative and a quantitative comparison of a selection of unstructured spline constructions, in particular the D-Patch, Almost-$C^1$, Analysis-Suitable $G^1$ and the Approximate $C^1$ constructions. Using this comparison, we aim to provide insight into the selection of methods for practical problems, as well as directions for future research. In the qualitative comparison, the properties of each method are evaluated and compared. In the quantitative comparison, a selection of numerical examples is used to highlight different advantages and disadvantages of each method. In the latter, comparison with weak coupling methods such as Nitsche's method or penalty methods is made as well. In brief, it is concluded that the Approximate $C^1$ and Analysis-Suitable $G^1$ converge optimally in the analysis of a bi-harmonic problem, without the need of special refinement procedures. Furthermore, these methods provide accurate stress fields. On the other hand, the Almost-$C^1$ and D-Patch provide relatively easy construction on complex geometries. The Almost-$C^1$ method does not have limitations on the valence of boundary vertices, unlike the D-Patch, but is only applicable to biquadratic local bases. Following from these conclusions, future research directions are proposed, for example towards making the Approximate $C^1$ and Analysis-Suitable $G^1$ applicable to more complex geometries

An isogeometric finite element formulation for phase transitions on deforming surfaces

This paper presents a general theory and isogeometric finite element implementation for studying mass conserving phase transitions on deforming surfaces. The mathematical problem is governed by two coupled fourth-order nonlinear partial differential equations (PDEs) that live on an evolving two-dimensional manifold. For the phase transitions, the PDE is the Cahn-Hilliard equation for curved surfaces, which can be derived from surface mass balance in the framework of irreversible thermodynamics. For the surface deformation, the PDE is the (vector-valued) Kirchhoff-Love thin shell equation. Both PDEs can be efficiently discretized using $C^1$-continuous interpolations without derivative degrees-of-freedom (dofs). Structured NURBS and unstructured spline spaces with pointwise $C^1$-continuity are utilized for these interpolations. The resulting finite element formulation is discretized in time by the generalized-$\alpha$ scheme with adaptive time-stepping, and it is fully linearized within a monolithic Newton-Raphson approach. A curvilinear surface parameterization is used throughout the formulation to admit general surface shapes and deformations. The behavior of the coupled system is illustrated by several numerical examples exhibiting phase transitions on deforming spheres, tori and double-tori.Comment: fixed typos, extended literature review, added clarifying notes to the text, added supplementary movie file

Isogeometric iFEM analysis of thin shell structures

Shape sensing is one of most crucial components of typical structural health monitoring systems and has become a promising technology for future large-scale engineering structures to achieve significant improvement in their safety, reliability, and affordability. The inverse finite element method (iFEM) is an innovative shape-sensing technique that was introduced to perform three-dimensional displacement reconstruction of structures using in situ surface strain measurements. Moreover, isogeometric analysis (IGA) presents smooth function spaces such as non-uniform rational basis splines (NURBS), to numerically solve a number of engineering problems, and recently received a great deal of attention from both academy and industry. In this study, we propose a novel “isogeometric iFEM approach” for the shape sensing of thin and curved shell structures, through coupling the NURBS-based IGA together with the iFEM methodology. The main aim is to represent exact computational geometry, simplify mesh refinement, use smooth basis/shape functions, and allocate a lower number of strain sensors for shape sensing. For numerical implementation, a rotation-free isogeometric inverse-shell element (isogeometric Kirchhoff–Love inverse-shell element (iKLS)) is developed by utilizing the kinematics of the Kirchhoff–Love shell theory in convected curvilinear coordinates. Therefore, the isogeometric iFEM methodology presented herein minimizes a weighted-least-squares functional that uses membrane and bending section strains, consistent with the classical shell theory. Various validation and demonstration cases are presented, including Scordelis–Lo roof, pinched hemisphere, and partly clamped hyperbolic paraboloid. Finally, the effect of sensor locations, number of sensors, and the discretization of the geometry on solution accuracy is examined and the high accuracy and practical aspects of isogeometric iFEM analysis for linear/nonlinear shape sensing of curved shells are clearly demonstrated

CAD-integrierte Isogeometrische Analyse und Entwurf leichter Tragwerke

Isogeometric methods are extended for the parametric design process of complex lightweight structures. Three novel methods for the coupling of different structural elements are proposed: rotational coupling, implicit geometry description, and frictionless sliding contact. Moreover, the necessary steps for the integration of the numerical analysis, including pre- and post-processing, in CAD are investigated. It is possible to base several different analyses on each other in order to parametrically represent a construction process with multiple steps.Die isogeometrischen Methoden werden zur Anwendung im parametrischen Entwurfsprozess von komplexen Leichtbaustrukturen erweitert. Hierzu werden drei neue Methoden zur Kopplung unterschiedlicher Strukturelemente vorgeschlagen: Rotationskopplung, implizite Geometriebeschreibung und reibungsfreier Gleitkontakt. Ferner werden die nötigen Schritte zur Einbindung von Pre- und Postprocessing für numerische Simulationen in CAD untersucht. Mehrere unterschiedliche Analysen können auf einander folgen und werden verlinkt, um den Aufbauprozess in mehreren Schritten vollparametrisch abzubilden

Isogeometric analysis of ice accretion on wind turbine blades

For wind turbines operating in cold weather conditions, ice accretion is an established issue that remains an obstacle in effective turbine operation. While the aerodynamic performance of wind turbine blades with ice accretion has received considerable research attention, few studies have investigated the structural impact of blade ice accretion. This work proposes an adaptable projection-based method to superimpose complex ice configurations onto a baseline structure. The proposed approach provides an efficient methodology to include ice accretion in the high fidelity isogeometric shell analysis of a realistic wind turbine blade. Linear vibration and nonlinear deflection analyses of the blade are performed for various ice configurations to demonstrate the impact of different ice accretion distributions on structural performance. These analyses indicate decreases in the blade natural frequencies and deflection under icing conditions. Such ice-induced changes clearly reveal the need for structural design consideration for turbines operating under icing conditions

Automated shape and thickness optimization for non-matching isogeometric shells using free-form deformation

Isogeometric analysis (IGA) has emerged as a promising approach in the field of structural optimization, benefiting from the seamless integration between the computer-aided design (CAD) geometry and the analysis model by employing non-uniform rational B-splines (NURBS) as basis functions. However, structural optimization for real-world CAD geometries consisting of multiple non-matching NURBS patches remains a challenging task. In this work, we propose a unified formulation for shape and thickness optimization of separately-parametrized shell structures by adopting the free-form deformation (FFD) technique, so that continuity with respect to design variables is preserved at patch intersections during optimization. Shell patches are modeled with isogeometric Kirchhoff--Love theory and coupled using a penalty-based method in the analysis. We use Lagrange extraction to link the control points associated with the B-spline FFD block and shell patches, and we perform IGA using the same extraction matrices by taking advantage of existing finite element assembly procedures in the FEniCS partial differential equation (PDE) solution library. Moreover, we enable automated analytical derivative computation by leveraging advanced code generation in FEniCS, thereby facilitating efficient gradient-based optimization algorithms. The framework is validated using a collection of benchmark problems, demonstrating its applications to shape and thickness optimization of aircraft wings with complex shell layouts

Manifold-based isogeometric analysis basis functions with prescribed sharp features

We introduce manifold-based basis functions for isogeometric analysis of surfaces with arbitrary smoothness, prescribed $C^0$ continuous creases and boundaries. The utility of the manifold-based surface construction techniques in isogeometric analysis was demonstrated in Majeed and Cirak (CMAME, 2017). The respective basis functions are derived by combining differential-geometric manifold techniques with conformal parametrisations and the partition of unity method. The connectivity of a given unstructured quadrilateral control mesh in $\mathbb R^3$ is used to define a set of overlapping charts. Each vertex with its attached elements is assigned a corresponding conformally parametrised planar chart domain in $\mathbb R^2$, so that a quadrilateral element is present on four different charts. On the collection of unconnected chart domains, the partition of unity method is used for approximation. The transition functions required for navigating between the chart domains are composed out of conformal maps. The necessary smooth partition of unity, or blending, functions for the charts are assembled from tensor-product B-spline pieces and require in contrast to earlier constructions no normalisation. Creases are introduced across user tagged edges of the control mesh. Planar chart domains that include creased edges or are adjacent to the domain boundary require special local polynomial approximants. Three different types of chart domain geometries are necessary to consider boundaries and arbitrary number and arrangement of creases. The new chart domain geometries are chosen so that it becomes trivial to establish local polynomial approximants that are always $C^0$ continuous across tagged edges. The derived non-rational manifold-based basis functions are particularly well suited for isogeometric analysis of Kirchhoff-Love thin shells with kinks