Geometric modeling of some engineering GBT-Bézier surfaces with shape parameters and their applications

This study is based on some C1, C2, and C3 continuous computer-based surfaces that are modeled by using generalized blended trigonometric Bézier (shortly, GBT-Bézier) curves with shape parameters. Initially, generalized blended trigonometric Bernstein-like (shortly, GBTB) basis functions with two shape parameters are derived in explicit expression which satisfied the basic geometric features of the traditional Bernstein basis functions. Moreover, the GBT-Bézier curves and tensor product GBT-Bézier surfaces with two shape parameters are also presented. All geometric features of the proposed GBT-Bézier curves and surfaces are similar to the traditional Bézier curves and surfaces, but the shape-adjustment is the additional feature that the traditional Bézier curves and surfaces do not hold. Finally, a class of some complex computer-based engineering surfaces via GBT-Bézier curves with shape parameters is presented. In addition, two adjacent GBT-Bézier surfaces segments are connected by higher C2 and C3 continuity constraints than the existing only C1 shape adjustable Bézier surfaces. Some practical examples are provided to show the efficiency of the proposed scheme and to prove it as another powerful way for the construction and modeling of various complex composite computer-based engineering surfaces using higher-order continuities.


Introduction
Bézier curves and surfaces techniques are very popular in computer technology, computer graphics (CG), and computer-aided geometric design (CAGD) due to their shape control parameters. In computer-aided manufacturing (CAM) and computer-aided design (CAD), Bézier curves and surfaces are robust tools for constructing free form curves and surfaces. Bézier curves and surfaces have an abundance of appealing in the fields of engineering, science, communications, and technology particularly in networks, animation, railway route, highway design, environment design, CAD system, and many other disciplines. The GBTB functions are constructed with only two shape parameters. 2 Only simple SG-Bézier general cylinder has been constructed.
Two GBT-Bézier general cylinders can be connected by higher C 2 continuity smoothness. 3 Two SG-Bézier ruled surface patches have been connected by only C 1 continuity.
Two GBT-Bézier ruled surfaces can be connected by higher C k (k = 2, 3) continuity constraints. 4 Two adjacent SG-Bézier swung surface patches have been connected by only C 1 continuity.
Two adjacent GBT-Bézier swung surfaces can be connected by higher C 2 continuity conditions. 5 Composition of two SG-Bézier swept surfaces by only C 0 continuity conditions.
Two adjacent GBT-Bézier swept surfaces can be connected by higher C 1 and C 2 continuity smoothness. 6 No set of algorithms has been constructed. A set of algorithms can be constructed for explaining how to enforce these engineering surfaces in practice.
Anyhow, due to the fixed shape and position relative to the control polygon [1], the traditional Bézier curves and surfaces still have some deficiencies. Practical applications of Bézier curves and surfaces in the area of geometric modeling in engineering are restricted due to these deficiencies. With the help of shape control parameters into Bézier scheme, a remarkable study has been done [2][3][4][5][6][7] to tackle the problems in adjusting and controlling the shapes of Bézier curves and surfaces. Using a recursive approach, Yan and Liang [8] defined the Bézier curve and rectangular Bézier with single shape control parameters which are based on a new kind of polynomial basis functions. The presence of shape control parameter in [8] enhance the control over the shape and position of the proposed curves and surfaces. Based on the basis functions proposed by Yan and Liang, Hu et al. joined the two Bézier-like curve and surface segments with G 2 continuity constrains in [9]. Hu et al. [10] presented a new efficient technique for quickly designing generalized Bézier rotation surfaces using shape-adjustable generalized Bézier curves along various shape control parameters. Some advantages over existing scheme [11] are presented in Table 1.
However, the traditional Bézier curves and surfaces have another limitation due to their polynomial representation because the polynomial functions have low smoothness. Thus many scholars have tried to resolve this issue in a non-polynomial function space. Considerable research has been accomplished during the last few years with the help of trigonometric functions or the blending of polynomial and trigonometric functions for the description of curves and surfaces. Trigonometric B-splines were firstly presented by Schoenberg in [12], and the iterative relationship of random order trigonometric B-splines was settled in [13]. The trigonometric polynomials play a key role in the areas of CAGD, medicine, and electronics [14,15]. In recent times, trigonometric polynomials have also acquired remarkable attention in the area of geometric modeling in engineering. A new technique for the construction of quartic quasi-Bézier rotation surfaces was presented by Hu in [16] using various shape control parameters. Yan [17] expressed cubic trigonometric nonuniform spline curves and surfaces. Using appropriate shape control parameters, Sharma suggested quartic trigonometric, quasi-quartic trigonometric, and a class of Bézier-type cubic trigonometric curves and surfaces sequentially in [18,19], and [20].
To solve some problems in the construction of symmetric revolutionary curves and surfaces, a new approach was defined by BiBi et al. [21] using generalized hybrid trigonomet-ric Bézier curve (shortly, GHT-Bézier) with shape control parameters. Some free-form complex curves with restriction of parametric continuity were also constructed by using these GHT-Bézier curves. Furthermore, to illustrate the efficiency of the scheme, the authors also presented some composite symmetric curves and surfaces using symmetric formulas and parametric continuity conditions. In 2020, Maqsood et al. defined generalized trigonometric Bézier (shortly, GT-Bézier) curves and surfaces with some shape control parameters and also discussed their continuity constraints and some applications of geometric modeling in engineering [22]. A new approach about G 3 continuity of GHT-Bézier curves is proposed by BiBi et al. in [23]. Modeling of free form complex figures and sketching is also presented in this literature by C k and G k (k ≤ 3) continuity.
The GBT-Bézier curves based on GBTB functions with two shape parameters will be constructed in this research work. Furthermore, taking these GBT-Bézier curves, six types of engineering surfaces like GBT-Bézier cylinders, GBT-Bézier bilinear surfaces, GBT-Bézier ruled surfaces, GBT-Bézier swung surfaces, GBT-Bézier swept surfaces, and GBT-Bézier rotation surfaces will be formulated.
As an alternative technique of representing curves and surfaces, these proposed curves and surfaces not only have the valuable features of Bézier curves and surfaces but also allow an efficient shape modification technique. As an alternative technique of representing curves and surfaces in Bézier approach, these GBT-Bézier curves and surfaces not only have the valuable features of the traditional Bézier curves and surfaces but will also prove an efficient shape modification technique in the area of computer vision, manufacturing industry, computer animation, and multimedia technology.
Some technical contributions are made in this study which are as follows: • Construction of a tensor product GBT-Bézier surface by a new set of GBTB functions with two shape parameters. • Construction of some computer-based engineering surfaces using GBT-Bézier with shape parameters. • The complex computer-based engineering surfaces using GBT-Bézier patches are composed by C k (k = 1, 2, 3) continuity conditions. • Algorithms of all complex GBT-Bézier engineering surfaces are constructed for easy understanding of interested readers. This research work is organized as follows: In Sect. 2, basic definitions and properties of GBTB functions and their corresponding GBT-Bézier curves and surfaces are given. On the bases of proposed GBT-Bézier curves along with their shape control parameters, six different engineering surfaces are designed in Sect. 3. Concluding remarks on this research are given in Sect. 4.

