Kirchhoff-Love Shell Research Articles

Magneto-Thermoelastic Principal Parameter Resonance of a Functionally Graded Cylindrical Shell with Axial Tension

The principal parameter resonance of a ferromagnetic functionally graded (FG) cylindrical shell under the action of axial time-varying tension in magnetic and temperature fields is investigated. The temperature dependence of physical parameters for functionally graded materials (FGMs) is considered. Meanwhile, the tension bending coupling effect is eliminated by introducing the physical neutral surface. The kinetic and strain energies are gained with the Kirchhoff–Love shell theory. Based on the nonlinear magnetization characteristics of ferromagnetic materials, the electromagnetic force acting on the shell is calculated. The nonlinear vibration equations are obtained through Hamilton’s principle. Afterward, the vibration equations are discretized and solved by Galerkin and multiscale methods, respectively. The stability criterion for steady-state motion is established utilizing Lyapunov stability theory. After example analysis, the effects of magnetic field intensity, temperature and power law index on the static deflection are elucidated. Subsequently, the impacts of these parameters, as well as axial tension, on the amplitude-frequency characteristics, resonance amplitude, and multiple solution regions are discussed explicitly. Results indicate that the stiffness can be enhanced due to the generation of static deflection. The amplitude decreases with increasing magnetic field intensity, temperature, and power law index. When the magnetic field intensity surpasses a threshold, the resonance phenomenon disappears.

Stochastic Dynamic Buckling Analysis of Cylindrical Shell Structures Based on Isogeometric Analysis

In this paper, we extend our previous work on the dynamic buckling analysis of isogeometric shell structures to the stochastic situation where an isogeometric deterministic dynamic buckling analysis method is combined with spectral-based stochastic modeling of geometric imperfections. To be specific, a modified Generalized-α time integration scheme combined with a nonlinear isogeometric Kirchhoff–Love shell element is used to simulate the buckling and post-buckling problems of cylindrical shell structures. Additionally, geometric imperfections are constructed based on NURBS surface fitting, which can be naturally incorporated into the isogeometric analysis framework due to its seamless CAD/CAE integration feature. For stochastic analysis, the method of separation is adopted to model the stochastic geometric imperfections of cylindrical shells based on a set of measurements. We tested the accuracy and convergence properties of the proposed method with a cylindrical shell example, where measured geometric imperfections were incorporated. The ABAQUS reference solutions are also presented to demonstrate the superiority of the inherited smooth and high-order continuous properties of the isogeometric approach. For stochastic dynamic buckling analysis, we evaluated the buckling load variability and reliability functions of the cylindrical shell with 500 samples generated based on seven nominally identical shells reported in the geometric imperfection data bank. It is noted that the buckling load variability in the cylindrical shell obtained with static nonlinear analysis is also presented to show the differences between dynamic and static buckling analysis.

Shape optimization of non-matching isogeometric shells with moving intersections

While shape optimization using isogeometric shells exhibits appealing features by integrating design geometries and analysis models, challenges arise when addressing computer-aided design (CAD) geometries comprised of multiple non-uniform rational B-splines (NURBS) patches, which are common in practice. The intractability stems from surface intersections within these CAD models. In this paper, we develop an approach for shape optimization of non-matching isogeometric shells incorporating intersection movement. Separately parametrized NURBS surfaces are modeled using Kirchhoff–Love shell theory and coupled using a penalty-based formulation. The optimization scheme allows shell patches to move without preserving relative location with other members during the shape optimization. This flexibility is achieved through an implicit state function, and analytical sensitivities are derived for the relative movement of shell patches. The introduction of differentiable intersections expands the design space and overcomes challenges associated with large mesh distortion, particularly when optimal shapes involve significant movement of patch intersections in physical space. Throughout optimization iterations, all members within the shell structures maintain the NURBS geometry representation, enabling efficient integration of analysis and design models. The optimization approach leverages the multilevel design concept by selecting a refined model for accurate analysis from a coarse design model while maintaining the same geometry. We adopt several example problems to verify the effectiveness of the proposed scheme and demonstrate its applicability to the optimization of the internal stiffeners of an aircraft wing.

Adaptive methods with [formula omitted] splines for multi-patch surfaces and shells

