Multilevel refinable triangular PSP-splines (Tri-PSPS)

A multi-level spline technique known as partial shape preserving splines (PSPS) (Li and Tian, 2011) has recently been developed for the design of piecewise polynomial freeform geometric surfaces, where the basis functions of the PSPS can be directly built from an arbitrary set of polygons that partitions a giving parametric domain. This paper addresses a special type of PSPS, the triangular PSPS (Tri-PSPS), where all spline basis functions are constructed from a set of triangles. Compared with other triangular spline techniques, Tri-PSPS have several distinctive features. Firstly, for each given triangle, the corresponding spline basis function for any required degree of smoothness can be expressed in closed-form and directly written out in full explicitly as piecewise bivariate polynomials. Secondly, Tri-PSPS are an additive triangular spline technique, where the spline function built from a given triangle can be replaced with a set of refined spline functions built on a set of smaller triangles that partition the initial given triangle. In addition, Tri-PSPS are a multilevel spline technique, Tri-PSPS surfaces can be designed to have a continuously varying levels of detail, achieved simply by specifying a proper value for the smoothing parameter introduced in the spline functions. In terms of practical implementation, Tri-PSPS are a parallel computing friendly spline scheme, which can be easily implemented on modern programmable GPUs or on high performance computer clusters, since each of the basis functions of Tri-PSPS can be directly computed independent of each other in parallel

Isogeometric FEM-BEM coupled structural-acoustic analysis of shells using subdivision surfaces

We introduce a coupled finite and boundary element formulation for acoustic scattering analysis over thin shell structures. A triangular Loop subdivision surface discretisation is used for both geometry and analysis fields. The Kirchhoff-Love shell equation is discretised with the finite element method and the Helmholtz equation for the acoustic field with the boundary element method. The use of the boundary element formulation allows the elegant handling of infinite domains and precludes the need for volumetric meshing. In the present work the subdivision control meshes for the shell displacements and the acoustic pressures have the same resolution. The corresponding smooth subdivision basis functions have the $C^1$ continuity property required for the Kirchhoff-Love formulation and are highly efficient for the acoustic field computations. We validate the proposed isogeometric formulation through a closed-form solution of acoustic scattering over a thin shell sphere. Furthermore, we demonstrate the ability of the proposed approach to handle complex geometries with arbitrary topology that provides an integrated isogeometric design and analysis workflow for coupled structural-acoustic analysis of shells

Multivariate Splines and Algebraic Geometry

Multivariate splines are effective tools in numerical analysis and approximation theory. Despite an extensive literature on the subject, there remain open questions in finding their dimension, constructing local bases, and determining their approximation power. Much of what is currently known was developed by numerical analysts, using classical methods, in particular the so-called Bernstein-B´ezier techniques. Due to their many interesting structural properties, splines have become of keen interest to researchers in commutative and homological algebra and algebraic geometry. Unfortunately, these communities have not collaborated much. The purpose of the half-size workshop is to intensify the interaction between the different groups by bringing them together. This could lead to essential breakthroughs on several of the above problems

Adaptive isogeometric analysis with hierarchical box splines

Isogeometric analysis is a recently developed framework based on finite element analysis, where the simple building blocks in geometry and solution space are replaced by more complex and geometrically-oriented compounds. Box splines are an established tool to model complex geometry, and form an intermediate approach between classical tensor-product B-splines and splines over triangulations. Local refinement can be achieved by considering hierarchically nested sequences of box spline spaces. Since box splines do not offer special elements to impose boundary conditions for the numerical solution of partial differential equations (PDEs), we discuss a weak treatment of such boundary conditions. Along the domain boundary, an appropriate domain strip is introduced to enforce the boundary conditions in a weak sense. The thickness of the strip is adaptively defined in order to avoid unnecessary computations. Numerical examples show the optimal convergence rate of box splines and their hierarchical variants for the solution of PDEs

Characterization of bivariate hierarchical quartic box splines on a three-directional grid

International audienceWe consider the adaptive refinement of bivariate quartic C 2-smooth box spline spaces on the three-directional (type-I) grid G. The polynomial segments of these box splines belong to a certain subspace of the space of quar-tic polynomials, which will be called the space of special quartics. Given a bounded domain Ω ⊂ R 2 and finite sequence (G ℓ) ℓ=0,...,N of dyadically refined grids, we obtain a hierarchical grid by selecting mutually disjoint cells from all levels such that their union covers the entire domain. Using a suitable selection procedure allows to define a basis spanning the hierarchical box spline space. The paper derives a characterization of this space. Under certain mild assumptions on the hierarchical grid, the hierarchical spline space is shown to contain all C 2-smooth functions whose restrictions to the cells of the hierarchical grid are special quartic polynomials. Thus, in this case we can give an affirmative answer to the completeness questions for the hierarchical box spline basis

Constructing IGA-suitable planar parameterization from complex CAD boundary by domain partition and global/local optimization

In this paper, we propose a general framework for constructing IGA-suitable planar B-spline parameterizations from given complex CAD boundaries consisting of a set of B-spline curves. Instead of forming the computational domain by a simple boundary, planar domains with high genus and more complex boundary curves are considered. Firstly, some pre-processing operations including B\'ezier extraction and subdivision are performed on each boundary curve in order to generate a high-quality planar parameterization; then a robust planar domain partition framework is proposed to construct high-quality patch-meshing results with few singularities from the discrete boundary formed by connecting the end points of the resulting boundary segments. After the topology information generation of quadrilateral decomposition, the optimal placement of interior B\'ezier curves corresponding to the interior edges of the quadrangulation is constructed by a global optimization method to achieve a patch-partition with high quality. Finally, after the imposition of C1=G1-continuity constraints on the interface of neighboring B\'ezier patches with respect to each quad in the quadrangulation, the high-quality B\'ezier patch parameterization is obtained by a C1-constrained local optimization method to achieve uniform and orthogonal iso-parametric structures while keeping the continuity conditions between patches. The efficiency and robustness of the proposed method are demonstrated by several examples which are compared to results obtained by the skeleton-based parameterization approach

Injective hierarchical free-form deformations using THB-splines

The free-form deformation (FFD) method deforms geometry in n-dimensional space by employing an n-variate function to deform (parts of) the ambient space. The original method pioneered by Sederberg and Parry in 1986 uses trivariate tensor-product Bernstein polynomials in R3 and is controlled as a Bézier volume. We propose an extension based on truncated hierarchical B-splines (THB-splines). This offers hierarchical and local refinability, an efficient implementation due to reduced supports of THB-splines, and intuitive control point hiding during FFD interaction. Additionally, we address the issue of fold-overs by efficiently checking the injectivity of the hierarchical deformation in real-time