Generalized blended trigonometric Bernstein-like basis functions
In this section, the generalized blended trigonometric Bernstein-like basis functions (GBTB) are defined.
The GBT-Bézier curves take over most advantageous features of the traditional Bézier curves, including symmetry, convex hull property, variation diminishing property, geometric invariance, shape adjustable property, and the terminal properties F ). The terminal properties of GBT-Bézier curves indicate that the GBT-Bézier curve segments interpolate to the both end points of their convex hull, and the values of their derivatives at both end points on GBT-Bézier curves will be easily modified by adjusting the shape parameters μ, ν in their respective value range, which brings a significant ease to smooth joining [24].
The tensor product GBT-Bézier surfaces shares all the properties with classical tensor product Bézier surfaces except for the shape adjustable property which is superior to the properties of classical tensor product Bézier surfaces.

Figures 2 and 3 represent a GBT-Bézier vase diagram and a bi-cubic GBT-Bézier surface
with different values of shape parameters.

GBT-Bézier cylinder with shape parameters
be a GBT-Bézier curve of degree m and V be a 3D unit vector, the representation of the GBT-Bézier cylinder T(z, z 1 , μ, ν) is obtained by sweeping F m (z, μ, ν) a distance α (α > 0) along V [25]. For sweep direction representing the parameter by z 1 , clearly the GBT-Bézier cylinder should meet the following requirements:  (8). Thus, in the perspective of the aforementioned conditions, the expression of the GBT-Bézier cylinder can be described as follows:

