• Tiada Hasil Ditemukan

TRIGONOMETRIC BEZIER TRANSITION CURVES AND APPLICATIONS

N/A
N/A
Protected

Academic year: 2022

Share "TRIGONOMETRIC BEZIER TRANSITION CURVES AND APPLICATIONS"

Copied!
185
0
0

Tekspenuh

(1)

TRIGONOMETRIC BEZIER TRANSITION CURVES AND APPLICATIONS

MD YUSHALIFY MISRO

UNIVERSITI SAINS MALAYSIA

2017

(2)

TRIGONOMETRIC BEZIER TRANSITION CURVES AND APPLICATIONS

by

MD YUSHALIFY MISRO

Thesis submitted in fulfilment of the requirements for the degree of

Doctor of Philosphy

August 2017

(3)

ACKNOWLEDGEMENT

In the name of Allah, the Most Gracious and Most Merciful. All praises to Allah and peace be upon His beloved prophet MUHAMMAD s.a.w. I am grateful for all the strength and blessings He bestowed me in order for me to complete my research. I was blessed to be backed up by wonderful people throughout this entire research. I would like to express my gratitude to all of them. My deepest gratitude goes to my supervisor, Dr. Ahmad Lutfi Amri Ramli for his guidance, moral injections, and encouragement at every stage of this research, especially during my falls. It has always been an honor to work with him. I am thankful for his commitments and contributions of ideas. I also appreciate the freedom he gave for me to make my own decision regarding this research. My sincere admiration goes to Prof. Jamaludin Md Ali, senior professor in Universiti Sains Malaysia (USM) for his endless encouragement since my degree years until now. His enthusiasm in this research area motivates me to follow his path.

My special thanks goes to Prof Hailiza Kamarulhaili for allowing me to pursue my doctoral studies and providing me a chance to work in this healthy research environ- ment. I would also like to acknowledge the Ministry of Higher Education Malaysia for conferring me a sufficient scholarship throughout this research period. No words can I express how indebted I am to my parents, Hamtiah Md Suadi and my late father, Misro Kasiran who were always there to support me even while battling stage 4 of cancer.

Their endless love, prayers, and priceless time have made me who I am today. Also, thanks to all my siblings and friends who stayed during the thick and thin. Our memo- ries will always be treasured. Lastly, thanks to the non academic staff of the school all for their help and cooperation.

(4)

TABLE OF CONTENTS

Acknowledgement . . . ii

Table of Contents . . . iii

List of Tables . . . vii

List of Figures . . . viii

List of Abbreviations . . . xv

List of Symbols . . . xvi

Abstrak . . . xix

Abstract . . . xx

CHAPTER 1 – INTRODUCTION 1.1 Introduction . . . 1

1.2 Motivation . . . 5

1.3 Objectives . . . 6

1.4 Outline of Thesis . . . 7

1.5 Problem Statements . . . 8

1.6 Scope and Limitations . . . 8

CHAPTER 2 – BACKGROUND AND LITERATURE REVIEW 2.1 Parametric Curves . . . 9

2.1.1 Bézier curve . . . 10

2.1.2 Trigonometric Bézier curve. . . 10

2.2 Spiral Curves . . . 13

(5)

2.3 Design of Highways . . . 15

2.3.1 Design Consistency . . . 16

2.3.2 Design Speed . . . 17

CHAPTER 3 – TRIGONOMETRIC BÉZIER CURVE 3.1 Introduction . . . 20

3.2 Cubic Trigonometric Bézier Curve . . . 21

3.2.1 Cubic Trigonometric Bézier Basis Function . . . 21

3.2.2 Properties of Cubic Trigonometric Bézier Curve . . . 23

3.3 Quintic Trigonometric Bézier Curve . . . 24

3.3.1 Quintic Trigonometric Bézier Basis Function . . . 24

3.3.2 Properties of Quintic Trigonometric Bézier Curve . . . 27

3.4 Curvature Distribution of Quintic Trigonometric Bézier Curve . . . 29

3.5 Composition of Two Quintic Trigonometric Bézier Curves . . . 30

3.5.1 Shape Control of the Quintic Trigonometric Bézier Curve . . . 37

3.6 Summary . . . 38

CHAPTER 4 – CUBIC TRIGONOMETRIC BÉZIER TRANSITION CURVE 4.1 Introduction . . . 40

4.1.1 Notation and Convention . . . 42

4.2 Straight Line to a Circle Spiral Curve . . . 43

4.3 Straight Line to Straight Line Transition Curve. . . 51

4.4 S-shaped Transition Curve Composed by Two Cubic Trigonometric Bézier Spiral . . . 54

4.5 S-shaped Transition Curve with Pair of Shape Parameters . . . 56

(6)

4.5.1 S-shaped Transition Curve with Same Radii of Circles . . . 57

4.5.2 S-shaped Transition Curve with Different Radii of Circles . . . 61

4.6 C-shaped Transition Curve Composed by Two Cubic Trigonometric Bézier Spiral . . . 65

4.7 C-shaped Transition Curve with Pair of Shape Parameters. . . 66

4.7.1 C-shaped Transition Curve with Same Radii of Circles. . . 66

4.7.2 C-shaped Transition Curve with Different Radii of Circles . . . 71

4.8 Circle in Circle with C-shaped Curve . . . 74

4.9 Summary . . . 76

CHAPTER 5 – QUINTIC TRIGONOMETRIC BÉZIER TRANSITION CURVE 5.1 Introduction . . . 79

5.1.1 Notation and Convention . . . 82

5.2 Straight Line to a Circle Spiral Curve . . . 83

5.3 Straight Line to Straight Line Transition Curve. . . 85

5.4 S-shaped Transition Curve Composed by Two Quintic Trigonometric Bézier Spirals . . . 88

5.5 S-shape Transition Curve with Pair of Shape Parameters . . . 90

5.5.1 S-shaped Transition Curve with Same Radii of Circles . . . 94

5.5.2 S-shaped Transition Curve with Different Radii of Circles . . . 96

5.6 C-shaped Transition Curve Composed by Two Quintic Trigonometric Bézier Spirals . . . 99

5.7 C-shaped Transition Curve with Shape Parameters . . . 102

5.7.1 C-shaped Transition Curve with Same Radii of Circles. . . 105

5.7.2 C-shaped Transition Curve with Different Radii of Circles . . . 108

5.8 Circle in Circle with C-shaped Curve . . . 111

(7)

5.9 Summary . . . 115

CHAPTER 6 – TRANSITION CURVES OF HIGHWAY DESIGN 6.1 Introduction . . . 116

6.2 Classical Theory . . . 121

6.2.1 Radial Forces . . . 121

6.2.2 Design Speed . . . 122

6.3 Research Area . . . 124

6.4 Speed Estimation . . . 126

6.4.1 Heuristic Method . . . 126

6.4.2 Cubic Bézier Curve . . . 129

6.4.3 Quintic Bézier Curve . . . 129

6.4.4 Quintic Trigonometric Bézier Curve . . . 133

6.5 Summary . . . 148

CHAPTER 7 – CONCLUSION 7.1 Introduction . . . 150