We introduce an adaptive isogeometric method for multi-patch surfaces and Kirchhoff–Love shell structures with hierarchical splines characterized by C1 continuity across patches. We extend the construction of smooth hierarchical splines from the multi-patch planar setting to analysis suitable G1 surfaces. The adaptive scheme to solve fourth order partial differential equations is presented in a general framework before showing its application for the numerical solution of the bilaplacian and the Kirchhoff–Love model problems. A selection of numerical examples illustrates the performance of hierarchical adaptivity on different multi-patch surface configurations.

Computationally-efficient locking-free isogeometric discretizations of geometrically nonlinear Kirchhoff–Love shells

Discretizations based on the Bubnov-Galerkin method and the isoparametric concept suffer from membrane locking when applied to Kirchhoff–Love shell formulations. Membrane locking causes not only smaller displacements than expected, but also large-amplitude spurious oscillations of the membrane forces. Continuous-assumed-strain (CAS) elements were originally introduced to remove membrane locking in quadratic NURBS-based discretizations of linear plane curved Kirchhoff rods (Casquero et al., CMAME, 2022). In this work, we propose CASs and CASns elements to overcome membrane locking in quadratic NURBS-based discretizations of geometrically nonlinear Kirchhoff–Love shells. CASs and CASns elements are interpolation-based assumed-strain locking treatments. The assumed strains have C0 continuity across element boundaries and different components of the membrane strains are interpolated at different interpolation points. CASs elements use the assumed strains to obtain both the physical strains and the virtual strains, which results in a global tangent matrix which is a symmetric matrix. CASns elements use the assumed strains to obtain only the physical strains, which results in a global tangent matrix which is a non-symmetric matrix. To the best of the authors’ knowledge, CASs and CASns elements are the first assumed-strain treatments to effectively overcome membrane locking in quadratic NURBS-based discretizations of geometrically nonlinear Kirchhoff–Love shells while satisfying the following important characteristics for computational efficiency: (1) No additional unknowns are added, (2) No additional systems of algebraic equations need to be solved, (3) The same elements are used to approximate the displacements and the assumed strains, (4) No additional matrix operations such as matrix inversions or matrix multiplications are needed to obtain the stiffness matrix, and (5) The nonzero pattern of the stiffness matrix is preserved. Analogously to the interpolation-based assumed-strain locking treatments for Lagrange polynomials that are widely used in commercial FEA software, the implementation of CASs and CASns elements only requires to modify the subroutine that computes the element residual vector and the element tangent matrix. The benchmark problems show that CASs and CASns elements, using either 2 × 2 or 3 × 3 Gauss–Legendre quadrature points per element, are effective locking treatments since these element types result in more accurate displacements for coarse meshes and excise the spurious oscillations of the membrane forces.

Wideband Vibro-Acoustic Coupling Investigation in Three Dimensions Using Order-Reduced Isogeometric Finite Element/Boundary Element Method

This study introduces an innovative model-order reduction (MOR) technique that integrates boundary element and finite element methodologies, streamlining the analysis of wideband vibro-acoustic interactions within aquatic and aerial environments. The external acoustic phenomena are efficiently simulated via the boundary element method (BEM), while the finite element method (FEM) adeptly captures the dynamics of vibrating thin-walled structures. Furthermore, the integration of isogeometric analysis within the finite element/boundary element framework ensures geometric integrity and maintains high-order continuity for Kirchhoff–Love shell models, all without the intermediary step of meshing. Foundational to our reduced-order model is the application of the second-order Arnoldi method coupled with Taylor expansions, effectively eliminating the frequency dependence of system matrices. The proposed technique significantly enhances the computational efficiency of wideband vibro-acoustic coupling analyses, as demonstrated through numerical simulations.

Study on fracture of hyperelastic Kirchhoff-Love plates and shells by phase field method

Thin walled structures such as plates and shells are widely used in many engineering fields. To Predict its fracture behavior is of great significance for integrity design and strength evaluation of engineering structures. Numerical simulation of the fracture behavior of hyperelastic plates and shells is a challenge due to complex kinematic description, hyperelastic constitutive relationship, geometric nonlinearity and the degradation on elastic parameter caused by fracture damage. Combining Kirchhoff Love (K-L) shell theory with the fracture phase field method, and numerically discretizing the first and second order partial derivatives of displacement field and phase field by using T-splines and meeting the requirements of K-L plate and shell theory for the C1 continuity of the shape function, a model for the isogeometric analysis numerical formulation of the phase field fracture in hyperelastic K-L plates and shells is established. The fracture failure behavior of hyperelastic K-L plates and shells under the uniform load and displacement load is simulated, and the effect of the Gaussian curvature on the fracture behavior of hyperelastic K-L shells is studied. The simulation results show that the present numerical scheme can effectively capture the complex crack propagation path of plates and shells under the uniform load, and the displacement field can effectively reflect the crack distribution of materials. The thin shell with negative Gaussian curvature shows the excellent fracture performance under the internal pressure, and can withstand the greater internal pressure.

