Generalized Developable Cubic Trigonometric B é zier Surfaces

: This paper introduces a new approach for the fabrication of generalized developable cubic trigonometric B é zier (GDCT-B é zier) surfaces with shape parameters to address the fundamental issue of local surface shape adjustment. The GDCT-B é zier surfaces are made by means of GDCT-B é zier-basis-function-based control planes and alter their shape by modifying the shape parameter value. The GDCT-B é zier surfaces are designed by maintaining the classic B é zier surface characteristics when the shape parameters take on different values. In addition, the terms are deﬁned for creating a geodesic interpolating surface for the GDCT-B é zier surface. The conditions appropriate and suitable for G 1 , Farin–Boehm G 2 , and G 2 Beta continuity in two adjacent GDCT-B é zier surfaces are also created. Finally, a few important aspects of the newly formed surfaces and the inﬂuence of the shape parameters are discussed. The modeling example shows that the proposed approach succeeds and can also signiﬁcantly improve the capability of solving problems in design engineering.


Introduction
The developable structure, as a kind of special and meaningful ruled structure, may be expanded onto a plane rather than being stretched or broken. Because of this particular feature, developable surfaces are of great importance in geometric design [1,2] and have broad applications in the manufacture of materials such as automotive components, aircraft skins, ship hull pipes, shoes, and clothing [3][4][5][6][7]. The specification of the surface to be produced is therefore of great significance for the metal plate and sheet metal industries. There are three basic types of producing surfaces: a conical surface, a tangent surface, and a cylindrical surface. Of course, a complex mixture of such surfaces may also be created.
A variety of experiments on developable surfaces have been carried out by researchers in the construction of the developable surfaces. The material mainly involves the design and precise description of the developable surface [8], quasi-developable surface [9], fitting [10], geometrical restriction [11,12], and smooth splicing [13]. The numerical representation advances the use of developable surfaces in architecture and engineering that can be classified into two groups: one is the representation of geometric points and the other is the representation of geometric lines and planes [11,[14][15][16][17]. Aumann [14] suggested sufficient conditions for the construction of a developable Bézier surface with two boundary curves confined to parallel planes and deduced a basic criterion. The general characterization (weights and control nets) of the developable rational Bézier surface was introduced by Lang and Roschel [15]. The biggest downside to this approach is that the characterizing equations are not linear. Nevertheless, another model of geometry proposed by Pottmann and Farin [16] is projective geometry, which describes developable surfaces as a projective space curve and avoids the non-linearity of the characterizing equations. However, to research developable surfaces, there is another much more suitable model called Laguerre geometry. Laguerre geometry is typically a geometric analysis of oriented circles and lines in two dimensions, or of higher dimensional orientated spheres and hyperplanes where the transformations of Laguerre reflect bijective maps that retain oriented contact [18,19]. Since a developable surface conserves tangent planes along each of its rulings, it is transformed into a curve in the Laguerre geometry [20]. Moreover, a study of using Laguerre geometry for envelopes of rotational cones and plane-based higher order interaction between revolution cones and a given reference surface that operates on a point model of a group of aligned planes is provided in [20]. However, a new dimension of Laguerre geometry is its function in geometrical design work. It turns out that if one uses the Laguerre geometry, it is surprisingly simple to overcome many geometric architecture issues [21].
Bodduluri and Ravani [17] created a new geometric representation for forming surfaces using the duality that occurs between the point and the plane. Meanwhile, [11,22] have designed the surface to interpolate curvature and geodesic curves, respectively. The above methods are almost limited to only the surfaces of the Bézier or B-spline. It is therefore uncontrollable, and due to its limited degrees of freedom, it is difficult to alter the structure of the surface. These drawbacks make it very difficult to follow the requirements of practical engineering. Of course, the use of a rational Bézier surface can be considered, but a higher power of rational Bézier can lead to various shortcomings as stated in [23,24]. In rational curves and surfaces, each control point has a weight and the collection of that weight is generally uncertain. That added weight freedom will create greater nuisance than genuine assistance [25]. Even transcendental curves like the helix and cycloid cannot be accommodated by the rational model. In addition, the repeated differentiation of the rational curve produces a higher degree of curve.
Some researchers have implemented curves and surfaces with shape parameters to address these deficiencies. Refs. [26,27] constructed a QTB-spline and a rational cubic trigonometric (CT)-Bézier curve with adjusting parameters respectively. Majeed and Qayyum [28] constructed a new rational cubic trigonometric B-spline where this proposed curve applies with various applicability and flexibility by using different weights and shape parameters. In [29], the authors have derived C 3 and C 5 continuity of the cubic trigonometric B-spline curves with uniform and non-uniform knots. Bashir et al. [30] provided a quadratic trigonometric (QT)-Bézier curve with two shape parameters. Yang et al. [31] explored the extension of the Bézier quartic trigonometric curve of G 2 and C 3 continuity. The quintic trigonometric Bézier curve with two dynamic shape parameters was developed by Misro et al. [32], which was then used to build five transition curve templates [33] and to estimate the maximum speed in highway designs [34]. Ammad and Misro [35] constructed a biquintic trigonometric Bézier surface with four shape parameters. Hu et al. [36,37] proposed a method for constructing generalized developable Bézier and H-Bézier surfaces. Recently, Bibi et al. [38] modeled symmetric revolutionary and rotational surfaces using a hybrid trigonometric Bézier surface that was later applied in engineering applications [39].
This paper attempts to resolve the issue of shape handling in developing surfaces using the cubic trigonometric Bézier basis functions and to analyze their favorable properties by extending the discussion to curves and surfaces. By enveloping developable and tangent curves of the spine, we build developable surfaces using a generalized developable cubic trigonometric Bézier (GDCT-Bézier) basis function. However, the shape of the surface is being managed by the design variables, making the design process of a particular product more manageable.
The rest of the paper is structured as follows. The definition of the GDCT-Bézier curve and properties are set out in Section 2. A new technique for developing GDCT-Bézier surfaces is presented in Section 3. The continuity conditions for the GDCT-Bézier surfaces are derived in Section 4. Practical examples are given in Section 5. Finally, Section 6 provides a brief conclusion.

Cubic Trigonometric Bézier Basis Functions
The definition of a cubic trigonometric Bézier (CT-Bézier) basis function is given as follows [40].
The CT-Bézier basis functions for two randomly chosen real values of shape parameters γ 1 and γ 2 are shown in Figure 1. The proof of Theorem 1 is provided in [40]. in Section 3. The continuity conditions for the GDCT-Bézier surfaces are derived in Section 4. Practical examples are given in Section 5. Finally, Section 6 provides a brief conclusion.

Cubic Trigonometric Bézier Basis Functions
The definition of a cubic trigonometric Bézier (CT-Bézier) basis function is given as follows [40].
The CT-Bézier basis functions for two randomly chosen real values of shape parameters γ 1 and γ 2 are shown in Figure 1. The proof of Theorem 1 is provided in [40].
By using the interpretation of the basic functions, the CT-Bézier curves have the following properties.
(a) Boundary properties: The CT-Bézier curve has a shape that is independent of the selection of the coordination, i.e., Equation (2) satisfies the following two equations: where q is an arbitrary vector in R 2 or R 3 and T is an arbitrary d × d matrix, d = 2 or 3. (d) Convex hull property: The entire segment of the CT-Bézier curve must lie inside the control polygon.

Dual Generation of Single-Parameter Family of Planes
Duality is a concept adapted from projective geometry. This concept is useful to construct developable surfaces. As per the duality theory among points and planes, the single parameter family of control points of a curve is dual to a single parameter family of planes. The single parameter family of planes {∏ u } can easily be reached if the control points R i in a CT-Bézier curve are treated as a control plane. Thus by Equation (2), the single parameter family of plane {∏ u } can be defined as where R i (i = 0, . . . , 3) are control points of {∏ u }, u is the family of parameter, and γ 1 , γ 2 are shape control parameters of {∏ u }. In a three-dimensional projective space, if the coordinates of the control plane are (3) can be written as where,

Generalized Enveloping Developable CT-Bézier Surface
Based on the properties and definition of a developable surface, it is comprehended that the envelope of a single-parameter family of planes is a developable surface. The plane corresponding to every value of u in Equation (2) can be defined as Equation (6) when differentiated will yield, According to the third definition from [41], the intersecting line between plane Equations (6) and (7) is the generator of a developable surface, which is symbolized by L(u; γ 1 , γ 2 ) and can be computed in terms of its Plücker coordinates L(u; γ 1 , and Assume φ(u) to be the nearest point to the origin that lies on the two planes, so it satisfies By solving this system of equations, we get Therefore, the line L(u; γ 1 , γ 2 ) can be represented as when the parameter u family range is [0, 1] and the line L(u; γ 1 , γ 2 ) produces a developable surface with shape parameter γ 1 , γ 2 . It is apparent from the above calculation process that a developable surface can be created when the control planes are provided.

Generalized Spine Curve Developable CT-Bézier surface
This section includes one more approach in constructing developable surfaces. With regard to Definition 3 in [41], the three successive planes in the family of a plane {∏ u } intersect at a characteristic point and can be achieved by intersecting Equations (6) and (7), and the second derivative of Equation (6) produces the third plane.
Significantly, the characteristic argument can be reached by solving the following set of equations: Hence, the coordinate of the characteristic point ρ(u) is given by When parameter u varies in the range [0, 1], ρ(u) can generate a spine curve which can be represented in parametric form such as below:

Developable Surface Interpolating Geodesic CT-Bézier Curve with Parameters
A geodesic is the shortest possible line between two points on a sphere or other curved surface. This section involves the construction of developable surface through a given cubic trigonometric Bézier curve, where the curve is the geodesic of the surface. Proof. Suppose that c(u) is a curve on the surface T(m, u) and φ is an arbitrary point on the curve c(u), N(m 0 , u) is the normal vector at point φ, ζ(u) is the principal normal, θ is the angle between ζ(u) and N(m 0 , u) and k g is the geodesic curvature. Then, according to the definitions of geodesic curve and k g k g = ±k sin θ = 0.
Since k = 0, the θ = 0 or θ = π then ζ(u) is ||N(m 0 , u). This proves the theorem. Theorem 3. When a CT-Bézier curve with shape parameters γ 1 , γ 2 is given, a developable surface that is geodesic to the developable surface must exist through it.
Proof. Given a CT-Bézier curve c(u; γ 1 , γ 2 ), we have the respective ruled surface When a curve c(u; γ 1 , γ 2 ) is geodesic to Theorem 2, we have where ζ(u; γ 1 , γ 2 ) and N(0, u; γ 1 , γ 2 ) are the principal normal and normal vector of the curve. For simplicity we omit the variable and parameters. Because ζ(u; γ 1 , γ 2 ) is parallel to N(0, u; γ 1 , γ 2 ) and c is perpendicular to N, the vector c 1 must be in the plane spanned by c and c ×c . Let According to the developable conditions, we have that is, Which interpolate the given curve c(u) and c(u) is a geodesic of the developable surface, where µ 1 and µ 2 satisfy the Equation (14).
In a special case, if we choose l2(t) = 1.

Analysis Properties of the GDCT-Bézier Surface
According to the single parameter family of planes {∏ u }, the 1st and 2nd derivatives of Equation (3) are Therefore, the planes corresponding to u = 0 and u = 1 in {∏ u } have the following characteristics: and The first Equations of (17) and (18) shows that the first and last planes corresponding to u = 0 and u = 1 in {∏ u } are given by the designer as the first and last control planes and they are both tangential to the GDCT-Bézier surface along its generator, at u = 0 and u = 1. In addition, the generator of the GDCT-Bézier surface at u = 0 is the intersection between two planes, represented by the first two Equations of (17), called the starting generator. Therefore, the generator L(0; γ 1 , γ 2 ) is the intersection of two planes R 0 and π 2 (2 + γ 1 )(R 1 − R 0 ) and can be written as where r 0 = (a 0 , b 0 , c 0 ), r 1 = (a 1 , b 1 , c 1 ), and X = (x, y, z). R 1 can be accomplished by the linear combination of two planes in Equation (19) and the generator L(0; γ 1 , γ 2 ) is the intersection of R 0 and R 1 . Furthermore, the generator L(0; γ 1 , γ 2 ) dual to the line that connects the points R 0 and R 1 , which is tangential to GDCT-Bézier curve defined by R 0 , R 1 , R 2 , R 3 . Similarly, the ending generator L 1;γ 1 ,γ 2 corresponding to u = 1 is the intersection of R 2 and R 3 in Equation (18).
In addition, the Equations (17) and (18) can be written in matrix form as: and Since the intersection of the S, S , S planes is the spine curve of the GDCT-Bézier surface, the characteristic point and plane coordinates of R 0 , R 1 , R 2 will be achieved through the linear combination of S, S , S according to Equation (20). So, the intersection of plane coordinates R 0 , R 1 , R 2 is the spine curve GDCT-Bézier surface characteristic point on T(m, 0; γ 1 , γ 2 ). Similarly, in Equation (21), the control points R 1 , R 2 , R 3 define the characteristic point on T(m, 1; γ 1 , γ 2 ).

Continuity Conditions between GDCT-Bézier Surfaces
Many complex engineering surfaces are built in on CAD/CAM systems, but the construction of a complex surface with a single surface is highly challenging. Therefore, the promising features seem very important to ensure a smooth relationship between the adjacent GDCT-Bézier surfaces [13,42]. However, the G 1 , Farin-Boehm G 2 , and G 2 Beta continuity are currently the standards for continuity measurement between two adjacent developable surfaces [17,43,44].

The G 1 Continuity Conditions of GDCT-Bézier Surfaces
In accordance with the description and analysis for G 1 smooth continuity between two parametrical curves in [43], the G 1 continuity conditions in Equation (22) can be computed as follows. Theorem 4. The necessary and sufficient conditions for G 1 smooth continuity between two adjacent GDCT-Bézier surfaces in Equation (22) at the joint are where ϕ > 0 is a constant.
This implies that for the stability of G 1 , the two surfaces first require a shared control plane. Next, based on Equations (17) and (18), S 1 (1) and S 2 (0) can be given by Combining Equation (26) according to the second Equation of (24), we get Finally, from Equation (25), the Equation of (27) can be written as Thus, Equations (25) and (28) are the G 1 continuity conditions of the two adjacent GDCT-Bézier surfaces.

Farin-Boehm G 2 Continuity Conditions of GDCT-Bézier Surfaces
Theorem 5. The conditions needed and necessary for Farin-Boehm G 2 continuity between the two GDCT-Bézier surfaces in Equation (22) are Proof. In order to obtain the Farin-Boehm G 2 continuity, the two surfaces need to satisfy the following Based on Equations (17) and (18), we have and Substitute Equations (31) and (32) into Equation (30), the conclusion in Equation (29) can be obtained.

G 2 Beta Continuity Conditions of GDCT-Bézier Surfaces
Theorem 6. The conditions needed and necessary for G 2 Beta continuity continuity between the two GDCT-Bézier surfaces in Equation (22) are where ϕ > 0, and ψ is an arbitrary constant.

Design Examples of GDCT-Bézier Surface
In this section, some examples are presented to demonstrate the construction of cubic trigonometric surfaces by means of control planes. Here, the plane coordinates (a i , b i , c i , d i ) of each R i (i = 0, 1, 2, 3) control plane can be defined by its center point (a i , b i , c i ) and d i = a 2 i + b 2 i + c 2 i is the distance. It should be noted that this sort of representation is not valid where the origin is the center point.
Distinctly, the envelop and spine curve family of developable surfaces, as well as piecewise-developable surfaces with some smooth continuity, are being generated by the proposed method used in this paper. The resulting shapes are shown in Figures 2-8. These surfaces can be viewed from different angles and the designer can change the shape manually by adjusting their shape parameter values without re-determining the control planes, which solves the problems of the methods in [1,[45][46][47][48][49][50]. Moreover, the GDCT-Bézier surfaces reserve many properties of the traditional Bézier surfaces, especially when the shape parameters have values equal to 1.

Examples of Enveloping GDCT-Bézier Surfaces
In this example, the construction of the cubic trigonometric enveloping developable Bézier surface is shown. Assume the coordinates of the control planes are Here, we presume that in designing the enveloping developable surfaces, the center points of four control planes are coplanar and each has an equal distance to the origin. When γ 1 , γ 2 assumes

Examples of Enveloping GDCT-Bézier Surfaces
In this example, the construction of the cubic trigonometric enveloping developable Bézier surface is shown. Assume the coordinates of the control planes are Here, we presume that in designing the enveloping developable surfaces, the center points of four control planes are coplanar and each has an equal distance to the origin. When γ 1 , γ 2 assumes different values, a family of enveloping GDCT-Bézier surfaces with different shapes can be constructed under the conditions of the specified control planes (see Figures 2 and 3).

•
When modifying the value of γ 1 and keeping γ 2 unchanged, the T(1; γ 1 , γ 2 ) generator retains the same length and location. However, the length of the generator T(0; γ 1 , γ 2 ) becomes longer when we increase the value of γ 1 , but its position remains unchanged.

•
If the value of shape parameter γ 1 is constant and γ 2 is adjusted, the generator T(0; γ 1 , γ 2 ) retains the same length and position. However, the position of generator T(1; γ 1 , γ 2 ) also remains unchanged, but its length increases when we increase the value of γ 2 .

Examples of Spine Curve GDCT-Bézier Surfaces
In this example, the construction of the spine curve GDCT-Bézier surface is demonstrated. Assume the coordinates of the control planes are In designing the developable surface of the spine curve, we assume that the four control planes have their center points in four different quadrants and that the distance to the origin is not equivalent from these center points.
In Figures 4 and 5, we can see the impact of shape parameters γ 1 , γ 2 on the spine curve developable surface as follows: • When modifying the value of γ 1 , and keeping γ 2 unchanged, the T(0; γ 1 , γ 2 ) generator retains the same length and location. However, the length of the generator T(1; γ 1 , γ 2 ) becomes shorter when we increase the value of γ 1 , but its position remains unchanged.
• If the value of shape parameter γ 1 is constant and γ 2 is adjusted, the generator T(1; γ 1 , γ 2 ) retains the same length. However, the length of the generator T(0; γ 1 , γ 2 ) becomes shorter when we increase the value of γ 2 , but its position remains unchanged.

Example of Developable Surface Interpolating Geodesic Cubic Trigonometric Bézier Curve with Parameters
In order to illustrate the effectiveness of the method constructed in Section 3.4, we design a developable surface with a special case when µ 2 = 1. The resulting developable surface interpolates the curve created by control points b 0 = [3, 3,3], b 1 = [4, 3,5], b 2 = [5, 3,5], b 3 = [6, 3,3]. Although the curve is geodesic to the developable surface, note that the shape parameters have an influence on the developable surface as shown in Figures 6  and 7.

Conclusions
Based on the CT-Bézier basis function with two shape parameters, this article has presented a novel approach for constructing a class of GDCT-Bézier function with local control for enveloping and spine developable surfaces. Furthermore, the conditions and evaluation of some interesting properties of the GDCT-Bézier surface to create geodesic interpolating surfaces are successfully derived. The shapes of the GDCT-Bézier surfaces can be easily modified due to the inclusion of a shape parameter without changing the control planes. This suggests that the Bézier trigonometric surfaces have increased surface design efficiency with added versatility. This developable surface inherits the properties of the classical developable Bézier surface. In addition, in order to handle complex shapes in engineering and architectural design, we have also derived sufficient and necessary continuity conditions between two adjacent GDCT-Bézier surfaces. Theoretical discussion and design examples, however, suggest that the GDCT-Bézier surfaces and their smooth continuity are not only easy and convenient to integrate but also provide a great deal of flexibility in the design of developing surfaces. Thus, they are prominently useful in constructing complex developable surfaces in CAD/CAM with different degrees of smoothness. However, our future aim is to investigate the approximation and interpolation of cubic trigonometric Bézier surfaces by exploring special algorithms for developable surfaces.