7.2 Research Conclusion . . . 150

7.3 Future Work and Recommendation . . . 154

REFERENCES. . . 155 APPENDICES

LIST OF PUBLICATIONS

(8)

LIST OF TABLES

Page Table 3.1 Curvature value comparison for blue and red curve 36

Table 3.2 Curvature value for second curve 36

Table 6.1 Speed Approximation Based on Manual Attempt 128 Table 6.2 Speed approximation on 10 control points of Tun Sardon

Road

131

Table 6.3 Curvature values, Velocity Designs, Color Gradients for 1st curve of planar quintic trigonometric Bézier

134

Table 6.4 Curvature values, Velocity Designs, Color Gradients for 2nd curve of planar quintic trigonometric Bézier

135

Table 6.5 Curvature values, Velocity Designs, Color Gradients for 3rd curve of planar quintic trigonometric Bézier

136

Table 6.6 Curvature values, Velocity Designs, Color Gradients for 4th curve of planar quintic trigonometric Bézier

137

Table 6.7 Curvature values, Velocity Designs, Color Gradients for 5th curve of planar quintic trigonometric Bézier

138

Table 6.8 Curvature values, Velocity Designs, Color Gradients for 6th curve of planar quintic trigonometric Bézier

139

Table 6.9 Curvature values, Velocity Designs, Color Gradients for 7th curve of planar quintic trigonometric Bézier

140

Table 6.10 Curvature values, Velocity Designs, Color Gradients for 8th curve of planar quintic trigonometric Bézier

141

Table 6.11 Curvature values, Velocity Designs, Color Gradients for 9th curve of planar quintic trigonometric Bézier

142

Table 6.12 Comparison of speed estimation of quintic Bézier curves and quintic trigonometric Bézier curves

146

(9)

LIST OF FIGURES

Page Figure 3.1 Cubic trigonometric Bézier basis function with different val-

ues of shape parameters

22

Figure 3.2 Quintic trigonometric Bézier basis function with different values of shape parameters

27

Figure 3.3 A curve withp=q(left) and its curvature distribution (right) 29 Figure 3.4 A curve with fixed p(left) and its curvature distribution

(right)

29

Figure 3.5 A curve with fixedq(left) and its curvature distribution (right)

29

Figure 3.6 Composite quintic trigonometric Bézier curve 35 Figure 3.7 Clover design using quintic trigonometric Bézier composite

curve

38

Figure 3.7(a) Red clover 38

Figure 3.7(b) Blue clover 38

Figure 3.7(c) Brown clover 38

Figure 3.7(d) Purple clover 38

Figure 3.7(e) Green clover 38

Figure 4.1 J-shaped transition curve 44

Figure 4.2 Curvature profile of J-shaped transition curve (left) and its curvature derivative (right)

44

Figure 4.3 J-shaped trigonometric spiral curve 49

Figure 4.4 J-shaped trigonometric spiral curves with shape parameters p=1,q=0 (red-dashed curve) and p=−1,q=0 (blue- dotted curve)

49

Figure 4.5 Curvature profile with shape parameters p=1,q=0 (red- dashed curve) and p=−1,q=0 (blue-dotted curve)

50

(10)

Figure 4.6 J-shaped trigonometric spiral curve with different values of

θ 50

Figure 4.7 Curvature profile of J-shaped spiral with differentθ values 51 Figure 4.8 Straight line to straight line spiral curve 53 Figure 4.9 Straight line to straight line spiral trigonometric Bézier curve

with the shape parameters p=2,q=0 (red-dashed) and p=1,q=0 (blue-dotted)

53

Figure 4.10 S-shaped curve composed by two cubic trigonometric Bézier spiral

55

Figure 4.11 An S-shaped transition curve with same radii (left) and its curvature profile (right)

58

Figure 4.12 S-shaped transition curves of same radii with shape parame- tersp=1,q=1 for black-lined,p=0,q=0 for red-dashed and p=−1,q=−1 for blue-dotted (left) and its curvature profile (right)

59

Figure 4.13 Curvature derivatives of S-shaped transition curves of same radii

60

Figure 4.14 S-shaped transition curves of same radii with shape param- eters p=0,q=1 for red-lined and p=1,q=0 for blue- dotted (left) and its curvature profile (right)

61

Figure 4.15 An S-shaped transition curve of different radii (left) and its curvature profile (right)

61

Figure 4.16 S-shaped transition curves in different radii with shape pa- rametersp=1,q=1 for black-lined curve,p=0,q=0 for red-dashed curve and p=1,q=1 for blue-dotted curve (left) and their curvature profiles (right)

62

Figure 4.17 Curvature derivatives of S-shaped transition curves of differ- ent radii

63

Figure 4.18 S-shaped transition curves with pairs of shape parameters p=1,q=0 for red curve and p=0,q=1 for blue-dotted curve (left) and their curvature profiles (right)

64

Figure 4.19 C-shaped curve composed by two cubic trigonometric Bézier spirals

65

(11)

Figure 4.20 A C-shaped transition curve of same radii (left) and its cur- vature profile (right)

67

Figure 4.21 C-shaped transition curves with shape parametersp=1,q= 1 for blue-dotted curve, p=0,q=0 for red-dashed curve andp=1,q=1 for black line curve (left) and their cur- vature profiles (right)

69

Figure 4.22 Curvature derivatives of C-shaped transition curves of same radii

70

Figure 4.23 C-shaped transition curves with shape parametersp=1,q= 0 for blue-dotted curve and p=0,q=1 for red-lined curve (left) and its curvature profile (right)

71

Figure 4.24 A C-shaped transition curve of different radii (left) and its curvature profile (right)

71

Figure 4.25 C-shaped transition curves with shape parametersp=1,q= 1 for blue-dotted, p = 0,q = 0 for red-dashed, and p =

1,q=1 for black-lined curve (left) and their curvature profiles (right)

72

Figure 4.26 Curvature derivatives of S-shaped transition curves of differ- ent radii

73

Figure 4.27 C-shaped transition curves with shape parametersp=0,q= 1 for red-dotted andp=−1,q=1 for blue-dashed (left) and its curvature profile (right)

73

Figure 4.28 A C-shaped trigonometric Bézier spiral curve (left) and its curvature profile (right)

74

Figure 4.29 A C-shaped trigonometric Bézier transition curve (left) and its curvature profile (right)

75

Figure 4.30 A C-shaped trigonometric Bézier transition curve with same pairs of shape parameters (left) and their curvature profiles (right)

76

Figure 5.1 A J-shaped spiral quintic trigonometric Bézier curve 84 Figure 5.2 Curvature profile of J-shaped spiral quintic trigonometric

Bézier curve (left) and its curvature derivative (right)

85

Figure 5.3 V-shaped composed of two PH quintic trigonometric Bézier spiral curves

87

(12)

Figure 5.4 V-shaped composed of two spiral curves without control points

87

Figure 5.5 S-shaped spiral curve with(m0,m1) = (3.293,1.293) 89 Figure 5.6 S-shaped spiral curve withSred= (m0,m1) = (2.707,1.707),

