CN103646150B - Continuous transition curve constitution method based on quintic Bezier curve - Google Patents

Continuous transition curve constitution method based on quintic Bezier curve Download PDF

Info

Publication number
CN103646150B
CN103646150B CN201310716983.4A CN201310716983A CN103646150B CN 103646150 B CN103646150 B CN 103646150B CN 201310716983 A CN201310716983 A CN 201310716983A CN 103646150 B CN103646150 B CN 103646150B
Authority
CN
China
Prior art keywords
prime
curve
omega
bezier
beta
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201310716983.4A
Other languages
Chinese (zh)
Other versions
CN103646150A (en
Inventor
宁涛
陈志同
席平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201310716983.4A priority Critical patent/CN103646150B/en
Publication of CN103646150A publication Critical patent/CN103646150A/en
Application granted granted Critical
Publication of CN103646150B publication Critical patent/CN103646150B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Generation (AREA)

Abstract

The invention discloses a G3 continuous transition curve constitution method based on a quintic Bezier curve. The method comprises the following steps: step 1, working out the first-order steering vector and the second-order steering vector at the end point of a primitive curve; step 2, enabling the primitive curve and a transition curve to satisfy a G2 continuity condition at the connection point; step 3, on the basis of the step 2, enabling the primitive curve and the transition curve to satisfy G3 continuity at the connection point; step 4, inversely calculating the reference mark of a quintic Bezier transition curve; step 5, normally calculating the quintic Bezier transition curve according to the control peak in the step 4. According to the invention, a Bezier curve is adopt to construct the transition curve, so that the advantages of controlling the shape of the curve easily and visually and providing convenience for interactive design can be achieved; besides, the few the parameters are, the few the orders are, so that the calculation frequency is low, and the memory space occupied in a computer is relatively small; the invention can be applied to construction of transition curves in CAD modeling.

Description

A kind of g3 continuous transition curve constructing method based on 5 bezier curves
Technical field
The present invention relates to a kind of g3 continuous transition curve constructing method based on 5 bezier curves, belong to cad technology neck Domain.
Background technology
G3 is continuously sometimes very useful in cad modeling, such as: the fillet curve of blade front and rear edge, iphone phone housing, vapour Car design etc..G3 continuously belongs to curvature variation continuously, and this continuous rank not only has g0, g1, g2 continuous level another characteristic Outside, at junction point, curvature variation is also continuous, and this makes Curvature varying more smooth, the surface of this continuous rank There is the visual effect more smooth than g2.Some cad softwares have the construction ability of g3 full curve curved surface, such as autodesk at present Alias2011essentials constructs g3 full curve with 7 order polynomials.But 7 order polynomial degree of curve are high, curve easily occurs " vibration " is so difficult model is interacted with formula design, and curve exponent number is higher it is desirable to the number of times calculating is more, and storage accounts for Also bigger with space.For above-mentioned problem present invention utilizes three order derivatives of two ends restrictive curve are it is proposed that use 5 Secondary bezier curve construction g3 continuous transition curve method.
Content of the invention
1st, purpose: it is an object of the present invention to provide a kind of g3 continuous transition curve constructing method based on 5 bezier curves, Improve the quality of easement curve moulding in cad modeling process.
2nd, technical scheme: the purpose of the present invention is achieved through the following technical solutions.
A kind of g3 continuous transition curve constructing method based on 5 bezier curves, it comprises the following steps:
Step one, obtain single order at end points for the primitive curve and lead resultant second order and lead arrow;
Step 2, primitive curve and easement curve is made to meet the g2 condition of continuity at junction point;
Step 3, on the basis of step 2, make primitive curve and easement curve meet g3 at junction point continuous;
Step 4, according to 5 bezier easement curve control point of step 3 inverse;
Step 5, just calculating 5 bezier easement curves according to the control vertex in step 4.
Wherein, " primitive curve " described in step one refers to need to carry out two parameter curves of transition connection, " single order Lead resultant second order and lead arrow " refer to the derivation to general parameterss for the parameter curve;
Wherein, " the g2 condition of continuity " described in step 2 refers to that two groups of curves meet g0, g1 continuous strip at junction point Under part, public curvature is had to swear (size and Orientation);G0 continuously refers to an end points of a curve and the one of another curve End points contacts;G1 continuously refer to two curves junction meet g0 continuous under the conditions of, cut arrow direction identical.
Wherein, " g3 is continuous " described in step 3 refers to 5 bezier curves and primitive curve in two junction points The derivative of the curvature at place is continuous;
Wherein, " 5 bezier easement curve control point of inverse " described in step 4 refer to drawing 5 bezier mistakes First control vertex is calculated, 5 times bezier easement curve control point has 6 before crossing curve;
Wherein, " just calculating 5 bezier easement curves " described in step 5 refers to according to the control top obtaining in step 4 Point draws the expression formula of 5 bezier easement curves.
3rd, advantage and effect
The invention provides with the continuous easement curve method of 5 bezier curve construction g3, using bezier curve structure Make easement curve and there is the shape being easy to intuitively controlling curve, the advantage facilitating interaction design, and number of parameters is fewer, rank Number is lower, and calculation times are also fewer, and memory space is less in a computer, and the present invention can be used for easement curve in cad modeling Construction.
Brief description
Fig. 1 combs schematic diagram for curvature of curve
Fig. 2 is for the present invention with 5 bezier curve construction g3 continuous transition curve synoptic diagrams
Fig. 3 is easement curve one, second dervative schematic diagram at junction point
Fig. 4 is FB(flow block) of the present invention.
In figure code name, symbol description are as follows:
G0, g1, g2, g3 curve g0, g1, g2, g3 are continuous
p0, p1, p2, p3, p4, p55 bezier curve controlled summits
r0Curve 0 and the junction point of easement curve
r1Curve 1 and the junction point of easement curve
R ' easement curve first derivative
R ' ' easement curve second dervative
Specific embodiment
See Fig. 1 Fig. 4, a kind of g3 continuous transition curve constructing method based on 5 bezier curves of the present invention, the method Specific implementation step is as follows:
Step one: obtain the single order in junction for the primitive curve and lead resultant second order and lead arrow,Represent respectively The one of one, two derivatives of curve 0 and easement curve junction point and curve 1 and easement curve junction point, second dervative, " point " represents To primitive curve derivation, " slash " represents to easement curve derivation;
Step 2: curve 0,1 and easement curve meet g2 continuously, and the expression formula meeting this condition is:r′′=λv12v2, wherein ω, λ are arbitrary coefficient;Can make in this stepWherein, α, β are unit vector, and have α ⊥ β, v2For?Vertical Projection vector on direction;
Step 3: according to curvature formulationsIf easily proving curve r (t) at end points, second dervative It is changed to r ' and r ' ', then curvature at end points for the r (t) is constant.By adjust ω make the mould length of first derivative take reasonable value (if Parameter field is [0,1], then | r ' | takes curve arc long), suitably choose λ value and make curve reach curvature derivative and be continuouslyAt this Can set in step f ( t ) = ( k 2 ) ′ = ( k 2 ) · = ( r ′ × r ′ , r ′ × r ′ ) ( r ′ , r ′ ) 3 , And set d 0 = | r · | , d 2 = | v 2 | Can get following formula:
f ′ ( t ) = ( r ′ × r ′ , r ′ × r ′ ) ′ ( r ′ , r ′ ) 3 - 6 ( r ′ × r ′ , r ′ × r ′ ) ( r ′ , r ′ ) 2 ( r ′ , r ′ ) ( r ′ , r ′ ) 6 = 2 ( r ′ × r ′ , r ′ × r ′ ′ ) ′ ( r ′ , r ′ ) - 6 ( r ′ × r ′ , r ′ × r ′ ) ( r ′ , r ′ ) ( r ′ , r ′ ) 4 = f ′ ( t ) = 2 ( r ′ ′ , β ) d 0 d 2 - 6 λωd 2 2 ω 2 d 0 5
Step 4: 5 times bezier curve representation formula is:bj,nT () is bernstein Base,
So, have at endpoint curve: p (0)=p0, p (1)=p5, according to even Continuous condition has r0=p0, r1=p5.One, two, three ranks at 5 bezier curve two-end-points lead arrow:
p0′=5(p1-p0) (2)
p5′=5(p5-p4) (3)
p0′′=20(p2-2p1+p0) (4)
p5′′=20(p5-2p4+p3) (5)
p0′′′=60(p3-3p2+3p1-p0) (6)
p′0、p′5、p″0、p″5、p″′0、p″′5(7)
Can be obtained by equation (1), (2), (4), (6):
- ω 0 d 02 d 00 λ 0 + ( v 11 , β 0 ) λ 1 = 3 ω 0 2 d 02 - ω 1 2 ( v 12 , β 0 ) - 20 ( - p 0 + 2 p 4 - p 5 , β 0 ) + ω 0 1 d 00 4 k 0 k · 0 3 d 02 - - - ( 8 )
Can be obtained by equation (1), (3), (5), (7):
( v 01 , β 1 ) λ 0 + ω 1 d 12 d 10 λ 1 = - ω 0 2 ( v 02 , β 1 ) + 3 ω 1 2 d 12 + 20 ( p 0 - 2 p 1 + p 5 , β 1 ) - ω 1 2 d 10 4 k 1 k · 1 3 d 12 - - - ( 9 )
Wherein, ωi、λi、kiβi, vi1、vi2、di0、di2Represent respectively curve i ω when being connected with easement curve, λ, k,β、v1、v2、d0、d2Value (i=0,1), λ can be solved by equation (8), (9)0、λ1
Step 5: the λ being drawn by step 40Obtain p with equation (2), (4)2, λ1Obtain p with equation (3), (5)3, All of control point all obtains afterwards, and then obtains 5 bezier easement curves of two groups of curves.
Present invention utilizes three order derivatives of two ends restrictive curve, only need to cutting along junction point in cad modeling process Control point p is moved in arrow direction1And p4, thus determine ω value, this movement can be designed to Interactive control, facilitate user Intuitively change the shape of easement curve, after ω determines, two other control point is obtained by step of the present invention, thus 5 times Bezier easement curve also determines that.