Proposition 3.1 The translational invariance property is used to define GBT-Bézier curves in
where α > 0 is a certain real constant. Taking F m (z, μ, ν) as directrix of the GBT-Bézier cylinder, expression (9) is called the GBT-Bézier cylinder.

Figure 3 Effect of shape parameters on bi-cubic GBT-Bézier surfaces
Moreover, the GBT-Bézier cylinder T cylinder is a surface which is generated by translating a small portion of the straight line from F m (0) to F m (0) + αV along an assigned curve F m (z, μ, ν). Anyhow, the cylinder in (9) leads to some complications in designing GBT-Bézier curves and surfaces as it is not a tensor product GBT-Bézier surface. To overcome this problem, equation (9) is transformed into the tensor product GBT-Bézier surface described in (7). It is observed that the mth-degree GBTB basis functions defined by (2) are only described for parameter m (m ≥ 2), but they are not given for m = 1. Thus, extending the definition of GBTB basis functions, we describe the first-degree GBTB basis functions Expression (10) shows that the GBTB basis functions of order one are the traditional Bernstein basis functions of order one. (9) can also be written in its tensor product form as follows:
Proof By taking control points Q k,0 = Q k and Q k,1 = Q k + αV = Q k,0 + αV and basis function given in (10), we have Therefore, GBT-Bézier cylinder in (9) can be exactly expressed in terms of a tensor product GBT-Bézier surface of degree (m, 1). Hence the proof of the result.
Since a GBT-Bézier curve is the directrix of a GBT-Bézier cylinder, it carries some remarkable benefits of the general Bézier cylinder. Furthermore, the GBT-Bézier cylinder (11) has an edge that by keeping the control mesh unchanged, the GBT-Bézier cylinder can be transformed into any shape merely by changing the shape parameters.
A closed GBT-Bézier cylinder is illustrated in Fig. 5 to show the impact of α and shape control parameters μ, ν on GBT-Bézier cylinder, whereas the C 2 continuity of two GBT-Bézier cylinders is illustrated in Fig. 6 with shape control parameters μ, μ 1 , ν, ν 1 . Figures 4, 5, and 6 present that the effect of shape parameters on GBT-Bézier cylinders (open or close) is the same as the effect on GBT-Bézier curves and surfaces.
However, the bi-linear surfaces described in (14) and (15) cannot be modified. This inconvenience of shape adjustment is due to the given four angular points. So, we need to generalize the classical bi-linear surface (14) into a GBT-Bézier bi-linear surface which is defined by taking a GBT-Bézier curve. Some properties of such surfaces are given as follows: 1. Angular points are interpolated by the surface; 2. Line segments are taken as the boundary curves of the surface; 3. T m,n (z, z 1 ) or T m,n (z, z 1 ) is a GBT-Bézier curve for fixed z or z 1 ; 4. The model of GBT-Bézier surface is used by the bi-linear surface. From the characteristics of GBT-Bézier curves and surfaces, the actual procedure for developing a GBT-Bézier bi-linear surface is described. The steps of the algorithm are as follows.
Modeling example of a third degree GBT-Bézier bi-linear surface is given below. By using control points and equation (15), the GBT-Bézier bi-linear surface can be designed as shown in Fig. 7 and Fig. 8. Bi-linear surfaces have sixteen control points, in which the exterior twelve control mesh points are coplanar. Furthermore, Fig. 8 exhibits the graphs of the GBT-Bézier bi-linear surfaces that have four non-coplanar boundary curves.