Sblue = (m0,m1) = (1.8,1) and Sblack = (m0,m1) = (3.293,1.293)

90

Figure 5.7 Quintic trigonometric Bézier S-shaped transition curve on same radii with shape parameters p=1,q=1

92

Figure 5.8 Curvature profile of S-shaped transition curve on same radii (left) and its curvature derivative (right)

92

Figure 5.9 Quintic trigonometric Bézier S-shaped transition curves on two different radii of circles

93

Figure 5.10 Curvature profile of quintic trigonometric Bézier S-shaped transition curve on two different radii (left) and its curvature derivative (right)

93

Figure 5.11 Quintic trigonometric S-shaped transition curves using dif- ferent parameterm

94

Figure 5.12 Curvature profile of quintic trigonometric S-shaped transi- tion curves using different parameterm(left) and their cur- vature derivatives (right)

95

Figure 5.13 Quintic trigonometric Bézier S-shaped transition curves with u0=2.5 (blue-dashed curve),u0=2.0 (black-dotted curve) andu0=1.5 (red-lined curve)

95

Figure 5.14 Curvature profile of quintic trigonometric S-shaped transi- tion curves using different parameteru0(left) and their cur- vature derivatives (right)

96

Figure 5.15 Quintic trigonometric Bézier S-shaped transition curves with m=2.2 (blue-dashed curve), m=1.7 (black-dotted curve) andm=1.2 (red-lined curve)

97

Figure 5.16 Curvature profile of PH quintic trigonometric S-shaped tran- sition curves using different parameterm(left) and their cur- vature derivative (right)

97

Figure 5.17 PH Quintic trigonometric S-shaped transition curve using different parameteru0

98

(13)

Figure 5.18 Curvature profile of PH quintic trigonometric S-shaped tran- sition curve using different parameteru0(left) and its curva- ture derivative (right)

99

Figure 5.19 C-shaped spiral curve withC= (m0,m1) = (2.707,1.707) 100 Figure 5.20 C-shaped spiral curve withCred= (m0,m1) = (2.707,1.707),

Cblue = (m0,m1) = (1.8,1) and Cblack = (m0,m1) = (3.293,1.293)

101

Figure 5.21 C-shaped transition curve on same radii 103

Figure 5.22 Curvature profile of C-shaped transition curve on same radii (left) and its curvature derivative (right)

104

Figure 5.23 Quintic trigonometric Bézier C-shaped transition curve with shape parameters p=1,q=1

104

Figure 5.24 Curvature profiles of quintic trigonometric Bézier C-shaped transition curve and its curvature derivative (right)

105

Figure 5.25 Quintic trigonometric Bézier C-shaped transition curve with m=1.8 (blue-dashed curve), m=1.5 (black-dotted curve) andm=1.2 (red-lined curve)

106

Figure 5.26 Curvature profiles of quintic trigonometric C-shaped transi- tion curves using different parameterm(left) and their cur- vature derivatives (right)

107

Figure 5.27 Quintic trigonometric C-shaped transition curve using dif-

ferent values ofθ 107

Figure 5.28 Curvature profiles of quintic trigonometric C-shaped transi- tion curves using different parameterθ (left) and their cur- vature derivative (right)

108

Figure 5.29 Quintic trigonometric C-shaped transition curves using dif- ferent values ofm

108

Figure 5.30 Curvature profiles of quintic trigonometric C-shaped transi- tion curves using different parametersm(left) and their cur- vature derivatives (right)

109

Figure 5.31 Quintic trigonometric Bézier C-shaped transition curve with θ =55 degree (blue-dashed curve), θ =50 degree (black- dotted curve) andθ =45 degree (red-lined curve)

109

(14)

Figure 5.32 Curvature profile of quintic trigonometric C-shaped transi- tion curves using different parameterθ (left) and their cur- vature derivatives (right)

110

Figure 5.33 Circle in circle spiral curve 111

Figure 5.34 Curvature profile of circle in circle spiral curve (left) and its curvature derivative (right)

112

Figure 5.35 Circle in circle transition curve 113

Figure 5.36 Curvature profile of circle in circle quintic trigonometric Bézier transition curve (left) and its curvature derivative (right)

114

Figure 6.1 A wholly transition curve 118

Figure 6.2 A compound transition curve 118

Figure 6.3 Two types of clothoids (Retrieved from Song (2006)) 119 Figure 6.4 A reverse curve (Retrieved from Song (2006)) 120 Figure 6.5 A circle in circle curve (Retrieved from Song (2006)) 120

Figure 6.6 A transitional section of curve 122

Figure 6.7 Superevelation in highway design 124

Figure 6.8 Topography Tun Sardon Road (Retrieved from Veloviewer) 124 Figure 6.9 3D Tun Sardon Road Elevation (Retrieved from Veloviewer) 125 Figure 6.10 Tun Sardon Road Elevation (Retrieved from Veloviewer) 126

Figure 6.11 Research Area 127

Figure 6.12 Quintic Bézier Curve Mapping on Road 132

Figure 6.13 1st piecewise of quintic trigonometric Bézier curve 134 Figure 6.14 2nd piecewise of quintic trigonometric Bézier curve 135 Figure 6.15 3rd piecewise of quintic trigonometric Bézier curve 136 Figure 6.16 4th piecewise of quintic trigonometric Bézier curve 137 Figure 6.17 5th piecewise of quintic trigonometric Bézier curve 138

(15)

Figure 6.18 6th piecewise of quintic trigonometric Bézier curve 139 Figure 6.19 7th piecewise of quintic trigonometric Bézier curve 140 Figure 6.20 8th piecewise of quintic trigonometric Bézier curve 141 Figure 6.21 9th piecewise of quintic trigonometric Bézier curve 142 Figure 6.22 Interpolation of quintic trigonometric Bézier curve with

color gradient schemes corresponding to maximum speed es- timation

143

Figure 6.23 Interpolation of quintic trigonometric Bézier curve for all curve segments with shape parameters p=4 andq=1 for the first curve segment and alternately for consecutive curves

144

Figure 6.24 Interpolation of quintic trigonometric Bézier curve for all curve segments with shape parameters isp=1 andq=−4 for the first curve segment and alternately for consecutive curves

145

Figure 1 Appendix A1 162

Figure 2 Appendix A2 163

(16)

LIST OF ABBREVIATIONS

CAGD Computer Aided Geometric Design CAD Computer Aided Design

CAM Computer Aided Manufacturing USM Universiti Sains Malaysia PH Pythagorean Hodograph JKR Jabatan Kerja Raya

CNC Computer Numerical Control

(17)

LIST OF SYMBOLS

× cross product

· dot product

square root

a vector a parallel

θ angle in degree μ angle in degree α angle in degree β angle in degree

Pi control points of P Qi control points of Q Bi control points of B p shape parameter left q shape parameter right s free parameter

m free parameter n free parameter r radius of circle t variable

(18)

t0 knot point t1 knot point

T0 beginning unit tangent T1 ending unit tangent

C0 position parametric continuity C1 tangent parametric continuity C2 curvature parametric continuity G0 position geometric continuity G1 tangent geometric continuity G2 curvature geometric continuity

