Stable splitting of bivariate spline spaces by Bernstein-Bézier methods

We develop stable splitting of the minimal determining sets for the spaces of bivariate C1 splines on triangulations, including a modified Argyris space, Clough-Tocher, Powell-Sabin and quadrilateral macro-element spaces. This leads to the stable splitting of the corresponding bases as required in Böhmer's method for solving fully nonlinear elliptic PDEs on polygonal domains

Gaussian quadrature for C1C^1 cubic Clough-Tocher macro-triangles

A numerical integration rule for multivariate cubic polynomials over n-dimensional simplices was designed by Hammer and Stroud [14]. The quadrature rule requires n + 2 quadrature points: the barycentre of the simplex and n + 1 points that lie on the connecting lines between the barycentre and the vertices of the simplex. In the planar case, this particular rule belongs to a two-parameter family of quadrature rules that admit exact integration of bivariate polynomials of total degree three over triangles. We prove that this rule is exact for a larger space, namely the C1 cubic Clough-Tocher spline space over macro-triangles if and only if the split-point is the barycentre. This results into a factor of three reduction in the number of quadrature points needed to integrate the Clough-Tocher spline space exactly

Watertight conversion of trimmed CAD surfaces to Clough-Tocher splines

The boundary representations (B-reps) that are used to represent shape in Computer-Aided Design systems create unavoidable gaps at the face boundaries of a model. Although these inconsistencies can be kept below the scale that is important for visualisation and manufacture, they cause problems for many downstream tasks, making it difficult to use CAD models directly for simulation or advanced geometric analysis, for example. Motivated by this need for watertight models, we address the problem of converting B-rep models to a collection of cubic C1C1 Clough–Tocher splines. These splines allow a watertight join between B-rep faces, provide a homogeneous representation of shape, and also support local adaptivity. We perform a comparative study of the most prominent Clough–Tocher constructions and include some novel variants. Our criteria include visual fairness, invariance to affine reparameterisations, polynomial precision and approximation error. The constructions are tested on both synthetic data and CAD models that have been triangulated. Our results show that no construction is optimal in every scenario, with surface quality depending heavily on the triangulation and parameterisation that are used.This research was supported by the Engineering and Physical Sciences Research Council through Grant EP/K503757/1.This is the final version. It was first published by Elsevier at http://www.sciencedirect.com/science/article/pii/S0167839615000795

B-spline-like bases for C2C^2 cubics on the Powell-Sabin 12-split

For spaces of constant, linear, and quadratic splines of maximal smoothness on the Powell-Sabin 12-split of a triangle, the so-called S-bases were recently introduced. These are simplex spline bases with B-spline-like properties on the 12-split of a single triangle, which are tied together across triangles in a B\'ezier-like manner. In this paper we give a formal definition of an S-basis in terms of certain basic properties. We proceed to investigate the existence of S-bases for the aforementioned spaces and additionally the cubic case, resulting in an exhaustive list. From their nature as simplex splines, we derive simple differentiation and recurrence formulas to other S-bases. We establish a Marsden identity that gives rise to various quasi-interpolants and domain points forming an intuitive control net, in terms of which conditions for $C^0$-, $C^1$-, and $C^2$-smoothness are derived

On numerical quadrature for C1C^1 quadratic Powell-Sabin 6-split macro-triangles

The quadrature rule of Hammer and Stroud [16] for cubic polynomials has been shown to be exact for a larger space of functions, namely the $C^1$ cubic Clough-Tocher spline space over a macro-triangle if and only if the split-point is the barycentre of the macro-triangle [21]. We continue the study of quadrature rules for spline spaces over macro-triangles, now focusing on the case of $C^1$ quadratic Powell-Sabin 6-split macro-triangles. We show that the $3$-node Gaussian quadrature(s) for quadratics can be generalised to the $C^1$ quadratic Powell-Sabin 6-split spline space over a macro-triangle for a two-parameter family of inner split-points, not just the barycentre as in [21]. The choice of the inner split-point uniquely determines the positions of the edge split-points such that the whole spline space is integrated exactly by a corresponding polynomial quadrature. Consequently, the number of quadrature points needed to exactly integrate this special spline space reduces from twelve to three. For the inner split-point at the barycentre, we prove that the two 3-node quadratic polynomial quadratures of Hammer and Stroud exactly integrate also the $C^1$ quadratic Powell-Sabin spline space if and only if the edge split-points are at their respective edge midpoints. For other positions of the inner and edge split-points we provide numerical examples showing that three nodes suffice to integrate the space exactly, but a full classification and a closed-form solution in the generic case remain elusive

A C r Trivariate Macro-Element Based on the Alfeld Split of Tetrahedra

Abstract We construct trivariate macro-elements of class C r for any r ≥ 1 over the Alfeld refinement of any tetrahedral partition in R 3 . In our construction, the degree of polynomials used for these macro-elements is the lowest possible. We also give the dimension formula for the subspace of consisting of these macro-elements

Exact conversion from Bézier tetrahedra to Bézier hexahedra

International audienceModeling and computing of trivariate parametric volumes is an important research topic in the field of three-dimensional isogeo-metric analysis. In this paper, we propose two kinds of exact conversion approaches from Bézier tetrahedra to Bézier hexahedra with the same degree by reparametrization technique. In the first method, a Bézier tetrahedron is converted into a degenerate Bézier hexahedron, and in the second approach, a non-degenerate Bézier tetrahedron is converted into four non-degenerate Bézier hexahedra. For the proposed methods, explicit formulas are given to compute the control points of the resulting tensor-product Bézier hexahedra. Furthermore, in the second method, we prove that tetrahedral spline solids with C k-continuity can be converted into a set of tensor-product Bézier volumes with G k-continuity. The proposed methods can be used for the volumetric data exchange problems between different trivariate spline representations in CAD/CAE. Several experimental results are presented to show the effectiveness of the proposed methods