Construction of Regular Developable B é zier Patches

: The modeling of developable surfaces is considered a very important application in plat-metal-based industries. Relating to the purpose, this discussion aims to obtain some formulas for constructing the regular developable B é zier patches, in which each boundary curve must be laid in two parallel planes. The results as follows: We ﬁnd some formulas of the equation systems that are described by the constant, linear, and quadratic control parameters of the regular developable B é zier patches criteria. The new approach is numerically tested for constructing the regular developable B é zier patches, in which their boundary curves are deﬁned, respectively, by the combination of four, ﬁve, and six degrees.


Introduction
Some methods for designing the developable surfaces have been reported.Aumann [1] defined a developable Bézier patch with two boundary curves (directrices).They are restricted to lie in parallel planes, their projection in the x-y plane is a rectangle, and the tangent vectors to each ruling are parallel.Frey and Bindschadler [2] extended the work of Aumann by generalizing the degree of the directrices.Via linear and quadratic form of regularity condition for Bézier developable, they can generate developable Bézier strip patches that allow their projections in the x-y plane to be trapezoids.In general, the method is simple and easy to use in modeling.However, the method to find the Bézier control points of the equations system is still not stable and deterministic.Chu and Séquin [3] and Aumann [4] developed a technique to determine the geometric design of developable Bézier surfaces based on the de Casteljau algorithm.Chalfant [5] proposed a method to construct a B-spline developable surface bounded by two directrices that lie on a parallel plane.Zhao and Wang [6] presented an approach to design a developable surface in the algebraic form using the surface pencil passing through a given curve.They derived the necessary and sufficient conditions for a surface to be developable by representing the surface with the combination of the given curve, and the three vectors decomposed along the directions of the Serret-Frenet frame.Meanwhile, Al-Ghefari and Abdel [7] presented an approach for designing a developable surface, possessing a given curve as the geodesic of it by analyzing the necessary and sufficient conditions when the resulting developable surface is a cylinder, cone, or tangent surface.Based on results from the optimal control of linear systems with quadratic criteria, Park et al. [8] introduced a mathematical framework for constructing developable surfaces.From a set of rulings, they analyze the existence of base curves, such that the resulting surface is developable.Then, they proposed a class of objective functions that are based on various physical criteria for constructing the base curve.In the context of defining developable approximation surfaces, Elber [9] presented an approximation method that trims the surface using isolines and interpolates each trimmed piece by ruled patches.Sun and Fiume [10] introduced a new modeling system that can be used to create developable surfaces from scratch.They divided a surface to be modeled into several pieces, then defined the shape of each surface piece with a generalized cone.
After that, they constructed separately each of 3D surface piece using piecewise continuous curves of the cross-section of the cones, then stitched together the surface pieces to construct the desired developable surface.Meanwhile, Xu et al. [11] presented a novel approach to constructing minimal surface from a given boundary by applying quasi-harmonic Bézier approximation and quasi-harmonic mask methods.Then, Xu et al. [12] proposed a general framework for constructing IGA-suitable planar B-spline parameterizations from given complex CAD boundaries.This method needs some complex operations to obtain the type surface, e.g., an operation of Bézier extraction and subdivision that are performed on each boundary curve to generate a planar parameterization, a robust planar domain partition framework, quadrangulation, and a local optimization process.A different approach of the methods is presented by Bodduluri and Ravani [13], Pottmann and Farin [14], and Pottmann and Wallner [15].They constructed the developable surfaces in terms of plane geometry using the concept of duality between points and planes in a 3D projective space.Mathematically, this method is interesting, but it may be less practical in the design of 3D shapes.Hu et al. [16] designed the generalized developable H-Bézier surfaces using control planes with generalized H-Bézier basis functions.Based on the Bodduluri and Ravani idea, Hu et al. [17] developed this method to define the developable Bézier-like surfaces with Bernstein-like basis functions of the multiple shape parameters.The shapes of both type surfaces can be adjusted by altering the values of the shape parameters in that basis functions.Choosing some different values of these shape parameters will produce, respectively, a family of developable H-Bézier and Bézier-like surfaces, without changing the characteristics of Bézier surfaces.The introduced methods generally have some of the following restrictions: The characteristic equations for constructing developable surfaces cannot be solved easily, the planar boundary curves of low degree are more useful for designing developable surfaces, and it is difficult to determine the vertex for a cone and the edge of regression for a tangent surface.To avoid these disadvantages, we propose a new method for designing regular developable Bézier patches, based on the work of Frey and Bindschadler in terms of the projections of the patches in the x-y plane as trapezoids, and their boundary curves to be defined by the combination of four, five, and six degrees.
The paper is organized in the following sections: In the first section, we discuss condition and regularity of the developable Bézier patches, which are supported by two parallel planes.In the second, we talk about the construction of the regular developable Bézier patches models that are defined by two curves planes from the combination of four, five, and six degrees.Finally, the results will be summarized in the conclusion section.