fi cubic trigonometric Bezier basis function

gi quintic trigonometric Bezier basis function ω0 first circle

ω1 second circle

∞ infinity

∑ summation λ lambda π pi

κ(t) curvature

κ(t) curvature derivative

υ(t) numerator of curvature derivative

(19)

u(t) quadratic function with coefficientui v(t) quadratic function with coefficientvi z(t) parametric curve

h=P0P1 magnitude fromP1toP0 k=P3P2 magnitude fromP3toP2 Cm cosine basis function

(20)

LENGKUNG PERALIHAN BEZIER TRIGONOMETRI DAN APLIKASI

ABSTRAK

Tesis ini menghuraikan pembinaan lengkung Bézier kuintik trigonometri dengan dua parameter bentuk melalui perbincangan fungsi asas, sifat, dan komposisi dua se- gmen lengkung. Lima templat lengkung peralihan yang boleh lentur dijana menggu- nakan lengkung trigonometri kubik dan kuintik dengan kehadiran parameter bentuk.

Lima templat lengkung peralihan diperoleh dengan menghubungkan garis lurus kepa- da bulatan, bulatan kepada bulatan dengan peralihan C, bulatan kepada bulatan dengan peralihan S, garis lurus kepada garis lurus, dan bulatan kepada bulatan yang mem- punyai bulatan kecil di dalam bulatan besar. Dua daripada templat tersebut membin- cangkan lengkung peralihan berbentuk S dan C menggunakan satu lengkung peralihan atau yang terbentuk daripada dua lingkaran Bézier. Tesis ini turut membincangkan aplikasi dalam reka bentuk lebuh raya dengan membuat anggaran halaju jalan raya menggunakan maklumat kelengkungan. Perbandingan antara lengkung Bézier kuintik dan lengkung Bézier kuintik trigonometri terhadap set data sebenar turut dibincangkan.

(21)

TRIGONOMETRIC BEZIER TRANSITION CURVES AND APPLICATIONS

ABSTRACT

This thesis describes the development of a new quintic trigonometric Bézier curve with two shape parameters by discussing its basis function, properties, and composition of two curve segments in the early part of this thesis. Five flexible templates of transi- tion curves are generated using both cubic and a newly proposed quintic trigonometric Bézier curves with the presence of shape parameters. The five templates of transition curves are achieved by joining a straight line to a circle, circle to circle with a broken back C transition, circle to circle with an S transition, straight line to straight line, and circle to circle where small circle lies inside the big circle. Two of the templates will discuss on S-shaped and C-shaped transition curve using single transition curve and two Bézier spiral segments. This thesis also includes works on the application in high- way design by approximating maximum speed estimation using curvature information.

Comparisons between quintic Bézier curve and quintic trigonometric Bézier curve on real data sets are also presented.

(22)

CHAPTER 1 INTRODUCTION

1.1 Introduction

Computers have become an indispensable tool in modelling and simulation. So- phisticated geometric shapes and physical objects are modelled in the domain of com- plex physical environments. As computational power drastically increased, users and applications also demand for a simpler method with less computational time. Ap- proaches for object modelling and designing involves non-physical methods. Non- physical method can be defined as individual or groups of control point are manually adjusted for shape editing and designing (Gibson and Mirtich, 1997).

Computer Aided Geometric Design (CAGD) is a discipline dealing with compu- tational geometry. Some applications of CAGD are construction and manipulation of free-form curves and surfaces (Farin, 2002a). The major breakthrough in CAGD were the theory of Bézier curves and surfaces, as well as B-spline method. B-spline actually originated from the word ‘spline’. Spline was not actually a mathematical term but was previously a draftsmen’s tool that consists of a flexible strip of metal or wood and several heavy pieces which hold or anchored the flexible strips for drawing ship (Farin, 2002b).

In 1959, De Casteljau had already developed curves similar to Bézier indepen- dently, but his work were never published. Therefore, when Pierre Étienne Bézier from a French car manufacturer, Renault, presented a parametric polynomial curve

(23)

representation technique in 1962, the whole theory of polynomial curves and surfaces in Bernstein form is now using Bézier’s name (Farin, 2002b). Bézier and B-spline are two important designing techniques with roots in the late 1950’s. Both of these techniques are used for interactive curve and surface design (Faux and Pratt, 1979;

Hoschek et al., 1993; Mathews, 1986).

Bézier and B-spline are mostly used for design purposes that will be incorporated and integrated into Computer Aided Design (CAD) or Computer Aided Manufacturing (CAM). Bézier is mostly useful in CAD/CAM to construct and manipulate free-form curves and surfaces to fabricate products such as automotive, aerospace and shipbuild- ing. Desired shape features play a pivotal role in constructing product of curves and surfaces. In some cases, it is difficult to distinguish between Bézier and B-spline be- cause any Bézier curve of arbitrary degree can be converted into B-spline and any B-spline can be transformed into one or more Bézier curves(Prautzsch et al., 2013;

Romani and Sabin, 2004).

In terms of providing a local control on the shape of the curve, B-spline is a bet- ter choice compared to Bézier, where Bézier offers a global shape control. B-spline curve is more flexible and a better piecewise than Bézier curves. Therefore, B-spline is a powerful generalization of Bézier curves (Racine, 2014). Although, with all the advantages of B-spline curves, it does have a shortcoming. B-spline are more com- plex than Bézier curve (Parekh et al., 2006) and thus industrial designers would prefer to use Bézier in terms of practicality. Several problems in applied mathematics and engineering would require splines to be one of their solutions.

(24)

This includes CAGD (Hoschek et al., 1993), function approximation (Nuernberger, 1989), optimal control (Schumaker, 2007), data fitting (Dierckx, 1995) and many more. Transition curve is useful in various fields, especially CAGD. It is also an important part in generating round corners or smooth transition between two circles.

Properties of transition curves which have gradual increasing or decreasing curvature between circles plays an important role in designing highways and railways where a transition path is needed to lead into one part of the road to another such as from one straight line to circle, circle to circle, and circle to straight line.

In road and rail design, passenger’s comfort and safety is another important as- pect to be considered. Transition curve may ease up the transition, prevent accidents caused by sharp changes in direction, and reduces shaking of vehicles in operation.

It also contributes in smoother, stable and safer transition. Baass (1984) proposed a design technique by implementing single ordinary Bézier function of continuity with five template of transition curve.

• Type 1. Straight line to circle with a J-shaped transition curve (Habib and Sakai, 2005b; Walton et al., 2003)

• Type 2. Circle to circle with a broken back C-shaped transition (Habib and Sakai, 2007a; Yahaya et al., 2008)

• Type 3. Circle to circle with an S-shaped transition (Habib and Sakai, 2007a;

Walton and Meek, 1999; Yahaya et al., 2009)

• Type 4. Straight line to straight line (Habib and Sakai, 2005b; Walton et al., 2003)

(25)

• Type 5. Circle in circle with C-shaped transition (Walton and Meek, 1996; Ya- haya et al., 2010)