GBT-Bézier ruled surface along shape parameters
Let b k (z) (k = 1, 2) be any two space curves (B-spline curves, Bézier curves, or NURBS curves), then in z 1 ruled direction, a ruled surface T(z, z 1 ) [25] is created if a linear interpolation is implemented among iso-parametric points of the two curves b 1 (z) and b 2 (z). Also, a ruled surface has an astonishing property that, for fixed z, T(z, z 1 ) is a line segment joining points b 1 (z) and b 2 (z) (line segments are called straight generatrices, whereas the curves b 1 (z) and b 2 (z) are named as guide lines, for short GLs). This portion particularly demonstrates how the ruled surface with shape parameters is expressed using tensor product GBT-Bézier surface design. The necessary condition to express the ruled surface into tensor product form is that both should have the same degree. Let be the two GLs of the ruled surface of degree m, where f l k,m (z), k = 0, 1, . . . , m, l = 0, 1, are the GBTB functions described by (2). The GBT-Bézier ruled surfaces T ruled (z, z 1 ) should fulfil the following two conditions according to its definition [25]: (1) For fixed value of z, the curve T ruled (z, z 1 ) is a line segment.
(2) The guild lines fulfil the T ruled (z, 0) = F m,0 (z, μ, ν) and T ruled (z, 1) = F m,1 (z, μ 1 , ν 1 ). In accordance with the definition of GBT-Bézier surfaces, the equation of the GBT-Bézier ruled surfaces T ruled (z, z 1 ) that meet the aforementioned requirement is where μ, μ l , ν, ν l are the shape parameters; Q l,k (l = 0, 1; k = 0, 1, . . . , m) are the mesh points defined by (17) and f k,m (z) and f l,1 (z 1 ) are the GBTB basis functions given in (2) and (10), respectively. Equation (18) represents the tensor product GBT-Bézier surface in the form of a GBT-Bézier modified ruled surface. The GBT-Bézier ruled surfaces defined in (18) take over all characteristics and benefits of GBT-Bézier surfaces because these surfaces are expressed in the form of tensor product GBT-Bézier surfaces and the GLs are two GBT-Bézier curves that have the same degree. The NURBS modeled ruled surface defined in [25] has some limitations, such like complex mathematics, ambiguous weight elements, the annoyance of finding the surface order and the node vector [26]. On the other hand, the required surface has some benefits including straightforward computation, obvious implication of the shape parameters, and the free choice of vertexes.

Proposition 3.3 It is worth recalling that the straight generatrices of a GBT-Bézier ruled
surfaces are the straight lines that are obtained by joining iso-parametric points, instead of joining the same length of an arc of the two GLs. Generally, joining the same length of an arc of any two GLs gives us a different ruled surface.

Proposition 3.4 The technique of splicing GBT-Bézier surfaces can be utilized to generate all types of complicated combined GBT-Bézier ruled surfaces for the reason that the GBT-
Bézier ruled surface is a special case of tensor product GBT-Bézier surface. In this technique, the two adjacent GBT-Bézier ruled surfaces satisfy C 1 , C 2 , and C 3 continuity conditions because both GLs of the ruled surfaces satisfy these conditions (see Fig. 10 and Fig. 11).
The algorithm for designing a GBT-Bézier ruled surface is as follows.

Draw the GBT-Bézier ruled surface in Wolfram Mathematica software with suitable
values of shape parameters. 6. Two adjacent GBT-Bézier ruled surfaces can be connected by C 1 , C 2 , and C 3 continuity conditions given in Theorem 2.
Example 3.3 Let F 4,l (z, μ 1 , ν 1 ) (l = 0, 1) be two GBT-Bézier curves of degree 4 in a 3dimensional space, then from equation (18) the equation of the GBT-Bézier ruled surfaces formed by these GBT-Bézier curves as the GLs is given by are the control points of GLs. Figure 9 illustrates the GBT-Bézier ruled surfaces determined by (19) with various shape parameters.
Two adjacent GBT-Bézier ruled surfaces satisfying the C 1 and C 2 continuity in z direction can be presented in Figs. 10 and 11 with third and fourth degree GBT-Bézier curves taken as GLs of the right and left GBT-Bézier ruled surface, respectively. Figure 12 shows the C 3 continuity of two adjacent GBT-Bézier ruled surfaces that have forth degree GBT-Bézier curves as GLs in z direction.

