isogeometric analysis using t splines we explore t splines a generalization of nurbs enabling local refinement as a basis for isogeometric analysis we review t splines as a surface design methodology and then develop it for engineering analysis applications we test t splines on some elementary two dimensional and three dimensional fluid and structural analysis problems and attain good results in all cases we summarize the current status of t splines their limitations and future possibilities

variational and projection methods for the volume constraint in finite deformation elasto plasticity this paper focuses on the treatment of volume constraints which in the context of elasto plasticity typically arise as a result of assuming volume preserving plastic flow projection methods based on the modification of the discrete gradient operator b often proposed on an ad hoc basis are systematically obtained in the variational context furnished by a three field hu washizu principle the fully nonlinear formulation proposed here is based on a local multiplicative split of the deformation gradient into volume preserving and dilatational parts without relying on rate forms of the weak form of momentum balance this approach fits naturally in a formulation of plasticity based on the multiplicative decomposition of the deformation gradient and enables one to exactly enforce the condition of volume preserving plastic flow within the framework proposed in this paper rate forms and incrementally objective algorithms are entirely bypassed

immersed b spline i spline finite element method for geometrically complex domains a novel b spline based immersed finite element method is introduced for the computation of geometrically and topologically complex problems the geometry description and the finite element analysis rely on a block structured logically cartesian mesh which encloses the domain of interest a signed distance function is used for representing the domain on the cartesian mesh whereby the domain boundary is the zeroth level set of the signed distance function away from the domain boundaries the standard b spline basis functions are used for the finite element interpolation close to domain boundaries a new approach has been developed for modifying the b spline basis functions so that they locally interpolate the dirichlet boundary conditions the efficiency and robustness of the proposed approach is demonstrated with a number of one two and three dimensional linear boundary value problems

isogeometric collocation methods with generalized b splines we introduce isogeometric collocation methods based on generalized b splines and we analyze their performance through numerical examples for univariate and multivariate scalar and vector valued problems in particular advection diffusion and linear elasticity model problems are addressed the resulting method combines the favorable properties of isogeometric collocation and the geometrical and analytical problem oriented advantages of generalized b splines

subdivision surfaces a new paradigm for thin shell finite element analysis we develop a new paradigm for thin shell finite element analysis based on the use of subdivision surfaces for i describing the geometry of the shell in its undeformed configuration and ii generating smooth interpolated displacement fields possessing bounded energy within the strict framework of the kirchhoff love theory of thin shells the particular subdivision strategy adopted here is loop s scheme with extensions such as required to account for creases and displacement boundary conditions the displacement fields obtained by subdivision are h2 and consequently have a finite kirchhoff love energy the resulting finite elements contain three nodes and element integrals are computed by a one point quadrature the displacement field of the shell is interpolated from nodal displacements only in particular no nodal rotations are used in the interpolation the interpolation scheme induced by subdivision is non local i e the displacement field over one element depend on the nodal displacements of the element nodes and all nodes of immediately neighbouring elements however the use of subdivision surfaces ensures that all the local displacement fields thus constructed combine conformingly to define one single limit surface numerical tests including the belytschko et al 10 obstacle course of benchmark problems demonstrate the high accuracy and optimal convergence of the method

a hierarchical approach to adaptive local refinement in isogeometric analysis abstract adaptive local refinement is one of the key issues in isogeometric analysis in this article we present an adaptive local refinement technique for isogeometric analysis based on extensions of hierarchical b splines we investigate the theoretical properties of the spline space to ensure fundamental properties like linear independence and partition of unity furthermore we use concepts well established in finite element analysis to fully integrate hierarchical spline spaces into the isogeometric setting this also allows us to access a posteriori error estimation techniques numerical results for several different examples are given and they turn out to be very promising

isogeometric analysis and hierarchical refinement for higher order phase field models abstract while the interest in higher order models in physics and mechanics grows their numerical simulation still poses a challenge especially for arbitrary shaped three dimensional domains this contribution presents the mathematical framework as well as the application to different problems in the field of material science fracture mechanics and diffusion problems all models under consideration require at least c 1 continuity which prevents the application of standard finite element analysis and local mesh refinements introducing isogeometric analysis iga for the discretization in a finite element framework enables us to deal with these requirements moreover a general hierarchical refinement scheme based on a subdivision projection is presented here for one two and three dimensional simulations this technique allows to enhance the approximation space using finer splines on each level but preserves the partition of unity as well as the continuity properties of the original discretization using this mathematical framework the improved convergence of a kuramoto sivashinsky model a mesh adapted thermal diffusion simulation and computations of a priori unknown crack propagation in different fracture modes underline the versatility of the presented hierarchical refinement scheme

multiscale phenomena green s functions the dirichlet to neumann formulation subgrid scale models bubbles and the origins of stabilized methods abstract an approach is developed for deriving variational methods capable of representing multiscale phenomena the ideas are first illustrated on the exterior problem for the helmholtz equation this leads to the well known dirichlet to neumann formulation next a class of subgrid scale models is developed and the relationships to bubble function methods and stabilized methods are established it is shown that both the latter methods are approximate subgrid scale models the identification for stabilized methods leads to an analytical formula for the intrinsic time scale whose origins have been a mystery heretofore

multigrid methods for isogeometric analysis with thb splines we propose a geometric multigrid algorithm for the solution of the linear systems which arise when using hierarchical spline spaces for isogeometric discretizations of elliptic partial differential equations in particular we focus on the use of hierarchical b spline hb and truncated hierarchical b spline thb bases to implement adaptively refined isogeometric discretizations we describe an approach for constructing the nested sequences of hierarchical spline spaces which are required for our geometric multigrid solver furthermore we give an algorithm for the computation of the so called prolongation matrices which map coarser space functions to their finer space representations for both hb and thb spline bases in several numerical experiments we study two dimensional poisson problems with singularities and employ a posteriori error estimators for the adaptive refinement of the hierarchical spline spaces using standard multigrid smoothers we observe that using thb spline bases consistently confers a dramatic advantage in terms of iteration numbers over the use of hb spline bases in particular for higher spline degrees when used as a preconditioner for conjugate gradient iteration the proposed multigrid method exhibits practical iteration numbers for two dimensional problems up to the tested case of fourth degree splines