Previous works that discussed on the construction of S-shaped transition curve had been done using cubic Bézier (Habib and Sakai, 2003; Walton and Meek, 1999) and cubic Bézier-like (Yahaya et al., 2008) and cubic trigonometric Bézier (Abbas et al., 2011). Walton and Meek (1999) constructed C-shaped transition curve using cubic Bézier curve whereas Habib and Sakai (2003) constructed transition curve but using a different approach. Other authors who had constructed C-shaped transition curve used a cubic Bézier-like curve such as in Yahaya et al. (2008). They used three control points which produced two segments.

Throughout this thesis, transition curves will be researched extensively using cu- bic trigonometric Bézier curves and quintic trigonometric Bézier curves on all five templates of transition curves. Cubic Bézier curve is the lowest degree of polynomial which can be used as transition curve to construct or design in CAD/CAM, whereas quintic curve is the lowest degree of Pythagorean Hodograph curve that may have an inflection point, as required for S-shaped transition curve (Habib and Sakai, 2003).

Single cubic curve can be used as transition curve between two circles with the guarantee that S-shaped transition curve does not have interior curvature extremum and C-shaped transition curve does not have inflection point or singularities. Even though cubic curve is smoother and simpler, it is not always helpful since they may have unwanted inflection points and singularities (Habib and Sakai, 2007b). A cubic segment has the following undesirable features which are; its arc-length is the integral

(26)

part of the square root of polynomial of its parameter, its offset is neither polynomial nor a rational algebraic function of its parameter, and it may have more curvature extrema than necessary.

Pythagorean Hodograph curve does not suffer from the first two undesirable fea- tures that were mentioned before. Pythagorean originally came from the word Pythago- ras, where is often referred to a pure mathematician (Ratner, 2009). Pythagoras con- tributed to studies of areas, polygons, proportions, and also perfect numbers. Pythagorean Hodograph curves were developed by Farouki (2008). Hodograph is defined as deriva- tives ofz(t)is the square of a complex polynomial. PH quintic curve is the isomorphic shape to general cubic curve.

1.2 Motivation

Bézier technique is a powerful tool that has been frequently used to generate free form curves and surfaces. It is also used for approximating, interpolating, preserving shape, designing shape of curves. Trigonometric Bézier acts as an alternative useful tools rather than tedious conventional rational Bézier which both are used as mediums to control desired curves. Trigonometric offers more flexibility by executing the results in simpler ways than rational Bézier.

Trigonometric Bézier curve is a useful function that inherits nice properties similar to rational Bézier, except the shape parameters are better in terms of less complexity, less time consuming, and less tedious works compared to rational Bézier. Transition curves are being constructed using Bézier, Bézier-like, Log aesthetic curves and also trigonometric Bézier. In this thesis, further discussion will be done on five templates

(27)

of transition curve using cubic and quintic trigonometric Bézier curves. The reason for the use of cubic and quintic trigonometric Bézier curves is due to trigonometric Bézier curves itself needs to satisfy some geometrical properties such as partition of unity and symmetry, therefore in order to have a more flexible curve, we need to have two shape parameters in its basis function. So when quartic trigonometric Bézier curve is used, it can only have five control points and its very difficult to control the curve because the properties of Bézier curves require it to be in symmetrical form in order to satisfy the conditions of transition curves. Some numerical analysis will also be provided by applying the function on real data sets.

1.3 Objectives

This thesis deals with the development of a new function called quintic trigonomet- ric Bézier curve with two shape parameters and also the construction of five templates of transition curve by using cubic and quintic trigonometric Bézier curve as transition curves. The objectives of this research are:

1. To propose a new generalized quintic trigonometric Bézier curve with two shape parameters by discussing its basis function, properties and composition of two curve segments.

2. To generate five flexible templates of transition curves using cubic trigonometric Bézier curve and a newly proposed quintic trigonometric Bézier curve with the presence of two shape parameters.

3. Applications in highway design by approximating maximum speed estimation

(28)

1.4 Outline of Thesis

The description of the chapters included in this thesis are as follows:

Chapter 2 presents an overview of parametric curves including Bézier and trigono- metric Bézier. In this section, parametric and geometric continuity are discussed with the important use of curvatures. Some notations and terminology that are essential to understand the work, detailed literature reviews and related works are presented.

In Chapter 3, the cubic trigonometric Bézier function and its properties are dis- cussed. A new function of quintic trigonometric Bézier function with its properties respectively will also be proposed and demonstrated in this chapter.

Chapter 4 deals with the construction of five templates of transition curve using cubic trigonometric Bézier curve. We will also discuss further on two templates which are S-shaped and C-shaped that can be represented by two piece curves or single curve, where one is classified as spiral while another one as transition curve.

In Chapter 5, the construction of five templates of design techniques using quintic trigonometric Bézier curves will be presented. Quintic trigonometric Bézier curves are isomorphic shape with cubic trigonometric Bézier curve that possessed almost similar features to cubic curve, but with more attractive properties.

Application on real data set using quintic Bézier and quintic trigonometric Bézier along their comparisons will be demonstrated in Chapter 6. Some classical theory with some different approaches to calculate maximum speed estimation are discussed.

(29)

Finally, the final conclusion followed by some suggestions for future researches will be provided in Chapter 7.

1.5 Problem Statements

1. To obtain a flexible curve by using trigonometric Bézier without having to use weightage in its denominator in order to alter the curve.

2. Previous author (Abbas et al., 2010), use cubic trigonometric Bézier curve to construct only one template of transition curve which is Type III without fur- ther discussion on different pairs of shape parameter or different radii of circles.

Therefore, this research attempt to complete the gaps left by previous author.

3. The quintic trigonometric Bézier curve have never been applied on Pythagorean Hodograph function using the same approaches and the use of this function in constructing five templates of transition curve were never discussed in any pre- vious researches.

1.6 Scope and Limitations

1. This research is focused on two-dimensional space as a basic guideline to form a strong foundation before it can be extended into a three-dimensional space of transition curve using both cubic and quintic trigonometric Bézier curves.

2. For the purpose of this research and to answer research objectives, the applica- tion of this transition curve is focused only on highway designing and approxi- mating maximum speed estimation.

(30)

CHAPTER 2

BACKGROUND AND LITERATURE REVIEW

2.1 Parametric Curves

Based on the book A History of Curves and Surfaces by Farin (2002b), curves were originally constructed by draftsmen. Most of the time, they use circles as their basic tool. French curve is a set of templates that function as common tools when they had to draw some of the ’free forms’. Meanwhile, wooden curves were designed carefully and consist of pieces of conics and spirals. Other mechanical tools such as spline was also used. This spline is a flexible strips of wood that was held in place and shape by metal weights known as ducks. The mechanical counterpart to a mechanical spline is a spline curve, one of the most fundamental parametric curve forms.

Later in the late 1800s, Serret/Frenet documented about the differential geometry of parametric curves. J. Ferguson and D. MacLaren constructed curves for the design of wings for US aircraft company Boeing in the middle of 1950s. They formed combined curves by piecing cubic space curves together which were overall twice differentiable.

Ferguson also used the cubic Hermite form which defines a cubic in terms of two endpoint and two endpoint derivatives. The most fundamental parametric curve forms are the Bézier curves (Farin, 2002b).