FEM-BEM analysis of acoustic interaction with submerged thin-shell structures under seabed reflection conditions

The paper introduces a novel method for simulating acoustic–shell interaction subject to seabed reflection. The vibration of thin-shell structure is analyzed using the finite element method (FEM) with Kirchhoff–Love shell elements. Exterior acoustic fields associated with seabed reflection is simulated by the boundary element method (BEM) with half-space fundamental solutions. Then the FEM and BEM are coupled to simulate interaction between acoustics wave and shell vibration in the context of isogeometric analysis. Finally, fast multipole method is formulated for the coupled system to accelerate computation and reduce memory storage. Numerical examples are provided to demonstrate the accuracy of the algorithm.

Efficient large-deformation shell elements for the simulation of electromechanical coupling effects in incompressible dielectric elastomers

So-called electro-active polymers are not only capable of undergoing very large deformations, but they also exhibit different electromechanical coupling effects due to their dielectric or electrostrictive characteristics. In the present paper the sub-class of dielectric elastomers is studied; in particular, we focus on thin dielectric elastomer shells, which undergo large deformations under an applied electric field. We develop a framework for the modeling and subsequent efficient simulation of thin structures made from incompressible dielectric elastomers. A thermodynamically consistent phenomenological continuum model based on the free energy density is adapted imposing the kinematic assumptions of Kirchhoff-Love shells in combination with including the thickness deformation and the electric field as independent unknown fields. To make the formulation accessible to finite element simulations, existing non-linear elastic shell mixed finite elements are extended to the present hyperelastic electromechanically coupled formulation for dielectric elastomer shells. Computational results of the shell theory are compared to three-dimensional results validating the proposed modeling and numerical simulation methodology, and showing high accuracy already for coarse finite element discretizations of the shell.

Sensitivity to intensity and distribution of the temperature field in the nonlinear thermo-mechanical analysis of laminated glass plates

Glass laminates consist of stiff glass plies permanently shear-coupled by polymeric interposed layers. When an external temperature rise occurs, the interlayers undergo a dramatic stiffness decay. As a consequence, not only the sectional warping typical of alternating stiff/soft composites is intensified, but also the overall behavior may evolve counter-intuitively. When slender elements prone to geometric nonlinearities are involved, even small thermal variations in intensity or distribution may act as uncertainty factors, strongly affecting the outcome. This paper proposes an efficient, robust, and accurate numerical framework to perform the sensitivity analysis to thermo-mechanical actions in glass plates. A large deformation isogeometric Kirchhoff-Love shell model enriched with through-the-thickness warping is employed, together with a generalized arc-length method involving a suitable temperature parameter as an additional unknown, namely the thermal amplifier or a spatial distribution coefficient. Numerical experiments are presented to highlight the effects that even small temperature variations produce on the equilibrium paths and the influence of the stiffness loss in the interlayer on the structural behavior and the accuracy of the models.

An efficient active-stress electromechanical isogeometric shell model for muscular thin film simulations

We propose an isogeometric approach to model the deformation of active thin films using layered, nonlinear, Kirchhoff–Love shells. Isogeometric Collocation and Galerkin formulations are employed to discretize the electrophysiological and mechanical sub-problems, respectively, with the possibility to adopt different element and time-step sizes.Numerical tests illustrate the capabilities of the active-stress-based approach to effectively simulate the contraction of thin films in both quasi-static and dynamic conditions.

Folding of a growing hyperelastic sheet in a viscous fluid

Understanding the growth-induced deformation of soft materials in viscous environments is important for a variety of problems in nature and engineering. Here, we focus on the fluid–structure interaction of a hyperelastic sheet growing in an incompressible Newtonian fluid in the Stokes flow regime. We develop a computational framework for simulating this problem, where the isogeometric boundary integral method is used with the Kirchhoff–Love shell formulation and elastic–plastic decomposition of the deformation gradient tensor. We quantify the relative effects of the growth rate, the sheet bending rigidity, and the fluid viscosity on the fold formation and development of the growing sheet. Our results suggest that the viscous resistance to in-plane deformation promotes fold formation, whereas the viscous resistance to out-of-plane deformation suppresses fold development. We also investigate the effects of the thickness and aspect ratio of the rectangular sheet. Finally, we compare the growth- and prestrain-induced deformations to find a common behavior of sheets under viscous environments.