a stabilised immersed boundary method on hierarchical b spline grids for fluid rigid body interaction with solid solid contact abstract an accurate efficient and robust numerical scheme is presented for the simulation of the interaction between flexibly supported rigid bodies and incompressible fluid flow with topology changes and solid solid contact the solution of the incompressible navier stokes equations is approximated by employing a stabilised formulation on cartesian grids discretised with hierarchical b splines the solid is modelled as a rigid body and represented by linear segments along its boundary kinematic conditions along the fluid rigid body interface are enforced weakly using nitsche s method while ghost penalty operators are employed to avoid excessive ill conditioning of the system matrix arising from small cut cells a staggered scheme is used for resolving the coupled fluid rigid body interaction the contact between moving or moving and fixed solid bodies is modelled with lagrange multipliers the excellent performance and wide range of applicability of the proposed scheme are demonstrated in a number of benchmark tests as well as industrially relevant model problems the examples cover the galloping phenomena particulate flow hydraulic check valves and a model turbine

isogeometric analysis for compound b spline surfaces abstract in this work an isogeometric analysis method for compound b spline surfaces is proposed the main advantage of the proposed method is that it directly utilizes the compound cad surface model generated with surface boolean operations in cad systems to achieve the analysis without additional efforts such as representing the cad model with multiple tensor product patches for this an integration scheme for trimmed elements is presented first which can well handle the geometric complexities of cells within a trimmed element then based on b spline representations of the curves on surfaces the displacement continuity requirement along an intersection curve is realized by the degree of freedom dof coupling of the intersected surfaces in which the dof coupling constraint is enforced by minimizing the augmented potential energy with lagrange multiplier the stationarity of the augmented potential energy yields the equilibrium equations for the displacements of the compound surfaces the validity of the proposed integration scheme is demonstrated with the convergence study on the problem of a hole within an infinite plate and the effectiveness of the proposed isogeometric analysis method is illustrated with an elaborated 3d compound surface shell structure constructed by boolean combination of several surfaces

optimal and reduced quadrature rules for tensor product and hierarchically refined splines in isogeometric analysis abstract we continue the study initiated in hughes et al 2010 in search of optimal quadrature rules for tensor product and hierarchically refined splines in isogeometric analysis these rules are optimal in the sense that there exists no other quadrature rule that can exactly integrate the elements of the given spline space with fewer quadrature points we extend the algorithm presented in hughes et al 2010 with an improved starting guess which combined with arbitrary precision arithmetic results in the practical computation of quadrature rules for univariate non uniform splines up to any precision explicit constructions are provided in sixteen digits of accuracy for some of the most commonly used uniform spline spaces defined by open knot vectors we study the efficacy of the proposed rules in the context of full and reduced quadrature applied to two and three dimensional diffusion reaction problems using tensor product and hierarchically refined splines and prove a theorem rigorously establishing the stability and accuracy of the reduced rules

a subdivision based implementation of the hierarchical b spline finite element method abstract a novel technique is presented to facilitate the implementation of hierarchical b splines and their interfacing with conventional finite element implementations the discrete interpretation of the two scale relation as common in subdivision schemes is used to establish algebraic relations between the basis functions and their coefficients on different levels of the hierarchical b spline basis the subdivision projection technique introduced allows us first to compute all element matrices and vectors using a fixed number of same level basis functions their subsequent multiplication with subdivision matrices projects them during the assembly stage to the correct levels of the hierarchical b spline basis the proposed technique is applied to convergence studies of linear and geometrically nonlinear problems in one two and three space dimensions

an isogeometric design through analysis methodology based on adaptive hierarchical refinement of nurbs immersed boundary methods and t spline cad surfaces we explore hierarchical refinement of nurbs as a basis for adaptive isogeometric and immersed boundary analysis we use the principle of b spline subdivision to derive a local refinement procedure which combines full analysis suitability of the basis with straightforward implementation in tree data structures and simple generalization to higher dimensions we test hierarchical refinement of nurbs for some elementary fluid and structural analysis problems in two and three dimensions and attain good results in all cases using the b spline version of the finite cell method we illustrate the potential of immersed boundary methods as a seamless isogeometric design through analysis procedure for complex engineering parts defined by t spline cad surfaces specifically a ship propeller and an automobile wheel we show that hierarchical refinement considerably increases the flexibility of this approach by adaptively resolving local features

aspects of the formulation and finite element implementation of large strain isotropic elasticity the paper presents a formulation of isotropic large strain elasticity and addresses some computational aspects of its finite element implementation on the theoretical side an eulerian setting of isotropic elasticity is discussed exclusively in terms of the finger tensor as a strain measure noval aspects are a direct representation of the eulerian elastic moduli in terms of the finger tensor and their rigorous decomposition into decoupled volumetric and isochoric contributions based on a multiplicative split of the finger tensor into spherical and unimodular parts the isochoric stress response is formulated in terms of the eigenvalues of the unimodular part of the finger tensor a constitutive algorithm for the computation of the stresses and tangent moduli for plane problems is developed and applied to a model problem of rubber elasticity on the computational side the implementation of the constitutive model in three possible finite element formulations is discussed after pointing out algorithmic techniques for the treatment of incompressible elasticity several numerical simulations are presented which show the performance of the proposed constitutive algorithm and the convergence behaviour of the different finite element fomulations for compressible and incompressible elasticity

a stabilized mixed finite element method for finite elasticity formulation for linear displacement and pressure interpolation abstract a stabilized mixed finite element method for finite elasticity is presented the method circumvents the fulfillment of the ladyzenskaya babuska brezzi condition by adding mesh dependent terms which are functions of the residuals of the euler lagrange equations to the usual galerkin method the weak form and the linearized weak form are presented in terms of the reference and current configuration numerical experiments using a tetrahedral element with linear shape functions for the displacements and for the pressure show that the method successfully yields a stabilized element

