Subdivision based isogeometric analysis for geometric flows

We present a new isogeometric analysis (IGA) approach based on extended Loop subdivision scheme for solving various geometric flows defined on subdivision surfaces. The studied flows include the second-order, fourth-order, and sixth-order geometric flows, such as averaged mean curvature flow, constant mean curvature flow, and minimal mean-curvature-variation flow, which are generally derived by minimizing the associate energy functionals with $L^2$-gradient flow respectively. The geometric flows are discretized by means of subdivision based IGA, where the finite element space is formulated by the limit form of the extended Loop subdivision for different initial control meshes. The basis functions, consisting of quartic box-splines corresponding to each subdivided control mesh, are utilized to represent the geometry exactly. For the cases of the evolution of open surfaces with any shape boundary, high-order continuous boundary conditions derived from the mixed variational forms of the geometric flows should be implemented to be consistent with the isogeometric concept. For time discretization, we adopt an adaptive semi-implicit Euler scheme. By several numerical experiments, we study the convergence behaviors of the proposed approach for solving the geometric flows with high-order boundary conditions. Moreover, the numerical results also show the accuracy and efficiency of the proposed method.     

Extended isogeometric analysis based on PHT‐splines for crack propagation near inclusions

Adaptive local refinement is one of the main issues for isogeometric analysis (IGA). In this paper, an adaptive extended IGA (XIGA) approach based on polynomial splines over hierarchical T‐meshes (PHT‐splines) for modeling crack propagation is presented. The PHT‐splines overcome certain limitations of nonuniform rational B‐splines–based formulations; in particular, they make local refinements feasible. To drive the adaptive mesh refinement, we present a recovery‐based error estimator for the proposed method. The method is based on the XIGA method, in which discontinuous enrichment functions are added to the IGA approximation and this method does not require remeshing as the cracks grow. In addition, crack propagation is modeled by successive linear extensions that are determined by the stress intensity factors under linear elastic fracture mechanics. The proposed method has been used to analyze numerical examples, and the stress intensity factors results were compared with reference results. The findings demonstrate the accuracy and efficiency of the proposed method.     

Adaptive isogeometric analysis for phase-field modeling of anisotropic brittle fracture

The surface energy a phase-field approach to brittle fracture in anisotropic materials is also anisotropic and gives rise to second-order gradients in the phase field entering the energy functional. This necessitates C¹ continuity of the basis functions which are used to interpolate the phase field. The basis functions which are employed in isogeometric analysis (IGA), such as nonuniform rational B-splines and T-splines naturally possess a higher order continuity and are therefore ideally suited for phase-field models which are equipped with an anisotropic surface energy. Moreover, the high accuracy of spline discretizations, also relative to their computational demand, significantly reduces the fineness of the required discretization. This holds a fortiori if adaptivity is included. Herein, we present two adaptive refinement schemes in IGA, namely, adaptive local refinement and adaptive hierarchical refinement, for phase-field simulations of anisotropic brittle fracture. The refinement is carried out using a subdivision operator and exploits the Bézier extraction operator. Illustrative examples are included, which show that the method can simulate highly complex crack patterns such as zigzag crack propagation. An excellent agreement is obtained between the solutions from global refinement and adaptive refinement, with a reasonable reduction of the computational effort when using adaptivity.     

Weight-adaptive isogeometric analysis for solving elastodynamic problems based on space-time discretization approach

In this paper, a novel adaptive isogeometric analysis (IGA) is introduced and its application in the numerical solution of two-dimensional elastodynamic problems based on the space-time discretization (STD) approach is studied. In the STD approach, the time is considered as an additional dimension and is discretized the same as the spatial domain. The weights of control points play the main role in the proposed method. In the conventional IGA, the same set of weights is used in the modeling of geometric and solution spaces. The idea is to define two groups of weights: geometric and solution weights. Geometric weights are known and can be determined based on the position of control points, but the solution weights are considered to be unknown and can be determined using a proper strategy so that the accuracy of the solution is optimized. This strategy is based on the minimization of an error function. The results obtained from the proposed method are compared with those obtained from the conventional IGA.     

An adaptive isogeometric analysis meshfree collocation method for elasticity and frictional contact problems

A collocation method has been recently developed as a powerful alternative to Galerkin's method in the context of isogeometric analysis, characterized by significantly reduced computational cost, but still guaranteeing higher-order convergence rates. In this work, we propose a novel adaptive isogeometric analysis meshfree collocation (IGAM-C) for the two-dimensional (2D) elasticity and frictional contact problems. The concept of the IGAM-C method is based upon the correspondence between the isogeometric collocation and reproducing kernel meshfree approach, which facilitates the robust mesh adaptivity in isogeometric collocation. The proposed method reconciles collocation at the Greville points via the discretization of the strong form of the equilibrium equations. The adaptive refinement in collocation is guided by the gradient-based error estimator. Moreover, the resolution of the nonlinear equations governing the contact problem is derived from a strong form to avoid the disadvantages of numerical integration. Numerical examples are presented to demonstrate the effectiveness, robustness, and straightforward implementation of the present method for adaptive analysis.     

On the use of adaptive refinement in isogeometric digitalimage correlation

A novel global digital image correlation method was developed using adaptive refinement of isogeometric shape functions. Non‐uniform rational B‐spline shape functions are used because of their flexibility and versatility, which enable them to capture a wide range of kinematics. The goal of this work was to explore the full potential of isogeometric shape functions for digital image correlation (DIC). This is reached by combining a global DIC method with an adaptive refinement algorithm: adaptive isogeometric GDIC. The shape functions are automatically adjusted to be able to describe the kinematics of the sought displacement field with an optimized number of degrees of freedom. This results in an accurate method without the need of making problem‐specific choices regarding the structure of the shape functions, which makes the method less user input dependent than regular global DIC methods, while keeping the number of degrees of freedom limited to realize optimum regularization of the ill‐posed DIC problem. The method's accuracy is demonstrated by a virtual experiment with a predefined, highly localized displacement field. Real experiments with a complex sample geometry demonstrate the effectiveness in practice. Copyright © 2015 John Wiley & Sons, Ltd.     

A large deformation frictional contact formulation using NURBS‐based isogeometric analysis

This paper focuses on the application of NURBS‐based isogeometric analysis to Coulomb frictional contact problems between deformable bodies, in the context of large deformations. A mortar‐based approach is presented to treat the contact constraints, whereby the discretization of the continuum is performed with arbitrary order NURBS, as well as C⁰‐continuous Lagrange polynomial elements for comparison purposes. The numerical examples show that the proposed contact formulation in conjunction with the NURBS discretization delivers accurate and robust predictions. Results of lower quality are obtained from the Lagrange discretization, as well as from a different contact formulation based on the enforcement of the contact constraints at every integration point on the contact surface. Copyright © 2011 John Wiley & Sons, Ltd.     

Development of stochastic isogeometric analysis (SIGA) method for uncertainty in shape

In this paper, a new method is proposed that extend the classical deterministic isogeometric analysis (IGA) into a probabilistic analytical framework in order to evaluate the uncertainty in shape and aim to investigate a possible extension of IGA in the field of computational stochastic mechanics. Stochastic IGA (SIGA) method for uncertainty in shape is developed by employing the geometric characteristics of the non-uniform rational basis spline and the probability characteristics of polynomial chaos expansions (PCE). The proposed method can accurately and freely evaluate problems of uncertainty in shape caused by deformation of the structural model. Additionally, we use the intrusive formulation approach to incorporate PCE into the IGA framework, and the C++ programming language to implement this analysis procedure. To verify the validity and applicability of the proposed method, two numerical examples are presented. The validity and accuracy of the results are assessed by comparing them to the results obtained by Monte Carlo simulation based on the IGA algorithm.     

Extended isogeometric analysis for cohesive fracture

The objective of this study is to present an extended isogeometric formulation for cohesive fracture. The approach exploits the higher order interelement continuity property of nonuniform rational B-splines (NURBS), in particular the higher accuracy that results for the stress prediction, which yields an improved estimate for the direction of crack propagation compared to customary Lagrangian interpolations. Shifting is used to ensure compatibility with the surrounding discretization, where, different from extended finite element methods, the affected elements stretch over several rows perpendicular to the crack path. To avoid fine meshes around the crack tip in case of cohesive fracture, a blending function is used in the extension direction of the crack path. To comply with standard finite element data structures, Bézier extraction is used. The absence of the Kronecker-delta property in the higher order interpolations of isogeometric analysis impedes the enrichment scheme and compatibility enforcement. These issues are studied comprehensively at the hand of several examples, while crack propagation analyses show the viability of the approach.     

Multiresolution topology optimization using isogeometric analysis

The paper introduces a novel multiresolution scheme to topology optimization in the framework of the isogeometric analysis. A new variable parameter space is added to implement multiresolution topology optimization based on the Solid Isotropic Material with Penalization approach. Design density variables defined in the variable space are used to approximate the element analysis density by the bivariate B‐spline basis functions, which are easily obtained using k‐refinement strategy in the isogeometric analysis. While the nonuniform rational B‐spline basis functions are used to exactly describe geometric domains and approximate unknown solutions in finite element analysis. By applying a refined sensitivity filter, optimized designs include highly discrete solutions in terms of solid and void materials without using any black and white projection filters. The Method of Moving Asymptotes is used to solve the optimization problem. Various benchmark test problems including plane stress, compliant mechanism inverter, and 2‐dimensional heat conduction are examined to demonstrate the effectiveness and robustness of the present method.     

An isogeometric analysis approach to gradient damage models

Continuum damage formulations are commonly used for the simulation of diffuse fracture processes. Implicit gradient damage models are employed to avoid the spurious mesh dependencies associated with local continuum damage models. The C⁰‐continuity of traditional finite elements has hindered the study of higher order gradient damage approximations. In this contribution we use isogeometric finite elements, which allow for the construction of higher order continuous basis functions on complex domains. We study the suitability of isogeometric finite elements for the discretization of higher order gradient damage approximations. Copyright © 2011 John Wiley & Sons, Ltd.     

An adaptive IgA-BEM with hierarchical B-splines based on quasi-interpolation quadrature schemes

The isogeometric formulation of the boundary element method (IgA-BEM) is investigated within the adaptivity framework. Suitable weighted quadrature rules to evaluate integrals appearing in the Galerkin BEM formulation of 2D Laplace model problems are introduced. The proposed quadrature schemes are based on a spline quasi-interpolation (QI) operator and properly framed in the hierarchical setting. The local nature of the QI perfectly fits with hierarchical spline constructions and leads to an efficient and accurate numerical scheme. An automatic adaptive refinement strategy is driven by a residual-based error estimator. Numerical examples show that the optimal convergence rate of the Galerkin solution is recovered by the proposed adaptive method.     

X‐FEM in isogeometric analysis for linear fracture mechanics

The extended finite element method (X‐FEM) has proven to be an accurate, robust method for solving problems in fracture mechanics. X‐FEM has typically been used with elements using linear basis functions, although some work has been performed using quadratics. In the current work, the X‐FEM formulation is incorporated into isogeometric analysis to obtain solutions with higher order convergence rates for problems in linear fracture mechanics. In comparison with X‐FEM with conventional finite elements of equal degree, the NURBS‐based isogeometric analysis gives equal asymptotic convergence rates and equal accuracy with fewer degrees of freedom (DOF). Results for linear through quartic NURBS basis functions are presented for a multiplicity of one or a multiplicity equal the degree. Copyright © 2011 John Wiley & Sons, Ltd.     

A framework for implementation of RVE‐based multiscale models in computational homogenization using isogeometric analysis

This study presents an isogeometric framework for incorporating representative volume element–based multiscale models into computational homogenization. First‐order finite deformation homogenization theory is derived within the framework of the method of multiscale virtual power, and Lagrange multipliers are used to illustrate the effects of considering different kinematical constraints. Using a Lagrange multiplier approach in the numerical implementation of the discrete system naturally leads to a consolidated treatment of the commonly employed representative volume element boundary conditions. Implementation of finite deformation computational strain‐driven, stress‐driven, and mixed homogenization is detailed in the context of isogeometric analysis (IGA), and performance is compared to standard finite element analysis. As finite deformations are considered, a numerical multiscale stability analysis procedure is also detailed for use with IGA. Unique implementation aspects that arise when computational homogenization is performed using IGA are discussed, and the developed framework is applied to a complex curved microstructure representing an architectured material.     

On the numerical integration of isogeometric interface elements

Zero‐thickness interface elements are commonly used in computational mechanics to model material interfaces or to introduce discontinuities. The latter class requires the existence of a non‐compliant interface prior to the onset of fracture initiation. This is accomplished by assigning a high dummy stiffness to the interface prior to cracking. This dummy stiffness is known to introduce oscillations in the traction profile when using Gauss quadrature for the interface elements, but these oscillations are removed when resorting to a Newton‐Cotes integration scheme 1. The traction oscillations are aggravated for interface elements that use B‐splines or non‐uniform rational B‐splines as basis functions (isogeometric interface elements), and worse, do not disappear when using Newton‐Cotes quadrature. An analysis is presented of this phenomenon, including eigenvalue analyses, and it appears that the use of lumped integration (at the control points) is the only way to avoid the oscillations in isogeometric interface elements. New findings have also been obtained for standard interface elements, for example that oscillations occur in the relative displacements at the interface irrespective of the value of the dummy stiffness. Copyright © 2015 John Wiley & Sons, Ltd.     

On the treatments of heterogeneities and periodic boundary conditions for isogeometric homogenization analysis

The treatments of heterogeneities and periodic boundary conditions are explored to properly perform isogeometric analysis (IGA) based on NURBS basis functions in solving homogenization problems for heterogeneous media with omni‐directional periodicity and composite plates with in‐plane periodicity. Because the treatment of the combination of different materials in IGA models is not trivial especially for periodicity constraints, the first priority is to clearly specify points at issue in the numerical modeling, or equivalently mesh generation, for IG homogenization analysis (IGHA). The most awkward, but important issue is how to generate patches for NURBS representation of the geometry of a rectangular parallelepiped unit cell to realize appropriate deformations in consideration of the convex‐hull property of IGA. The issue arises from the introduction of overlapped control points located at angular points in the heterogeneous unit cell, which must satisfy multiple point constraint (MPC) conditions associated with periodic boundary conditions (PBCs). Although two measures may be conceivable, we suggest the use of multiple patches along with double MPC that imposes PBCs and the continuity conditions between different patches simultaneously. Several numerical examples of numerical material and plate tests are presented to demonstrate the validity of the proposed strategy of IG modeling for IGHA. Copyright © 2016 John Wiley & Sons, Ltd.     

A generalized finite element formulation for arbitrary basis functions: From isogeometric analysis to XFEM

Many of the formulations of current research interest, including iosogeometric methods and the extended finite element method, use nontraditional basis functions. Some, such as subdivision surfaces, may not have convenient analytical representations. The concept of an element, if appropriate at all, no longer coincides with the traditional definition. Developing a new software for each new class of basis functions is a large research burden, especially, if the problems involve large deformations, non‐linear materials, and contact. The objective of this paper is to present a method that separates as much as possible the generation and evaluation of the basis functions from the analysis, resulting in a formulation that can be implemented within the traditional structure of a finite element program but that permits the use of arbitrary sets of basis functions that are defined only through the input file. Elements ranging from a traditional linear four‐node tetrahedron through a higher‐order element combining XFEM and isogeometric analysis may be specified entirely through an input file without any additional programming. Examples of this framework to applications with Lagrange elements, isogeometric elements, and XFEM basis functions for fracture are presented. Copyright © 2010 John Wiley & Sons, Ltd.