(31)

2.1.1 Bézier curve

Pierre Étienne Bézier’s initial idea was to represent a basic curve as the intersec- tion of two elliptical cylinders. Later, Bézier moved to polynomial formulations of this inital concept and also extended it to higher degrees. Bézier’s work were widely published and caught A.R Forrest’s attention. His works on Bézier curves helped to popularize and make Bézier curves influential. He then realized that Bézier curves can be expressed in terms of Bernstein polynomial. Several Bernstein basis properties con- trol the behaviour of Bézier curves like symmetry, recursion, non-negativity, partition of unity, variation diminishing property, derivatives and so (Farin, 2002b).

Han et al. (2008) have achieved similar results by extending the classical Bernstein basis function that inherits most properties of Bézier curve. The basis function was introduced with n adjustable shape parameters λ to control the shape of the curve without changing the control polygon. From the view point of basis function, Bézier curves were also investigated using trigonometric polynomials.

2.1.2 Trigonometric Bézier curve

Bézier curves had been investigated and expanded using trigonometric polynomials (Fitter et al., 2014). Parametric representation through trigonometric polynomial had been emphasized for its benefit in offering local control over the shape through shape parameters and its continuity conditions are being examined consistently under various degrees like quadratic trigonometric Bézier curve (Bashir et al., 2012), cubic trigono- metric Bézier curve (Han et al., 2009), quartic trigonometric Bézier curve (Dube and Sharma, 2013), quasi-quintic trigonometric Bézier curves (Bashir et al., 2013).

(32)

A recent development in CAGD is the introduction of trigonometric spline based on trigonometric polynomials and a blending of algebraic and trigonometric polyno- mials. Schoenberg (1964) was the pioneer in presenting trigonometric spline whose components were piecewise trigonometric polynomials. This introduction of trigono- metric splines is a new beginning for researchers. Expanding previous works, Lyche and Winther (1979) established a three term recurrence relation for trigonometric B- splines of arbitrary order, analogous to that for polynomial B-splines. The method to differentiate trigonometric B-splines and a trigonometric analog of Marsden’s identity were also given (Hoschek et al., 1993).

The trigonometric B-splines studied by these authors did not sum up to unity as op- posed to polynomial B-splines and thus did not lie in the convex hull of control points.

To get the better of this deficiency, Koch et al. (1995) introduced control curves for trigonometric splines instead of control polygon and proved two important properties, which are convex hull property and variation diminishing property for these splines.

The study also includes knot insertion algorithm. It was concluded that by adding extra knots in trigonometric splines, it will lead the control curves in converging the splines.

Certain CAGD application involve curves with such complex shapes that cannot be represented by single curve segments (Buss, 2003). One solution is degree rais- ing that would add flexibility to the curve. However, the processing effort for curve manipulation is considerably increased. In addition, a curve with higher degrees may induce some numerical noise in computations. To cope with this, a curve with complex shape is often represented by stitching together a number of low degree curves (called segments) while meeting certain specifications of continuity at the joints.

(33)

Such curve is termed as a piecewise curve. A curve that is made of several Bézier curves is called a composite Bézier curve (Farin, 2002a). In order to create a good composite Bézier curve, smoothness is a very important element. Smoothness is a primitive requirement in generating curves. The order of smoothness is application dependent. For example, in architectural design, it is sufficient that the curve segments only bear position continuity. It is not enough to just match the curves exclusively at their adjoining points.

A smooth joint is possible only if some extra conditions are imposed at those points.

Thus, parametric and geometric continuity conditions are used in order to get a smooth joint of two or more pieces of a single curve. Parametric continuity is accomplished by matching the parametric derivatives of adjacent segments of curves at their common boundary (Farin, 2002a). The first nth-order parametric continuity where the curve segments meet are known asCn. The three most commonly used levels of parametric continuity are:

C0stands for position continuity. The adjoining curve segments are merely con- nected at their respective endpoints.

C1suggests that the first order of parametric derivatives of successive segments at the shared points to be equal. Therefore, the two curve segments have the same tangent at the common point, both in magnitude and direction.

C2 means that both first order and second order parametric derivatives of two curve segments are same at the intersection point. Thus,C2denotes the continu- ity of acceleration vectors.

(34)

A relaxed form of nth-order parametric continuity, Cn has been developed and dubbed nth-order geometric continuity or Gn (Barsky and DeRose, 1989). Geomet- ric continuity depends upon the geometry of the curve (Meek and Walton, 2002).

G0is equivalent toC0.

G1 requires that the two curve segments and their unit tangent vectors are con- tinuous at the joint.

G2means that the two adjacent curves areG1and in addition their signed curva- tures at the joining point are continuous.

Curvature formulas for parametrically defined curves and surfaces are well-known both in the classical literature on Differential Geometry (Spivak, 1975; Struik, 2012) and in the contemporary literature on Geometric Modeling (Farin, 2002a; Hoschek et al., 1993). For planar curves, curvature has several equivalent definitions by Gold- man (2005). The curvature of a curve is the measure of its deviation from a straight line in a neighborhood of a given point, and the curvature becomes greater as this deviation becomes greater.

2.2 Spiral Curves

Spiral curves are generally used to provide a gradual change in curvature from a straight section of road to a curved section (Sipes and Sipes, 2013). Spiral was first used in the late 1800s, and its use peaked in the design of the parkways of the 1930s (Myers, 2001). The use of spirals was first documented in the late 1600s in

"Sino Loria," by James Bernoulli (Higgins, 1922). Spirals were rediscovered in 1874

(35)

by Cornu (Kimia et al., 2003) and used in optics and later spirals began to replace parabolic curves in easing transitions for railroads.

A transition curve, or easement curve, as it is sometimes called, is a curve of vary- ing radius used to connect circular curves with tangents for the purpose of avoiding the shock and disagreeable lurch of trains, due to the instant change of direction and also to the sudden change from level to inclined track (Talbot, 1904). The term spiral is interchangeable with easement or transition curve. Spirals are curves used to transition between a circular curve with a specific radius and degree of curvature and a straight tangent. Spiral provides a gradual transition from moving in a straight line to moving in a curve around a point.

A spiral curve is a geometric feature that can be added to a regular circular curve which the radius and sharpness of a spiral curve are increased uniformly along its length (Myers, 2001). Spiral is the most important characteristic in the development of railroad and highway by allowing a safe transition from straight to curved section of track. The use of a spiral is to make roads or tracks to follow the same form that the vehicle naturally takes. The length and degree of curvature of a spiral curve are based on the anticipated speed of traffic and the sharpness of the circular curve that the spiral must meet. In cars, you don’t go directly from going straight to fully turning. There is a transition area where you slowly turn the steering wheel.

Lateral acceleration is slowly increased as the spiral is entered, or it is slowly de- creased as the spiral is exited. Spiral curve enable railroad vehicles to proceed into simple curve without derailing. There are a few advantages of spiral curves. Firstly,

(36)

it provides a natural path such that the lateral force increase and decrease gradually as the vehicle enters and leaves the circular curve. Secondly, the length of the transition curve provides a suitable location for the superelevation runoff. Next, the spiral curve facilitates the transition in width where the travelled way is widened on circular curve.

