
Designing complex surfaces is one of the major problems in industries such as the automotive, shipbuilding and aerospace industries. To solve this problem, continuity conditions between surfaces are applied to construct the complex surfaces. The geometric and parametric continuities are the two metrics that usually have been used in connecting surfaces. However, the conventional geometric and parametric continuities have significant limitations. The existing continuity conditions only allow the two surfaces to be joined at the end of the boundary point. Therefore, if the designers want to connect at any arbitrary line of the first surface, the designers must use the subdivision method to splice the surfaces. Nevertheless, this method is tedious and involves a high computational cost, especially when dealing with a higher degree order of surfaces. Thus, this paper presents fractional continuity of degree two (or F2) for generalized fractional Bézier surfaces. The fractional parameter embedded in the generalized fractional Bézier basis functions will solve the mentioned limitation by introducing fractional continuity. The generalized fractional Bézier surface also has excellent shape parameters that can alter the shape of the surface without changing the control points. Thus, the shape parameters enable the control of the shape flexibility of the surfaces, while fractional parameters enable the control of the adjustability of the surfaces' size. The F2 continuity for generalized fractional Bézier surfaces can become an easier and faster alternative to the subdivision method. Therefore, the fractional continuity for generalized fractional Bézier surfaces will be a good tool to generate complex surfaces due to its flexibility and adjustability of shape and fractional parameters.
Citation: Syed Ahmad Aidil Adha Said Mad Zain, Md Yushalify Misro. A novel technique on flexibility and adjustability of generalized fractional Bézier surface patch[J]. AIMS Mathematics, 2023, 8(1): 550-589. doi: 10.3934/math.2023026
[1] | Salwa Syazwani Mahzir, Md Yushalify Misro, Kenjiro T. Miura . Preserving monotone or convex data using quintic trigonometric Bézier curves. AIMS Mathematics, 2024, 9(3): 5971-5994. doi: 10.3934/math.2024292 |
[2] | Qi Xie, Yiting Huang . A class of generalized quadratic B-splines with local controlling functions. AIMS Mathematics, 2023, 8(10): 23472-23499. doi: 10.3934/math.20231193 |
[3] | Samia BiBi, Md Yushalify Misro, Muhammad Abbas . Smooth path planning via cubic GHT-Bézier spiral curves based on shortest distance, bending energy and curvature variation energy. AIMS Mathematics, 2021, 6(8): 8625-8641. doi: 10.3934/math.2021501 |
[4] | Joseph Lifton, Tong Liu, John McBride . Non-linear least squares fitting of Bézier surfaces to unstructured point clouds. AIMS Mathematics, 2021, 6(4): 3142-3159. doi: 10.3934/math.2021190 |
[5] | Xuanyi Zhao, Jinggai Li, Shiqi He, Chungang Zhu . Geometric conditions for injectivity of 3D Bézier volumes. AIMS Mathematics, 2021, 6(11): 11974-11988. doi: 10.3934/math.2021694 |
[6] | Mohammad Ayman-Mursaleen, Md. Nasiruzzaman, Nadeem Rao, Mohammad Dilshad, Kottakkaran Sooppy Nisar . Approximation by the modified $ \lambda $-Bernstein-polynomial in terms of basis function. AIMS Mathematics, 2024, 9(2): 4409-4426. doi: 10.3934/math.2024217 |
[7] | Yujun Li . Characteristics of planar sextic indirect-PH curves. AIMS Mathematics, 2024, 9(1): 2215-2231. doi: 10.3934/math.2024110 |
[8] | Lanyin Sun, Kunkun Pang . Numerical solution of unsteady elastic equations with C-Bézier basis functions. AIMS Mathematics, 2024, 9(1): 702-722. doi: 10.3934/math.2024036 |
[9] | Guorong Zhou, Qing-Bo Cai . Bivariate $ \lambda $-Bernstein operators on triangular domain. AIMS Mathematics, 2024, 9(6): 14405-14424. doi: 10.3934/math.2024700 |
[10] | Muhsin Incesu . LS (3)-equivalence conditions of control points and application to spatial Bézier curves and surfaces. AIMS Mathematics, 2020, 5(2): 1216-1246. doi: 10.3934/math.2020084 |
Designing complex surfaces is one of the major problems in industries such as the automotive, shipbuilding and aerospace industries. To solve this problem, continuity conditions between surfaces are applied to construct the complex surfaces. The geometric and parametric continuities are the two metrics that usually have been used in connecting surfaces. However, the conventional geometric and parametric continuities have significant limitations. The existing continuity conditions only allow the two surfaces to be joined at the end of the boundary point. Therefore, if the designers want to connect at any arbitrary line of the first surface, the designers must use the subdivision method to splice the surfaces. Nevertheless, this method is tedious and involves a high computational cost, especially when dealing with a higher degree order of surfaces. Thus, this paper presents fractional continuity of degree two (or F2) for generalized fractional Bézier surfaces. The fractional parameter embedded in the generalized fractional Bézier basis functions will solve the mentioned limitation by introducing fractional continuity. The generalized fractional Bézier surface also has excellent shape parameters that can alter the shape of the surface without changing the control points. Thus, the shape parameters enable the control of the shape flexibility of the surfaces, while fractional parameters enable the control of the adjustability of the surfaces' size. The F2 continuity for generalized fractional Bézier surfaces can become an easier and faster alternative to the subdivision method. Therefore, the fractional continuity for generalized fractional Bézier surfaces will be a good tool to generate complex surfaces due to its flexibility and adjustability of shape and fractional parameters.
In Computer Aided Geometric Design (CAGD) and Computer-Aided Manufacturing, the development of curves and surfaces is crucial. Researchers studied the representation of free-form curves and surfaces, applying them in numerous fields, such as engineering, manufacturing and designing. In real-world problems, most curves and surfaces are complex to model. Therefore, having a tool to create flexible and adjustable curves and surfaces becomes necessary to model any complex shapes. However, finding a representation of desired complex curves and surfaces is difficult and impractical [2,3]. To tackle the problem of modeling the complex shapes, the curves and surfaces are broken down into simpler curves and patches, joining them to form the complex curves or surfaces. To join such curves or surfaces, a continuity concept is applied to the adjacent curves or surfaces.
Generally, the higher the order of continuity, the smoother the curves or surfaces are connected. Therefore, connecting the curves or surfaces is important with a high degree of continuity. Two standard metrics have been used in connecting these curves and surfaces: parametric and geometric continuity. Parametric continuity of order r, or Cr continuity, is the simplest form of the continuity concept. Nonetheless, parametric continuity has its restrictions. For example, when the surfaces are connected by C1 continuity, they still need to possess a common tangent at their boundary points. Hence, parametric continuity cannot be the only exact standard method in constructing smooth curves or surfaces [4]. As a result, researchers have developed an upgraded version of parametric continuity called geometric continuity, or Gr continuity. Gr continuity is the less restrictive form where scale factors are embedded in the continuity. These scale factors overcome a common tangent for the curves or surfaces.
For shape designing and geometric representation in CAGD, parametric curves and surfaces are useful tools. Classical Bézier curves and surfaces are one of many parametric curves or surfaces. The Bézier method has become one of the famous methods in modeling curves and surfaces due to the simple formulation and excellent geometric properties [5]. However, there is a constraint in the classical Bézier method; the shape of curves or surfaces cannot be changed without changing the control points. To overcome this constraint, researchers have developed aesthetic Bézier curves and surfaces that still preserve the excellent features of classical Bézier curves and surfaces but with added flexibility and adjustability. Rational Bézier curve is one of the aesthetic curves that was introduced earlier. The weight factors in rational Bézier curves and surfaces allow shape modification without altering the control points [6]. However, compared to the classical Bézier curves and surfaces, the rational Bézier curves and surfaces have considerably more difficult computation, convoluted integrals and repeated differentiation, since rational functions define them [7].
To overcome the limitations of the rational Bézier, scholars have created several Bézier curves and surfaces with shape parameters to preserve the benefits of the Bézier model and increase the shape adjustability of curves and surfaces. With three shape parameters, [7] created enhanced classical Bézier curves and surfaces dubbed the shape-adjustable generalized Bézier (or SG-Bézier for short). For example, [8] presented a class of quasi-quintic trigonometric Bézier curves with two shape parameters. Moreover, [9] produced quintic trigonometric Bézier basis functions with two shape parameters. The generalized Hybrid-Trigonometric Bézier (GHT-Bézier) basis functions with embedding exponential and trigonometric functions along with shape parameters have been formulated by [10]. On the other hand, [11] proposed new cubic hyperbolic Bézier basis functions derived from hyperbolic functions.
Researchers have used these aesthetic Bézier curves to develop complex curves and surfaces using continuity. GC1 continuity conditions between adjacent rectangular and triangular Bézier surface patches have been developed by [12]. On the other hand, G1 continuity of piecewise Bézier surfaces has been discussed in [13]. Meanwhile, [14] discussed G2 continuity for conditions for generalized Bézier-like surfaces with multiple shape parameters. Surface construction using C2 continuity for quintic trigonometric Bézier curves is proposed in [15]. Furthermore, [16] proposed G2 continuity for the biquintic trigonometric surface and modeled a few engineering surfaces. In addition, SG-Bézier surfaces have been modeled up until G2 continuity in [17], while [18] constructed the GHT-Bézier surfaces up until G2 continuity. Continuity also has been used in the modeling of the developable surface, while developable λ-Bézier surfaces have been discussed with their geometric design and continuity conditions up until G2 in [19]. Apart from that, [20] introduced the generalized developable cubic trigonometric surface up to G2 continuity. Moreover, [21] developed a generalized developable hybrid trigonometric surface up until G3 continuity.
Cr and Gr continuity enable modeling complicated curves and surfaces by joining several simple curves and surfaces to form the desired curves and surfaces. Nevertheless, the existing continuity has significant limitations, in which the curves or surfaces can only be connected at the endpoints or boundary points. For example, if the designers want to connect the midpoint of the first curve to the second curve, it can only be done by splicing the first curve using the subdivision method first. Then, after the curve splicing, the new curve will join the second curve using continuity. The subdivision method is good for splicing the curve or surfaces, but it has a high computational time when involving a higher degree curve.
Recently, [22] has developed an improved version of geometric continuity dubbed fractional continuity or Fr continuity and connects the curves up until F2. Here, the fractional continuity enables the curves or surfaces to be connected at the endpoints or boundary points and any arbitrary point on the first curve and arbitrary line (according to u or v direction) of the first surface as well. The fractional parameter embedded in fractional continuity is the key to the development of fractional continuity. This paper will discuss the F2 continuity for the generalized fractional Bézier surfaces. The fractional parameter will become a key point in the fractional continuity. By varying the fractional parameter, the common boundary between two surfaces can be adjusted along the first surface while preserving the degree of continuity. Shape parameters are also included in the generalized fractional Bézier surface. Hence, shape parameters will be added as an additional tool for the designer to change the shape of the surface.
The work is organized as follows. The generalized fractional Bézier basis functions are defined in Section 2. Section 3 discussed F2 continuity conditions between three consecutive curves. In Section 4, the definition and the properties of the generalized fractional Bézier tensor product surface are explained. Next, in Section 5, the F2 continuity conditions for the generalized fractional Bézier surfaces are discussed. The procedures and some examples for the F2 continuity will be shown in Section 6. Last but not least, in Section 7, some conclusions and recommendations for future work will be addressed.
Generalized fractional Bézier basis functions
Definition 2.1 (Riemann-Liouville fractional integral for f(t)=t). Let Re(w)>0 and f(t)=t. Then, for t>0, the Riemann-Liouville fractional integral of f(t)=t of order w is as follows:
D−wt(t)=1Γ(w)∫t0(t−x)w−1(x)dx=1Γ(w+2)tw+1. | (2.1) |
The general form and derivation of the definition can be seen in [1].
Definition 2.2 (Generalized fractional Bézier basis functions). For t∈[0,1], ai∈R and w≥0, the following function is defined as the generalized fractional Bézier curve basis function of degree n with n shape parameters:
ˉfi,n(t)=fi,n(t)(1+ain−i+1(1−D−wt(t))−ai+1i+1D−wt(t)),−(n−i+1)<ai<i,a0=an+1=0,i=0,1,...,n, | (2.2) |
where fi,n(t)=(ni)(1−D−wt(t))n−iD−wt(t)i and D−wt(t)=1Γ(w+2)tw+1.
For n≥2, the generalized fractional Bézier basis functions can be defined recursively as follows:
ˉfi,n(t)=(1−D−wt(t))ˉfi,n−1(t)+D−wt(t)ˉfi−1,n−1(t). | (2.3) |
Theorem 2.1 (Generalized fractional Bézier basis function). The generalized fractional basis functions with shape parameters have the following properties, such as degeneracy, non-negativity, a partition of unity, symmetry and linear independence.
Proof. The proof for degeneracy, non-negativity, a partition of unity, and symmetry can be seen in [22]. The linear independence property will be proven as follows:
Linear independence: ∑ni=0ˉfi,n(t)ci=0⟺ci=0,i=0,1,…,n;n≥1. The sufficient condition is clear, and the principle of mathematical induction will be used for the necessary condition.
Base step: n=1.
1∑i=0ˉfi,1(t)ci=0, | (2.4) |
c0ˉf0,1(t)+c1ˉf1,1(t)=0. | (2.5) |
Substituting Eq (2.2) into the above equation yields:
c0(1−D−wt(t))(1−a1D−wt(t))+c1D−wt(t)(1+a1(1−D−wt(t))=0, | (2.6) |
c0(1−D−wt(t)−a1D−wt(t)+a1(D−wt(t))2)+c1(D−wt(t)+a1D−wt(t)−a1(D−wt(t))2)=0, | (2.7) |
c0+((−c0+c1)+(c0+c1)a1)D−wt(t)+(c0+c1)a1(D−wt(t))2=0. | (2.8) |
By comparing coefficients, the following is obtained:
{c0=0,(−c0+c1)+(c0+c1)a1=0,(c0+c1)a1=0. | (2.9) |
By simplifying the above equations, c0=c1=0 is obtained. Hence, the base case is true. Now, the inductive step will be proven:
Inductive step: n=k is true; show for n=k+1.
Assume that ∑ki=0ˉfi,k(t)ci=0⟺ci=0,i=0,1,…,k;k≥1 is true. Consider
k+1∑i=0ˉfi,k+1(t)ci=0. | (2.10) |
Substituting Eq (2.3) in Eq (2.10) and rearranging it, the following is obtained:
(1−D−wt(t))k+1∑i=0ˉfi,k(t)ci+D−wt(t)k+1∑i=0ˉfi−1,k(t)ci=0. | (2.11) |
Since D−wt(t) is an arbitrary value on interval t∈[0,1],
k+1∑i=0ˉfi,k(t)ci=k∑i=0ˉfi,k(t)ci+ˉfk+1,k(t)ck+1=k∑i=0ˉfi,k(t)ci=0, | (2.12) |
k+1∑i=0ˉfi−1,k(t)ci=ˉf−1,k(t)c0+k+1∑i=1ˉfi−1,k(t)ci=k+1∑i=1ˉfi−1,k(t)ci=0. | (2.13) |
Note that ˉfi,k(t)=0 for i=−1 or i>k. Hence, from hypothesis, we obtain ci=0,i=0,1,…,k for Eq (2.12) and ci=0,i=1,2,…,k+1 for Eq (2.13). By the principle of mathematical induction, linear independence is satisfied. It is clear that the generalized fractional Bézier basis functions span any polynomial equation with at most degree n(w+1)+1. Therefore, the generalized fractional Bézier basis functions are a basis for polynomial space with at most degree n(w+1)+1.
Figure 1 illustrates the quadratic fractional Bézier basis functions with multiple spectrum values of shape and fractional parameters. Figure 1 also indicates the same as the classical Bernstein basis functions, especially when the values of the shape and fractional parameters are zero.
Definition 3.1 (Generalized fractional Bézier curve). The generalized fractional Bézier curve of n degree with n shape parameters is defined as follows:
C(t;w,a1,a2,...,an)=n∑i=0ˉfi,n(t;w,a1,a2,...,an)Pi,t∈[0,1], | (3.1) |
where Pi for i=0,1,...,n is the set of control points in R2.
Theorem 3.1 (Properties of generalized fractional Bézier curve). The generalized fractional Bézier curve has the following properties:
(1) Endpoint terminal. The generalized fractional Bézier curves are always interpolated their first control points (t=0) regardless the values of w. For the last control point, the curve will interpolated the last control point (t=1) when w=0. However, for w>0, the endpoint for t=1 can still be calculated using the given equation below:
{C(0,w,a1,a2,...,an)=P0,C(1,w,a1,a2,...,an)=∑ni=0ˉfi,n(1;w,a1,a2,...,an)Pi,C(1,0,a1,a2,...,an)=Pn. | (3.2) |
(2) Endpoint tangent. The endpoint tangent for the generalized fractional Bézier curve at t=0 is independent of the w. However for t=1, the endpoint tangent depends on w. The equations for endpoint tangent are given as follows:
{C′(0,w,a1,a2,...,an)=(n+a1)(P1−P0),C′(1,w,a1,a2,...,an)=ddt(∑ni=0ˉfi,n(t;w,a1,a2,...,an)Pi)t=1,C′(1,0,a1,a2,...,an)=(n−an)(Pn−Pn−1). | (3.3) |
(3) Shape adjustable property. The numbers of shape parameters of the generalized fractional Bézier curve depend on the curve's degree. These shape parameters allow local control of the curve's shape.
(4) Fractional curve adjustable property. One fractional parameter is associated with the generalized fractional Bézier curve. The fractional parameter allows the control of the curve's length adjustability along the u direction. Please refer Section 3.2 for more details.
Proof. The rest of the properties and their proofs can be seen in [22].
The fractional parameter provides the curve adjustable property, as in Theorem 3.1. In Eqs (3.2) and (3.3), it is clear that for t=0, the fractional parameter does not affect both equations. However, for t=1, both equations depend on the fractional parameter. This implies that at t=1 and w>0, the curve does not interpolate at the last control point, Pn. However, the new endpoint (position/coordinate) can be calculated using the equation below:
C(1,w,a1,a2,...,an)=n∑i=0ˉfi,n(1;w,a1,a2,...,an)Pi. | (3.4) |
Example 3.1. Figure 2 shows the curve adjustability of the fractional parameter, w1 for the cubic fractional Bézier curve, C1(t;w1). In this figure, the shape parameters used are (a1,a2,a3)=(−1,0,1), and the control points are represented by the red points. Each sub-figure has an arrow that is pointing to the endpoints of the curve with their respective value. The arrow is pointing towards endpoint at t=1 with arbitrary value of w1. Note that the interval of t remains on the interval [0,1] regardless of the value of w1.
Let three consecutive generalized fractional Bézier curves be as follows:
{C1(t;w1,a1,a2,...,an)=∑ni=0ˉfi,n(t;v1,a1,a2,...,an)Pi,t∈[0,1],n≥2,C2(t;w2,b1,b2,...,bm)=∑mj=0ˉfj,m(t;w2,b1,b2,...,bm)Qj,t∈[0,1],m≥2,C3(t;w3,c1,c2,...,cm)=∑lk=0ˉfk,l(t;w3,c1,c2,...,cl)Rk,t∈[0,1],l≥2, | (3.5) |
where
(1) ˉfi,n, ˉfj,m and ˉfk,l are the generalized fractional Bézier basis functions of degree n, m and l, respectively.
(2) Pi(i=0,1,...,n), Qj(j=0,1,...,m) and Rk(k=0,1,...,l) are the control points.
(3) a1,a2,…,an, b1,b2,…,bm and c1,c2,…,cl are the shape parameters.
(4) w1, w2 and w3 are the fractional parameters.
Definition 3.2 (Fractional continuity for three consecutive generalized fractional Bézier curves, Fr). Suppose three curves C1(t;w1), C2(t;w2) and C3(t;w3) are Fr continuous and will connect to each other consecutively, i.e., C1 connects to C2, and C2 connects to C3. Fr continuity between three consecutive curves can be achieved if the following conditions are satisfied:
{C1(1;w1)=C2(0;w2=0)C′1(1;w1)=ϕ1C′2(0;w2=0)C″1(1;w1)=ϕ21C″2(0;w2=0)+ϕ2C′3(0;w3=0)⋮C(r)1(1;w1)=ϕr1C(r)2(0;w2=0)+ϕr−12C(r−1)2(0;w2=0)+...+ϕrC′2(0;w2=0). | (3.6) |
{C2(1;w2)=C3(0;w3=0)C′2(1;w2)=χ1C′3(0;w3=0)C″2(1;w2)=χ21C″3(0;w3=0)+χ2C′3(0;w3=0)⋮C(r)2(1;w2)=χr1C(r)3(0;w3=0)+χr−12C(r−1)3(0;w3=0)+...+χrC′3(0;w3=0). | (3.7) |
The range of the scale factor, ϕi,χi for i=0,1,...,r depends on the degree of continuity, r. Generally, for Fr continuity, the range for the scale factors are ϕ1,χ1>0 and ϕi,χi∈R for i=2,3,...,r.
Theorem 3.2 (Fractional continuity for generalized fractional Bézier curve, Fr). Consider three generalized fractional Bézier curves as in Eq (3.5), and the necessary and sufficient conditions for fractional continuity at the joint points are given as follows.
(1) F0 continuity:
Q0=n∑iˉfi,n(1;w1,a1,a2,...,an)Pi, | (3.8) |
R0=m∑jˉfj,n(1;w2,b1,b2,...,bm)Qi. | (3.9) |
(2) F1 continuity:
{Q0=∑niˉfi,n(1;w1,a1,a2,...,an)Pi,Q1=1ϕ1(m+b1)(ddt(∑niˉfi,n(t;w1,a1,a2,...,an)Pi)t=1)+Q0. | (3.10) |
{R0=∑mjˉfj,m(1;w2,b1,b2,...,bm)Qi,R1=1χ1(l+c1)(ddt(∑mjˉfj,m(t;w2,b1,b2,...,bm)Qi)t=1)+R0. | (3.11) |
(3) F2 continuity:
{Q0=∑niˉfi,n(1;w1,a1,a2,...,an)Pi,Q1=1ϕ1(m+b1)(ddt(∑niˉfi,n(t;w1,a1,a2,...,an)Pi)t=1)+Q0,Q2=1mϕ21(b2+m−1)(d2dt2(∑niˉfi,n(t;w1,a1,a2,...,an)Pi)t=1+(m((1−m)ϕ21)+ϕ2)+(−2mϕ21+ϕ2)b1)Q0+((2mϕ21−ϕ2)b1+m((−1+m)2ϕ21−ϕ2+ϕ21b2)Q1). | (3.12) |
{R0=∑mjˉfj,n(1;w2,b1,b2,...,bm)Qi,R1=1χ1(l+c1)(ddt(∑mjˉfj,n(t;w2,b1,b2,...,bm)Qi)t=1)+R0,R2=1lχ21(c2+l−1)(d2dt2(∑niˉfj,m(t;w2,b1,b2,...,bm)Qi)t=1+(l((1−l)χ21)+χ2)+(−2lχ21+χ2)c1)R0+((2lχ21−χ2)c1+l((−1+l)2χ21−χ2+χ21l2)R1). | (3.13) |
Proof. The theorem can be proven as follows:
By using the endpoint terminal and tangent properties in Theorem 3.1, Theorem 3.2 can be derived. From Eqs (3.2) and (3.3), the endpoint terminal and tangent for t=0 are independent of w. This implies that for the right-hand side of the equation of Definition 3.2, the value of the fractional parameter can be set to zero to simplify the equations. Generally, if the curves C1,C2…,Cn are connected consecutively, the fractional parameters w1,w2,…,wn−1 (fractional parameters of the first curve, C1 until (n−1)th curve, Cn−1) are needed to satisfy the continuity conditions, while the last fractional parameter, wn, can be used to control the curve adjustability of the last curve, Cn.
(1) F0 continuity condition is obtained by solving Q0 and R0 in C1(1;w1)=C2(0;w2=0) and C2(1;w2)=C3(0;w3=0), respectively.
(2) F0 continuity must be achieved first. F1 continuity condition is obtained by solving for Q1 and R1 in C′1(1;w1)=ϕ1C′2(0;w2=0) and C′2(1;w2)=χ1C′3(0;w3=0), respectively.
(3) F1 continuity must be satisfied first. F2 continuity condition is achieved by solving for Q2 and R2 in C″1(1;w1)=ϕ21C″2(0;w2=0)+ϕ2C′2(0;w2=0) and C″2(1;w2)=χ21C″3(0;w3=0)+χ2C′3(0;w3=0), respectively.
Example 3.2. Figures 3 and 4 depict the adjustment of the common point via w1 and the curve adjustability of second curve via w2, respectively. Blue and orange curves with degree three are C1 and C2, respectively. The control points and shape parameters are the same as in Example 3.1. Q0, Q1 and Q2 can be calculated by Theorem 3.2. Other control point, shape parameters and scale factors are Q3=(3,−5), (b1,b2,b3)=(0.5,0,0.5), ϕ1=0.75 and ϕ2=−0.5, respectively. The values of Q0, Q1 and Q2 for each sub-figure in Figure 3 are
(1) Figure 3a:
Q0=(3,0),Q1=(3.7619,−1.52381),Q2=(3.52205,−4.59965). |
(2) Figure 3b:
Q0=(2.35078,0.884407),Q1=(3.69117,−0.306394),Q2=(5.77607,−6.54795). |
(3) Figure 3c:
Q0=(1.92252,1.16194),Q1=(3.33641,0.568314),Q2=(6.50264,−4.93599). |
(4) Figure 3d:
Q0=(1.5,1.25),Q1=(2.83333,1.25),Q2=(6.56173,−2.30556). |
From Figure 4, it is obvious that w2 did not affect the control points of C2. This means that the common point or joined point between C1 and C2 is independent of w2.
Example 3.3. Figure 5 shows three consecutive cubic generalized fractional Bézier curves with F2 continuity. Blue, orange and green curves are C1, C2 and C3, respectively. The control points and shape parameters are the same as in Example 3.2. Q0, Q1, Q2, R0, R1 and R2 can be obtained by Theorem 3.2. Other control point, shape parameters and scale factors are R3=(0,−2), (c1,c2,c3)=(−1.5,0,1.5), χ1=1.25 and χ2=−0, respectively. The values of Q0, Q1, Q2, R0, R1 and R2 for each sub-figure in Figure 5 are
(1) Figure 5a:
Q0=(3,0),Q1=(3.7619,−1.52381),Q2=(3.52205,−4.59965),R0=(3,−5),R1=(2.02551,−5.74733),R2=(2.1651,−3.78947). |
(2) Figure 5b:
Q0=(2.35078,0.884407),Q1=(3.69117,−0.306394),Q2=(5.77607,−6.54795),R0=(3,−5),R1=(−2.182,−2.11049),R2=(−3.59036,1.92539). |
(3) Figure 5c:
Q0=(3,0),Q1=(3.7619,−1.52381),Q2=(3.52205,−4.59965),R0=(3.19263,−4.74017),R1=(2.33428,−6.48591),R2=(2.32423,−4.16135). |
(4) Figure 5d:
Q0=(2.35078,0.884407),Q1=(3.69117,−0.306394),Q2=(5.77607,−6.54795),R0=(3.89875,−5.30721),R1=(0.525349,−5.29025),R2=(−1.91411,−0.038501). |
(5) Figure 5e:
Q0=(2.35078,0.884407),Q1=(3.69117,−0.306394),Q2=(5.77607,−6.54795),R0=(4.4075,−4.99986),R1=(3.06888,−7.704),R2=(0.413905,−2.68839). |
(6) Figure 5f:
Q0=(1.5,1.25),Q1=(2.83333,1.25),Q2=(6.56173,−2.30556),R0=(4.64757,−2.76389),R1=(3.61053,−7.91349),R2=(−0.442755,−6.46024). |
Figure 6 illustrates the curve adjustability for the third curve via w3. Based on Figure 6, clearly w3 did not affect the control points of C3. This means that the common point or joined point between C2 and C3 is independent of w3.
In Figure 5b, when the value of the fractional parameter w1 varies, it will cause the control points Q0,Q1,Q2,R1, and R2 to change their positions to maintain F2 continuity between the curves C1 and C2 without affecting the convex hull of the curve, C1, and the control points Q3 and R0.
In Figure 5c, when the fractional parameter w2 differs, the coordinates of the control points R0,R1, and R2 are moved to maintain the F2 continuity between the curves C2 and C3 without changing the convex hull of the curves C1 and C2 and the control point R3.
Figure 5d–5f illustrates that the convex hulls of the curves C2 and C3 are shifted when the values of the fractional parameters w1 and w2 are altered. Figure 6 depicts that changing the value of the fractional parameter w3 does not affect the convex hull of all curves, but the fractional parameter w3 can be used to adjust the curve adjustability of the curve C3.
Definition 4.1 (Generalized fractional Bézier tensor product surface). Suppose Pi,j, where i=0,1,…,m and j=0,1,…,n are the control points. The generalized fractional Bézier tensor product surface is defined as follows:
S(u,v;w1,w2,a1,a2,…,am,b1,b2,…,bn)=m∑i=0n∑j=0ˉfi,m(u;w1,a1,a2,…,am)ˉfi,n(v;w2,b1,b2,…,bn)Pi,j,(u,v)∈[0,1]×[0,1], | (4.1) |
where ai and bj are the shape parameters for i=1,…,m and j=1,…,n, while w1 and w2 are the fractional parameters. Moreover, ˉfi,m(u;w1,a1,a2,…,am) and ˉfi,n(v;w2,b1,b2,…,bn) are the generalized fractional Bézier basis functions. Also, note that a0=am+1=b0=bn+1=0. Hence, the generalized fractional Bézier tensor product surface has m+n shape parameters with an additional two fractional parameters.
Remark 1. The tensor product generalized fractional Bézier surfaces (see Eq (4.1)) inherited almost all the properties of the classical tensor product Bézier surfaces, such as geometric and affine invariance, convex hull, endpoint interpolation (when the fractional parameters are zero), degeneracy and symmetry properties.
Remark 2. The generalized fractional Bézier tensor product surface has two additional properties: local shape flexibility and surface patch adjustability. Shape parameters give flexibility in changing the shape of the surface locally without changing the control points. The fractional parameters enable the adjustability of the surface patch. By varying the values of the fractional parameters, the fraction of the surface can be generated.
A classical Bézier tensor product surface has a lack of flexibility in changing shape. This is because the surface can be changed only by altering its control points. Moving the control points can cause the whole surface to change. Some of the modeling surfaces preferred only the local part of the surface to change. Hence, the implementation of shape parameters enables the shape to change locally or globally by simply changing the spectrum value of the shape parameters.
The generalized fractional Bézier tensor product surface has excellent local dynamic shape parameters. This shape parameter can be used to change the shape of the surface. This section will demonstrate the geometric effect of the shape parameters on the surface.
Proposition 4.1. On the premise that the generalized fractional Bézier surface's control points and fractional parameters remain unchanged, the following holds:
(1) Based on Eq (4.1), the number of shape parameters is m+n. The corner points, P0,0, Pm,0, P0,n and Pm,n will be associated with two shape parameters. Control points Pk,0, P0,l, Pk,n, and Pm,l where k=1,2,…,m−1 and l=1,2,…,n−1 corresponded to three shape parameters. The rest of the control points will have four shape parameters. For a clear visual, Figure 7 shows each control point in the control mesh corresponding to their respective shape parameters.
(2) Increasing the value of ai causes the surface to move closer to the control point Pi,j but further from Pi−1,j. Conversely, if the value of ai decreases, the surface will move further from the control point Pi,j. For bj, increasing the value will cause the surface to move closer to the control point Pi,j but further from Pi,j−1. Meanwhile, decreasing the value of bi causes the surface to move further from the control point Pi,j.
(3) Hence, to generate the surface that is as close as possible to the control point Pi,j, the shape parameters need to be closer to the following values:
(a) sup(ai),sup(ai−1),…,sup(a1).
(b) sup(bj),sup(bj−1),…,sup(b1).
(c) inf(ai+1),inf(ai+2),…,inf(am).
(d) inf(bj+1),min(bj+2),…,inf(bn).
Conversely, to generate the shape of a surface that is as far as possible from the control point of Pi,j, the shape parameters need to be closer to the following values:
(a) inf(ai),inf(ai−1),…,inf(a1).
(b) inf(bj),inf(bj−1),…,inf(b1).
(c) sup(ai+1),sup(ai+2),…,sup(am).
(d) sup(bj+1),sup(bj+2),…,sup(bn).
Remark 3. The terms sup and inf used in Proposition 4.1 are based on Definition 2.2. Since the range of the shape parameters is open range, the value of shape parameters can only be close to the extreme values. Therefore, sup and inf are used to indicate that the shape parameter values can only approach the extreme values.
Example 4.1. Figure 8 indicates the biquartic fractional Bézier surface with zero shape parameters. This example shows how to control the shape of the surface to get as close as possible to the control point P2,2. Here, the shape parameters are (a1,a2,a3,a4,b1,b2,b3,b4).
Example 4.2. The shape-changing process of the biquartic fractional Bézier surface to get as far as possible from point P2,2 is illustrated in Figure 8. The control points are the same as in Example 4.1. At the same time, the fractional parameters are zero.
From Figures 8 and 9, designers can manipulate these shape parameters to adjust the shape of the surface to their liking. To change the shape locally, the designer can manipulate one of the shape parameters. To change the shape globally, the designers need to change ai (i=1,…,m) or change bj (j=1,…,n).
The generalized fractional Bézier surface has another kind of parameter called the fractional parameter. This fractional parameter gives a surface patch adjustability property to the surface. The surface patch adjustability property means the surface can adjust its length in u and v directions according to the respective fractional parameter.
Any degree of the surfaces will have two fractional parameters, w1 and w2. Here, w1 is used to adjust the surface patch in u direction while w2 controls the length of surface generated in the v direction. By observing the geometric behavior of each fractional parameter, it is possible to derive the fractional continuity either in u, v or both directions (will be discussed further in Section 5).
Example 4.3. Figures 10 and 11 illustrate the surface patch adjustability using the fractional parameters. The control points and the value of the shape parameters for Figures 10 and 11 are the same as in Figures 8d and 9d, respectively.
Remark 4. There will be several figures of surfaces that have been created with reduced opacity color from this part till the end of the section. These surfaces with a lower opacity color depict the constructed surfaces when fractional parameters are set to zero. This type of modeling makes a comparison between surfaces with zero fractional value parameters and surfaces with non-zero fractional value parameters easier.
Figures 10 and 11 depict that w1 is used to control the surface length in the u direction while w2 is for controlling the surface length in the v direction.
Previously, if the designers wanted to connect the second surface to the first surface's mid-line, the first surface must be spliced using the subdivision approach and then combined to the second surface. This method is time-consuming and expensive to compute, especially with high surface degrees. A new sort of continuity needs to be developed to address these restrictions.
Fractional continuity is an upgraded version of geometric continuity, where its definitions and theorems for curves up to degree two are discussed meticulously in [22]. By referring to Figures 10–12 in [23], the behavior of the fractional continuity can be seen clearly with the help of visualization of the curvature plot and comb. The same paper also shows that the degree of continuity is still maintained regardless of the position of the common point and the values of fractional parameters. Any part of the first curve can be connected to the second curve as the method of cut and combine is introduced in [24]. Using the cut and combine technique in [24], it is possible to construct outline of the real-life objects and animations with low curvature profile.
Although parametric and geometric continuities only allow two surfaces to be joined at the boundary ends, fractional continuity, on the other hand, allows two surfaces to be connected at any arbitrary line along the first surface in either u or v directions. The fractional parameter plays an essential role in making the fractional continuity feasible. This is due to the surface patch adjustability from the fractional parameter that makes the surfaces can be joined along any line of the first surface either in u or v direction. Hence, the fractional continuity is a good substitution for the subdivision method since the same result can be obtained by simply varying the fractional parameter.
Remark 5. Note that it is recommended for the readers to go through the definitions, theorems and analysis of the fractional continuity in [22,23,24].
Before deriving the F2 continuity conditions, two adjacent generalized fractional Bézier surfaces must be defined first to facilitate the proving and derivation of the following theorems. Suppose there are two adjacent generalized fractional Bézier surfaces as follows:
{S1(u,v;w1,1,w1,2,a1,1,a1,2,…,a1,m1,b1,1,b1,2,…,b1,n1)=∑m1i=0∑n1j=0ˉfi,m1(u;w1,1,a1,1,a1,2,…,a1,m1)ˉfi,n1(v;w1,2,b1,1,b1,2,…,b1,n1)Pi,j,S2(u,v;w2,1,w2,2,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2)=∑m2i=0∑n2j=0ˉfi,m2(u;w2,1,a2,1,a2,2,…,a2,m2)ˉfi,n2(v;w2,2,b2,1,b2,2,…,b2,n2)Qi,j. | (5.1) |
The first and second surfaces will be denoted as S1(u,v;w1,1,w1,2) and S2(u,v;w2,1,w2,2), respectively, to make the notation simpler for derivation and proving since the original notation is too long.
Definition 5.1 (Fractional continuity for generalized fractional Bézier surface, Fr). Consider two generalized fractional Bézier surfaces S1(u,v;w1,1,w1,2) and S2(u,v;w2,1,w2,2) on (u,v)∈[0,1]×[0,1]. The two surfaces are Fr continuous in the u direction if the following condition is satisfied:
S1(u,1;0,w1,2)=S2(u,0;0,0),∂∂vS1(u,1;0,w1,2)=ϕ∂∂vS2(u,0;0,0),∂2∂v2S1(u,1;0,w1,2)=ϕ2∂2∂v2S2(u,0;0,0),⋮∂r∂vrS1(u,1;0,w1,2)=ϕr∂r∂vrS2(u,0;0,0). |
For the v direction, the following condition must be satisfied:
S1(1,v;w1,1,0)=S2(0,v;0,0),∂∂uS1(1,v;w1,1,0)=ϕ∂∂uS2(0,v;0,0),∂2∂u2S1(1,v;w1,1,0)=ϕ2∂2∂u2S2(0,v;0,0),⋮∂r∂urS1(1,v;w1,1,0)=ϕr∂r∂urS2(0,v;0,0). |
For u and v directions, the following condition must be satisfied:
S1(u,1;0,w1,2)=S2(0,v;0,0),∂∂vS1(u,1;0,w1,2)=ϕ∂∂uS2(0,v;0,0),∂2∂v2S1(u,1;0,w1,2)=ϕ2∂2∂u2S2(0,v;0,0),⋮∂r∂vrS1(u,0;0,w1,2)=ϕr∂r∂urS2(0,v;0,0). |
Generally, for Fr continuity, ϕ>0 and ϕ∈R.
Theorem 5.1. If the two adjacent generalized fractional Bézier surfaces S1(u,v;w1,1,w1,2) and S2(u,v;w2,1,w2,2) satisfy all the following conditions:
{w1,1=w2,1,m1=m2,a1,k=a2,l,(k=0,1,…,m1;l=0,1,…,m2),Fori=0,1,…,m1,Qi,0=∑n1jˉfj,n1(1;w1,2,b1,1,b1,2,…,b1,n1)Pi,j,Qi,1=1ϕ(n2+b2,1)(ddv(∑n1jˉfj,n1(v;w1,2,b1,1,b1,2,…,b1,n1)Pi,j)v=1)+Qi,0,Qi,2=1n2ϕ2(b2,2+n2−1)(d2dv2(∑n1jˉfj,n1(v;w1,2,b1,1,b1,2,…,b1,n1)Pi,j)v=1+(n2((1−n2)ϕ2))+(−2n2ϕ2)b2,1)Qi,0+((2n2ϕ2)b2,1+n2((−1+n2)2ϕ2+ϕ2b2,2))Qi,1), | (5.2) |
then F2 continuity will be achieved for the u direction with ϕ>0.
Proof. To obtain F2 continuity, F1 continuity at the joints must be accomplished first where for every point on the same boundary, the two surfaces must have a shared tangent plane. Similarly, to obtain F1 continuity, the F0 continuity condition must be satisfy first by having the same common boundary of the two surfaces. Note that w1,1, w2,1 and w2,2 are set to zero to simplify the equation. This is due to the fractional surface adjustability property since w1,2 is the only fractional parameter responsible for changing the boundary point for S1(u,1) in the u direction. Hence, it can be written in the equation form as
S1(u,1;0,w1,2,a1,1,a1,2,…,a1,m1,b1,1,b1,2,…,b1,n1)=S2(u,0;0,0,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2). |
Using the boundary properties of the surfaces in Theorem 2.1, the previous equation can be written as
m1∑iˉfi,m1(u;w1,2,b1,1,b1,2,…,b1,n1)Pi,n1=m2∑iˉfi,m2(u;0,b2,1,b2,2,…,b2,n1)Qi,0. | (5.3) |
Since the basis functions in (5.3) are linearly independent, it can be simplified by comparing the coefficients as follows:
{w1,1=w2,1m1=m2,a1,k=a2,l(k=0,1,…,m1;l=0,1,…,m2),For i=0,1,…,m1,Qi,0=∑n1jˉfj,n1(1;w1,2,b1,1,b1,2,…,b1,n1)Pi,j. | (5.4) |
By F1 continuity definition, the two surfaces must have a common tangent at any point at their common boundary. Thus, the following equation needs to be satisfied:
∂∂vS1(u,1;0,w1,2,a1,1,a1,2,…,a1,m1,b1,1,b1,2,…,b1,n1)×∂∂uS1(u,1;0,w1,2,a1,1,a1,2,…,a1,m1,b1,1,b1,2,…,b1,n1)=ϕ(u)∂∂vS2(u,0;0,0,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2)×∂∂uS2(u,0;0,0,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2), | (5.5) |
where ϕ(u) is the scaling factor between the normal vectors at the joint and ϕ(u)>0. To simplify the calculation, the Faux method will be used, in which Eq (5.5) can be simplified as follows:
∂∂vS1(u,1;0,w1,2,a1,1,a1,2,…,a1,m1,b1,1,b1,2,…,b1,n1)=ϕ∂∂vS2(u,0;0,0,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2), | (5.6) |
where ϕ is a positive real constant. Here, Eq (5.6) shows that the tangent vector of the cross-border at their common boundary should be continuous.
Using the boundary tangent property in the v direction and substituting in the previous equation, it can be expressed as follows:
ddv(m1∑iˉfi,m1(v;w1,2,b1,1,b1,2,…,b1,n1)Pi,n1)v=1=ϕddv(m2∑iˉfi,m2(v;0,b2,1,b2,2,…,b2,n1)Qi,0)v=1. | (5.7) |
Since the fractional parameters w1,1, w2,1 and w2,2 are independent in Eq (5.7), it can be further simplified to
ddv(m1∑iˉfi,m1(1;w1,2,b1,1,b1,2,…,b1,n1)Pi,n1)v=1=ϕm2∑i(Qi,1−Qi,0)(n2+b2,1). | (5.8) |
Combining the results from Eqs (5.4) and (5.8) can be simplified by taking Qi,1 as a subject, which yields
Qi,1=1ϕ(n2+b2,1)(ddv(n1∑jˉfj,n1(v;w1,2,b1,1,b1,2,…,b1,n1)Pi,j)v=1)+Qi,0,(i=0,1,…,m1). | (5.9) |
In pursuance of F2 smooth continuity, the two surfaces need to hold the same normal curvature at any point on the common boundary; hence, they need to satisfy
∂2∂v2S1(u,1,0,w1,2,a1,1,a1,2,…,a1,m1,b1,1,b1,2,…,b1,n1),=ϕ2∂2∂v2S2(u,0;0,0,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2),+2ϕχ(u)∂2∂u∂vS2(u,0;0,0,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2),+χ2(u)∂2∂u∂vS2(u,0;0,0,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2),+ψ∂∂vS2(u,0;0,0,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2),+ω(u)∂∂uS2(u,0;0,0,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2), | (5.10) |
where χ(u) and ω(u) are linear functions of u, while ψ and ϕ are the arbitrary constants. In practical applications, χ(u)=ω(u)=ψ=0 are set to simplify the previous equation, which can be written as follows:
∂2∂v2S1(u,1;0,w1,2,a1,1,a1,2,…,a1,m1,b1,1,b1,2,…,b1,n1)=ϕ2∂2∂v2S2(u,0;0,0,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2). | (5.11) |
The second-order derivatives of the generalized fractional Bézier basis functions (i=0,1,…,n;n≤2) at terminal point u=0 are given by
¯f″i,n(0;w,a1,a2,…,an)={n(2a1+n−1),i=0,−n(2a1+a2+2n−2),i=1,−n(a2+n−1),i=2,0,i=3,4,…,n, | (5.12) |
¯f″i,n(1;w,a1,a2,…,an)=d2du2(ˉfi,n(u;w,a1,a2,…,an))u=1. | (5.13) |
Equation (5.12) shows that the fractional parameter, w, is independent for any derivative at terminal points u=0. Thus, w=0 is set to simplify the calculation. However, for terminal point u=1, the second derivative becomes the linear combination of function u [22].
Thus, from Eqs (5.12) and (5.13), we have
{∂2∂v2S1(u,1;0,w1,2,a1,1,a1,2,…,a1,m1,b1,1,b1,2,…,b1,n1),=d2dv2(n1∑jˉfj,n1(u;w1,2,b1,1,b1,2,…,b1,n1)Pi,j),∂2∂v2S2(u,0;0,0,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2)=n2∑i(n2(2b2,1+n2−1)Qi,0,−n2(2b2,1+b2,2+2n2−2)Qi,1−n2(b2,2+n2−1)Qi,2). | (5.14) |
Substituting Eq (5.14) in Eq (5.11) yields
d2dv2(n1∑jˉfj,n1(v;w1,2,b1,1,b1,2,…,b1,n1)Pi,j)v=1=ϕ2n2∑i(n2(2b2,1+n2−1)Qi,0−n2(2b2,1+b2,2+2n2−2)Qi,1−n2(b2,2+n2−1)Qi,2). | (5.15) |
Finally, combining the F0 and F1 continuity conditions of Eqs (5.4) and (5.9) with Eq (5.15), F2 continuity conditions can be obtained as follows:
Qi,2=1n2ϕ2(b2,2+n2−1)(d2dv2(n1∑jˉfj,n1(v;w1,2,b1,1,b1,2,…,b1,n1)Pi,j)v=1+(n2((1−n2)ϕ2))+(−2n2ϕ2)b2,1)Qi,0+((2n2ϕ2)b2,1+n2((−1+n2)2ϕ2+ϕ2b2,2))Qi,1). | (5.16) |
To conclude, if the S1(u,v;w1,1,w1,2) and S2(u,v;w2,1,w2,2) satisfy Eqs (5.4), (5.9) and (5.15), the surface will be connected with F2 continuity in the u direction at the joint. Hence, Theorem 5.1 is proven. To achieve F1 continuity in the u direction at the joint, the two surfaces need to satisfy only Eqs (5.4) and (5.9).
Theorem 5.2. If the two adjacent generalized fractional Bézier surfaces S1(u,v;w1,1,w1,2) and S2(u,v;w2,1,w2,2) satisfy all the following conditions:
{w1,1=w2,2,m1=n2,a1,k=b2,l,(k=0,1,…,m1;l=0,1,…,n2),Fori=j=0,1,…,m1,andϕ>0,Q0,j=∑n1jˉfj,n1(1;w1,2,b1,1,b1,2,…,b1,n1)Pi,j,Q1,j=1ϕ(m2+a2,1)(ddv(∑n1jˉfj,n1(v;w1,2,b1,1,b1,2,…,b1,n1)Pi,j)v=1)+Qi,0,Q2,j=1m2ϕ2(a2,2+m2−1)(d2dv2(∑n1jˉfj,n1(v;w1,2,b1,1,b1,2,…,b1,n1)Pi,j)v=1+(m2((1−m2)ϕ2))+(−2m2ϕ2)a2,1)Q0,j+((2m2ϕ2)a2,1+m2((−1+m2)2ϕ2+ϕ2a2,2))Q1,j), | (5.17) |
then F2 continuity will be achieved for u and v directions.
Proof. Suppose that S1(u,v;w1,1,w1,2) in the u direction wants to be connected at F2 continuity with S2(u,v;w2,1,w2,2) in the v direction. The derivation is similar to Theorem 5.1, where the two surfaces must achieve F1 continuity conditions in the u and v directions. Hence,
{w1,1=w2,2,m1=n2,a1,k=b2,l,(k=0,1,…,m1;l=0,1,…,n2),For i=j=0,1,…,m1, and ϕ>0,Q0,j=∑n1jˉfj,n1(1;w1,2,b1,1,b1,2,…,b1,n1)Pi,j,Q1,j=1ϕ(m2+a2,1)(ddv(∑n1jˉfj,n1(v;w1,2,b1,1,b1,2,…,b1,n1)Pi,j)v=1)+Qi,0. | (5.18) |
When F1 continuity conditions have been fulfilled, the two surfaces need to get the same normal curvature at any point on their common boundary. This can be written in the equation as follows:
∂2∂v2S1(u,1;0,w1,2,a1,1,a1,2,…,a1,m1,b1,1,b1,2,…,b1,n1)=ϕ2∂2∂u2S2(0,v;0,0,a2,1,a2,2,…,a2,m2,b2,1,b2,2,…,b2,n2), | (5.19) |
where ϕ is similar in Eq (5.18).
Using the same method in Theorem 5.1, i.e., using the second-order cross-border tangent vector in Eq (5.19), the following equation can be obtained:
d2dv2(n1∑jˉfj,n1(v;w1,2,b1,1,b1,2,…,b1,n1)Pi,j)v=1=ϕ2m2∑i(m2(2a2,1+m2−1)Q0,j−m2(2a2,1+a2,2+2m2−2)Q1,j−m2(a2,2+m2−1)Q2,j). | (5.20) |
Since a1,k=b2,l for k=0,1,…,m1 and l=0,1,…,n2, it can be simplified further as follows:
Q2,j=1m2ϕ2(a2,2+m2−1)(d2dv2(n1∑jˉfj,n1(v;w1,2,b1,1,b1,2,…,b1,n1)Pi,j)v=1+(m2((1−m2)ϕ2))+(−2m2ϕ2)a2,1)Q0,j+((2m2ϕ2)a2,1+m2((−1+m2)2ϕ2+ϕ2a2,2))Q1,j). | (5.21) |
To sum up, when the surfaces S1(u,v;w1,1,w1,2) and S2(u,v;w2,1,w2,2) satisfy the conditions in Eqs (5.19) and (5.21), then the two surfaces achieve the F2 continuity in the direction of u and v at the joint. Thus, Theorem 5.2 is proven.
Theorem 5.3. If the two adjacent generalized fractional Bézier surfaces S1(u,v;w1,1,w1,2) and S2(u,v;w2,1,w2,2) satisfy all the following conditions:
{w1,2=w2,2,n1=n2,b1,k=b2,l,(k=0,1,…,n1;l=0,1,…,n2),Forj=0,1,…,n1,Q0,j=∑m1iˉfi,m1(1;w1,1,a1,1,a1,2,…,a1,n1)Pi,j,Q1,j=1ϕ(m2+a2,1)(ddu(∑m1iˉfi,m1(u;w1,1,a1,1,a1,2,…,a1,m1)Pi,j)v=1)+Q0,j,Q2,j=1m2ϕ2(a2,2+m2−1)(d2du2(∑m1iˉfi,m1(u;w1,1,a1,1,a1,2,…,a1,m1)Pi,j)u=1+(m2((1−m2)ϕ2))+(−2m2ϕ2)a2,1)Q0,j+((2m2ϕ2)a2,1+m2((−1+m2)2ϕ2+ϕ2a2,2))Q1,j), | (5.22) |
then F2 continuity will be achieved for v direction with ϕ>0.
Proof. The proof is similar to Theorem 5.1, implying that it will not be covered. Note that the fractional parameter w1,1 will be the only fractional parameter contributing to F2 continuity in the v direction. Hence, the other fractional parameters can be set to zero to simplify the calculations in the proof.
The F2 continuity can be reverted to G2 continuity by simply setting all the fractional parameters to zero.
Fractional continuity for generalized fractional Bézier surfaces is a key in constructing numerous complex surfaces. In addition, shape flexibility and length adjustability are the main features of versatile surface designing. This section will show procedures for F2 continuity of two generalized fractional Bézier surfaces. On account of Theorem 5.1, the procedures are as follows:
(1) The following requirement must be fulfilled first: the order m1, n1 of S1(u,v;w1,1,w1,2) and its control points, Pi,j as well as shape parameters a1,i and b1,j (i=0,1,…,m1;j=0,1,…,n1) and fractional parameters w1,1 and w1,2 are set to any arbitrary values within their respective ranges.
(2) Suppose w1,1=w2,1, m1=m2, a1,k=a2,l (k=0,1,…,m1;l=0,1,…,m2) and Qi,0 (i=0,1,…,m1); this enables S1(u,v;w1,1,w1,2) and S2(u,v;w2,1,w2,2) to achieve F0 continuity.
(3) Set the values of shape parameters b2,j (j=0,1,…,n2), the constant ϕ>0, the fractional parameter w2,2 and the order of n2 for the S2(u,v;w2,1,w2,2). Then, calculate for Qi,1 (i=0,1,…,m1) of S2(u,v;w2,1,w2,2).
(4) After calculating Qi,1, then proceed for Qi,2 (i=0,1,…,m1) calculation.
(5) The remaining n2−2 control points of S2(u,v;w2,1,w2,2) are up to the designers' discretion. Hence, F2 continuity between two generalized fractional Bézier surfaces in the u direction is achieved.
The procedures above can be repeated to connect multiples of generalized fractional Bézier surface with F2 continuity.
To demonstrate the F2 continuity, some examples will be shown. In this section, F0 and F1 continuity will be shown first to compare different values of fractional parameters. Then, for F2 continuity, some examples are given with different values of fractional parameters, shape parameters and scale factors.
Remark 6. The remaining control points for Examples 6.1 – 6.6 are up to the designer's discretion to choose, due to the fact that the remaining control points do not involve as a constraint in continuity conditions.
The fractional parameter plays an essential role in fractional continuity. By varying the fractional parameter, the designer can control the boundary line between two surfaces while maintaining continuity.
Example 6.1. Figure 12 shows the F0 continuity in the u direction between two biquartic fractional Bézier surfaces with different values of fractional parameters. w1,2 is for controlling the connection line between the two surfaces in the u direction. Meanwhile, w2,2 is for controlling the length of the second surface in the u direction. Moreover, w1,1=w2,1 are for controlling the length of both surfaces in v direction.
Example 6.2. Figures 13 and 14 depict the F1 continuity and F2 continuity in the u direction between two biquartic fractional Bézier surfaces with different values of fractional parameters, respectively. The scale factor, ϕ=1.25, is used in this example, and the control points are the same as in Example 6.1.
The shape parameter gives an advantage to the aesthetic curve in terms of flexibility. This is due to the shape parameter enabling the shape of the surface to change without changing the control points. The shape parameter also allows the shape of the surface to change locally. To change shape locally, the designers need to change one of the shape parameters. To change shape globally, the designers need to change ai (i=0,1,…,m) or change bj (j=0,1,…,n).
Proposition 4.1 gives some insight into controlling the shape using shape parameters. The control points and convex hull will guide how to control the surface shape at a specific point.
Example 6.3. The F2 continuity between two biquartic fractional Bézier surfaces in u direction with the variation of shape parameters is depicted in Figure 15. Here, w1,1=w2,1=w2,2=0, w1,2=0.75 and ϕ=1.25 are chosen. The shape parameters used for each sub-figure are
a. (a1,1,a1,2,a1,3,a1,4,b1,1,b1,2,b1,3,b1,4,b1,1,b1,2,b1,3,b1,4)=(0,0,0,0,0,0,0,0,0,0,0,0).
b.
(a1,1,a1,2,a1,3,a1,4,b1,1,b1,2,b1,3,b1,4,b1,1,b1,2,b1,3,b1,4)=(−2,−1,1,2,0.5,1.5,−1.5,−0.5,0,−2,2,0). |
c.
(a1,1,a1,2,a1,3,a1,4,b1,1,b1,2,b1,3,b1,4,b1,1,b1,2,b1,3,b1,4)=(0.5,0.5,−0.5,−0.5,−1,−1.5,1.5,0.5,−1,−1,0.5,1). |
d.
(a1,1,a1,2,a1,3,a1,4,b1,1,b1,2,b1,3,b1,4,b1,1,b1,2,b1,3,b1,4)=(−1,−2,2,1,0.5,1.5,−1.5,−0.5,0.5,1.5,−1.5,0.5). |
Example 6.4. Figure 16 illustrates the F2 continuity in the u direction between two biquartic fractional Bézier surfaces with different scale factor values, ϕ. The control points, shape and fractional parameters are similar to Figure 15a.
The F2 continuity conditions can also be applied to generate three consecutive surfaces either in u or v directions. The third generalized fractional Bézier tensor surface will be defined as follows:
S3(u,v;w3,1,w3,2,a3,1,a3,2,…,a3,m3,b3,1,b3,2,…,b3,n3)=m3∑i=0n3∑j=0ˉfi,m3(u;w3,1,a3,1,a3,2,…,a3,m3)ˉfi,n3(v;w3,2,b3,2,…,b3,n3)Ri,j. | (6.1) |
Example 6.5. Figure 17 depicts the Möbius strip generated by three consecutive surfaces with F2 continuity in u direction with different fractional parameters. The green surface is the third surface defined by Eq (6.1). The F2 continuity is between magenta and orange surfaces with orange and green surfaces. Green and magenta surfaces are connected by simple C0/G0 continuity; hence the value of w3,2=0 is set and cannot be changed, or else they will be disconnected. w1,1=w2,1=w3,1=w3,2=0 are set.
In this section, the application of F2 continuity will be shown in modeling of one of the engineering surfaces known as the ruled surface. Assume four generalized fractional Bézier space curves given as follows:
{x1,1(u;w1,1,a1,1,a1,2,...,a1,n1)=∑n1j=0P1,jˉfj,n1(u),for u∈[0,1],x1,2(u;w1,2,b1,1,b1,2,...,b1,m1)=∑m1j=0P2,jˉfj,m1(u),for u∈[0,1],x2,1(u;w2,1,a2,1,a2,2,...,a2,n2)=∑n2j=0Q1,jˉfj,n2(u),for u∈[0,1],x2,2(u;w2,2,b2,1,b2,2,...,b2,m2)=∑m2j=0Q2,jˉfj,m2(u),for u∈[0,1]. | (6.2) |
By using Eq (6.2), two ruled surfaces will be defined as follows:
Definition 6.1 (Two generalized fractional Bézier ruled surface). Suppose x1,1(u;w1,1,a1,1,a1,2,...,a1,m1), x1,2(u;w1,2,b1,1,b1,2,...,b1,n1), x2,1(u;w2,1,a2,1,a2,2,...,a2,m2) and x2,2(u;w2,2,b2,1,b2,2,...,b2,n2) with m1, n1, m2 and n2 degree are the generalized fractional Bézier space curves, respectively. Two ruled surfaces are defined as follows:
Sr1(u,v;w1,1,w1,2,a1,1,a1,2,...,a1,m1,b1,1,b1,2,...,b1,n1)=(1−v)x1,1(u;w1,1,a1,1,a1,2,...,a1,m1)+vx1,2(u;w1,2,b1,1,b1,2,...,b1,n1), | (6.3) |
Sr2(u,v;w2,1,w2,2,a2,1,a2,2,...,a2,m2,b2,1,b2,2,...,b2,n2)=(1−v)x2,1(u;w2,1,a2,1,a2,2,...,a2,m2)+vx2,2(u;w2,2,b2,1,b2,2,...,b2,n2). | (6.4) |
Example 6.6. Figure 18 shows the F2 continuity between two bicubic fractional Bézier ruled surfaces with a variation of fractional parameters. The control points, shape parameters and scale factors are as follows:
P1,0=(1,1,1)P1,1=(2,1.5,−1)P1,2=(3,0,1)P1,3=(4,1.5,1)P2,0=(1.5,1,1)P2,1=(2.5,7,1.25)P2,2=(3.5,7,1,5)P2,3=(4.5,8,1)Q1,3=(8,2,−1)Q2,3=(8,8,−1.5)a1,1=−0.5a1,2=0a1,3=0.5b1,1=−1.5b1,2=0b1,3=1.5a2,1=−0.5a2,2=0a2,3=0.5b2,1=−2b2,2=0b2,3=2ϕ1,1=0.75ϕ1,2=1.25ϕ2,1=0.5ϕ2,2=−0.5 |
The control points Qi,j for i=1,2 and j=0,1,2 are obtained by F2 continuity conditions as in [23]. The red, blue, green and pink curves are x1,1(u), x1,2(u), x2,1(u) and x2,2(u), respectively. In order to connect two ruled surfaces with F2 continuity, x1,1(u) and x1,2(u) must be connected with F2 continuity with x2,1(u) and x2,2(u), respectively.
This paper discusses the F2 smooth continuity conditions for generalized fractional Bézier surfaces. The benefit of this work is that fractional continuity overcomes the conventional parametric and geometric continuity in the aspect of connecting surfaces at any arbitrary boundary along the first surface. Hence, fractional continuity simplifies two different processes (subdivision method and continuity conditions) into one process (continuity conditions) by simply varying the fractional parameters. The generalized fractional Bézier basis functions have complex formulation compared to the classical Bernstein basis functions, which becomes our limitation of study. Future research recommends expanding the basis functions in constructing developable surfaces and optimizing the shape and fractional parameters in modeling smart surface manufacturing for the industry.
This research is supported by Universiti Sains Malaysia under Short Term Grant (Khas) (304/PMATHS/6315587) and the School of Mathematical Sciences, Universiti Sains Malaysia. In addition, the authors are very grateful to the anonymous referees for their valuable suggestions.
The authors declare that they have no competing interests in this paper.
[1] | K. Miller, B. Ross, An introduction to the fractional calculus and fractional differential equations, Wiley, New York, 1993. |
[2] | G. E. Farin, G. Farin, Curves and surfaces for CAGD: A practical guide, Morgan Kaufmann, USA, 2002. |
[3] | D. M. Yip-Hoi, Teaching surface modeling to CAD/CAM technologists, 2011 ASEE Annual Conference Exposition, Vancouver, BC, 2011. |
[4] | F. Shi, Computer aided geometric design and non-uniform rational B-spline, Higher Education Press, Beijing, 2001. |
[5] | H. Prautzsch, W. Boehm, M. Paluszny, Bézier and B-spline techniques, Springer Berlin, Heidelberg, 2002. https://doi.org/10.1007/978-3-662-04919-8 |
[6] |
E. Mainar, J. M. Peña, J. Sánchez-Reyes, Shape preserving alternatives to the rational Bézier model, Comput. Aided Geom. D., 18 (2001), 37–60. https://doi.org/10.1016/S0167-8396(01)00011-5 doi: 10.1016/S0167-8396(01)00011-5
![]() |
[7] |
G. Hu, J. Wu, X. Qin, A novel extension of the Bézier model and its applications to surface modeling, Adv. Eng. Softw., 125 (2018), 27–54. https://doi.org/10.1016/j.advengsoft.2018.09.002 doi: 10.1016/j.advengsoft.2018.09.002
![]() |
[8] |
U. Bashir, M. Abbas, M. N. H. Awang, J. M. Ali, A class of quasi-quintic trigonometric Bézier curve with two shape parameters, Sci. Asia S, 39 (2013), 11–15. https://doi.org/10.2306/scienceasia1513-1874.2013.39S.011 doi: 10.2306/scienceasia1513-1874.2013.39S.011
![]() |
[9] |
M. Y. Misro, A. Ramli, J. M. Ali, Quintic trigonometric Bézier curve with two shape parameters, Sains Malays., 46 (2017), 825–831. https://doi.org/10.17576/jsm-2017-4605-17 doi: 10.17576/jsm-2017-4605-17
![]() |
[10] |
S. BiBi, M. Abbas, M. Y. Misro, G. Hu, A novel approach of hybrid trigonometric Bézier curve to the modeling of symmetric revolutionary curves and symmetric rotation surfaces, IEEE Access, 7 (2019), 165779–165792. https://doi.org/10.1109/ACCESS.2019.2953496 doi: 10.1109/ACCESS.2019.2953496
![]() |
[11] |
F. Li, G. Hu, M. Abbas, K. T. Miura, The generalized H-Bézier model: Geometric continuity conditions and applications to curve and surface modeling, Mathematics, 8 (2020), 924. https://doi.org/10.3390/math8060924 doi: 10.3390/math8060924
![]() |
[12] |
D. Liu, J. Hoschek, GC1 continuity conditions between adjacent rectangular and triangular Bézier surface patches, Comput. Aided D., 21 (1989), 194–200. https://doi.org/10.1016/0010-4485(89)90044-4 doi: 10.1016/0010-4485(89)90044-4
![]() |
[13] |
W. Du, F. J. Schmitt, On the G1 continuity of piecewise Bézier surfaces: A review with new results, Comput. Aided Design, 22 (1990), 556–573. https://doi.org/10.1016/0010-4485(90)90041-A doi: 10.1016/0010-4485(90)90041-A
![]() |
[14] |
G. Hu, H. Cao, X. Wang, X. Qin, G2 continuity conditions for generalized Bézier-like surfaces with multiple shape parameters, J. Inequal. Appl., 2017 (2017), 1–17. https://doi.org/10.1186/s13660-017-1524-7 doi: 10.1186/s13660-017-1524-7
![]() |
[15] |
N. H. M. Ismail, M. Y. Misro, Surface construction using continuous trigonometric Bézier curve, AIP Conf. Proc., 2266 (2020), 040012. https://doi.org/10.1063/5.0018101 doi: 10.1063/5.0018101
![]() |
[16] |
M. Ammad, M. Y. Misro, Construction of local shape adjustable surfaces using quintic trigonometric Bézier curve, Symmetry, 12 (2020), 1205. https://doi.org/10.3390/sym12081205 doi: 10.3390/sym12081205
![]() |
[17] |
G. Hu, C. Bo, G. Wei, X. Qin, Shape-adjustable generalized Bezier surfaces: Construction and it is geometric continuity conditions, Appl. Math. Comput., 378 (2020), 125215. https://doi.org/10.1016/j.amc.2020.125215 doi: 10.1016/j.amc.2020.125215
![]() |
[18] |
S. Bibi, M. Abbas, M. Y. Misro, A. Majeed, T. Nazir, Construction of generalized hybrid trigonometric Bézier surfaces with shape parameters and their applications, J. Math. Imaging Vis., 2021, 1–25. https://doi.org/10.1007/s10851-021-01046-y doi: 10.1007/s10851-021-01046-y
![]() |
[19] |
G. Hu, H. Cao, X. Qin, X. Wang, Geometric design and continuity conditions of developable λ-Bézier surfaces, Adv. Eng. Softw., 114 (2017), 235–245. https://doi.org/10.1016/j.advengsoft.2017.07.009 doi: 10.1016/j.advengsoft.2017.07.009
![]() |
[20] |
M. Ammad, M. Y. Misro, M. Abbas, A. Majeed, Generalized developable cubic trigonometric Bézier surfaces, Mathematics, 9 (2021), 283. https://doi.org/10.3390/math9030283 doi: 10.3390/math9030283
![]() |
[21] |
S. BiBi, M. Y. Misro, M. Abbas, A. Majeed, T. Nazir, G3 Shape adjustable GHT-Bézier developable surfaces and their applications, Mathematics, 9 (2021). https://doi.org/10.3390/math9192350 doi: 10.3390/math9192350
![]() |
[22] |
S. A. A. A. Said Mad Zain, M. Y. Misro, K. T. Miura, Generalized fractional Bézier curve with shape parameters, Mathematics, 9 (2021), 2141. https://doi.org/10.3390/math9172141 doi: 10.3390/math9172141
![]() |
[23] |
S. A. A. A. Said Mad Zain, M. Y. Misro, Shape analysis and fairness metric of generalized fractional Bézier curve, Comput. Appl. Math., 41 (2022), 1–24. https://doi.org/10.1007/s40314-022-01983-3 doi: 10.1007/s40314-022-01983-3
![]() |
[24] |
S. A. A. A. Said Mad Zain, M. Y. Misro, K. T. Miura, Curve fitting using generalized fractional Bézier curve, Comput. Aided Design Appl., 20 (2023), 350–363. https://doi.org/10.14733/cadaps.2023.350-363 doi: 10.14733/cadaps.2023.350-363
![]() |
1. | Moavia Ameer, Muhammad Abbas, Madiha Shafiq, Tahir Nazir, Asnake Birhanu, Sameer Sheshrao Gajghate, Generalized Bézier-like model and its applications to curve and surface modeling, 2024, 19, 1932-6203, e0303340, 10.1371/journal.pone.0303340 | |
2. | Syed Ahmad Aidil Adha Said Mad Zain, Md Yushalify Misro, Kenjiro T. Miura, Enhancing flexibility and control in κ-curve using fractional Bézier curves, 2024, 89, 11100168, 71, 10.1016/j.aej.2024.01.047 | |
3. | Mohamad Ekram Nordin, Md Yushalify Misro, Modified progressive iterative approximation techniques for interoperable conversion of QT-Bézier and rational Bézier curve, 2025, 16, 20904479, 103256, 10.1016/j.asej.2024.103256 |