Claims (1)

1. a kind of g3 continuous transition curve constructing method based on 5 bezier curves it is characterised in that: it includes following step Rapid:
Step one: obtain the single order in junction for the primitive curve and lead resultant second order and lead arrow,Represent curve 0 respectively With the one of one, two derivatives of easement curve junction point and curve 1 and easement curve junction point, second dervative, " point " represents to former Beginning curve derivation, " slash " represents to easement curve derivation;
Step 2: curve 0,1 and easement curve meet g2 continuously, and the expression formula meeting this condition is:R "=λ v12v2, wherein ω, λ are arbitrary coefficient;Make in this step Wherein, α, β are unit vector, and have α ⊥ β, v2For?Vertical direction on projection Vector;
Step 3: according to curvature formulationsIf easily proving that curve r (t) at end points, second dervative are changed to R ' and r ", then curvature at end points for the r (t) is constant;The mould length of first derivative is made to take reasonable value by adjusting ω, if parameter field For [0,1], then | r ' | takes curve arc long, chooses λ value and makes curve reach curvature derivative and is continuouslyThis step 3 setsAnd setObtain following formula:
f ′ ( t ) = ( r ′ × r ′ ′ , r ′ × r ′ ′ ) ′ ( r ′ , r ′ ) 3 - 6 ( r ′ × r ′ ′ , r ′ × r ′ ′ ) ′ ( r ′ , r ′ ) 2 ( r ′ , r ′ ′ ) ( r ′ , r ′ ) 6 = 2 ( r ′ × r ′ ′ , r ′ × r ′ ′ ′ ) ′ ( r ′ , r ′ ) - 6 ( r ′ × r ′ ′ , r ′ × r ′ ′ ) ( r ′ , r ′ ′ ) ( r ′ , r ′ ) 4 = f ′ ( t ) = 2 ( r ′ ′ ′ , β ) d 0 d 2 - 6 λωd 2 2 ω 2 d 0 5 - - - ( 1 )
Step 4: 5 times bezier curve representation formula is:bj,nT () is bernstein base,So, have at endpoint curve: p (0)=p0, p (1)=p5, according to the condition of continuity There is r0=p0, r1=p5;One, two, three ranks at 5 bezier curve two-end-points lead arrow:
p0'=5 (p1-p0) (2)
p5'=5 (p5-p4) (3)
p0"=20 (p2-2p1+p0) (4)
p5"=20 (p5-2p4+p3) (5)
p0" '=60 (p3-3p2+3p1-p0) (6)
p′0、p′5、p″0、p″5、p″′0、p″′5(7)
Obtained by equation (1), (2), (4), (6):
- ω 0 d 02 d 00 λ 0 + ( v 11 , β 0 ) λ 1 = 3 ω 0 2 d 02 - ω 1 2 ( v 12 , β 0 ) - 20 ( - p 0 + 2 p 4 - p 5 , β 0 ) + ω 0 2 d 00 4 k 0 k · 0 3 d 02 - - - ( 8 )
Obtained by equation (1), (3), (5), (7):
( v 01 , β 1 ) λ 0 + ω 1 d 12 d 10 λ 1 = - ω 0 2 ( v 02 , β 1 ) + 3 ω 1 2 d 12 + 20 ( p 0 - 2 p 1 + p 5 , β 1 ) - ω 1 2 d 10 4 k 1 k · 1 3 d 12 - - - ( 9 )
Wherein, ωi、λi、kiβi, vi1、vi2、di0、di2Represent respectively curve i ω when being connected with easement curve, λ, k,β、 v1、v2、d0、d2Value, i=0,1, solve λ by equation (8), (9)0、λ1
Step 5: the λ being drawn by step 40Obtain p with equation (2), (4)2, λ1Obtain p with equation (3), (5)3, finally own Control point all obtain, and then obtain 5 bezier easement curves of two groups of curves.
CN201310716983.4A 2013-12-23 2013-12-23 Continuous transition curve constitution method based on quintic Bezier curve Expired - Fee Related CN103646150B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310716983.4A CN103646150B (en) 2013-12-23 2013-12-23 Continuous transition curve constitution method based on quintic Bezier curve

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310716983.4A CN103646150B (en) 2013-12-23 2013-12-23 Continuous transition curve constitution method based on quintic Bezier curve

Publications (2)

Publication Number Publication Date
CN103646150A CN103646150A (en) 2014-03-19
CN103646150B true CN103646150B (en) 2017-01-25

Family

ID=50251363

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310716983.4A Expired - Fee Related CN103646150B (en) 2013-12-23 2013-12-23 Continuous transition curve constitution method based on quintic Bezier curve

Country Status (1)

Country Link
CN (1) CN103646150B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3136192A1 (en) * 2015-08-24 2017-03-01 Siemens Aktiengesellschaft Control method for the movement of a tool and control device
WO2017113416A1 (en) * 2015-12-31 2017-07-06 深圳配天智能技术研究院有限公司 Smooth transition method for processing trajectories and processing device
CN107980109B (en) * 2017-01-04 2021-05-07 深圳配天智能技术研究院有限公司 Robot motion trajectory planning method and related device
CN109521731B (en) * 2017-09-19 2021-07-30 沈阳高精数控智能技术股份有限公司 G2 continuous Bezier tool path smoothing algorithm based on tolerance zone
CN108153245A (en) * 2017-12-26 2018-06-12 深圳市汇川技术股份有限公司 Smooth trajectory forwarding method and system
CN110749332B (en) * 2019-10-18 2021-04-20 广州文远知行科技有限公司 Curvature optimization method and device of RS curve, computer equipment and storage medium
CN112146667B (en) * 2020-09-29 2022-10-14 广州小鹏自动驾驶科技有限公司 Method and device for generating vehicle transition track

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001184514A (en) * 1999-10-12 2001-07-06 Sanei Giken:Kk Cad system
CN102831283A (en) * 2012-09-19 2012-12-19 河海大学常州校区 Complicated product model construction method based on surface feature

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001184514A (en) * 1999-10-12 2001-07-06 Sanei Giken:Kk Cad system
CN102831283A (en) * 2012-09-19 2012-12-19 河海大学常州校区 Complicated product model construction method based on surface feature

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
几何连续性约束下的自由曲面过渡问题研究;施侃乐;《中国博士学位论文全文数据库 信息科技辑》;20130715(第7期);第I138-45页 *
高阶几何连续过渡曲面的造型研究与应用;石雄辉;《中国优秀博硕士学位论文全文数据库 (硕士) 信息科技辑》;20040315(第1期);第I138-627页 *

Also Published As

Publication number Publication date
CN103646150A (en) 2014-03-19

Similar Documents

Publication Publication Date Title
CN103646150B (en) Continuous transition curve constitution method based on quintic Bezier curve
US7196702B1 (en) Geometric design and modeling system using control geometry
US8836701B1 (en) Surface patch techniques for computational geometry
CN108062453B (en) Flood high-efficiency simulation and high-fidelity visual dynamic display method for water conservancy information system
CN105183405A (en) 3D printing method for user-defined surface hollow model
CN102521863B (en) Three-dimensional fluid scalar vector uniform dynamic showing method based on particle system
CN101493954B (en) Three-dimensional modelling approach based on framework sketch drafting
CN101908235A (en) Novel method for modeling B spline surface
CN101482979A (en) Fairing-optimized CAD method for NURBS space curve curvature continuous split joint
CN104715508B (en) A kind of mesh modeling method and device
CN104570928A (en) Method for numerical control machining and path planning on mesh surface based on conformal parameterization
CN103559740B (en) A kind of prism mess generation method realizing interactive operation
US20240193313A1 (en) Surface patch techniques for computational geometry
JP2015001789A (en) Curved surface generation device, curved surface generation program, and curved surface generation method
CN102903138A (en) Shape-considered two-dimensional digital character skeleton operation method
CN104750995B (en) The computational methods of multiloop soil direct-buried cable conductor temperature
CN101908234B (en) User controlled method for generating highly structured triangular meshes
CN101299278A (en) Method for split joint of product external form space curve based on extension
CN104166158B (en) A kind of mutual method for tracing of geological data non-homogeneous spacing isopleth
CN101546351A (en) Geometric parameterization modeling method for optimizing variable-complexity shape
CN102930586A (en) Interactive geometry deformation method based on linear rotation invariant differential coordinates
CN113158307A (en) Method for dividing free-form surface nonsingular building grids
CN104091023A (en) Method for applying standards in electromagnetic field edge element method
CN110135083A (en) A kind of three-dimensional pipe network automated construction method of included space attribute
CN104463960A (en) Mean quadratic error metrics edge contraction simplifying method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170125

Termination date: 20201223