GBT-Bézier swung surface along shape parameters
A swung surface can be generated from the surface of revolution [25]. Using a tensor product GBT-Bézier surface design with shape parameters, a swung surface is described in this portion. Let be the two GBT-Bézier curves described by (3), where Q k,0 = (Q X k,0 , 0, Q Z k,0 ) and Q l,1 = (Q X l,1 , Q Y l,1 , 0) are their control points, and F m (z, μ, ν) and F n (z 1 , μ 1 , ν 1 ) are known as profile curves (PCs) and trajectory curves (TCs), respectively. They are described in the XOZ plane and XOY plane, respectively. Vector form of the these curves and their corresponding nonzero coordinate functions F X m (z), F Z m (z) and F X n (z 1 ), F Y n (z 1 ) are used to construct a and the GBT-Bézier swung surface is defined by [25,26] T swung (z, where λ (λ > 0) is an arbitrary scaling factor.

Proposition 3.5
The GBT-Bézier swung surface expressed by (22) has the following geometric characteristics: if the profile curve F m (z) is swung about z-axis and re-scale in the X and Y directions, the swung surface T swung (z, z 1 ) is achieved in Fig. 13 with different values of scaling factor λ. 1. If the z line is fixed, the curves T swung (z, z 1 ) which have identical trajectory curve F n (z 1 ) but are scaled as λF X m (z) in X and Y directions are produced. 2. The fixed z 1 line produces curves Z swung (z, z 1 ) with the same profile curve F m (z), and the connection between the curves T swung (z, z 1 ) and the PCs F m (z) is: the resulting curves T swung (z, z 1 ) are achieved by rotating the PCs F m (z) into the plane with the vector (F X n (z 1 ), F Y n (z 1 ), 0), and use the scale factor λ|F n (z 1 )| to scale the rotated curve along X and Y coordinates while Z coordinate stays unscaled.
Thus, a tensor product GBT-Bézier surface design of degree (n, m) is used to express the swung surfaces (22).
The swung surface described by Theorem 4 is called GBT-Bézier swung surface. It appears that the GBT-Bézier swung surface carried all the characteristics and benefits of GBT-Bézier surface. The shape parameters are used to modify GBT-Bézier swung surface by keeping control net points fixed. For smooth connection among adjacent GBT-Bézier swung surfaces like previous surfaces, if the PCs or the TCs are combined by C 1 or C 2 continuity for GBT-Bézier curve, then the resulting GBT-Bézier swung surface is also a C 1 or C 2 continuous surface in the direction of profile or trajectory curve. It is worth noting that if the swung surface T swung (z, z 1 ) is either open or close in z and z 1 directions, then the corresponding profile curve F m (z) and the trajectory curve F n (z 1 ) both will be either open or closed. Additionally, if the scale factor λ = 1 and the trajectory curve F n (z 1 ) circle with a unit radius with a center at the origin, then the GBT-Bézier swung surface is a surface of revolution. However, the surface of revolution no longer is a tensor product GBT-Bézier swung surface which is described by Theorem 4.
The detailed algorithm for designing a GBT-Bézier swung surface is described as follows.

Algorithm 4
1. Take a profile curve in XOZ and a trajectory curve in the XOY plane in the form of GBT-Bézier curves. 2. Define a GBT-Bézier swung surface by using equation (23) It can be observed from Fig. 13 that the scaling of swung surfaces in both directions is handled by scale factor λ keeping shape parameters μ and ν same and fixed. Figures 14  and 15 exhibit the graphs of GBT-Bézier swung surfaces along the same control points of the PCs and TCs as taken in Fig. 13, but the values of shape parameters are different and scaling factor λ = 0.15 is fixed. Figures 16, 17, and 18 depict the C 1 and C 2 continuity connection of two adjacent GBT-Bézier swung surfaces for various values of shape control parameters and scaling factor λ.

GBT-Bézier swept surfaces with shape parameters
In this section, we especially tackle the problem of constructing a surface by sweeping a section curve (SC) around an arbitrary TC [11]. The GBT-Bézier curves of order m and n are taken as SCs and TCs in a 3D space, respectively. These two curves are given by Generally, a swept surface is given by [25] T swept (z, z 1 ; μ, ν, μ 1 , ν 1 ) = F n (z 1 ; μ 1 , ν 1 ) + S(z 1 )F m (z; μ, ν), where S(z 1 ) is a 3 × 3 matrix integrating scaling and rotation of the SC F m (z). In this study, the SC F m (z) and the TCs F n (z 1 ) both are GBT-Bézier curves, but the type and shape might Practically, swept surfaces can be generated by the following two categories: • M(z 1 ) is an identity matrix, which for all z 1 and F m (z) is just translated by F n (z 1 ).
• M(z 1 ) is not an identity matrix. The second type is very complicated and debatable, which can be referred to literature [25]. Particularly, we consider the first case of swept surfaces which can be defined by T swept (z, z 1 ; μ, ν, μ 1 , ν 1 ) = F n (z 1 ) + F m (z). (29) where Q k,l = Q k,m + Q l,n (k = 0, 1, . . . , m; l = 0, 1, . . . , n) are the control net points. Thus, the GBT-Bézier swept surface described by (29) can be exactly expressed as tensor product GBT-Bézier surface design.   Figure 19 and Fig. 20 depict the composite curves of two contiguous GBT-Bézier swept surfaces with different shape parameters which can be connected by C 1 and C 2 continuity constraints, respectively.

GBT-Bézier rotation surfaces with shape parameters
Geometric shapes of many objects are rotation surfaces in nature. Rotation surface is a type of special and common geometric profile in the field of product modeling design and has a broad demand in many fields such as aeronautic, architectural engineering, industrial designing, machining, and computer graphics. In CAD/CAM technology, how we can obtain 3D mathematical designs of rotation surfaces efficiently and swiftly is an important research issue. This portion is especially focused on how to develop a C 1 or C 2 continuous 3D mathematical model of GBT-Bézier rotation surface by taking a combined GBT-Bézier curve. Rotation surfaces are created by revolving a curve along a straight line in space with a constant angle. When the GBT-Bézier curves F m (z; μ, ν) (0 < z ≤ 1) are considered as the generating line and φ (0 < φ ≤ 2π ) as a fixed angle around the rotation axis, then a rotation surface R rotation can be constructed in the XOZ plane. Let us imagine that the generating line is a GBT-Bézier curve F m (z; μ, ν) (0 < z ≤ 1) in the XOZ plane and has Z-axis as the rotation axis, then the expression for rotation surface R rotation is described  The graphs of GBT-Bézier rotation surfaces using control points (35) are presented in Fig. 21 with different shape parameters μ, ν and with similar control polygon of generating lines.

Conclusions
This research work is mainly focused on the construction of engineering surfaces in Bézier approach. For this purpose, utilizing GBT-Bézier curves proposed in [24], six different types of tensor product GBT-Bézier engineering surfaces including GBT-Bézier cylinders, GBT-Bézier bilinear surfaces, GBT-Bézier ruled surfaces, GBT-Bézier swung surfaces, GBT-Bézier swept, and GBT-Bézier rotation surfaces are constructed. All geometric features of proposed GBT-Bézier surfaces are similar to the traditional Bézier, but the shape-adjustment feature is an additional feature that is not present in the traditional Bézier surface approach. The detailed algorithms for designing these GBT-Bézier engi- neering surfaces are also a part of this study. Furthermore, the construction of composite GBT-Bézier surfaces with C 1 , C 2 , and C 3 continuous connection is also presented. These special engineering surfaces not only have some parameters for their shape control but can also be expressed in terms of tensor product GBT-Bézier surface design (excluding rotation surfaces).
Some modeling samples of the proposed GBT-Bézier surfaces are presented here to illustrate that the designed GBT-Bézier surfaces can approach their convex hull better than the traditional Bézier surfaces approaches. Moreover, in the field of geometric modeling, these GBT-Bézier surfaces are very easy to implement. This study is meaningful and con-