Lastly, the appearance of the highway is enhanced by the application of spiral curves.

2.3 Design of Highways

Highway is an important aspect when it comes to transportation. Transport plays an important role in our daily activities, such as for economic purposes where raw ma- terials, fuels, food, and manufactured goods are moved from where there are surpluses to where there are shortages. There are various modes of transportation such as air, rail, road, water, cable, pipeline, and space. Basically, it operates using three main ways which are land, water, and air. There are certain factors that needed to be con- sidered when choosing a transport, such as speed, comfort, cost, flexibility, regularity, and safety. Considering the speed, air transports are the quickest mode, but its also the costliest of them all.

A research on safety criteria for highways’ curve design had been done (Anderson et al., 1999; Lamm et al., 1995). The research was conducted to evaluate the validity of design criteria for horizontal highway curves. The evaluation was specifically con- cerned with the design equation, assumed levels of tire-pavement side friction capabil- ity, safe side friction factors, maximum degree of curvature, maximum superelevation, and design factors of safety. Besides that, Ivan et al. (2009) had done their research on designing roads that guide drivers to choose safer speeds. Their report describes

(37)

an investigation on how actual vehicle running speeds are associated with the occur- rence and severity of motor vehicles crashed in conjunction with roadway and roadside characteristics.

In Malaysia, guideline on geometric design of roads had been published by Jabatan Kerja Raya (1986) (JKR) and Hamzah et al. (1992). These guidelines were applied to all new construction and improvements of roads for vehicular traffic . In JKR’s module, they categorized the area of roads into two categories, which are rural and urban areas.

Then, the specific value of superelevation rate was applied to those areas in calculating the speed for drivers to operate their vehicles.

For practical purposes in establishing the design criteria for the roads, a maximum superelevation rate of 0.10 is used for roads in rural areas and 0.06 for roads in urban areas. The roads in rural areas are further classified into five categories by function namely expressway, highway, primary road, secondary road and minor road, and into four categories in urban areas, namely, expressway, arterial, collector and local street.

Lower design speeds are usually adopted for urban roads to take into account the nature of traffic and adjoining land use. Thus, every categorized area of the road should be connected and consistently designed.

2.3.1 Design Consistency

Consistency of the design can be quantified by examining speed measures, align- ment measures and workload measures. Consistency of design can be defined as the conformance of a road’s geometric and operational features with driver expectancy (Wooldridge, 2003). Design consistency is a crucial element in designing roads to

(38)

avoid situations where drivers are surprised by unexpected features of the road. Thus, it is important to construct roads with good continuity and less sharp corners to im- prove safety.

Researches were put on the table over the years to develop an efficient methodol- ogy. The primary purpose was to remove the element of surprise from driving tasks.

Surprise elements can put drivers under pressure which lead them to making wrong decisions and possibly resulting in losing control of the vehicle (Wolhuter, 1994). Fur- ther researches was made to predict operating speeds on curves of a range of radii.

Various methods have been developed to quantify the consistency of designs in terms of variations in speed along the length of the road or drivers’ workload. It is now possible to assess the consistency of design objectives as opposed to reliance on sub- jective judgement. Good operating speed of the roadway can be assessed by its design consistency.

2.3.2 Design Speed

Design speed was the first measure to be applied to the achievement of the con- sistency of design but has been shown to have weaknesses. A methodology had been developed where consistency can be quantified (Lamm et al., 1999). Roads are de- signed according to a speed design which is constant for a given stretch of roadway.

Thus, a vehicle must be able to travel comfortably and safely at the length of a given stretch of road at the design speed regardless of bends. Design speed was defined as the maximum speed at which the road could be traversed safely when only the geome- try of the road dictated what this speed could be. A major weakness of design speed as

(39)

a measure of consistency is that it only defines the minimum geometric standards that could be adopted in the design of roads. Consistency could thus imply that the entire road should have geometric standards that are close to minimum values.

This definition of consistency had been replaced by the statement that ‘Design speed is the speed selected to determine the various geometric design features of the road for design’ by (AASHTO, 2001). Design speed is thus no longer a speed per se but rather constitutes a grouping of geometric standards. Consistency of design dictates that it should be logical with regard to the topography, adjacent land used, and the functional classification of the road, all of which would have a bearing on expected operating speeds.

Vehicle running speed associates with consistency of design of the pathway. Two long tangents connected by a 60 km/h radius curve could quite correctly and yet totally misleadingly be described as having a design speed of 60 km/h. The next development was to use the operating speed from 85 per cent of maximum design speed, as the basis for achieving consistency of design (Fitzpatrick et al., 2003). This new lower speed was adopted for design, albeit briefly. It was then found that, under conditions of long tangents and gentle radii of curvature, operating speed could actually be higher than design speed by anything up to 10 or even 20 km/h. In consequence, operating speed regained respectability as a measure of design consistency.

A successful design will lead to a road that provides safety, accessibility, mobility, and convenience to the movement of people and goods with minimum side effects. A road can never be guaranteed totally safe, but it can be designed to have safety features.

(40)

Usually, crash victims would say that the factors that contributed to the accident were black spots rather than individual error. Black spots occur on roads where every ele- ment conforms to standard and suggest that many drivers are making similar mistakes at the same places (Wolhuter, 1994).

(41)

CHAPTER 3

TRIGONOMETRIC BÉZIER CURVE

3.1 Introduction

Curves are commonly used in industry and engineering field (Han et al., 2009), while surfaces are mostly being used by medical practitioners to aid them during med- ical operational process (Iscan and Steyn, 2013). Uses of curves such as B-spline and Bézier are based on their basis that span the entire element. Over the years, other than polynomial space, many basis had been presented in new spaces.

Peña (1997) constructed a basis Cm =span{1,cost,....,cosmt} where m∈ Z+. Curves in the space span of{1,cost,sint,....,cosmt,sinmt}have been developed by Zhang (1996). Ahmad et al. (2014) constructed a new generalization basis called Bézier-like curve. However, Mainar et al. (2001) also constructed α-basis of quar- tic Bernstein. The other approach by Chen and Wang (2003), used integral as their basis called C-Bézier.

Most of these bases do not enclose free form curves of higher-order polynomial.

Recently, the degree of trigonometric Bézier curve expanded from quadratic by Bashir et al. (2012), cubic trigonometric Bézier curve with shape parameter by (Han et al., 2009), quartic (Dube and Sharma, 2013; Zhu et al., 2012), quasi-quintic (Bashir et al., 2013), and quintic (Dube and Yadav, 2014).

In this chapter, we present a new basis of quintic trigonometric Bézier curve with

(42)

two shape parameters. This new basis provides another alternative curve compared to Dube and Yadav (2014). Dube and Yadav (2014) used four control points that mimic cubic trigonometric Bézier curves with a single shape parameter. Here, six control points with two shape parameters and some properties are discussed. Later, we will concentrate on how different parameters affect those types of curve.

3.2 Cubic Trigonometric Bézier Curve

