CN103903304A - Axle wire generating method of random quadrangular curved surface for product three-dimensional data processing - Google Patents

Abstract

The invention discloses an axle wire generating method of a random quadrangular curved surface for product three-dimensional data processing. Combination of a locus method and a movable Frenet frame is adopted, and a simple and quick axle wire generating algorithm is provided. A B-rep three-dimensional solid model which is frequently applied to engineering problems has specific border definition, information of top points, sides and faces of the model are specifically stored in a model file, the specific information can be used for simplifying an axle wire generating method in the traditional locus method, and therefore data size needed in generation of an axle wire is reduced. The specific border information in the B-rep model is used as conditions of the axle wire generating method, calculation of branch points is omitted through features of a bisector, and factors influencing the precision of the axle wire are reduced. Density of inserting points can be controlled to guarantee the precision of the axle wire.

Description

For the arbitrary quadrilateral curved surface axis generation method of product three-dimensional data processing

Technical field

The present invention relates to a kind of arbitrary quadrilateral curved surface axis generation method for the processing of product three-dimensional data, be particularly useful for the CAD/CAE technology such as model simplification, model reconstruction, finite element grid generation, solid modelling.Relating to Patent classificating number G06 calculates; Calculate; The 3D modeling of G06T17/00 for computerized mapping processed or produced to the view data that counting G06T is general.

Background technology

Model skeleton or center line are also in axis, conventionally be described as follows: in two-dimensional closed region, there is the disk of a radius variable at this intra-zone, keep this disk and zone boundary to have two point of contacts at least, then mobile this disk, after path that disk is walked covers All Ranges, the line that the track that the center of circle of disk stays forms is called axis, is called as middle axial plane in three-dimensional model.Axis is in each key link extensive application of CAD/CAE system, such as model reconstruction, model analysis, computer vision, solid modelling and in feature extraction of geometry designs etc.But axis is also restricted in some applications, because do not have at present fast a kind of and calculate accurately the algorithm of axis, particularly under curved boundary conditions, generate axis.

Blum proposed the concept of axis for the first time in 1967, the advantageous characteristic such as that axis has is unique, reversible, symmetrical, homeomorphic.It is very large that axis is affected by the continuity on border, and be also difficult at that time extract curved boundary accurately.Before 10 years, main discovery is the introducing of bisector, and it is a special case of axis, but the method needs cut operator, has increased the complexity of processing procedure.In recent years, Teixeira and Zucker have studied the relation between bent limit and axis curvature of a curve, tangent line and normal line vector three, have provided the differential equation of the curve axis motion of being described by mean curvature.Although there is the algorithm much generating about axis, generating stable and accurate axis is not also a pipe course.

The axis generating algorithm of commonly using in CAD/CAE system at present is roughly divided into two classes: method of loci and Voronoi figure method.Method of loci is the definition according to axis, allows the circle of a radius variable remain tangent with border, at least two, this region, and the left track of this round heart is exactly axis.But this algorithm need to judge whether distance and this point on each each border of some arrival are take-off point, have two positions with incision superius.This has increased data volume virtually, will inevitably affect the efficiency of calculating.Voronoi figure method be first by the border of given area point set by its discretize, generate Voronoi figure according to this discrete point set, what the summit of Voronoi figure formed is exactly axis.But this algorithm is more suitable for two dimensional model, need to increases extra alternative condition and just can be extended to three-dimensional model.

Summary of the invention

The present invention is directed to the proposition of above problem, and a kind of arbitrary quadrilateral curved surface axis generation method for the processing of product three-dimensional data of development has following steps:

-model decomposition step: quadrilateral surface model is decomposed into four limits; Non-conterminous both sides in four described limits are divided into one group, and two groups of limits that obtain are respectively as playing initial line and stopping limit and two binding sides;

-interpolation procedure: obtain the merging vector of two binding sides by calculating the vector of two described binding sides, using the direction of this merging vector as interpolation direction, insert multiple insertion points by step-length δ;

Step is moved in-insertion point: in the time that described multiple insertion points are all positioned at model inside, calculate the distance of each insertion point and described two binding sides; Use the mode of iteration, move described insertion point, make the distance approximately equal of this insertion point and described two binding sides;

-axis determining step: repeat above-mentioned steps; In the time of the equal approximately equal of distance of described insertion point and two binding sides, utilize described each insertion point structure B-spline curve, this curve approximation is the axis of quadrilateral surface model.

Step-length δ of the present invention is two distances between insertion point, utilizes step-length δ can control the quantity of insertion point.Step-length δ is less, and insertion point is more, and the axis precision generating is higher; Insertion point is fewer, and the calculated amount that generates axis is less.

Further, consider due to the difference of having chosen initial line and termination limit, may cause the volume coordinate of insertion point to be positioned at outside quadrilateral model.Therefore insertion point of the present invention is moved and is at least also had insertion point in step and judge and mobile step:

The distance d of-calculating insertion point P and two binding sides ₁and d ₂, and and d ₁and d ₂intersection point q with two binding sides ₁and q ₂, q ₁and q ₂between distance be d, d=||q ₁q ₂||; Obtain the direction vector of insertion point P according to Frenet frame formula ${\stackrel{\→}{e}}_1=q_1{q}_2/\left|\right|q_1{q}_2{\left|\right|}_2;$

-relatively d ₁, d ₂relation with d: if d ₁>d or d ₂>d, insertion point P is in model area outside;

If-d ₁>d ₂, put P along

direction displacement Δ=| (d ₁-d ₂) |/2; If d ₁<d ₂, put P along reverse direction displacement Δ=| (d ₁-d ₂) |/2; Insertion point P is moved to behind model inside, upgrade this point coordinate;

If-d ₁<d and d ₂<d, shows that a P, in model inside, does not temporarily adjust.

By above-mentioned decision process, can ensure that each insertion point is positioned at the inside of quadrilateral model, has ensured the carrying out of algorithm.

Preferably, described " use alternative manner, move described insertion point, make the distance approximately equal of this insertion point and described two binding sides " specifically comprises following steps:

The distance d of-calculating insertion point P and two binding sides ₁and d ₂, determining that insertion point is in the time that quadrilateral model is inner: obtain according to Frenet frame formula the direction vector that P position is adjusted

If-d ₁>d ₂, some P along

direction displacement Δ ₁;

If-d ₁<d ₂, some P along

reverse direction displacement Δ ₁;

-repeat above-mentioned steps, until d ₁/ d ₂=ε, stops computing; Wherein 0.99< ε <1.01; Displacement Δ ₁value be | (d ₁-d ₂) |/2 ⁿ, wherein n value is twice inconsistent number of times of continuous moving direction.

In described interpolation procedure, described " obtaining the merging vector of two binding sides by calculating the vector of two described binding sides " is specially:

Set A ₁and A ₂for two described binding sides; Press initial line to the direction that stops limit, obtained two binding side A ₁and A ₂vector

utilize vector merging formula (1) to calculate the interpolation direction vector of insertion point

$\stackrel{\&RightArrow;}{e}={\stackrel{\&RightArrow;}{e}}_{a1}+{\stackrel{\&RightArrow;}{e}}_{a2}---\left(1\right)$

Described " inserting multiple insertion points by step-length δ " is specially:

On an initial line, get mid point P ₁as starting point, according to the direction vector of interpolation

with step-length δ, utilize formula (2) to calculate the coordinate P of next insertion point ₂(x, y)

$\left\{\begin{array}{c}{P}_2.x=P_1.x+\stackrel{\&RightArrow;}{e}.x*\mathrm{\&delta;}\\ P_2.y=P_1.y+\stackrel{\&RightArrow;}{e}.y*\mathrm{\&delta;}\end{array}\right.---\left(2\right)$

If P ₂be less than step-length δ to the distance d that stops limit, represent that new insertion point approaches and stops limit, no longer needs new insertion point, so P ₂it is exactly last terminating point; Otherwise, continue to insert new insertion point P*.

Owing to having adopted technique scheme, the arbitrary quadrilateral curved surface axis generation method for the processing of product three-dimensional data provided by the invention, adopts method of loci in conjunction with mobile Frenet frame, has proposed a kind of simple and quick axis generating algorithm.In engineering problem, applying more is B-rep three-dimensional entity model, this model has clear and definite boundary definition, the information of the top points, edges, faces of model is clear and definite being stored in model file all, utilize these clear and definite information just can simplify traditional method of loci axis generation method, generate needed data volume thereby reduced axis.The present invention utilizes the condition of boundary information clear and definite in B-rep model as axis generation method, and utilizes the feature of bisector to eliminate the calculating to take-off point, reduces the factor that affects axis degree of accuracy.In order to ensure the degree of accuracy of axis, can adopt the density of insertion point to control.

Brief description of the drawings

For the technical scheme of clearer explanation embodiments of the invention or prior art, introduce simply the accompanying drawing of required use in embodiment or description of the Prior Art being done to one below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, do not paying under the prerequisite of creative work, can also obtain according to these accompanying drawings other accompanying drawing.

Fig. 1 is axis generating algorithm process flow diagram;

Fig. 2 is the each limit grouping of quadrilateral schematic diagram;

Fig. 3 is insertion point algorithm schematic diagram;

Fig. 4 is for adjusting insertion point to model schematic internal view;

Fig. 5 is for adjusting position, insertion point schematic diagram on axis;

Fig. 6 is that four irregular curved tetragonal axis, limit generate design sketch, and Fig. 6 a is original quadrilateral model; Fig. 6 b is the axis schematic diagram obtaining by algorithm of the present invention;

Fig. 7 is Fillet Feature face decomposing schematic representation;

Fig. 8 is that mechanical component model fillet is simplified the rear model of design sketch a master pattern b Fillet Feature simplification;

Fig. 9 utilizes axis parted pattern generating mesh a initial model b model four c limit, limit ab, cd to intersect for the axis that constraint condition generates axis e both direction for constraint condition generates d limit, axis bc, da, forms grid.

Embodiment

For making object, technical scheme and the advantage of embodiments of the invention clearer, below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is known to complete description:

As Figure 1-Figure 5: a kind of arbitrary quadrilateral curved surface axis generation method for the processing of product three-dimensional data, mainly comprises following step:

In model decomposition step, divide into groups shown in schematic diagram 2 in the each limit of quadrilateral.Every limit of quadrilateral is divided into two groups, A by " limit and limit are non-conterminous " principle ₁and A ₂be one group as binding side, limit B ₁and B ₂for another group, respectively as playing initial line and stopping limit.Select the main consideration on initial line peace treaty bundle limit to want the trend of the axis of selecting.

The described interpolation procedure of content of the present invention is as shown in Figure 3: insertion point algorithm schematic diagram.Set B ₁for playing initial line, B ₂for stopping limit, A ₁and A ₂for binding side.In the time specifically having selected initial line and stopped limit, the main trend of considering to want the target axis obtaining, as binding side, remaining two limits can be used as initial limit and stop limit with two limits of target axis approximate parallel (moving towards consistent).

Press initial line to the direction that stops limit, obtained the vector of two binding sides

utilizing vector to merge formula (1) calculates

be the interpolation direction vector of insertion point.

$\stackrel{\&RightArrow;}{e}={\stackrel{\&RightArrow;}{e}}_{a1}+{\stackrel{\&RightArrow;}{e}}_{a2}---\left(1\right)$

On an initial line, get mid point P ₁as starting point, according to the direction vector of interpolation

with step-length δ, utilize formula (2) to calculate the coordinate P of next insertion point ₂(x, y).

$\left\{\begin{array}{c}{P}_2.x=P_1.x+\stackrel{\&RightArrow;}{e}.x*\mathrm{\&delta;}\\ P_2.y=P_1.y+\stackrel{\&RightArrow;}{e}.y*\mathrm{\&delta;}\end{array}\right.---\left(2\right)$

If P ₂be less than step-length δ to the distance d that stops limit, represent that new insertion point approaches and stops limit, no longer needs new insertion point, so P ₂it is exactly last terminating point.Otherwise, continue to insert new insertion point P ^*.

In the present embodiment step-length δ be in quadrilateral model minor face 1/5th.

Further, consider due to the difference of having chosen initial line and termination limit, may cause the volume coordinate of insertion point to be positioned at outside quadrilateral model.Therefore insertion point of the present invention move in step, at least also have that insertion point is judged and mobile step as shown in Figure 4.Adjust insertion point to model schematic internal view.To acquired insertion point P, calculate the bee-line d of two binding sides of this some arrival ₁and d ₂, and the some q corresponding with bee-line ₁and q ₂(be positioned at binding side A ₁and A ₂on point), calculate direction vector according to Frenet frame formula.

${\stackrel{\&RightArrow;}{e}}_1=q_1{q}_2/\left|\right|q_1{q}_2{\left|\right|}_2;$

Judge that whether the position of insertion point P is in model area inside, if not in model inside, will adjust to model inside, judgement and method of adjustment are as follows:

According to a P to binding side A ₁and A ₂the corresponding point q of bee-line ₁and q ₂, obtain apart from d=||q ₁q ₂||.Comparison point P is to the bee-line d of binding side ₁, d ₂and the relation between d.

If d ₁>d or d ₂>d, puts P in model area outside.If d ₁>d ₂, put P along

direction move a certain distance Δ=| (d ₁-d ₂) |/2; If d ₁<d ₂, put P along

reverse direction move a certain distance Δ=| (d ₁-d ₂) |/2.Insertion point P is moved to behind model inside, upgrade this point coordinate;

If d ₁<d and d ₂<d, shows that a P is in model inside, does not temporarily need to adjust.

As Fig. 5: adjust position, each insertion point and make it drop on schematic diagram on axis.To each insertion point, adjust by following rule:

(a) if d ₁>d ₂, some P along

direction moves a certain distance Δ ₁;

(b) if d ₁<d ₂, some P along

reverse direction moves a certain distance Δ ₁;

(c) if d ₁=d ₂(approximately equal), upgrades some P coordinate figure.

Displacement Δ in the present invention ₁initial value is made as | (d ₁-d ₂) |/2.In the time that twice continuous moving direction is inconsistent, secondary displacement Δ ₁be reduced into | (d ₁-d ₂) |/4, adopt iterative algorithm repeatedly to adjust position a little until two constraint back gauges of some P arrival are equal.Displacement Δ ₁value be | (d ₁-d ₂) |/2 ⁿ, wherein n value is twice inconsistent number of times of continuous moving direction.Adopt alternative manner repeatedly to adjust position a little, the point inserting if make arrives binding side distance and equates there is certain difficulty completely, so approximately equal, i.e. d ₁/ d ₂=ε, ε approaches 1.For high-precision axis, it is 1 better that ε more approaches, but operation efficiency can reduce.

In the present invention, each model instance ε value is 0.99< ε <1.01.

As Fig. 6: four irregular curved tetragonal axis, limit generate design sketchs.In commercial production cad model, be regular compared with multi-model, but also have some irregular models.Adopt the high irregular boundary model of complexity to test axis generating algorithm.As shown in Figure 6 a, adopt the quadrilateral of the four irregular bendings in limit.Import after this model, application axis generating algorithm, effect is as shown in Figure 6 b.Although bilateral bending makes the direction of insertion point not be equal to the direction vector of binding side, through the adjustment to insertion point, insertion point can be dropped on axis.Recycling has obtained the some structure B-spline curve on axis, and this curve is the axis of bending quadrilateral curved surface, has met the symmetry of axis.

Utilizing axis to realize Fillet Feature simplifies.Patent " the Fillet Feature short-cut method based on fillet axis " (201310578463.1) discloses utilizes axis to cut apart Fillet Feature curved surface, the divided curved surface of reconstruct and corresponding seating surface, thus realize the method that Fillet Feature is simplified.Most of Fillet Feature are made up of four edges, two smooth limits and two common limits.With smooth limit, as constraint condition, common limit, as playing initial line and stopping limit, meets the axis generating algorithm that this patent proposes, just therefore can generate axis on this Fillet Feature curved surface.As shown in Figure 7, with smooth limit e ₄and e ₇for binding side, common limit e ₅and e ₆for playing initial line and stopping limit, generate axis e ₈fillet Feature curved surface is divided into two face F3 and F4; Two end points of axis are by limit e simultaneously ₅be divided into two sections of a ₁and a ₂, limit e ₆be divided into b ₁and b ₂two sections.Fig. 8 is mechanical component, is made up of 20 faces, wherein comprises multiple Fillet Feature, utilizes axis to cut apart reconstruct fillet surface, simplifies Fillet Feature, and the face number of model reduces to 11, and the model after simplification as shown in Figure 8 b.

Utilizing axis parted pattern to realize grid generates.In axis generation method, the constraint condition border having been generated as axis.Utilize axis parted pattern, the border of each subregion is also subject to the constraint on border, makes all subregion boundary similarity.The Jiang Gai region, axis generating in given area is divided into reciprocity two parts, if these two regions are continued to divide, just can obtain some approximately equalised subregions.Exchange the relation that binding side has been followed initial line, stopped limit, again carry out same division, obtained a series of sequential cells like this are the grid of generation.

Given closed region as shown in Fig. 9 a, is divided into quadrilateral: limit ab, bc, cd, da composition.Generate axis ef using limit ab, cd as constraint condition, by this Region Decomposition, form two regions, as shown in Fig. 9 b.Then judge decision condition: in each region, minor face is less than or equal to initial mesh border setting value.Dissatisfied decision condition is divided again, and as shown in Fig. 9 c, respectively using limit ae, fd and limit eb, cf as two axis of constraint condition regeneration, former region is divided into four sub regions, now meets decision condition and stops dividing.Exchange the relation that master mould binding side has been followed initial line, stopped limit, again carry out same division using limit bc, da as constraint condition, generate as three of Fig. 9 d axis.Twice different directions divided six axis that generate and intersected, and crossed node is grid node, and a series of continuous quadrilateral being formed by connecting by grid node is grid cell, as shown in Fig. 9 e.

The above; it is only preferably embodiment of the present invention; but protection scope of the present invention is not limited to this; any be familiar with those skilled in the art the present invention disclose technical scope in; be equal to replacement or changed according to technical scheme of the present invention and inventive concept thereof, within all should being encompassed in protection scope of the present invention.

Claims (5)

1. the arbitrary quadrilateral curved surface axis generation method for the processing of product three-dimensional data, has following steps:

-model decomposition step: quadrilateral surface model is decomposed into four limits; Non-conterminous both sides in four described limits are divided into one group, and two groups of limits that obtain are respectively as playing initial line and stopping limit and two binding sides;

-interpolation procedure: obtain the merging vector of two binding sides by calculating the vector of two described binding sides, using the direction of this merging vector as interpolation direction, insert multiple insertion points by step-length δ;

Step is moved in-insertion point: in the time that described multiple insertion points are all positioned at model inside, calculate the distance of each insertion point and described two binding sides; Use the mode of iteration, move described insertion point, make the distance approximately equal of this insertion point and described two binding sides;

-axis determining step: repeat above-mentioned steps; In the time of the equal approximately equal of distance of described insertion point and two binding sides, utilize described each insertion point structure B-spline curve, this curve approximation is the axis of quadrilateral surface model.

2. the arbitrary quadrilateral curved surface axis generation method for the processing of product three-dimensional data according to claim 1, is further characterized in that described insertion point moves and in step, at least also have insertion point and judge and mobile step:

The distance d of-calculating insertion point P and two binding sides ₁and d ₂, and and d ₁and d ₂intersection point q with two binding sides ₁and q ₂, q ₁and q ₂between distance be d, d=||q ₁q ₂||; Obtain according to Frenet frame formula the direction vector that P position is adjusted ${\stackrel{\&RightArrow;}{e}}_1=q_1{q}_2/\left|\right|q_1{q}_2{\left|\right|}_2;$

-relatively d ₁, d ₂relation with d: if d ₁>d or d ₂>d, insertion point P is in model area outside;

If-d ₁>d ₂, put P along

direction displacement Δ=| (d ₁-d ₂) |/2; If d ₁<d ₂, put P along

reverse direction displacement Δ=| (d ₁-d ₂) |/2; Insertion point P is moved to behind model inside, upgrade this point coordinate;

If-d ₁<d and d ₂<d, shows that a P, in model inside, does not temporarily adjust.

3. the arbitrary quadrilateral curved surface axis generation method for the processing of product three-dimensional data according to claim 1 and 2, be further characterized in that described " use alternative manner; mobile described insertion point, make the distance approximately equal of this insertion point and described two binding sides " specifically comprises following steps:

The distance d of-calculating insertion point P and two binding sides ₁and d ₂, determining that insertion point is in the time that quadrilateral model is inner: obtain according to Frenet frame formula the direction vector that P position is adjusted

If-d ₁>d ₂, some P along

direction displacement Δ ₁;

If-d ₁<d ₂, some P along

reverse direction displacement Δ ₁;

-repeat above-mentioned steps, until d ₁/ d ₂=ε, stops computing; Wherein 0.99< ε <1.01; Displacement Δ ₁value be | (d ₁-d ₂) |/2 ⁿ, wherein n value is twice inconsistent number of times of continuous moving direction.

4. the arbitrary quadrilateral curved surface axis generation method for the processing of product three-dimensional data according to claim 1, be further characterized in that in described interpolation procedure, described " obtaining the merging vector of two binding sides by calculating the vector of two described binding sides " is specially:

Set A ₁and A ₂for two described binding sides; Press initial line to the direction that stops limit, obtained two binding side A ₁and A ₂vector

utilize vector merging formula (1) to calculate the interpolation direction vector of insertion point

$\stackrel{\&RightArrow;}{e}={\stackrel{\&RightArrow;}{e}}_{a1}+{\stackrel{\&RightArrow;}{e}}_{a2}---\left(1\right)$

5. the arbitrary quadrilateral curved surface axis generation method for the processing of product three-dimensional data according to claim 4, is further characterized in that described " inserting multiple insertion points by step-length δ " is specially:

On an initial line, get mid point P ₁as starting point, according to the direction vector of interpolation with step-length δ, utilize formula (2) to calculate the coordinate P of next insertion point ₂(x, y)

$\left\{\begin{array}{c}{P}_2.x=P_1.x+\stackrel{\&RightArrow;}{e}.x*\mathrm{\&delta;}\\ P_2.y=P_1.y+\stackrel{\&RightArrow;}{e}.y*\mathrm{\&delta;}\end{array}\right.---\left(2\right)$

If P ₂be less than step-length δ to the distance d that stops limit, represent that new insertion point approaches and stops limit, no longer needs new insertion point, so P ₂it is exactly last terminating point; Otherwise, continue to insert new insertion point P*.

Images