Abstract:This work presents numerical techniques to enforce continuity constraints on multi-patch surfaces for three distinct problem classes. The first involves structural analysis of thin shells that are described by general Kirchhoff-Love kinematics. Their governing equation is a vector-valued, fourth-order, nonlinear, partial differential equation (PDE) that requires at least $C^1$-continuity within a displacement-based finite element formulation. The second class are surface phase separations modeled by a phase field. Their governing equation is the Cahn-Hilliard equation - a scalar, fourth-order, nonlinear PDE - that can be coupled to the thin shell PDE. The third class are brittle fracture processes modeled by a phase field approach. In this work, these are described by a scalar, fourth-order, nonlinear PDE that is similar to the Cahn-Hilliard equation and is also coupled to the thin shell PDE. Using a direct finite element discretization, the two phase field equations also require at least a $C^1$-continuous formulation. Isogeometric surface discretizations - often composed of multiple patches - thus require constraints that enforce the $C^1$-continuity of displacement and phase field. For this, two numerical strategies are presented: For this, two numerical strategies are presented: A Lagrange multiplier formulation and a penalty method. The curvilinear shell model including the geometrical constraints is taken from Duong et al. (2017) and it is extended to model the coupled phase field problems on thin shells of Zimmermann et al. (2019) and Paul et al. (2020) on multi-patches. Their accuracy and convergence are illustrated by several numerical examples considering deforming shells, phase separations on evolving surfaces, and dynamic brittle fracture of thin shells.

What problem does this paper attempt to address?

This paper aims to solve the consistency constraint problem of multi - block surfaces under the fourth - order deformation and phase - field model. Specifically, three different problem categories are focused on in the research: 1. **Thin - shell structure analysis**: Thin - shell structures described by Kirchhoff - Love dynamics. Its governing equation is a fourth - order nonlinear partial differential equation (PDE), which requires at least C1 continuity in the displacement - based finite - element formulation. 2. **Surface phase separation**: Surface phase separation described by the phase - field model. Its governing equation is the Cahn - Hilliard equation - a scalar fourth - order nonlinear PDE, which can be coupled with the thin - shell PDE. 3. **Brittle fracture process**: The brittle fracture process modeled by the phase - field method. These processes are described by a scalar fourth - order nonlinear PDE similar to the Cahn - Hilliard equation, and are also coupled with the thin - shell PDE. In order to maintain C1 continuity in the multi - block surface discretization, the paper proposes two numerical strategies: the Lagrange multiplier method and the penalty method. These methods are used to enforce C1 continuity of displacement and phase - field at the block interfaces to ensure accurate transfer of stress, bending moment, mass flux and damage gradient at the interfaces. ### Main contributions - **Propose new C1 continuity constraints**: Applicable to the phase - field on multi - block NURBS. - **Fully formulated in curvilinear coordinate systems**: Use these constraints to simulate coupling problems on multi - block NURBS. - **Excellent accuracy compared with single - block discretization**: Can describe the branching of cracks on folded surfaces in high - order phase - field models. - **Verify the proposed constraints through multiple nonlinear benchmark problems**. ### Structure overview - **Section 2**: Summarizes the thin - shell theory framework. - **Section 3**: Outlines the calculation formula of the Cahn - Hilliard phase - separation model. - **Section 4**: Outlines the phase - field model of brittle shells. - **Section 5**: Discusses the calculation formula of the required block - interface constraints. - **Section 6**: Presents several numerical examples, highlighting the impact of these constraints on the numerical solution. - **Section 7**: Summarizes the full text and looks forward to future work. ### Key techniques - **Geometric description**: Use the mapping between the parameter domain and the physical surface to describe the curved surface. - **Finite - element discretization**: Perform finite - element approximation of the main fields (deformation and phase - field). - **Weak - form discretization**: Derive the discretization formula of the mechanical weak - form. - **Material constitutive relations**: Derive the constitutive behavior of the material starting from the Helmholtz free - energy function. - **Phase - separation and brittle - fracture models**: Respectively introduce the discretized weak - form and constitutive relations of the Cahn - Hilliard model and the phase - field brittle - fracture model. - **Continuity constraints of block - interfaces**: Propose new C1 continuity constraints to ensure the continuity of multi - block surfaces at the interfaces. Through these methods, the paper provides an effective solution for solving the continuity problem of multi - block surfaces in complex engineering structures.

Isogeometric continuity constraints for multi-patch shells governed by fourth-order deformation and phase field models

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

An Interior Penalty coupling strategy for Isogeometric non-conformal Kirchhoff-Love shell patches

An isogeometric finite element formulation for boundary and shell viscoelasticity based on a multiplicative surface deformation split

An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS

Elasto-plastic large deformation analysis of multi-patch thin shells by isogeometric approach

An isogeometric finite element formulation for phase transitions on deforming surfaces

Large Deformation Frictional Contact Formulations for Isogeometric Kirchhoff–Love Shell

Smooth multi-patch scaled boundary isogeometric analysis for Kirchhoff-Love shells

An embedding-aware continuum thin shell formulation

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

Weak Boundary Condition Enforcement for Linear Kirchhoff-Love Shells: Formulation, Error Analysis, and Verification

Coupling of non-conforming trimmed isogeometric Kirchhoff-Love shells via a projected super-penalty approach

Nonlinear Deformation of Flexible Shallow Shells of Complex Shape Made of Materials with Different Resistance to Tension and Compression

Non-conforming multipatches for NURBS-based finite element analysis of higher-order phase-field models for brittle fracture

Internal-interfacial cracking interaction: Combined phase-field and discontinuous Galerkin/cohesive zone modeling

Variationally consistent isogeometric analysis of trimmed thin shells at finite deformations, based on the STEP exchange format

A general isogeometric finite element formulation for rotation-free shells with in-plane bending of embedded fibers

Phase-Field Modelling of Damage and Fracture—Convergence and Local Mesh Refinement

Coupled peridynamic model for geometrically nonlinear deformation and fracture analysis of thin shell structures

Isogeometric Bezier dual mortaring: The Kirchhoff-Love shell problem