reduced bezier element quadrature rules for quadratic and cubic splines in isogeometric analysis abstract we explore the use of various element based reduced quadrature strategies for bivariate and trivariate quadratic and cubic spline elements used in isogeometric analysis the rules studied encompass tensor product gauss and gauss lobatto rules and certain so called monomial rules that do no possess a tensor product structure the objective of the study is to determine quadrature strategies which enjoy the same accuracy and stability behavior as full gauss quadrature but with significantly fewer quadrature points several cases emerge that satisfy this objective and also demonstrate superior efficiency compared with standard c 0 continuous finite elements of the same order

initial data sets of minimal energy the total mass of a gravitational system the adm mass is considered as a functional of the initial data sets for sources and gravitation from the condition that such a functional has a minimum under a restricted class of variations tensorial equations for the gravitational initial data are obtained the solutions to these equations whose existence for a large class of source fields is asserted represent initial data sets for gravitation satisfying the constraint equations and may be considered as having no gravitational radiation

isogeometric analysis using lr b splines abstract the recently proposed locally refined b splines denoted lr b splines by dokken et al 2013 6 may have the potential to be a framework for isogeometric analysis to enable future interoperable computer aided design and finite element analysis in this paper we propose local refinement strategies for adaptive isogeometric analysis using lr b splines and investigate its performance by doing numerical tests on well known benchmark cases the theory behind lr b spline is not presented in full details but the main conceptual ingredients are explained and illustrated by a number of examples

the role of bezier extraction in adaptive isogeometric analysis local refinement and hierarchical refinement summary we present two adaptive refinement techniques namely adaptive local refinement and adaptive hierarchical refinement for isogeometric analysis an element wise point of view is adopted exploiting bezier extraction which facilitates the implementation of adaptive refinement in isogeometric analysis locally refined t splines lr t splines and hierarchical t splines are employed for the description of the geometry as well as for the approximation of the solution space in the analysis the refinement is carried out with the aid of a subdivision operator which is computed by again exploiting the bezier extraction operator the concept and algorithm of an element based adaptive isogeometric analysis are illustrated numerical examples are given to examine the accuracy the convergence and the condition number this article is protected by copyright all rights reserved

stabilized mixed finite element formulations for materially nonlinear partially saturated two phase media standard and enhanced bbar or eas low order finite elements applied to the problem of consolidation of two phase nonlinear continua do not satisfy the lbb condition when the same interpolation functions are used for both displacement and pore pressure fields strong spatial pressure oscillations are the main consequence of the lbb condition violation the class of direct stabilized methods widely used in the field of fluid mechanics is a powerful tool to circumvent violation of the aforementioned condition three stabilized formulations designed for the problem of consolidation of fully or partially saturated media are presented galerkin least squares gls with the least squares term construction based on the residuum of the local fluid mass conservation equation the pressure stabilized formulation fpl in which the rate of the pore pressure laplacian residual free in the incompressibility limit is added to the fluid mass conservation equation and finally a stabilized formulation in which the stabilization term is constructed based on the rate of the residuum of the local momentum equation an h convergence study and test problems for the three stabilized schemes are discussed

a posteriori error estimation and adaptivity for elastoplasticity using the reciprocal theorem we present a posteriori error estimators and adaptive methods for the finite element approximation of nonlinear problems and especially elastoplasticity the main characteristic of the proposed method is the introduction of duality techniques or in other notions the reciprocal theorem for error estimation at an equilibrium point the nonlinear boundary value problem and an additional linearized dual problem are considered the loading of the dual problem is specifically designed for capturing the influence of the errors of the entire domain to the considered variable our approach leads to easy computable refinement indicators for locally or integrally defined variables for instationary problems as elastoplasticity in a first step we neglect the errors due to time discretization and evaluate the error indicators within each time step for a stationary problem the versatility of the presented framework is demonstrated with numerical examples

on the completeness of hierarchical tensor product b splines abstract given a grid in r d consisting of d bi infinite sequences of hyperplanes possibly with multiplicities orthogonal to the d axes of the coordinate system we consider the spaces of tensor product spline functions of a given degree on a multi cell domain such a domain consists of finite set of cells which are defined by the grid a piecewise polynomial function belongs to the spline space if its polynomial pieces on adjacent cells have a contact according to the multiplicity of the hyperplanes in the grid we prove that the connected components of the associated set of tensor product b splines whose support intersects the multi cell domain form a basis of this spline space more precisely if the intersection of the support of a tensor product b spline with the multi cell domain consists of several connected components then each of these components contributes one basis function in order to establish the connection to earlier results we also present further details relating to the three dimensional case with single knots only a hierarchical b spline basis is defined by specifying nested hierarchies of spline spaces and multi cell domains we adapt the techniques from giannelli and juttler 2013 to the more general setting and prove the completeness of this basis in the sense that its span contains all piecewise polynomial functions on the hierarchical grid with the smoothness specified by the grid and the degrees under certain assumptions on the domain hierarchy finally we introduce a decoupled version of the hierarchical spline basis that allows to relax the assumptions on the domain hierarchy in certain situations such as quadratic tensor product splines the decoupled basis provides the completeness property for any choice of the domain hierarchy

goal adaptive isogeometric analysis with hierarchical splines abstract in this work a method of goal adaptive isogeometric analysis is proposed we combine goal oriented error estimation and adaptivity with hierarchical b splines for local h refinement the goal oriented error estimator is computed with a p refined discrete dual space which is adaptively refined alongside the primal space this discrete dual space is proven to be a strict superset of the primal space hierarchical refinements are introduced in marked regions that are formed as the union of chosen coarse level spline supports from the primal basis we present two ways of extracting localized refinement indicators suitable for the hierarchical refinement procedure one based on a partitioning of the dual weighted residual into contributions of basis function supports and one based on the combination of element indicators within a basis function support the proposed goal oriented adaptive strategy is exemplified for the poisson problem and a free surface flow problem numerical experiments on these problems show convergence of the adaptive method with optimal rates furthermore the corresponding goal oriented error estimators are shown to be accurate with effectivity indices in the range of 0 7 1 1