Mathematical Formulations of the Regular Developable Bézier Patches
The mathematical formulations of the regular developable patches in the algebraic form and Bézier expression are defined as follows: with f: I → R 3 and g: I → R 3 of class C n in I and o / ∈ g(I).
Here, (I, f) is called the directrix curve and g(u) in the lines form is called the generatrix.
Definition 2. The ruled surface S(u,v) = f(u) + v g(u) is developable, if all u ∈ I the vectors [g (u),f (u),g(u)] are coplanar, namely there exists unique real scalars δ(u) and σ(u) such that Definition 3. Developable Bézier surface is developable surface that can be constructed by two Bézier curves C 1 (u) and C 2 (u) in the form The Bézier curves C 1 (u) and C 2 (u) are defined by the form . Definition 2 implies that the ruled Bézier patch S(u,v) is developable, if the vectors g (u), C 1 (u) and g(u) are coplanar.This means that it exists two real scalars α(u) and σ(u)such that If C 1 (u) and C 2 (u) are to be conditioned respectively in two parallel planes [Ψ 1 ,Ψ 2 ], the developable condition (1) can be simplified as Meanwhile for its regular condition, Frey [2] and Kusno [18] state that the developable Bézier patch S(u,v) is defined by two curves [C 1 (u),C 2 (u)] respectively on two parallel planes [Ψ 1 ,Ψ 2 ] are regular, if the choice of the real scalar ρ(u) > 0.

Construction of the Regular Developable Bézier Patches
Mathematically, one aspect of the most exciting properties of developable surfaces is their ability to be laid flat on a plane without stretching, tearing, or deforming.So, studies about the developable surfaces are considered very important in application in plat-metal-based industries, including ship hulls, aircrafts industries, trains, and automobile parts, such as upholstery and body panels [2,8,9].For that purpose, each small surface section (patch) of those objects, generally, must be modeled by their two boundary curves C 1 (u) and C 2 (u) that respectively lie in the plane Ψ 1 parallel to the plane Ψ 2 .The problem that will be discussed here is how to construct the developable Bézier patches by using both curves.
Consider the regular developable Bézier patches that are defined by two curves C 1 (u) and C 2 (u) of degree n and m, respectively.We note it as [RDBP(n,m)] and the curves C 1 (u) and C 2 (u) are expressed as Then, we evaluate some different value of control parameter ρ(u) in Equation ( 2) that can affect and modify the shape of the modeled regular developable Bézier patches.

