CN108304686B - Adaptive dynamic evolution calculation method for rough machining process - Google Patents

Adaptive dynamic evolution calculation method for rough machining process Download PDF

Info

Publication number
CN108304686B
CN108304686B CN201810372957.7A CN201810372957A CN108304686B CN 108304686 B CN108304686 B CN 108304686B CN 201810372957 A CN201810372957 A CN 201810372957A CN 108304686 B CN108304686 B CN 108304686B
Authority
CN
China
Prior art keywords
region
characteristic
rough machining
residual
machining
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.)
Active
Application number
CN201810372957.7A
Other languages
Chinese (zh)
Other versions
CN108304686A (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.)
Changzhou Campus of Hohai University
Original Assignee
Changzhou Campus of Hohai 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 Changzhou Campus of Hohai University filed Critical Changzhou Campus of Hohai University
Priority to CN201810372957.7A priority Critical patent/CN108304686B/en
Publication of CN108304686A publication Critical patent/CN108304686A/en
Application granted granted Critical
Publication of CN108304686B publication Critical patent/CN108304686B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Numerical Control (AREA)

Abstract

The invention discloses a rough machining process self-adaptive dynamic evolution calculation method which is used for machining complex cavity parts. The type of the residual region is then identified based on the medial axis topology of the residual region. And secondly, establishing a dynamic evolution model of the rough machining area and the residual area of the part based on the initial process scheme according to the characteristic primary rough machining area and the residual area. And finally, providing a self-adaptive evolution mechanism of the rough machining process of the part driven by the process design intention so as to eliminate the problems of interference and under-cut in the rough machining process and ensure the connectivity of the rough machining process situation. The invention decides the self-adaptive evolution method of the rough machining process according to the types of the residues and the time of the residues, thereby making up the defects of the existing method from the geometric analysis, reducing the burden of process designers and improving the numerical control programming efficiency.

Description