an unfitted hp adaptive finite element method based on hierarchical b splines for interface problems of complex geometry generating finite element discretizations with direct interface parameterizations constitutes a considerable computational expense in case of complex interface geometries the paper at hand introduces a b spline finite element method which circumvents parameterization of interfaces and offers fast and easy meshing irrespective of the geometric complexity involved its core idea is the adaptive approximation of discontinuities by hierarchical grid refinement which adds several levels of local basis functions in the close vicinity of interfaces but unfitted to their exact location so that a simple regular grid of knot span elements can be maintained numerical experiments show that an hp refinement strategy which simultaneously increases the polynomial degree of the b spline basis and the levels of refinement around interfaces achieves exponential rates of convergence despite the presence of discontinuities it is also demonstrated that the hierarchical b spline fem can be used to transfer the recently introduced finite cell concept to geometrically nonlinear problems its computational performance imposition of unfitted boundary conditions and fast hierarchical grid generation are illustrated for a set of benchmark problems in one two and three dimensions and the advantages of the regular grid approach for complex geometries are demonstrated by the geometrically nonlinear simulation of a voxel based foam composite

imposing essential boundary conditions in mesh free methods imposing essential boundary conditions is a key issue in mesh free methods the mesh free interpolation does not verify the kronecker delta property and therefore the imposition of prescribed values is not as straightforward as for the finite element method the aim of this paper is to present a general overview on the existing techniques to enforce essential boundary conditions in galerkin based mesh free methods special attention is paid to the mesh free coupling with finite elements for the imposition of prescribed values and to methods based on a modification of the galerkin weak form particular examples are used to analyze and compare their performance in different situations

isogeometric finite element data structures based on bezier extraction of t splines we develop finite element data structures for t splines based on bezier extraction generalizing our previous work for nurbs as in traditional finite element analysis the extracted bezier elements are defined in terms of a fixed set of polynomial basis functions the so called bernstein basis the bezier elements may be processed in the same way as in a standard finite element computer program utilizing exactly the same data processing arrays in fact only the shape function subroutine needs to be modified while all other aspects of a finite element program remain the same a byproduct of the extraction process is the element extraction operator this operator localizes the topological and global smoothness information to the element level and represents a canonical treatment of t junctions referred to as hanging nodes in finite element analysis and a fundamental feature of t splines a detailed example is presented to illustrate the ideas copyright 2011 john wiley sons ltd

blended b spline construction on unstructured quadrilateral and hexahedral meshes with optimal convergence rates in isogeometric analysis abstract we present a novel blended b spline method to construct bicubic tricubic splines over unstructured quadrilateral and hexahedral meshes for isogeometric analysis c 1 and truncated c 2 b spline functions are used in regular elements whereas c 0 and truncated c 1 b spline functions are adopted in boundary elements and interior irregular elements around extraordinary edges vertices the truncation mechanism is employed for a seamless transition from irregular to regular elements the resulting smoothness of the blended construction is c 2 continuous everywhere except c 0 continuous around extraordinary edges and c 1 continuous across the interface between irregular and regular elements the blended b spline construction yields consistent parameterization during refinement and exhibits optimal convergence rates spline functions in the blended construction form a non negative partition of unity are linearly independent and support bezier extraction such that the construction can be used in existing finite element frameworks several examples provide numerical evidence of optimal convergence rates

a collocated c0 finite element method reduced quadrature perspective cost comparison with standard finite elements and explicit structural dynamics summary we demonstrate the potential of collocation methods for efficient higher order analysis on standard nodal finite element meshes we focus on a collocation method that is variationally consistent and geometrically flexible converges optimally embraces concepts of reduced quadrature and leads to symmetric stiffness and diagonal consistent mass matrices at the same time it minimizes the evaluation cost per quadrature point thus reducing formation and assembly effort significantly with respect to standard galerkin finite element methods we provide a detailed review of all components of the technology in the context of elastodynamics that is weighted residual formulation nodal basis functions on gauss lobatto quadrature points and symmetrization by averaging with the ultra weak formulation we quantify potential gains by comparing the computational efficiency of collocated and standard finite elements in terms of basic operation counts and timings our results show that collocation is significantly less expensive for problems dominated by the formation and assembly effort such as higher order elastostatic analysis furthermore we illustrate the potential of collocation for efficient higher order explicit dynamics throughout this work we advocate a straightforward implementation based on simple modifications of standard finite element codes we also point out the close connection to spectral element methods where many of the key ideas are already established copyright 2014 john wiley sons ltd

shape optimisation with multiresolution subdivision surfaces and immersed finite elements abstract we develop a new optimisation technique that combines multiresolution subdivision surfaces for boundary description with immersed finite elements for the discretisation of the primal and adjoint problems of optimisation similar to wavelets multiresolution surfaces represent the domain boundary using a coarse control mesh and a sequence of detail vectors based on the multiresolution decomposition efficient and fast algorithms are available for reconstructing control meshes of varying fineness during shape optimisation the vertex coordinates of control meshes are updated using the computed shape gradient information by virtue of the multiresolution editing semantics updating the coarse control mesh vertex coordinates leads to large scale geometry changes and conversely updating the fine control mesh coordinates leads to small scale geometry changes in our computations we start by optimising the coarsest control mesh and refine it each time the cost function reaches a minimum this approach effectively prevents the appearance of non physical boundary geometry oscillations and control mesh pathologies like inverted elements independent of the fineness of the control mesh used for optimisation on the immersed finite element grid the domain boundary is always represented with a relatively fine control mesh of fixed resolution with the immersed finite element method there is no need to maintain an analysis suitable domain mesh in some of the presented two and three dimensional elasticity examples the topology derivative is used for introducing new holes inside the domain the merging or removing of holes is not considered

