CN103646150B - Continuous transition curve constitution method based on quintic Bezier curve - Google Patents
Continuous transition curve constitution method based on quintic Bezier curve Download PDFInfo
- 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
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
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′′=λv1+ω2v2, 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 And set Can get following formula:
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):
Can be obtained by equation (1), (3), (5), (7):
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 "=λ v1+ω2v2, 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:
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):
Obtained by equation (1), (3), (5), (7):
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.
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)
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)
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 |
-
2013
- 2013-12-23 CN CN201310716983.4A patent/CN103646150B/en not_active Expired - Fee Related
Patent Citations (2)
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)
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 |