Adaptive dynamic evolution calculation method for rough machining process
Technical Field
The invention relates to a rough machining process self-adaptive dynamic evolution calculation method for complex cavity parts, and belongs to the field of feature-based numerical control process design and reuse in the manufacturing industry.
Background
In recent years, the content-based three-dimensional CAD model retrieval technology is widely researched, and abundant research results are obtained, so that designers can be helped to effectively find similar three-dimensional geometry, and a brand-new support means is provided for refining and reusing process data. The basic principle is as follows: through the relevance of the three-dimensional geometry and the numerical control process, the similar numerical control process is found by the similar three-dimensional geometry, then through the reusability of the numerical control process, the similar numerical control process is transplanted to the query three-dimensional geometry on the part to be manufactured, and the query three-dimensional geometry is adapted through the self-adaption of the similar numerical control process. However, the existing three-dimensional CAD model retrieval and numerical control process reuse technology still stays at the discovery level of similar numerical control processes, and is mainly used for designers to refer to, and the reuse process of the technology is very dependent on the experience and knowledge of the designers. In addition, even though the existing research results can automatically set the process parameters in the commercial CAM system (such as CATIA, UG, etc.), the designer needs to perform interactive iterative modification on the similar numerical control process according to the simulation results (including undercutting and interference, etc.) of each processing stage of the commercial CAM system based on the process design intention, and particularly in the rough processing stage of the part, so as to ensure that each processing stage is a complete and through process situation (no undercutting or undercutting exists), and thus a great amount of human-computer interaction is still needed. In order to realize the process data driven adaptive numerical control process design and further reduce the burden of process designers, the following key problems need to be solved: (1) how to efficiently and dynamically calculate a machining area and a residual area of the complex cavity characteristics according to a similar numerical control process; (2) how to adaptively change the similar numerical control process according to the process design intention of each processing stage of the part and the dynamic evolution calculation of the characteristic processing region and the residual region.
Disclosure of Invention
In order to overcome the defects, the invention provides a rough machining procedure self-adaptive dynamic evolution calculation method, which is oriented to complex cavity parts and adopts a rough machining procedure self-adaptive evolution mechanism to eliminate the problems of interference and under-cut in the rough machining process and ensure the connectivity of the rough machining process situation.
In order to achieve the purpose, the technical scheme of the invention is realized as follows:
a rough machining procedure self-adaptive dynamic evolution calculation method comprises the following steps:
(a) the complex contour of the feature is divided into a plurality of simple sub-contours through dynamic evolution calculation of a complex cavity rough machining region MR and a residual region UR guided by a feature center axis, and the machining region MR is obtained through the bias of the simple sub-contours;
(b) and a self-adaptive evolution mechanism of a part rough machining process driven by a process design intention is provided to eliminate the problems of interference and under-cut in the rough machining process and ensure the connectivity of the rough machining process situation.
Preferably, the dynamic evolution calculation of the rough machining region MR and the residual region UR of the complex cavity guided by the characteristic central axis includes the following specific steps:
1.1: introducing a local inscribed circle moving representation model, and calculating a primary rough machining region MR based on a characteristic central axis;
1.2: constructing an attribute adjacency graph according to the local inscribed circle moving representation model in the step 1.1, and classifying and calculating the primary rough machining residual region UR;
1.3: and (3) according to the processing region MR and the residual region UR calculated in the steps 1.1 and 1.2, realizing the construction of a dynamic evolution model of the characteristic rough processing region MR and the residual region UR.
Preferably, the specific calculation steps of the primary rough machining region MR based on the characteristic central axis in step 1.1 are as follows:
1.1.1: moving representation model r of each local inscribed circle of feature central axis MA (F) of traversing complex cavity feature Fi(s) calculating a local inscribed circle movement representation model ri(s) and the horizontal line (R + delta) of the radius of the inscribed arcr) Is the center of the inscribed arc and is according to ri(s) and R + deltarComparing the sizes of the characteristic axes MA (F) with the side e of the characteristic axis MA (F)iDivided into multiple inscribed arcs siDividing the machining area constrained by the machining allowance into a zone association area machinable M and a zone association area non-machinable U;
1.1.2: representing the model r by local inscribed circle movementi(s) and its class, construct ri(s) attribute adjacency graph G;
1.1.3: the machinable M of the interval related area is taken as constraint and the depth is takenR in step 1.1.2i(s) obtaining a plurality of connected subgraphs G from the attribute adjacency graph GiEach connectivity sub-graph GiCorresponding to a simple sub-contour, i.e. connected sub-graph GiThe number is consistent with the number of the simple sub-outlines;
1.1.4: according to the connected subgraph G obtained in 1.1.3iOffset distance delta by feature profilerAnd calculating to obtain a processing area MR required by the complex cavity characteristic F in a primary rough processing stage.
Preferably, the specific calculation steps of the primary rough machining residual region UR in step 1.2 are as follows:
1.2.1: and (5) traversing r in the depth direction 1.1.2 by taking the unmanageable U of the interval association region as a constrainti(s) obtaining a plurality of connected subgraphs G from the attribute adjacency graph GiEach connectivity sub-graph giCorresponding to a sub-process residual SURiI.e. process residue URP=∪SURi
1.2.2: due to each connectivity graph giDifferent topological structures are provided, and each topological structure corresponds to one type of process residues; thus, according to the connectivity sub-graph g in 1.2.1iThe topology of (1.2.1), the process residues in (1.2.1) are divided into fillet residues F-URPResidual B-UR in single bottle neckPResidual MB-UR in multiple bottle necksPAnd composite residual C-URP
1.2.3: the process residue UR in step 1.2.2PThe range is controlled by the radius R of the tool, when R approaches 0, only the residual UR remainsAIt is recorded as
Figure GDA0003263361460000041
Therefore, URPThe calculation of (d) translates into a machining region MR (0, delta) with a tool radius of 0r) With a machining region MR (R, delta) of tool radius Rr) Between Boolean operations, i.e. URP=MR(0,δr)-MR(R,δr);
1.2.4: according to the inscribed arc s in step 1.1.1iAnd connectivity sub-graph g in 1.2.1iThe corresponding relation between the vertexes of (1) and the establishment of the SUR of the sub-process residueiAnd giIn betweenCorrelation relationship by which and in 1.2.2 according to giThe process residue classification of the topological structure of (1) judges the SUR of the sub-process residueiType (c) of the cell.
Preferably, the representation method of the local inscribed circle movement representation model is as follows: giving a complex cavity characteristic F, wherein the central axis MA is composed of a plurality of central axis edges eiComposition, i.e. MA (F) { ei},1≤i≤nmWherein n ismCharacteristic middle shaft edge e of complex cavity characteristic FiQuantity, characteristic axial edge eiAt an arbitrary point v oni(s),s∈[smin,smax]Where s is a characteristic arc length parameter corresponding to a point on the axial edge, there being one or more points each having a length vi(s) is a circle center Oi(s) local inscribed circle LIC, when LIC is along eiWhile moving, its radius ri(s) changing continuously, i.e. using ri(s) characterization of eiBalance ri(s) is eiThe LIC moving representation model of (1) is recorded as
Figure GDA0003263361460000042
The characteristic central axis of the complex cavity characteristic F is represented as MA (F) ═ ri(s)},1≤i≤nm
Preferably, said ri(s) attribute adjacency graph G vertices correspond to
Figure GDA0003263361460000051
The attribute of the endpoint is taken as the attribute of the vertex; r isi(s) Attribute adjacency graph G characteristic Axis edge eiThe corresponding movement representation model is
Figure GDA0003263361460000052
Class of which is characteristic of medial axis edge eiWherein the category includes a machinable M section associated region and a non-machinable U section associated region.
Preferably, said step 1.2.2 is according to giThe topology of (2) implements residual region classification, which is expressed as formula (1):
Figure GDA0003263361460000053
wherein a isiDenotes giVertex v in (1)iDegree, | G in Gi| represents attribute adjacency graph giNumber of vertices of (1), T (g)i) Denotes giCorresponding to the type of process residue.
Preferably, the dynamic evolution model of step 1.3 is denoted as Mk,jDynamically evolving model F from multiple featuresi.Mk,jConstitution Fi.Mk,jDriven by parameters reflecting the design intention of the process and realized by the process of removing the margin between adjacent processes in a dynamic evolution model, i.e.
Figure GDA0003263361460000054
Fi.Mk,j=Fi.Mk,j-1-Fi.Mk,j
MR(Fi.Mk,j)=f(Fi.Mk,j-1,R,δr)
Fi.Mk,j=(MR,zt,zb,n)
UR(Fi.Mk,j)=UR(Fi.Mk,j-1)-MR(Fi.Mk,j)
Wherein the parameters comprise a tool T and a radial allowance deltarAxial margin deltaaN is FiIn the axial direction of the tool, MR (F)i.Mk,j) A processing area F (F) obtained by projecting the characteristic dynamic evolution model to an XY plane along the Z-axis direction in a local coordinate systemi.Mk,j-1,R,δr) Constructing a function for the effective machining area, R being the radius of the tool T, zt、zbAre respectively Fi.Mk,jTop and bottom surface constraint of UR (F)i.Mk,j) And representing a residual region obtained by projecting the characteristic dynamic evolution model to an XY plane along the Z-axis direction in a local coordinate system.
Preferably, the adaptive evolution mechanism of the part rough machining process driven by the process design intent is divided into the following two types:
the first method comprises the following steps: longitudinal adaptation, in one process step, the complex cavity features F of which are machined in sequence, and each machining operation only machines one feature, given any one machining operation opiSuppose opiFor feature FnOf the machining region MRniProcessing is carried out when feature FmAnd feature FnWhen the formula (2) is satisfied, MR in the processing region is requiredniFeature F from top to bottom before cuttingmThe layers of (2) are left to be cut off, and the formula (2) is as follows:
Figure GDA0003263361460000061
wherein, FmAnd FnRepresenting two features with a longitudinal dependency between them, i.e. Fm→FnAnd which are machined in the roughing stage using different tools, URmj(j-1, 2) represents FmTwo-layer residue of (1), MRniIs represented by FnI th of (1)thCutting the secondary machining area;
and the second method comprises the following steps: lateral adaptation, assuming t1Each feature of time FmIs URm1For any of its sub-residual SURsiIf SUR is notiThe type of the (D) is fillet residue, so that the treatment is not needed; if SURiIf the type of (d) is non-fillet residue, then according to its corresponding sub-graph giEach side eiSelecting a tool T satisfying the condition shown in formula (3):
Figure GDA0003263361460000062
where R is the radius of the tool T, where T is along giWith each edge moving without interference, T being a tool, by adding a new roughing step to cut off the SURi
Has the advantages that: the invention provides a rough machining process self-adaptive dynamic evolution calculation method, which is applied to the machining of complex cavity parts and has the following advantages:
(1) establishing a local inscribed circle movement representation model to calculate a rough machining area and a residual area of the cutter considering the radial allowance, and judging the type of the residual according to the middle shaft topological structure of the residual area;
(2) a self-adaptive evolution mechanism of the rough machining process is provided, and the connectivity of the rough machining process situation of the part is ensured through the longitudinal and transverse self-adaptation of the step correlation characteristics in the rough machining stage of the part, so that the process design intention of the rough machining of the part is realized.
(3) The self-adaptive evolution method of the rough machining process is decided according to the type of the residue and the time of the residue, and the process design intention is used as a drive, so that the defects of the existing method from the geometric analysis are overcome;
(4) the burden of process designers is reduced, and the numerical control programming efficiency is improved.
Drawings
FIG. 1 is a general block diagram of the method of the present invention;
FIG. 2 is a schematic diagram of the calculation of the primary rough machining area based on the characteristics of the middle axis transformation in the method of the present invention;
FIG. 3 is a schematic illustration of the classification of process residues according to the method of the present invention;
FIG. 4 is a schematic illustration of a calculation of a primary rough process residue according to the method of the present invention;
FIG. 5 is a schematic diagram of the longitudinal adaptation of the roughing process of the present invention.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application are clearly and completely described below, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
A rough machining procedure self-adaptive dynamic evolution calculation method comprises the following steps:
(a) the complex contour of the feature is divided into a plurality of simple sub-contours through dynamic evolution calculation of a complex cavity rough machining region MR and a residual region UR guided by a feature center axis, and the machining region MR is obtained through the bias of the simple sub-contours;
(b) and a self-adaptive evolution mechanism of a part rough machining process driven by a process design intention is provided to eliminate the problems of interference and under-cut in the rough machining process and ensure the connectivity of the rough machining process situation.
Preferably, the dynamic evolution calculation of the rough machining region MR and the residual region UR of the complex cavity guided by the characteristic central axis includes the following specific steps:
1.1: introducing a local inscribed circle moving representation model, and calculating a primary rough machining region MR based on a characteristic central axis;
1.2: constructing an attribute adjacency graph according to the local inscribed circle moving representation model in the step 1.1, and classifying and calculating the primary rough machining residual region UR;
1.3: and (3) according to the processing region MR and the residual region UR calculated in the steps 1.1 and 1.2, realizing the construction of a dynamic evolution model of the characteristic rough processing region MR and the residual region UR.
Preferably, the specific calculation steps of the primary rough machining region MR based on the characteristic central axis in step 1.1 are as follows:
1.1.1: moving representation model r of each local inscribed circle of feature central axis MA (F) of traversing complex cavity feature Fi(s) calculating a local inscribed circle movement representation model ri(s) and the horizontal line (R + delta) of the radius of the inscribed arcr) Is the center of the inscribed arc and is according to ri(s) and R + deltarComparing the sizes of the characteristic axes MA (F) with the side e of the characteristic axis MA (F)iDivided into multiple inscribed arcs siDividing the machining area constrained by the machining allowance into a zone association area machinable M and a zone association area non-machinable U;
1.1.2: according to local inscribed circle movementRepresentation model ri(s) and its class, construct ri(s) attribute adjacency graph G;
1.1.3: traversing r in step 1.1.2 in depth by taking machinable M of the interval associated region as constrainti(s) obtaining a plurality of connected subgraphs G from the attribute adjacency graph GiEach connectivity sub-graph GiCorresponding to a simple sub-contour, i.e. connected sub-graph GiThe number is consistent with the number of the simple sub-outlines;
1.1.4: according to the connected subgraph G obtained in 1.1.3iOffset distance delta by feature profilerAnd calculating to obtain a processing area MR required by the complex cavity characteristic F in a primary rough processing stage.
Preferably, the specific calculation steps of the primary rough machining residual region UR in step 1.2 are as follows:
1.2.1: and (5) traversing r in the depth direction 1.1.2 by taking the unmanageable U of the interval association region as a constrainti(s) obtaining a plurality of connected subgraphs G from the attribute adjacency graph GiEach connectivity sub-graph giCorresponding to a sub-process residual SURiI.e. process residue URP=∪SURi
1.2.2: due to each connectivity graph giDifferent topological structures are provided, and each topological structure corresponds to one type of process residues; thus, according to the connectivity sub-graph g in 1.2.1iThe topology of (1.2.1), the process residues in (1.2.1) are divided into fillet residues F-URPResidual B-UR in single bottle neckPResidual MB-UR in multiple bottle necksPAnd composite residual C-URP
1.2.3: the process residue UR in step 1.2.2PMainly due to the capability of the tool, the extent of which is governed by the radius R of the tool, UR being the greater RPThe more, and when R is smaller, URPThe less, and thus when R goes to 0, there is only a residual URAIt is recorded as
Figure GDA0003263361460000101
Therefore, URPThe calculation of (d) translates into a machining region MR (0, delta) with a tool radius of 0r) With a machining region MR (R, delta) of tool radius Rr) BetweenBy Boolean operation of, i.e. URP=MR(0,δr)-MR(R,δr);
1.2.4: according to the inscribed arc s in step 1.1.1iAnd connectivity sub-graph g in 1.2.1iThe corresponding relation between the vertexes of (1) and the establishment of the SUR of the sub-process residueiAnd giThe correlation between the two and the g in 1.2.2iThe process residue classification of the topological structure of (1) judges the SUR of the sub-process residueiType (c) of the cell.
Preferably, the representation method of the local inscribed circle movement representation model is as follows: giving a complex cavity characteristic F, wherein the central axis MA is composed of a plurality of central axis edges eiComposition, i.e. MA (F) { ei},1≤i≤nmWherein n ismCharacteristic middle shaft edge e of complex cavity characteristic FiQuantity, characteristic axial edge eiAt an arbitrary point v oni(s),s∈[smin,smax]Where s is a characteristic arc length parameter corresponding to a point on the axial edge, there being one or more points each having a length vi(s) is a circle center Oi(s) local inscribed circle LIC, when LIC is along eiWhile moving, its radius ri(s) changing continuously, i.e. using ri(s) characterization of eiBalance ri(s) is eiThe LIC moving representation model of (1) is recorded as
Figure GDA0003263361460000102
The characteristic central axis of the complex cavity characteristic F is represented as MA (F) ═ ri(s)},1≤i≤nm
Preferably, said ri(s) attribute adjacency graph G vertices correspond to
Figure GDA0003263361460000103
The attribute of the endpoint is taken as the attribute of the vertex; r isi(s) Attribute adjacency graph G characteristic Axis edge eiThe corresponding movement representation model is
Figure GDA0003263361460000111
Class of which is characteristic of medial axis edge eiWherein the categories includeThe section related area can be processed by M and the section related area can not be processed by U.
Preferably, said step 1.2.2 is according to giThe topology of (2) implements residual region classification, which is expressed as formula (1):
Figure GDA0003263361460000112
wherein a isiDenotes giVertex v in (1)iDegree, | G in Gi| represents attribute adjacency graph giNumber of vertices of (1), T (g)i) Denotes giCorresponding to the type of process residue.
From the formula (1), the fillet residue F-URPHas only one edge and contains a vertex with a degree of 1, indicating a F-URPPresent at concave corners of the feature profile; for bottleneck residual B-URPAnd MB-URPEach of the edges of (a) having a degree of each vertex of not less than 2, such that they are generated primarily from non-concave corner elongate regions of the feature profile; and compound residual C-URPConsists of fillet remnants and bottleneck remnants, so that at least one vertex with the degree of 1 and a vertex with the degree of not less than 2 exist, and the 2 vertexes need a vertex with the degree of not less than 3 to connect the 2 vertexes together, so that the 2 vertexes mainly exist in an elongated area adjacent to the concave corner.
Preferably, the dynamic evolution model of step 1.3 is denoted as Mk,jDynamically evolving model F from multiple featuresi.Mk,jConstitution Fi.Mk,jDriven by parameters reflecting the design intention of the process and realized by the process of removing the margin between adjacent processes in a dynamic evolution model, i.e.
Figure GDA0003263361460000121
Fi.Mk,j=Fi.Mk,j-1-Fi.Mk,j
MR(Fi.Mk,j)=f(Fi.Mk,j-1,R,δr)
Fi.Mk,j=(MR,zt,zb,n)
UR(Fi.Mk,j)=UR(Fi.Mk,j-1)-MR(Fi.Mk,j)
Wherein the parameters comprise a tool T and a radial allowance deltarAxial margin deltaaN is FiIn the axial direction of the tool, MR (F)i.Mk,j) A processing area F (F) obtained by projecting the characteristic dynamic evolution model to an XY plane along the Z-axis direction in a local coordinate systemi.Mk,j-1,R,δr) Constructing a function for the effective machining area, R being the radius of the tool T, zt、zbAre respectively Fi.Mk,jTop and bottom surface constraint of UR (F)i.Mk,j) And (3) representing a residual region obtained by projecting the characteristic dynamic evolution model to an XY plane along the Z-axis direction in a local coordinate system (Z ═ -n).
Preferably, the adaptive evolution mechanism of the part rough machining process driven by the process design intent is divided into the following two types:
the first method comprises the following steps: longitudinal adaptation, in one process step, the complex cavity features F of which are machined in sequence, and each machining operation only machines one feature, given any one machining operation opiSuppose opiFor feature FnOf the machining region MRniProcessing is carried out when feature FmAnd feature FnWhen the formula (2) is satisfied, MR in the processing region is requiredniFeature F from top to bottom before cuttingmThe layers of (2) are left to be cut off, and the formula (2) is as follows:
Figure GDA0003263361460000122
wherein, FmAnd FnRepresenting two features with a longitudinal dependency between them, i.e. Fm→FnAnd which are machined in the roughing stage using different tools, URmj(j=1,2) Is represented by FmTwo-layer residue of (1), MRniIs represented by FnI th of (1)thCutting the secondary machining area;
and the second method comprises the following steps: lateral adaptation, assuming t1Each feature of time FmIs URm1For any of its sub-residual SURsiIf SUR is notiThe type of the (D) is fillet residue, so that the treatment is not needed; if SURiIf the type of (d) is non-fillet residue, then according to its corresponding sub-graph giEach side eiSelecting a tool T satisfying the condition shown in formula (3):
Figure GDA0003263361460000131
where R is the radius of the tool T, where T is along giWith each edge moving without interference, T being a tool, by adding a new roughing step to cut off the SURi
FIG. 1 is a general block diagram of a rough machining process adaptive dynamic evolution calculation method for complex cavity parts, and the working flow of the invention is as follows: firstly, according to the corresponding relation between the central axis of the characteristic and the circle center of an inscribed circular arc of a processing area, a local inscribed circle movement representation model is established to guide the calculation of the characteristic primary rough processing area MR and a residual area UR; then identifying the type of the middle shaft topological structure according to the residual region; secondly, establishing a dynamic evolution model of the rough machining area and the residual area of the part based on the initial process scheme according to the characteristic primary rough machining area and the residual area; and finally, providing a self-adaptive evolution mechanism of the rough machining process of the part driven by the process design intention so as to eliminate the problems of interference and under-cut in the rough machining process and ensure the connectivity of the rough machining process situation.
FIG. 2 is a schematic diagram of MR calculation of a characteristic primary rough machining region based on medial axis transformation, and FIG. 2(a) is a characteristic FiCentral axis MA (white line). FIG. 2(b) is a middle axis edge e1When F represents the modeliWhen a rough machining (2 mm radial allowance) is carried out once by using a cutter with the radius of 13mm, e1Is divided into two intersection points v4(s1) And v5(s2) Divided into 3 sections, i.e. v3-v4,v4-v5,v5-v6Wherein v is3-v4And v5-v6Is located at a horizontal line l1Top, indicating that these regions are workable; and v is4-v5At l1Below, it is indicated that these areas are not machinable. Similarly, for other middle axle edge ei(i is 1,2, …,11), and the attribute of the corresponding side of each subinterval is determined. R as shown in FIG. 2c can be constructed from the segmented medial axis edgesiAnd the attribute connection graph G is formed, wherein the attribute of the black line segment is M, and the attribute of the white line segment is U. Obtaining 2 connected subgraphs G through depth traversal G1And G2. FIG. 2d is a schematic representation ofi(i 1,2) with the vertex with the intersection point as the attribute value as the center (e.g., v)2,v4,v3Etc.) build radius R + δrInscribed arc of (F)iIs divided into 2 simple sub-profiles L1And L2. At this time, simple sub-profile L1And L2At deltarIs an offset of the distance to obtain F i2 sub-machining regions C to be cut off at one rough machining stage1And C2I.e. MRi=C1∪C2As shown in fig. 2 e.
FIG. 3 is a schematic diagram showing the classification of UR residues in a process, wherein 4F-UR residues are contained in CaseP2B-URPWith 1C-URP
FIG. 4 shows feature FiSchematic diagram of the primary rough machining process residue calculation. In fig. 4a, the radius of the tool is 0, and 1 machining area C is obtained by the method of step 13(white line), FIG. 4b shows 2 sub-machining regions C with a tool radius R1And C2From equation (2), 7 sub-residual regions (closed black lines) can be calculated, as shown in fig. 4 c. Due to g4Vertex v of18Is an inscribed arc s1Is at the center of the circle of1Belong to SUR4Thus SUR4And g4Correlation, SUR is known from the formula (1)4Is of type B-URP
Fig. 5 is a schematic diagram of longitudinal adaptation of a rough machining process. Suppose F2Is absent, F1After longitudinal adaptation, if non-fillet residuals (e.g., SUR) are still present1) Then a suitable tool T is selected according to equation (3), at T1A secondary rough machining step is added at any moment to cut off the SUR1So that F1To achieve radial allowance (delta) at rough machining stager2) and axial margin (δ)a2) as a processing target.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Two modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims (1)