multi level bezier extraction for hierarchical local refinement of isogeometric analysis abstract one of the main topics of research on isogeometric analysis is local refinement among the various techniques currently studied and developed one of the most appealing referred to as hierarchical b splines consists of defining a suitable set of basis functions on different hierarchical levels this strategy can also be improved for example to recover partition of unity resorting to a truncation operation giving rise to the so called truncated hierarchical b splines despite its conceptual simplicity implementing the hierarchical definition of shape functions into an existing code can be rather involved in this work we present a simple way to bring the hierarchical isogeometric concept closer to a standard finite element formulation practically speaking the hierarchy of functions and knot spans is flattened into a sequence of elements being equipped with a standard single level basis in fact the proposed multi level extraction is a generalization of the classical bezier extraction and analogously offers a standard element structure to the hierarchical overlay of functions moreover this approach is suitable for an extension to non linear problems and for a parallel implementation the multi level extraction is presented as a general concept that can be applied to different kinds of refinements and basis functions finally few basic algorithms to compute the local multi level extraction operator for knot insertion on spline spaces are outlined and compared and some numerical examples are presented

subdivision stabilised immersed b spline finite elements for moving boundary flows an immersed finite element method is presented to compute flows with complex moving boundaries on a fixed cartesian grid the viscous incompressible fluid flow equations are discretized with b spline basis functions the two scale relation for b splines is used to implement an elegant and efficient technique to satisfy the lbb condition on non grid aligned fluid domains and at moving boundaries the boundary conditions are enforced with a consistent penalty method as originally proposed by nitsche in addition a special extrapolation technique is employed to prevent the loss of numerical stability in presence of arbitrarily small cut cells the versatility and accuracy of the proposed approach is demonstrated by means of convergence studies and comparisons with previous experimental and computational investigations

adaptive fem based nonrigid image registration using truncated hierarchical b splines we present an efficient approach of finite element method fem based nonrigid image registration in which the spatial transformation is constructed using truncated hierarchical b splines thb splines the image registration framework minimizes an energy functional using an fem based method and thus involves solving a large system of linear equations this framework is carried out on a set of successively refined grids however due to the increased number of control points during subdivision large linear systems are generated which are generally demanding to solve instead of using uniform subdivision an adaptive local refinement scheme is carried out only refining the areas of large change in deformation of the image by incorporating the key advantages of thb spline basis functions such as linear independence partition of unity and reduced overlap into the fem based framework we improve the matrix sparsity and computational efficiency the performance of the proposed method is demonstrated on 2d synthetic and medical images

truncated hierarchical tricubic c0 spline construction on unstructured hexahedral meshes for isogeometric analysis applications abstract we present a new method for truncated hierarchical tricubic spline th spline3d construction to enable adaptive isogeometric analysis on unstructured hexahedral meshes taking the input unstructured hexahedral mesh as the control mesh we first develop blending functions with the aid of tricubic bernstein polynomials this development is an extension of defining bicubic bernstein polynomial blending functions over quadrilateral meshes we further build the hierarchical structure and apply the truncation mechanism to the developed blending functions for highly localized refinement during the refinement of th spline3d high level blending functions are added in the solution space whereas certain low level ones are discarded or truncated depending on the high level subdomain the blending functions are piecewise polynomials that form a partition of unity their support overlapping is also reduced due to the truncation mechanism resulting in sparser stiffness matrices compared to the classical hierarchical refinement th spline3d supports bezier extraction such that it can be easily incorporated into existing finite element frameworks several examples are used to demonstrate the analysis suitability and efficiency of the proposed method

a fictitious domain distributed lagrange multiplier based fluid structure interaction scheme with hierarchical b spline grids abstract we present a numerical scheme for fluid structure interaction based on hierarchical b spline grids and fictitious domain distributed lagrange multipliers the incompressible navier stokes equations are solved over a cartesian grid discretised with b splines the fluid grid near the immersed solids is refined locally using hierarchical b splines the immersed solid is modelled as geometrically exact beam discretised with standard linear lagrange shape functions the kinematic constraint at the fluid solid interface is enforced with distributed lagrange multipliers the unconditionally stable and second order accurate generalised method is used for integration in time for both the fluid and solid domains a fully implicit and fully coupled solution scheme is developed by using the newton raphson method to solve the non linear system of equations obtained with galerkin weak formulation first the spatial and temporal convergence of the proposed scheme is assessed by studying steady and unsteady flow past a fixed cylinder then the scheme is applied to several benchmark problems to demonstrate the efficiency and robustness of the proposed scheme the results obtained with the present scheme are compared with the reference values

weak coupling for isogeometric analysis of non matching and trimmed multi patch geometries nitsche s method can be used as a coupling tool for non matching discretizations by weakly enforcing interface constraints we explore the use of weak coupling based on nitsche s method in the context of higher order and higher continuity b splines and nurbs we demonstrate that weakly coupled spline discretizations do not compromise the accuracy of isogeometric analysis we show that the combination of weak coupling with the finite cell method opens the door for a truly isogeometric treatment of trimmed b spline and nurbs geometries that eliminates the need for costly reparameterization procedures we test our methodology for several relevant technical problems in two and three dimensions such as gluing together trimmed multi patches and connecting non matching meshes that contain b spline basis functions and standard triangular finite elements the results demonstrate that the concept of nitsche based weak coupling in conjunction with the finite cell method has the potential to considerably increase the flexibility of the design through analysis process in isogeometric analysis