Fast parametric analysis of trimmed multi-patch isogeometric Kirchhoff-Love shells using a local reduced basis method

AbstractThis contribution presents a model order reduction framework for real-time efficient solution of trimmed, multi-patch isogeometric Kirchhoff-Love shells. In several scenarios, such as design and shape optimization, multiple simulations need to be performed for a given set of physical or geometrical parameters. This step can be computationally expensive in particular for real world, practical applications. We are interested in geometrical parameters and take advantage of the flexibility of splines in representing complex geometries. In this case, the operators are geometry-dependent and generally depend on the parameters in a non-affine way. Moreover, the solutions obtained from trimmed domains may vary highly with respect to different values of the parameters. Therefore, we employ a local reduced basis method based on clustering techniques and the Discrete Empirical Interpolation Method to construct affine approximations and efficient reduced order models. In addition, we discuss the application of the reduction strategy to parametric shape optimization. Finally, we demonstrate the performance of the proposed framework to parameterized Kirchhoff-Love shells through benchmark tests on trimmed, multi-patch meshes including a complex geometry. The proposed approach is accurate and achieves a significant reduction of the online computational cost in comparison to the standard reduced basis method.

Broadband topology optimization of three-dimensional structural-acoustic interaction with reduced order isogeometric FEM/BEM

This paper presents a model order reduction method to accelerate broadband topology optimization of structural-acoustic interaction systems by coupling Finite Element Methods and Boundary Element Methods. The finite element method is used for simulating thin-shell vibration and the boundary element method for exterior acoustic fields. Moreover, the finite element and boundary element methods are implemented in the context of isogeometric analysis, whereby the geometric accuracy and high order continuity of Kirchhoff-Love shells can be guaranteed and meantime no meshing is necessary. The topology optimization method takes continuous material interpolation functions in the density and bulk modulus, and adopts adjoint variable methods for sensitivity analysis. The reduced order model is constructed based on second-order Arnoldi algorithm combined with Taylor's expansions which eliminate the frequency dependence of the system matrices. Numerical results show that the proposed algorithm can significantly improve the efficiency of broadband topology optimization analysis.

A Coarsened-Shell-Based Cosserat Model for the Simulation of Hybrid Cables

The simulation of elastic slender objects like cables is essential for industrial applications in predicting elastic behaviors and life cycles. The Cosserat model and its variants are the dominant approaches due to their high efficiency and accuracy. However, these assume cables with homogeneous interiors and thus cannot simulate hybrid cables containing different materials. We address this by developing a novel coarsened-shell-based Cosserat (CSC) model. The CSC model constructs a material-aware elastic energy function along the cable’s cross-section to describe the global elastic behavior. The CSC model is specifically developed by carefully leveraging the strengths of three approaches: the Cosserat theory to model slender cables, the Kirchhoff–Love shell theory to model the cable’s cross-sectional energy, and numerical coarsening to reduce the degrees of freedom in the shell simulation via constructing a set of new types of material-aware shape/base functions. This allows the more accurate computation of the local and global deformations of hybrid cables, surpassing the classical Cosserat models in accuracy.

A hyperelastic extended Kirchhoff–Love shell model with out-of-plane normal stress: II. An isogeometric discretization method for incompressible materials

This is Part II of a multipart article on a hyperelastic extended Kirchhoff–Love shell model with out-of-plane normal stress. We introduce an isogeometric discretization method for incompressible materials and present test computations. Accounting for the out-of-plane normal stress distribution in the out-of-plane direction affects the accuracy in calculating the deformed-configuration out-of-plane position, and consequently the nonlinear response of the shell. The return is more than what we get from accounting for the out-of-plane deformation mapping. The traction acting on the shell can be specified on the upper and lower surfaces separately. With that, the model is now free from the “midsurface’ location in terms of specifying the traction. In dealing with incompressible materials, we start with an augmented formulation that includes the pressure as a Lagrange multiplier and then eliminate it by using the geometrical representation of the incompressibility constraint. The resulting model is an extended one, in the Kirchhoff–Love category in the degree-of-freedom count, and encompassing all other extensions in the isogeometric subcategory. We include ordered details as a recipe for making the implementation practical. The implementation has two components that will not be obvious but might be critical in boundary integration. The first one is related to the edge-surface moment created by the Kirchhoff–Love assumption. The second one is related to the pressure/traction integrations over all the surfaces of the finite-thickness geometry. The test computations are for dome-shaped inflation of a flat circular shell, rolling of a rectangular plate, pinching of a cylindrical shell, and uniform hydrostatic pressurization of the pinched cylindrical shell. We compute with neo-Hookean and Mooney–Rivlin material models. To understand the effect of the terms added in the extended model, we compare with models that exclude some of those terms.