Case: ρ(u) Constant
We set two Bézier curves C 1 (u) and C 2 (u) of the same degree n, and the scalar ρ(u) is positive constant, i.e., ρ(u) = α ∈ R + .From the condition (2), we have Because of i = 0, . . ., (n − 1), the polynomials for all i = 0, . . ., (n − 1).When we add those equations, we will find an equation of the Bézier polygon control points of the curves C 1 (u) and C 2 (u) as follows So, to construct a regular developable Bézier patch which is supported by two curves C 1 (u) and C 2 (u) of degree n and conditioned by ρ(u) positive constant, we note as [RDBP(n,n)] C , must verify the Equations ( 4) and ( 5), namely 1.
The two vectors parallel (q n − q 0 ) and (p n − p 0 ) must be in the same direction to determine α value; 2.
In relation to treating the geometric continuity order 1 between two adjacent joined patches, in Equation ( 4), it is necessary to set four boundary control points [p 0 ,q 0 ,p n ,q n ] and two control points [p 1 ,p n−1 ] of the Bézier curve C 1 (u) that are used as the fixed data points.Because of the practical reasons, we are interested in the construction of [RDBP(4,4)] C and [RDBP (6,6)] C based on those constrains.
When the conditions (p n − p 0 )//(q n − q 0 ) have been verified, the curves shape C 1 (u) and C 2 (u) that can define [RDBP(4, 4)] C are fully be affected by the determination of four Bézier control points [p 0 ,p 1 ,p 3 ,p 4 ] for the curve C 1 (u) and two Bézier control points [q 0 ,q 4 ], tangent vector (p 1 − p 0 ) and (p 4 − p 3 ) for the curve C 2 (u). Figure 1a shows that the Bézier control points q 1 and q 3 are defined, respectively, by [q 0 ,p 0 ,p 1 ] and [q 4 ,p 3 ,p 4 ].The Bézier control points q 2 and p 2 , respectively, depend on the control points [q 1 ,q 3 ] and [p 1 ,p 3 ].By changing the two tangent vectors direction (Figure 1c,d), we can modify the curve shape C 1 (u) and C 2 (u).On the other hand, the curves shape C 1 (u), and C 2 (u) of [RDBP (6,6)] C depend on some Bézier control points of both curves.Figure 1b shows that we can elect some Bézier control points of the curve C 1 (u) and C 2 (u) respectively in plan Ψ 1 and Ψ 2 to obtain the various forms of that curves.Figure 1g presents one result example of that method.To ensure that the patches are developable, Figure 1e,f,h show that (at all parameter values, u, certainly) its tangent vectors and generatrixes must be coplanar.Based on the result of the discussion, we can conclude that the construction method [RDBP(n,n)] C can facilitate to model the shape of developable patches.