on error estimates and adaptivity in elastoplastic solids applications to the numerical simulation of strain localization in classical and cosserat continua the a posteriori error estimates based on the post processing approach are introduced for elastoplastic solids the standard energy norm error estimate established for linear elliptic problems is generalized here to account for the presence of internal variables through the norm associated with the complementary free energy this is known to represent a natural metric for the class of elastoplastic problems of evolution in addition the intrinsic dissipation functional is utilized as a basis for a complementary a posteriori error estimates a posteriori error estimates and adaptive refinement techniques are applied to the finite element analysis of a strain localization problem as a model problem the constitutive equations describing a generalization of standard j2 elastoplasticity within the cosserat continuum are used to overcome serious limitations exhibited by classical continuum models in the post instability region the proposed a posteriori error estimates are appropriately modified to account for the cosserat continuum model and linked with adaptive techniques in order to simulate strain localization problems superior behaviour of the cosserat continuum model in comparison to the classical continuum model is demonstrated through the finite element simulation of the localization in a plane strain tensile test for an elastopiastic softening material resulting in convergent solutions with an h refinement and almost uniform error distribution in all considered error norms

mixed stabilized finite element methods in nonlinear solid mechanics part ii strain localization this paper deals with the question of strain localization associated with materials which exhibit softening due to tensile straining a standard local isotropic rankine damage model with strain softening is used as exemplary constitutive model both the irreducible and mixed forms of the problem are examined and stability and solvability conditions are discussed lack of uniqueness and convergence difficulties related to the strong material nonlinearities involved are also treated from this analysis the issue of local discretization error in the pre localization regime is deemed as the main difficulty to be overcome in the discrete problem focus is placed on low order finite elements with continuous strain and displacement fields triangular p1p1 and quadrilateral q1q1 although the presented approach is very general numerical examples show that the resulting procedure is remarkably robust it does not require the use of auxiliary tracking techniques and the results obtained do not suffer from spurious mesh bias dependence

a local solution approach for adaptive hierarchical refinement in isogeometric analysis abstract adaptively refining the state variables and fast solving the equations are two key issues in isogeometric analysis currently existing research more focuses on the first issue while the efficiency issue is not elaborated when the overall equations are solved at each refinement step this article develops an approach to promote the efficiency by locally solving the equations that govern the refined variables at each refinement step here hierarchical refinement of the state variables in nurbs is adopted and transformational matrix that transforms the basis functions from the current level to the refined level is utilized to handle the relations between the equations at different levels by dividing the coefficient matrix into blocks and estimating the error of state variables for different blocks we represent the global solution method in a hierarchical structure and separately construct and solve equations to obtain the solution for each level in the structure we validate the proposed method by some steady heat transfer and structural analysis problems in domains with different dimensionality we have compared the numerical result accuracy as well as the time cost of the proposed method with those of current available refinement approaches and it turns out to be very promising

on the similarities and differences between classical hierarchical truncated hierarchical and lr b splines abstract smooth spline functions such as b splines and nurbs are already an established technology in the field of computer aided design cad and have in recent years been given a lot of attention from the computer aided engineering cae community the advantages of local refinement are obvious for anyone working in either field and as such several approaches have been proposed among others we find the three strategies classical hierarchical b splines truncated hierarchical b splines and locally refined b splines we will in this paper present these three frameworks and highlight similarities and differences between them in particular we will look at the function space they span and the support of the basis functions we will then analyse the corresponding stiffness and mass matrices in terms of sparsity patterns and conditioning numbers we show that the basis functions in general do not span the same space and that conditioning numbers are comparable moreover we show that the weighting needed by the classical hierarchical basis to maintain partition of unity has significant implications on the conditioning numbers

a stabilised immersed boundary method on hierarchical b spline grids abstract in this work an immersed boundary finite element method is proposed which is based on a hierarchically refined cartesian b spline grid and employs the non symmetric and penalty free version of nitsche s method to enforce the boundary conditions the strategy allows for h and p refinement and employs a so called ghost penalty term to stabilise the cut cells an effective procedure based on hierarchical subdivision and sub cell merging which avoids excessive numbers of quadrature points is used for the integration of the cut cells a basic laplace problem is used to demonstrate the effectiveness of the cut cell stabilisation and of the penalty free nitsche method as well as their impact on accuracy the methodology is also applied to the incompressible navier stokes equations where the supg pspg stabilisation is employed simulations of the lid driven cavity flow and the flow around a cylinder at low reynolds number show the good performance of the methodology excessive ill conditioning of the system matrix is robustly avoided without jeopardising the accuracy at the immersed boundaries or in the field

truncated hierarchical catmull clark subdivision with local refinement in this paper we present a new method termed truncated hierarchical catmull clark subdivision thccs which generalizes truncated hierarchical b splines to control grids of arbitrary topology thccs basis functions satisfy partition of unity are linearly independent and are locally refinable thccs also preserves geometry during adaptive hh refinement and thus inherits the surface continuity of catmull clark subdivision namely c2c2 continuous everywhere except at the local region surrounding extraordinary nodes where the surface continuity is c1c1 adaptive isogeometric analysis is performed with thccs basis functions on a benchmark problem with extraordinary nodes local refinement on complex surfaces is also studied to show potential wide application of the proposed method

adaptive finite element methods in computational mechanics we present a general approach to adaptivity for finite element methods and give applications to linear elasticity non linear elasto plasticity and nonlinear conservation laws including numerical results

a posteriori analysis of a b spline based finite element method for the stationary quasi geostrophic equations of the ocean abstract we provide a posteriori error analysis of a b spline based finite element method for the streamfunction formulation of the stationary quasi geostrophic equations to model the large scale wind driven ocean circulation the analysis gives rise to a posteriori error estimator for more efficient and accurate modeling of a strong western boundary layer via adaptive local refinement of the computational geometries to verify the capability of the error estimator we use hierarchical b splines in conjunction with dorfler marking strategy numerical studies are performed for several benchmark examples on rectangular and l shape ocean geometries