Goal-adaptive Meshing of Isogeometric Kirchhoff–Love Shells

AbstractMesh adaptivity is a technique to provide detail in numerical solutions without the need to refine the mesh over the whole domain. Mesh adaptivity in isogeometric analysis can be driven by Truncated Hierarchical B-splines (THB-splines) which add degrees of freedom locally based on finer B-spline bases. Labeling of elements for refinement is typically done using residual-based error estimators. In this paper, an adaptive meshing workflow for isogeometric Kirchhoff–Love shell analysis is developed. This framework includes THB-splines, mesh admissibility for combined refinement and coarsening and the Dual-Weighted Residual (DWR) method for computing element-wise error contributions. The DWR can be used in several structural analysis problems, allowing the user to specify a goal quantity of interest which is used to mark elements and refine the mesh. This goal functional can involve, for example, displacements, stresses, eigenfrequencies etc. The proposed framework is evaluated through a set of different benchmark problems, including modal analysis, buckling analysis and non-linear snap-through and bifurcation problems, showing high accuracy of the DWR estimator and efficient allocation of degrees of freedom for advanced shell computations.

An interior penalty coupling strategy for isogeometric non-conformal Kirchhoff–Love shell patches

This work focuses on the coupling of trimmed shell patches using Isogeometric Analysis, based on higher continuity splines that seamlessly meet the C1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C^1$$\\end{document} requirement of Kirchhoff–Love-based discretizations. Weak enforcement of coupling conditions is achieved through the symmetric interior penalty method, where the fluxes are computed using their correct variationally consistent expression that was only recently proposed and is unprecedentedly adopted herein in the context of coupling conditions. The constitutive relationship accounts for generically laminated materials, although the proposed tests are conducted under the assumption of uniform thickness and lamination sequence. Numerical experiments assess the method for an isotropic and a laminated plate, as well as an isotropic hyperbolic paraboloid shell from the new shell obstacle course. The boundary conditions and domain force are chosen to reproduce manufactured analytical solutions, which are taken as reference to compute rigorous convergence curves in the L2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L^2$$\\end{document}, H1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H^1$$\\end{document}, and H2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H^2$$\\end{document} norms, that closely approach optimal ones predicted by theory. Additionally, we conduct a final test on a complex structure comprising five intersecting laminated cylindrical shells, whose geometry is directly imported from a STEP file. The results exhibit excellent agreement with those obtained through commercial software, showcasing the method’s potential for real-world industrial applications.

Large rotation isogeometric shell model for alternating stiff/soft curved laminates including warping and interlayer thickness change

The mechanics of laminates made up of elastic alternating stiff/soft layers is dominated by bending and membrane actions in the stiff layers, while transverse shear deformations concentrate in the soft interlayers producing significant zigzag warping effects. Additionally, curved geometry and large deformations can induce an interlayer thickness strain affecting the overall response of the laminate. This work presents a rotation-free Total Lagrangian large deformation curved shell model to accurately capture such a structural response. Each stiff layer is modeled as a Kirchhoff–Love shell. Without additional DOFs, each soft interlayer is modeled as a solid-shell by expressing, in finite kinematics, the displacement of top and bottom interfaces as a function of the mid-surface displacement of the upper and lower stiff layers it couples. The model is characterized by a sparse stiffness matrix compared to high-order theories and requires a small number of DOFs (3 per stiff layer) to reproduce the 3D solution. A NURBS discretization meets the high continuity of the weak form with a patch-wise reduced quadrature to avoid locking and enhance efficiency. All thickness integrals are pre-computed just once.

Full-Scale Isogeometric Topology Optimization of Cellular Structures Based on Kirchhoff–Love Shells

Full-Scale Isogeometric Topology Optimization of Cellular Structures Based on Kirchhoff–Love Shells