1. A rough machining procedure self-adaptive dynamic evolution calculation method is characterized by comprising the following steps: the method comprises the following steps:
(a) the complex contour of the feature is divided into a plurality of simple sub-contours through dynamic evolution calculation of a complex cavity rough machining region MR and a residual region UR guided by a feature center axis, and the machining region MR is obtained through the bias of the simple sub-contours;
the method comprises the following specific steps of calculating the dynamic evolution of the MR and the UR of the rough machining region of the complex cavity guided by the characteristic center shaft:
1.1: introducing a local inscribed circle moving representation model, and calculating a primary rough machining region MR based on a characteristic central axis;
the specific calculation steps of the primary rough machining region MR based on the characteristic central axis in the step 1.1 are as follows:
1.1.1: moving representation model r of each local inscribed circle of feature central axis MA (F) of traversing complex cavity feature Fi(s) calculating a local inscribed circle movement representation model ri(s) and inscribed arc radius horizontal line R + deltarIs the center of the inscribed arc and is according to ri(s) and R + deltarComparing the sizes of the characteristic axes MA (F) with the side e of the characteristic axis MA (F)iDivided into multiple inscribed arcs siDividing the machining area constrained by the machining allowance into a zone association area machinable M and a zone association area non-machinable U;
1.1.2: representing the model r by local inscribed circle movementi(s) and its class, construct ri(s) attribute adjacency graph G;
1.1.3: traversing r in step 1.1.2 in depth by taking machinable M of the interval associated region as constrainti(s) obtaining a plurality of connected subgraphs G from the attribute adjacency graph GiEach connectivity sub-graph GiCorresponding to a simple sub-contour, i.e. connected sub-graph GiThe number is consistent with the number of the simple sub-outlines;
1.1.4: according to the connected subgraph G obtained in 1.1.3iOffset distance delta by feature profilerCalculating to obtain a processing region MR required by the complex cavity characteristic F in a primary rough processing stage;
1.2: constructing an attribute adjacency graph according to the local inscribed circle moving representation model in the step 1.1, and classifying and calculating the primary rough machining residual region UR;
the specific calculation steps of the primary rough machining residual region UR in the step 1.2 are as follows:
1.2.1: and (5) traversing r in the depth direction 1.1.2 by taking the unmanageable U of the interval association region as a constrainti(s) obtaining a plurality of connected subgraphs G from the attribute adjacency graph GiEach connectivity sub-graph giCorresponding to a sub-process residual SURiI.e. process residue URP=∪SURi
1.2.2: due to each connectivity graph giDifferent topological structures are provided, and each topological structure corresponds to one type of process residues; thus, according to the connectivity sub-graph g in 1.2.1iThe topology of (1.2.1), the process residues in (1.2.1) are divided into fillet residues F-URPResidual B-UR in single bottle neckPMultiple bottle neck residueReserve MB-URPAnd composite residual C-URP
According to g in said step 1.2.2iThe topology of (2) implements residual region classification, which is expressed as formula (1):
Figure FDA0003248668360000021
wherein a isiDenotes giVertex v in (1)iDegree, | G in Gi| represents attribute adjacency graph giNumber of vertices of (1), T (g)i) Denotes giThe type of the corresponding process residue;
1.2.3: the process residue UR in step 1.2.2PThe range is controlled by the radius R of the tool, when R approaches 0, only the residual UR remainsAIt is recorded as
Figure FDA0003248668360000022
Therefore, URPThe calculation of (d) translates into a machining region MR (0, delta) with a tool radius of 0r) With a machining region MR (R, delta) of tool radius Rr) Between Boolean operations, i.e. URP=MR(0,δr)-MR(R,δr);
1.2.4: according to the inscribed arc s in step 1.1.1iAnd connectivity sub-graph g in 1.2.1iThe corresponding relation between the vertexes of (1) and the establishment of the SUR of the sub-process residueiAnd giThe correlation between the two and the g in 1.2.2iThe process residue classification of the topological structure of (1) judges the SUR of the sub-process residueiType of (d);
the expression method of the local inscribed circle moving expression model is as follows: giving a complex cavity characteristic F, wherein the central axis MA is composed of a plurality of central axis edges eiComposition, i.e. MA (F) { ei},1≤i≤nmWherein n ismCharacteristic middle shaft edge e of complex cavity characteristic FiQuantity, characteristic axial edge eiAt an arbitrary point v oni(s),s∈[smin,smax]Wherein s is a characteristic central axisThe arc length parameters of the point corresponding to the edge are all present by one vi(s) is a circle center Oi(s) local inscribed circle LIC, when LIC is along eiWhile moving, its radius ri(s) changing continuously, i.e. using ri(s) characterization of eiBalance ri(s) is eiThe LIC moving representation model of (1) is recorded as
Figure FDA0003248668360000031
The characteristic central axis of the complex cavity characteristic F is represented as MA (F) ═ ri(s)},1≤i≤nm
Said ri(s) attribute adjacency graph G vertices correspond to
Figure FDA0003248668360000032
The attribute of the endpoint is taken as the attribute of the vertex; r isi(s) Attribute adjacency graph G characteristic Axis edge eiThe corresponding movement representation model is
Figure FDA0003248668360000033
Figure FDA0003248668360000034
Class of which is characteristic of medial axis edge eiThe category comprises a machinable interval associated region M and an unmachined interval associated region U;
1.3: according to the processing region MR and the residual region UR calculated in the steps 1.1 and 1.2, realizing the construction of a dynamic evolution model of the characteristic rough processing region MR and the residual region UR;
the dynamic evolution model of step 1.3 is denoted as Mk,jDynamically evolving model F from multiple featuresi.Mk,jConstitution Fi.Mk,jDriven by parameters reflecting the design intention of the process and realized by the process of removing the margin between adjacent processes in a dynamic evolution model, i.e.
Figure FDA0003248668360000041
Fi.Mk,j=Fi.Mk,j-1-Fi.Mk,j
MR(Fi.Mk,j)=f(Fi.Mk,j-1,R,δr)
Fi.Mk,j=(MR,zt,zb,n)
UR(Fi.Mk,j)=UR(Fi.Mk,j-1)-MR(Fi.Mk,j)
Wherein the parameters comprise a tool T and a radial allowance deltarAxial margin deltaaN is FiIn the axial direction of the tool, MR (F)i.Mk,j) A processing area F (F) obtained by projecting the characteristic dynamic evolution model to an XY plane along the Z-axis direction in a local coordinate systemi.Mk,j-1,R,δr) Constructing a function for the effective machining area, R being the radius of the tool T, zt、zbAre respectively Fi.Mk,jTop and bottom surface constraint of UR (F)i.Mk,j) Representing a residual region obtained by projecting the characteristic dynamic evolution model to an XY plane along the Z-axis direction in a local coordinate system;
(b) a self-adaptive evolution mechanism of a part rough machining procedure driven by a process design intention is provided to eliminate the problems of interference and under-cut in the rough machining process and ensure the connectivity of the rough machining process situation;
the adaptive evolution mechanism of the rough machining process of the parts driven by the process design intention is divided into the following two mechanisms:
the first method comprises the following steps: longitudinal adaptation, in one process step, the complex cavity features F of which are machined in sequence, and each machining operation only machines one feature, given any one machining operation opiSuppose opiFor feature FnOf the machining region MRniProcessing is carried out when feature FmAnd feature FnWhen the formula (2) is satisfied, MR in the processing region is requiredniFeature F from top to bottom before cuttingmThe layers of (2) are left to be cut off, and the formula (2) is as follows:
Figure FDA0003248668360000051
wherein, FmAnd FnRepresenting two features with a longitudinal dependency between them, i.e. Fm→FnAnd which are machined in the roughing stage using different tools, URmj(j-1, 2) represents FmTwo-layer residue of (1), MRniIs represented by FnI th of (1)thThe cutting of the secondary machining area is carried out,
and the second method comprises the following steps: lateral adaptation, assuming t1Each feature of time FmIs URm1For any one of its sub-process residual SURiIf SUR is notiThe type of the (D) is fillet residue, so that the treatment is not needed; if SURiIf the type of (d) is non-fillet residue, then according to its corresponding sub-graph giEach side eiSelecting a tool T satisfying the condition shown in formula (3):
Figure FDA0003248668360000052
where R is the radius of the tool T, where T is along giWith each edge moving without interference, T being a tool, by adding a new roughing step to cut off the SURi
CN201810372957.7A 2018-04-24 2018-04-24 Adaptive dynamic evolution calculation method for rough machining process Active CN108304686B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810372957.7A CN108304686B (en) 2018-04-24 2018-04-24 Adaptive dynamic evolution calculation method for rough machining process

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810372957.7A CN108304686B (en) 2018-04-24 2018-04-24 Adaptive dynamic evolution calculation method for rough machining process