Cubic trigonometric Bézier curve with two shape parameters had been proposed by Han et al. (2009). The two shape parameters allow designers to change the curve according to their desired shape while keeping the control polygon unchanged. More- over, this unique function can also approximate ellipses and arc segments of circles.

3.2.1 Cubic Trigonometric Bézier Basis Function

Cubic trigonometric Bézier curve with two shape parameters from Han et al. (2009) is given by:

z(t) =

3

i=0

Pifi(t), t [0,1] (3.1)

wherePiare the control points and fiare the cubic trigonometric basis functions for all i=0,1,2,3 with values of shape parameters p,q∈[−2,1]. The fiare defined as:

f0(t) = (1sinπt

2)2(1−psinπt 2), f1(t) =sinπt

2 (1sinπt

2 )(2+p−psinπt 2 ), f2(t) =cosπt

2 (1cosπt

2 )(2+q−qcosπt 2 ), f3(t) = (1cosπt

2)2(1−qcosπt 2 ).

(3.2)

(43)

Figure 3.1 shows plots of cubic trigonometric Bézier basis function for two arbi- trarily selected real values of p and q, where p,q∈[−2,1]. In Figure 3.1, different value of shape parameters give different projections of the basis function as such red- dotted curve (p=2,q=2), blue-dotted curve(p=1,q=1), purple-dashed curve(p=0,q=0)whereas black-line curve(p=1,q=1).

Figure 3.1: Cubic trigonometric Bézier basis function with different values of shape parameters

Cubic trigonometric Bézier basis function in Equation (3.2) have the following properties:

i. Non-negativity: fi t

0, fori=0,1,2,3.

ii. Partition of unity:∑3i=0 fi(t)1.

iii. Symmetry: fi(t;p,q) = f3−i(1−t;p,q)for everyi=0,1,2,3.

The proof for properties of cubic trigonometric Bézier can be obtained in Han et al.

(2009).

(44)

3.2.2 Properties of Cubic Trigonometric Bézier Curve

Cubic trigonometric Bézier curve in Equation (3.1) has the following properties (Han et al., 2009):

(a) Endpoint terminal

The endpoint terminal for cubic trigonometric Bézier curve are provided in Equation (3.3) until Equation (3.5):

z(0) = P0 z(1) = P3

(3.3)

z(0) = π2(P1−P0)(2+p) z(1) = π2(P3−P2)(2+q)

(3.4)

z(0) = π22[(2p+1)(P1−P0) + (P2−P1)]

z(1) = π22[(2q+1)(P3−P2) + (P2−P1)]

(3.5)

(b) Convex hull

The entire cubic trigonometric Bézier curve segment must lie inside its convex hull spanned byP0,P1,P2,P3.

(c) Symmetry

P0,P1,P2,P3andP3,P2,P1,P0define the same trigonometric Bézier curve in different parametrization, i.e.,z(t;p,q:P0,P1,P2,P3) =z(1−t;p,q:P3,P2,P1,P0) where 0≤t≤1 and p,q∈[−2,1]if and only if p=q.

(45)

(d) Geometric invariance

Geometric invariance for cubic trigonometric Bézier curve are as follows:

z(t;p,q:P0+m,P1+m,P2+m,P3+m) =z(t;p,q:P0,P1,P2,P3) +m z(t;p,q:P0×n,P1×n,P2×n,P3×n) =z(t;p,q:P0,P1,P2,P3)×n

where 0≤t≤1, p,q∈[−2,1],mis an arbitrary vector inR2orR3, andnis an arbitraryd×dmatrix withd=3 or 4.

3.3 Quintic Trigonometric Bézier Curve

In this section, a new type of trigonometric function called quintic trigonometric Bézier curve with two shape parameters is proposed. This function consists of six basis functions and possesses some of the properties that will be discussed below.

3.3.1 Quintic Trigonometric Bézier Basis Function

Quintic trigonometric Bézier curve with two shape parameters are given as follows:

z(t) =i=0

5 Pigi(t), t[0,1] (3.6)

wherePiare the control points andgi are the quintic trigonometric basis functions for alli=0,1,2,3,4,5 with values of shape parameters p,q∈[−4,1]. Thegi are defined as:

(46)

g0(t) = (1sinπt

2 )4(1−psinπt 2 ), g1(t) =sinπt

2 (1sinπt

2)3(4+p−psinπt 2 ), g2(t) = (1sinπt

2 )2(1cosπt

2 )(8 sinπt

2 +3 cosπt 2 +9), g3(t) = (1cosπt

2 )2(1sinπt

2 )(8 cosπt

2 +3 sinπt 2 +9), g4(t) =cosπt

2 (1cosπt

2)3(4+q−qcosπt 2 ), g5(t) = (1cosπt

2 )4(1−qcosπt 2 ).

(3.7)

Quintic trigonometric Bézier basis function have the following properties:

i. Non-negativity:gi t

0, fori=0,1,2,3,4,5. ii. Partition of unity:∑5i=0gi(t)1

iii. Symmetry:gi(t;p,q) =g5−i(1−t;p,q), ifp=qfori=0,1,2,3,4,5. Proof.

i. Fort∈[0,1]andp,q∈[−4,1], then 1−psinπt2 0, 1−qcosπt2 0, 1sinπt2 0, 1cosπt2 0, sinπt2 0, cosπt2 0, 4+p−psinπt2 0, 4+q−qcosπt2 0.

Therefore,

g0= (1sinπt2)4(1−psinπ2t)0,

g1=sinπt2(1sinπt2)3(4+p−psinπt2)0,

g2= (1sinπt2)2(1cosπt2)(8 sinπt2 +3 cosπ2t+9)0, g3= (1cosπ2t)2(1sinπt2)(8 cosπt2 +3 sinπ2t+9)0, g4=cosπt2(1cosπt2)3(4+q−qcosπt2)0,

g5= (1cosπ2t)4(1−qcosπt2)0.

Rujukan

DOKUMEN BERKAITAN

If cubic trigonometric B-spline was a superior method than cubic B-spline for solving two-point boundary value problems and one-dimensional hyperbolic equation with trigonometric

Figure 8 Cubic Timmer basis functions (except inner points) Based on Goodman and Said [8], the value of each control point for B ѐ zier and Ball triangular patch is based

Furthermore, G 2 continuity at joints when composing a curve from the segments using cubic Bezier and PH quintic spiral seem as a smooth transition, but it is not adequate

Therefore, this study is to design the S and C-shaped transition curves using Said-Ball cubic curve based on the third and fifth cases of circle to circle templates with

Integrating Spur Gear Teeth Design and its Analysis with G 2 Parametric Bézier-Like Cubic Transition and Spiral Curves, Global Engineers &.. Technologists

A 4×4 full factorial design with the principal experimental factors as pod storage (0, 3, 7 and 10 days) and roasting time (0, 15, 30 and 45 minutes) at 120 o C were used

(2015), leading to the inefficiency of bell pepper in inhibiting microbial growth due to level of capsaicin is too low to perform its antibacterial function. and

In this paper, we constructed the series of fuzzified fuzzy Bezier curve based on the modeling uncertainty data points through Bezier curve function and also the