divergence conforming discretization for stokes problem on locally refined meshes using lr b splines abstract to solve the incompressible flow problems using isogeometric analysis the div compatible spline spaces were originally introduced by buffa et al 2011 and later developed by evans 2011 in this paper we extend the div compatible spline spaces with local refinement capability using locally refined lr b splines over rectangular domains we argue that the spline spaces generated on locally refined meshes will satisfy compatibility provided they span the entire function spaces as governed by mourrain 2014 dimension formula we will in this work use the structured refined lr b splines as introduced by johannessen et al 2014 further we consider these div compatible lr b spline spaces to approximate the velocity and pressure fields in mixed discretization for stokes problem and a set of standard benchmark tests are performed to show the stability efficiency and the conservation properties of the discrete velocity fields in adaptive isogeometric analysis

truncated t splines fundamentals and methods in this paper we present truncated t splines as a new type of t splines suitable for both geometric design and analysis supporting highly localized refinement truncated t spline basis functions are piece wise polynomials that are linearly independent and form a partition of unity refinement of truncated t splines produces nested spline spaces furthermore we study truncated t splines and local refinement on the general domain 2 manifold with extraordinary points in the t mesh g1g1 continuity is attained around extraordinary points by properly capping quartic bezier patches where a constrained optimization problem is solved in the end we study benchmark problems using truncated t splines in the context of isogeometric analysis we also apply truncated t splines to complex geometries to show the smooth surfaces and simulation results under local refinement

efficient quadrature for nurbs based isogeometric analysis abstract we initiate the study of efficient quadrature rules for nurbs based isogeometric analysis a rule of thumb emerges the half point rule indicating that optimal rules involve a number of points roughly equal to half the number of degrees of freedom or equivalently half the number of basis functions of the space under consideration the half point rule is independent of the polynomial order of the basis efficient rules require taking into account the precise smoothness of basis functions across element boundaries several rules of practical interest are obtained and a numerical procedure for determining efficient rules is presented we compare the cost of quadrature for typical situations arising in structural mechanics and fluid dynamics the new rules represent improvements over those used previously in isogeometric analysis

stabilized finite element approximation of transient incompressible flows using orthogonal subscales abstract in this paper we present a stabilized finite element method to solve the transient navier stokes equations based on the decomposition of the unknowns into resolvable and subgrid scales the latter are approximately accounted for so as to end up with a stable finite element problem which in particular allows to deal with convection dominated flows and the use of equal velocity pressure interpolations three main issues are addressed the first is a method to estimate the behavior of the stabilization parameters based on a fourier analysis of the problem for the subscales secondly the way to deal with transient problems discretized using a finite difference scheme is discussed finally the treatment of the nonlinear term is also analyzed a very important feature of this work is that the subgrid scales are taken as orthogonal to the finite element space in the transient case this simplifies considerably the numerical scheme

subdivision shells with exact boundary control and non manifold geometry we introduce several new extensions to subdivision shells that provide an improved level of shape control over shell boundaries and facilitate the analysis of shells with non smooth and non manifold joints to this end extended subdivision schemes are used that enable to relax the continuity of the limit surface along prescribed crease edges and to create surfaces with prescribed limit positions and normals furthermore shells with boundaries in the form of conic sections such as circles or parabolas are represented with rational subdivision schemes which are defined in analogy to rational b splines in terms of implementation the difference between the introduced and conventional subdivision schemes is restricted to the use of modified subdivision stencils close to the mentioned geometric features hence the resulting subdivision surface is in most parts of the domain identical to standard smooth subdivision surfaces the particular subdivision scheme used in this paper constitutes an improved version of the original loop s scheme and is as such based on triangular meshes as in the original subdivision shells surfaces created with the modified scheme are used for interpolating the reference and deformed shell configurations at the integration points the subdivision surface is evaluated using a newly developed discrete parameterization approach in the resulting finite elements the only degrees of freedom are the mid surface displacements of the nodes and additional lagrange parameters for enforcing normal constraints the versatility of the newly developed elements is demonstrated with a number of geometrically nonlinear shell examples copyright 2011 john wiley sons ltd

fundamental aspects of shape optimization in the context of isogeometric analysis we develop a mathematical foundation for shape optimization problems under state equation constraints where both state and control are discretized by b splines or nurbs in other words we use isogeometric analysis iga for solving the partial differential equation and a nodal approach to change domains where control points take the place of nodes and where thus a quite general class of functions for representing optimal shapes and their boundaries becomes available the minimization problem is solved by a gradient descent method where the shape gradient will be defined in isogeometric terms this gradient is obtained following two schemes optimize first discretize then and reversely discretize first optimize then we show that for isogeometric analysis the two schemes yield the same discrete system moreover we also formulate shape optimization with respect to nurbs in the optimize first ansatz which amounts to finding optimal control points and weights simultaneously numerical tests illustrate the theory

isogeometric collocation cost comparison with galerkin methods and extension to adaptive hierarchical nurbs discretizations abstract we compare isogeometric collocation with isogeometric galerkin and standard c 0 finite element methods with respect to the cost of forming the matrix and residual vector the cost of direct and iterative solvers the accuracy versus degrees of freedom and the accuracy versus computing time on this basis we show that isogeometric collocation has the potential to increase the computational efficiency of isogeometric analysis and to outperform both isogeometric galerkin and standard c 0 finite element methods when a specified level of accuracy is to be achieved with minimum computational cost we then explore an adaptive isogeometric collocation method that is based on local hierarchical refinement of nurbs basis functions and collocation points derived from the corresponding multi level greville abscissae we introduce the concept of weighted collocation that can be consistently developed from the weighted residual form and the two scale relation of b splines using weighted collocation in the transition regions between hierarchical levels we are able to reliably handle coincident collocation points that naturally occur for multi level greville abscissae the resulting method combines the favorable properties of isogeometric collocation and hierarchical refinement in terms of computational efficiency local adaptivity robustness and straightforward implementation which we illustrate by numerical examples in one two and three dimensions