Case: ρ(u) Linear
Let two Bézier curves, C 1 (u) and C 2 (u), respectively, of degree n and n + 1, in two parallel planes [Ψ 1 ,Ψ 2 ].The scalar ρ(u) = α 0 (1 − u) + α 1 u is a real positive.We have the identity [19] (1 where Therefore, from Equations ( 2) and ( 8), we have Because of i = 0, . . ., (n − 1) the polynomials B n−1 i (u) are not zero, this mean that for all i = 0, . . ., n.When we add those equations, we will find an equation of control points of both curves and the centroid point of the Bézier polygon C 1 (u) as follows and p = 1 n ∑ n i=0 p i .Thus, to construct the regular developable Bézier patches with the two Bézier curves C 1 (u) and C 2 (u) of degree n and n + 1 that lie, respectively, in two parallel planes [Ψ 1 ,Ψ 2 ] and the parameter ρ(u) is positive linear, that is noted [RDBP(n,n+1)] L , can be summarized as follows: 1.
Because α 0 and α 1 must be positive, the position of the vector (q n − q 0 ) in Equation ( 10) must lie on the interior of the angle between two vectors (p − p 0 ) and (p n − p).If the control points [p 0 ,p n ,q 0 ,q n+1 ] are fixed, then the selection of the centroid point p must satisfy the conditions that the direction of the vector (p − p 0 ) and (p n − p) are in the different sides to the vector (q n − q 0 ).

2.
The value of α 0 and α 1 in Equation ( 10) can be calculated by To find the unique solution of the system (9), we use the strategy as follows.It has to determine (n + 2) control points of (2n + 3) control points in the system, such that the number (n + 1) calculated control points are the same as the number of these equations, i.e., (n + 1) equations.Also, the determinant of its coefficients matrix of the system must be different from zero.In relation to treating the geometric continuity order 1 between two adjacent joined patches, it is necessary that the points [p 0 ,p 1 ,p n−1 ,p n ,q 0 ,q n+1 ] are included in the list of the (n + 2) fixed points.
with λ∈ (0, 1) and δ, γ > 0. In the other hand, for τ > 0, we can choose p in other form (Figure 2c) This mean that if αo ≠ α1, then the unique solution of Equations ( 14  In case of the construction [RDBP (5,6)]L, if we determine the control points [p0, p1, p4, p5,q0, q3, q6] are shown in the black points of Figure 4a, then from the equations (10) we will  Let [RDBP(4,5)] L and from Equation ( 9), for i = 0, . . ., 4 we find the control points Equation (10) in the form (q 5 − q 0 ) = α 0 (p − p 0 ) + α 1 (p 4 − p).(12) In order to obtain an unique solution of the equation system, if the points [p 0 ,p 1 ,p 3 ,p 4 ,q 0 ,q 5 ] are fixed (Figure 2a) by constraint (1), we can determine the centroid point p of Bézier polygon of the curve C 1 (u), such that the direction of the vectors (p − p 0 ) and (p 4 − p) are in the different sides to the vector (q 5 − q 0 ).Clearly, in Figure 2b,d,e, let in plane Ψ 1 , the vectors (b − p 0 ) = (p 4 − c) = (q 5 − q 0 ) and (p 4 − p 0 ) = (d − p 4 ) = (p 0 − a).So, to facilitate the modeling patch, one of our options is to put the centroid point p in position with λ ∈ (0, 1) and δ, γ > 0. In the other hand, for τ > 0, we can choose p in other form (Figure 2c) with respectively the unity vectors that are perpendicular to (b − a) and (d − c) in Ψ 1 .After that, we calculate α 0 and α 1 in Equation ( 12) that have been formulated in Equation (11).Moreover, from the equation system (9) we have five equations to calculate the five control points [q 1 ,q 4 ,p 2 ,q 2 ,q 3 ] and we will find −5 q 2 + 5 The application examples of the method are shown in Figure 2d-f.The Equation ( 14) have a unique solution if its determinant of the co-efficient matrix is different from zero, i.e., This mean that if α o = α 1 , then the unique solution of Equations ( 14) are In Figure 3a,b, we show the effect of the centroid points p of the [RDBP(4, 5)] L when they are positioned in above and below of (p 4 − p 0 ).Meanwhile, Figure 3c,d present that along parameter values u of the boundary curves C 1 (u) and C 2 (u), tangent vectors and the generatrixes of the patches are coplanar.
The control points setting to model the shape of the curve C 1 (u) and C 2 (u) of construction [RDBP(n,n+1)] L are principally like the construction [RDBP(n,n)] C .However, the condition of two fixed vectors (p n − p o ) and (q n+1 − q 0 ) of the forth boundary control points [p o ,p n ,q 0 ,q n+1 ] must not be parallel.Also, we can change the centroid point P in position to the left, right, up, or down to modify the surface shape of the wanted patches as shown in the examples Figure 2d-f, Figure 3a,b, and Figure 4b-d.To ensure that the patches are developable, Figure 3c,d presents the tangent vectors and generatrixes of the patches that must be coplanar (at all parameter values, u, certainly).
Because α 0 and α 1 must be positive, the election α 1 in Equation ( 19) must verify the expression (20) or the vector [(q n+2 − q 0 ) − α 1 v] in Equation ( 19) must lie on the interior of the angle between two vectors (p − p 0 ) and (p n − p).
To find the exact solution of the system (18), we have to determine (n + 2) control points of (2n + 4) control points in the system, such that the number (n + 2) calculated control points are the same as the number of these Equation (18).The determinant of its coefficients matrix is different from zero and the control points [p 0 ,p 1 ,p n ,p n−1 ,q 0 ,q n+2 ] must be included in the list of the (n + 2) fixed control points.
Let [RDBP(4,6)] Q , and from the Equation ( 18) for i = 0, . . .,5, we find the control points and the centroid point Equation ( 19) in the form The points [p 0 ,p 1 ,p 3 ,p 4 ,q 0 ,q 6 ] are fixed and the elected value α 1 verify Equation (20).Then, by constraint (1), we can determine the centroid point p in a position likely in the case of ρ(u) linear, i.e., expressing of the formulae (13).The value α 0 and α 1 in Equation ( 21) can be determined by the form (11).After that, we can calculate the control points [q 1 ,q 5 ,p 2 ,q 2 ,q 3 ,q 4 ] in the equation system (18) of the six equations and case of its determinant of coefficients matrix is different from zero, i.e., α 2 = α 0 , we get where Based on the discussion, results were obtained as follows: The construction method [RDBP(n,n + 2)] Q .has many advantages for modeling the shape of regular developable patches.First, it can facilitate the determination of some control points of both curves C 1 (u) and C 2 (u) to model the curves.Second, by moving the centroid point position p and changing the parameter value α 1 , it will modify the various forms of the patches.Figure 5a,b show the patches shape changing when the position p is moved down and Figure 5c,d present the patches shape changing when the value α 1 = (0.6).w and α 1 = (−0.8).w.
[RDBP(n,n+2)]Q.has many advantages for modeling the shape of regular developable patches.First, it can facilitate the determination of some control points of both curves C1(u) and C2(u) to model the curves.Second, by moving the centroid point position p and changing the parameter value α1, it will modify the various forms of the patches.Figure 5a

Conclusions
We developed the method to design a regular developable Bézier patch through two given Bézier curves from the combination of four, five, and six degrees laid respectively in two parallel planes.Using the determination of some Bézier control points of the curves and the parameter fixation of the developable condition, we can control and model various shapes of the Bézier developable patches as follows: The shape [RDBP(n,n)]C with ρ(u) positive constant is entirely determined by the conditions of the vectors (pn − p0)//(qn − q0) and some fixed Bézier control points for Bézier curves C1(u) and C2(u).The shape modeling [RDBP(n,n+1)]L with ρ(u) positive linear depends on some fixed Bézier control points and the determination of centroid point position p of the Bézier polygon of curve C1(u).Finally, the shape of patches

Conclusions
We developed the method to design a regular developable Bézier patch through two given Bézier curves from the combination of four, five, and six degrees laid respectively in two parallel planes.Using the determination of some Bézier control points of the curves and the parameter fixation of the developable condition, we can control and model various shapes of the Bézier developable patches as follows: The shape [RDBP(n,n)] C with ρ(u) positive constant is entirely determined by the conditions of the vectors (p n − p 0 )//(q n − q 0 ) and some fixed Bézier control points for Bézier curves C 1 (u) and C 2 (u).The shape modeling [RDBP(n,n + 1)] L with ρ(u) positive linear depends on some fixed Bézier control points and the determination of centroid point position p of the Bézier polygon of curve C 1 (u).Finally, the shape of patches [RDBP(n,n + 2)] Q with ρ(u) positive quadratic can be modeled by the fixation of some control points of the Bézier curves C 1 (u) and C 2 (u), the modification position of the centroid point p, and the given parameter value α 1 .
Using two curves that lie, respectively, in the plane parallel and the developable surface, criteria can be designed to determine the regular developable Bézier patches of four, five, and six degrees.In further studies, discussions should focus on how to construct the developable patches in which the two boundaries are space curves and can be adjusted by some shape parameters.