IGA. Robin Bouclier

Чтение книги онлайн.

Figure 1.5. Examples of B-spline basis functions of different degrees. The same interior knots were used (0.15, 0.50 and 0.75)

For illustration purposes, we display several examples of B-spline basis functions in Figure 1.5. We can observe the influence of the degree. Many properties of these B-spline basis functions are given in pioneering works (Cohen et al. 1980, 2001; Piegl and Tiller 1997; Farin 2002). In particular, an interesting feature of B-splines is their higher degree of regularity. Indeed, a B-spline function of order p can reach a Cp−¹ regularity at knot ξi₁, if this one is single. There can be more than one knot at a given location of the parameter space. More precisely, if m is the multiplicity of a given knot, the functions are Cp−m continuous at that location, which is in contrast to standard FE, where only a C⁰ regularity is encountered on the element boundaries. Figure 1.6 highlights this point: the second derivatives of the cubic basis functions are discontinuous at the double knot, while they remain continuous at single knots.

Figure 1.6. Differentiability of cubic B-spline functions in the case of a uniform knot-vector (interior knots: 2, 4 and 6) (left), and in the case of a double knot (interior knots: 2.4, 2.4 and 6) (right)

Last but not least, it may be noted that the computation of the kth derivatives of a B-spline basis can be efficiently performed, similarly to that for building the basis functions themselves using [1.13]. Indeed, the derivatives of B-spline basis functions are represented in terms of B-spline lower-degree bases:

      [1.14]

      Again, we refer to Piegl and Tiller (1997) for efficient algorithms regarding these calculations.

NOTE.– Looking more closely at the differentiability properties of B-Spline functions, it may be noted that the B-spline technology encompasses the standard FE technology based on Lagrange polynomials. Indeed, if we seek to recover the subspace of standard Lagrange polynomials with the B-spline formalism, we only need to take a knot-vector where the multiplicity of each (interior) knot is equal to p. This would generate piecewise polynomials of degree p that are C⁰ at the element boundaries, which is equivalent to taking standard Lagrange polynomials to generate the approximation subspace (see Bouclier and Passieux (2022) for more information regarding the relation between FEM and IGA and, for example, Schillinger et al. (2016b) and Tirvaudey (2019)). Figure 1.5 (top left) clearly illustrates this statement for the case where p = 1: the same basis functions as standard FEM are constructed. As a result, we will obviously consider splines of degree p ≥ 2 and, as much as possible, interior knots of multiplicity 1, in order to benefit from the increased smoothness of this technology.

      1.2.2.3. NURBS curves

      The NURBS basis functions, denoted by in the following, are defined from a B-spline basis of the same degree . In order to do so, it is necessary to introduce the weights wi for the control points xi. More precisely, the piecewise rational functions read as follows:

[1.15]

      Then, NURBS curves are defined similarly to those in [1.9]:

      [1.16]

      Consequently, a NURBS curve takes the same inputs as a B-spline curve,

i.e. a degree p, a knot-vector Ξ₁ and control points Pi, plus an additional set of n₁ weights wi. The weights give more shape control to the designer. By changing their values, a curve can be modified, as shown in Figure 1.7, for example. In fact, the positions of the control points and the values of the associated weights can be adjusted in order to build conical sections exactly (circles, ellipses, etc.; see, for example, Cottrell et al. (2007, 2009)). Given equation [1.15] (and verifying that the B-spline functions satisfy the partition of unity), it may be noted that if all weights are equal, the NURBS entity turns out to be a B-spline entity.

Figure 1.7. An example of a NURBS curve: influence of the weights on the basis functions and on the final curve with comparison to the B-spline case (the weights СКАЧАТЬ