Publications (2)

Publication Number Publication Date
CN108304686A CN108304686A (en) 2018-07-20
CN108304686B true CN108304686B (en) 2021-12-10

Family

ID=62846302

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810372957.7A Active CN108304686B (en) 2018-04-24 2018-04-24 Adaptive dynamic evolution calculation method for rough machining process

Country Status (1)

Country Link
CN (1) CN108304686B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110457735B (en) * 2019-06-11 2023-07-04 上海航天精密机械研究所 Coarse machining unit calculation method for complex groove cavity characteristics
CN112099434A (en) * 2020-09-21 2020-12-18 无锡职业技术学院 Contour numerical control milling method and system
CN113468800A (en) * 2021-05-27 2021-10-01 河海大学 Complex cavity cutter combination optimization method based on genetic algorithm
CN113591235A (en) * 2021-06-29 2021-11-02 河海大学 Multi-level structured numerical control process model construction method based on dynamic processing characteristics
CN113962511A (en) * 2021-09-07 2022-01-21 河海大学 Process data structured representation method fusing process design intents

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003006283A2 (en) * 2001-07-13 2003-01-23 Knorr-Bremse Systeme für Nutzfahrzeuge GmbH Method and device for detecting residual pressure in a braking system, in particular in a braking cylinder
CN101916085A (en) * 2010-08-05 2010-12-15 华南理工大学 Control system and method of glue-injecting motor and glue-smelting motor of electric injection molding machine
CN102831283A (en) * 2012-09-19 2012-12-19 河海大学常州校区 Complicated product model construction method based on surface feature
CN102922013A (en) * 2012-10-25 2013-02-13 南京航空航天大学 Cavity-characteristic-based high-efficiency rough machining method of aircraft structural part
CN104008181A (en) * 2014-06-09 2014-08-27 中国电子科技集团公司第十四研究所 A retrieval method of similar numerical control technics of electronic parts based on characters of a three-dimensional model
CN104281098A (en) * 2014-10-27 2015-01-14 南京航空航天大学 Modeling method for dynamic machining features of complex curved surface
EP3242177A1 (en) * 2016-04-29 2017-11-08 Rolls-Royce plc Adaptive repair method for aerofoil blades

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003006283A2 (en) * 2001-07-13 2003-01-23 Knorr-Bremse Systeme für Nutzfahrzeuge GmbH Method and device for detecting residual pressure in a braking system, in particular in a braking cylinder
CN101916085A (en) * 2010-08-05 2010-12-15 华南理工大学 Control system and method of glue-injecting motor and glue-smelting motor of electric injection molding machine
CN102831283A (en) * 2012-09-19 2012-12-19 河海大学常州校区 Complicated product model construction method based on surface feature
CN102922013A (en) * 2012-10-25 2013-02-13 南京航空航天大学 Cavity-characteristic-based high-efficiency rough machining method of aircraft structural part
CN104008181A (en) * 2014-06-09 2014-08-27 中国电子科技集团公司第十四研究所 A retrieval method of similar numerical control technics of electronic parts based on characters of a three-dimensional model
CN104281098A (en) * 2014-10-27 2015-01-14 南京航空航天大学 Modeling method for dynamic machining features of complex curved surface
EP3242177A1 (en) * 2016-04-29 2017-11-08 Rolls-Royce plc Adaptive repair method for aerofoil blades

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Computation of the medial axis of planar domains based on saddle point programming;Lixin Gao 等;《Computer-Aided Design》;20110831;第43卷(第8期);第1020-1037页 *
Local feature extraction and matching partial objects;Dmitriy Bespalov 等;《Computer-Aided Design》;20060930;第38卷(第9期);第979-988页 *
基于中轴变换的槽腔特征相似性评价方法;徐昌鸿 等;《计算机集成制造系统》;20140325;第20卷(第12期);第3066-3074页 *
基于层切法的模具自适应数控粗加工研究;刘伟军 等;《工具技术》;20050428;第21-23、44页 *
数控铣床加工模具内腔的几个问题;靳玉成;《黑龙江科技信息》;20071008;第28页 *
相似特征数控工艺的重用评价方法;张曦卯 等;《工具技术》;20140720;第48卷(第7期);第20-24页 *