stabilized finite element method for the transient navier stokes equations based on a pressure gradient projection abstract in this paper we present a stabilized finite element formulation for the transient incompressible navier stokes equations the main idea is to introduce as a new unknown of the problem the projection of the pressure gradient onto the velocity space and to add to the incompresibility equation the difference between the laplacian of the pressure and the divergence of this new vector field this leads to a pressure stabilization effect that allows the use of equal interpolation for both velocities and pressures in the case of the transient equations we consider the possibility of treating the pressure gradient projection either implicitly or explicity in the first case the number of unknowns of the problem is substantially increased with respect to the standard galerkin formulation nevertheless iterative techniques may be used in order to uncouple the calculation of the pressure gradient projection from the rest of unknowns velocity and pressure when this vector field is treated explicitly the increment of computational cost of the stabilized formulation with respect to the galerkin method is very low we provide a stability estimate for the case of the simple backward euler time integration scheme for both the implicit and the explicit treatment of the pressure gradient projection

a stabilized formulation for incompressible elasticity using linear displacement and pressure interpolations in this paper a stabilized finite element method to deal with incompressibility in solid mechanics is presented a mixed formulation involving pressure and displacement fields is used and a continuous linear interpolation is considered for both fields to overcome the babuska brezzi condition a stabilization technique based on the orthogonal sub scale method is introduced the main advantage of the method is the possibility of using linear triangular or tetrahedral finite elements which are easy to generate for real industrial applications results are compared with standard galerkin and q1p0 mixed formulations for nearly incompressible problems in the context of linear elasticity

stabilization of incompressibility and convection through orthogonal sub scales in finite element methods two apparently different forms of dealing with the numerical instability due to the incompressibility constraint of the stokes problem are analyzed in this paper the first of them is the stabilization through the pressure gradient projection which consists of adding a certain least squares form of the difference between the pressure gradient and its l2 projection onto the discrete velocity space in the variational equations of the problem the second is a sub grid scale method whose stabilization effect is very similar to that of the galerkin least squares gls method for the stokes problem it is shown here that the first method can also be recast in the framework of sub grid scale methods with a particular choice for the space of sub scales this leads to a new stabilization procedure whose applicability to stabilize convection is also studied in this paper

imposing dirichlet boundary conditions with nitsche s method and spline based finite elements a key challenge while employing non interpolatory basis functions in finite element methods is the robust imposition of dirichlet boundary conditions the current work studies the weak enforcement of such conditions for b spline basis functions with application to both second and fourth order problems this is achieved using concepts borrowed from nitsche s method which is a stabilized method for imposing constraints on surfaces conditions for the stability of the system of equations are derived for each class of problem stability parameters in the nitsche weak form are then evaluated by solving a local generalized eigenvalue problem at the dirichlet boundary the approach is designed to work equally well when the grid used to build the splines conforms to the physical boundary of interest as well as to the more general case when it does not through several numerical examples the approach is shown to yield optimal rates of convergence copyright 2010 john wiley sons ltd

bezier extraction and adaptive refinement of truncated hierarchical nurbs abstract this contribution presents bezier extraction of truncated hierarchical b splines and the application of the approach to adaptive isogeometric analysis the developed procedures allow for the implementation of hierarchical b splines and nurbs without the need for an explicit truncation of the basis moreover standard procedures of adaptive finite element analysis for error estimation and marking of elements are directly applicable due to the strict use of an element viewpoint starting from a multi level nested mesh that results from uniform h refinement standard bezier extraction is applied to active elements that contribute to the hierarchical approximation this results in a multi level system of equations without communication between individual hierarchy levels a hierarchical subdivision operator is developed to recover this communication by transforming the multi level system of equations into a hierarchical system of equations it is demonstrated that this approach implicitly defines the truncated hierarchical basis in terms of a simple matrix multiplication in this way the implementation effort is reduced to a minimum as shape function routines and bezier extraction procedures remain unchanged compared to standard isogeometric analysis the convergence and the computational efficiency of the approach are examined in three different demonstration problems of heat conduction linear elasticity and the phase field modelling of brittle fracture

triangles and tetrahedra in explicit dynamic codes for solids explicit dynamic codes which are used currently for the study of plastic deformations in impact or with some modification for metal forming suffer two serious limitations first only quadrilateral or hexahedral linear elements can be used thus limiting the possibilities of adaptive refinement and adaptive meshing second even with the use of such elements special devices such as reduced integration must be introduced to avoid locking and reduce costs these necessitate complex hour glass control mending type procedures the main difficulties are those due to the need of treatment of almost incompressible deformation modes recently similar difficulties have been overcome in the context of fluid dynamics soil dynamics and we show here how the processes introduced there can be adopted effectively to the present problem thus allowing an almost unrestricted choice of element interpolations 1998 john wiley sons ltd

a class of mixed assumed strain methods and the method of incompatible modes a three field mixed formulation in terms of displacements stresses and an enhanced strain field is presented which encompasses as a particular case the classical method of incompatible modes within this frame work incompatible elements arise as particular compatible mixed approximations of the enhanced strain field the conditions that the stress interpolation contain piece wise constant functions and be l2 ortho gonal to the enhanced strain interpolation ensure satisfaction of the patch test and allow the elimination of the stress field from the formulation the preceding conditions are formulated in a form particularly convenient for element design as an illustration of the methodology three new elements are developed and shown to exhibit good performance a plane 3d elastic plastic quad an axisymmetric element and a thick plate bending quad the formulation described herein is suitable for non linear analysis

extended truncated hierarchical catmull clark subdivision abstract in this paper we present an extended truncated hierarchical catmull clark subdivision ethccs method which improves the efficiency of local refinement in truncated hierarchical catmull clark subdivision thccs based on the stam s method we first build a set of basis functions over arbitrary quadrilateral meshes and apply them to thccs then a new basis function insertion scheme is developed with the aid of the truncation mechanism which refines one ring neighboring elements rather than two ring neighborhoods therefore ethccs significantly improves the efficiency of local refinement compared with thccs as demonstrated by one benchmark problem and several complex models moreover ethccs is also proved to preserve the input geometry and produce nested spaces