Also Published As

Publication number Publication date
CN108304686A (en) 2018-07-20

Similar Documents

Publication Publication Date Title
CN108304686B (en) Adaptive dynamic evolution calculation method for rough machining process
US4837703A (en) Method for generating tool path
CN105354353A (en) MBD model based processing feature identification and modeling method
US7110849B2 (en) Horizontally-structured CAD/CAM modeling for virtual fixture and tooling processes
CN107562015B (en) Process geometric model construction method based on numerical control machining programming
CN102436217B (en) Method for reconstructing web processing drive geometry of slot characteristic of aircraft structure
Kailash et al. A volume decomposition approach to machining feature extraction of casting and forging components
Cai et al. Freeform machining feature recognition with manufacturability analysis
Liu et al. Extracting the turning volume and features for a mill/turn part with multiple extreme faces
CN110412938B (en) Multi-cavity structural part processing step sequence decision method based on transition feature simplification
Flutter et al. A machining strategy for toolmaking
CN110457735B (en) Coarse machining unit calculation method for complex groove cavity characteristics
US20020133803A1 (en) Enhancement to horizontally-structured CAD/CAM modeling
CN113962511A (en) Process data structured representation method fusing process design intents
Zhou et al. An automated CNC programming approach to machining pocket with complex islands and boundaries by using multiple cutters in hybrid tool path patterns
CN107301310B (en) Method for constructing three-dimensional process template of mechanical part
Huang et al. Multilevel structured NC machining process model based on dynamic machining feature for process reuse
Han et al. A cutter selection method for 2 1/2-axis trochoidal milling of the pocket based on optimal skeleton
Yıldız et al. Development of a feature based CAM system for rotational parts
Chan et al. A high-efficiency rough milling strategy for mould core machining
Rui et al. An effective adaptive dynamic evolution computing approach of roughing process for part with complex pockets
Dai et al. A Persistent Naming Discrimination Method Based on the Sweeping Direction
JP3745526B2 (en) Roughing intermediate shape model generator for NC data creation
CN110653401B (en) Cavity partition plunge milling cutter path planning method based on maximum cutting wrap angle constraint
Tolouei-Rad et al. A hybrid approach to automatic generation of NC programs

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant