CN111061218A - Complex curved surface machining cutting force prediction method and device based on ACIS - Google Patents
Complex curved surface machining cutting force prediction method and device based on ACIS Download PDFInfo
- Publication number
- CN111061218A CN111061218A CN201911413493.0A CN201911413493A CN111061218A CN 111061218 A CN111061218 A CN 111061218A CN 201911413493 A CN201911413493 A CN 201911413493A CN 111061218 A CN111061218 A CN 111061218A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- cutter
- boundary
- cutting edge
- force
- 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.)
- Granted
Links
- 238000005520 cutting process Methods 0.000 title claims abstract description 122
- 238000000034 method Methods 0.000 title claims abstract description 44
- 241001502050 Acis Species 0.000 title claims abstract description 26
- 238000000627 alternating current impedance spectroscopy Methods 0.000 title claims abstract description 25
- 238000003754 machining Methods 0.000 title claims abstract description 20
- 238000007654 immersion Methods 0.000 claims abstract description 39
- 239000007787 solid Substances 0.000 claims abstract description 7
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 6
- 230000008569 process Effects 0.000 claims description 18
- 238000006243 chemical reaction Methods 0.000 claims description 14
- 238000003801 milling Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 239000002826 coolant Substances 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011438 discrete method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/19—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/35—Nc in input of data, input till input file format
- G05B2219/35349—Display part, programmed locus and tool path, traject, dynamic locus
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Numerical Control (AREA)
Abstract
The invention belongs to the field of computer numerical control machining, and discloses a complex curved surface machining cutting force prediction method and device based on ACIS. The method comprises the following steps: adopting ACIS geometric modeling to construct a solid model of the envelope surface of the workpiece and the cutter; obtaining an engagement area and an engagement boundary between the workpiece and the tool enveloping body by using a Boolean intersection algorithm on the solid model; converting the meshing boundary under the global coordinate system into a meshing boundary under a cutter coordinate system; dispersing the meshing boundary of the cutter and the workpiece in the z-axis direction of a cutter coordinate system, and obtaining the radial immersion angle range of each infinitesimal position obtained by dispersion; the cutting force of each element is calculated according to the radial immersion angle, and then the elements in the meshing boundary are summed to obtain the cutting force of the whole cutter on the workpiece. The method can accurately obtain the meshing boundary and realize the prediction of the cutting force of the complex curved surface.
Description
Technical Field
The invention belongs to the field of computer numerical control machining, and particularly relates to a complex curved surface machining cutting force prediction method and device based on ACIS.
Background
Over the past century, companies and research institutes have made tremendous efforts to improve the performance of machining. By optimizing the tool geometry, cutting material, coolant and coating, tool life and machining performance are improved. However, there is no mature solution for predicting cutting force, especially for predicting cutting force of complex curved surface machining under a multi-axis machine tool.
Cutting forces directly affect the shape error and surface quality of the machined part. Once the cutting forces can be simulated, shape errors and vibrations caused by tool and workpiece deflections, as well as problems affecting surface quality, can be effectively predicted. The cutting force depends on the material of the tool and the workpiece, the geometry of the tool and the cutting conditions. In a great deal of previous research, the cutting force calculation of a smooth simple plane can be calculated by using an analytical method, a very mature theoretical result is formed, and commercial cutting force prediction software is introduced. However, for complex curved surfaces, due to the irregularity of the complex curved surfaces, it is difficult to obtain an accurate engagement boundary between the tool and the workpiece, so that the analytic method has not been well applicable in the calculation of the cutting force of the complex curved surfaces.
Disclosure of Invention
Aiming at the defects or improvement requirements of the prior art, the invention provides a complex curved surface machining cutting force prediction method and equipment based on ACIS, and aims to obtain a cutting boundary, a meshing area and a meshing boundary, perform discretization differential processing and calculate instantaneous cutting force by modeling the meshing condition of a tool and a workpiece in any contact state through three-dimensional model simulation by virtue of the modeling advantage of an ACIS engine, so as to realize the prediction of the complex curved surface cutting force.
To achieve the above object, according to an aspect of the present invention, there is provided an ACIS-based cutting force prediction method for complex curved surface machining, comprising the steps of:
(1) adopting ACIS geometric modeling to construct a solid model of the envelope surface of the workpiece and the cutter;
(2) on the established entity model of the workpiece and the tool enveloping surface, obtaining a meshing area between the workpiece and the tool enveloping body by using a Boolean intersection algorithm, and extracting the boundary of the meshing area so as to obtain the meshing boundary of the tool and the workpiece under a global coordinate system;
(3) converting the meshing boundary under the global coordinate system into a meshing boundary under a cutter coordinate system, wherein the cutter coordinate system takes a cutter location point as an original point, a cutter axis vector as a z-axis, and a rotating shaft of a cutter forward inclination angle and a rotating shaft of a cutter side inclination angle as an x-axis and a y-axis respectively;
(4) discretizing the meshing boundary of the cutter and the workpiece in the z-axis direction of the cutter coordinate system, and obtaining the radial immersion angle range of each infinitesimal position obtained by discretization;
(5) the cutting force of each element is calculated according to the radial immersion angle, and then the elements in the meshing boundary are summed to obtain the cutting force of the whole cutter on the workpiece.
Further, taking a world coordinate system XYZ as a global coordinate system, and modeling in the step (1) under the world coordinate system XYZ;
setting a process coordinate system FCN and a tool coordinate system xyz in the step (3), wherein:
the process coordinate system FCN takes the cutter location point as an origin, and an F axis, a C axis and an N axis are respectively parallel to an X axis, a Y axis and a Z axis of the world coordinate system;
the tool coordinate system xyz is obtained by rotating the process coordinate system FCN by the same angle around the axis C according to the size of the tool rake angle and then rotating by the same angle around the axis F according to the size of the tool side rake angle.
Further, the coordinate system is converted into an expression:
wherein T is a transformation matrix, l is a cutter rake angle, and T is a cutter roll angle.
Further, in step (4), the method for solving the radial immersion angle range is as follows:
setting a infinitesimal step length dz according to the meshing boundary outline under a tool coordinate system, decomposing a meshing region into a plurality of arc-shaped infinitesimals parallel to an xy plane along a z axis according to upper and lower boundaries zLow and zHigh of the meshing region in the z direction, and determining the z-axis coordinate z of the ith infinitesimaliThe expression of (a) is as follows:
zi=zLOW+i×dz (3)
in the formula,and i is an integer, zLow is the lower boundary of the meshing boundary, and zHigh is the upper boundary of the meshing boundary;
obtaining cut-in boundary and cut-out boundary on the meshing boundary contour through simulation, and obtaining coordinate z of ith infinitesimaliAfter, in ziMaking and cutting into the boundary at the location andcutting the intersection point of the boundary to obtain the corresponding radial immersion angle at the heightThe upper and lower bounds of (1) are startAngle and endAngle, respectively.
Further, in the step (5),
in the case of each instantaneous cut, the cutting force is only generated if the cutting edge is in the region of engagement, so that for each cutting edge at each infinitesimal position ziRadial immersion angle ofIf it is judged thatAt [ startAngle, endAngle ]]In the range, the following calculation formula of infinitesimal force is used for solving the tangential, radial and axial forces of the cutting edge;
wherein, Ktc、Krc、Kac、Kte、Kre、KaeIs the coefficient of cut; the subscript j indicates the jth cutting edge on the tool; subscripts t, r, a are the t, r and a axes of a cutting edge coordinate system tra established in tangential, radial and axial directions of the cutting edge;is the cutting thickness h of the jth cutting edgejWith respect to z-coordinate ziAnd radial immersion angleA function of (a);
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle ofInstantaneous tangential force in time;
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle ofInstantaneous radial force in time;
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle ofInstantaneous axial force in time;
then, the following formula is used for conversion, and tangential, radial and axial forces of the cutting edge are converted into forces under a tool coordinate system;
wherein, TxyzThe matrix is a conversion matrix of a cutting edge coordinate system and a cutter coordinate system during forward milling or reverse milling;
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle ofInstantaneous tangential force in time;
is the j-th cutting edge under the coordinate system of the toolAt infinitesimal position ziHas a radial immersion angle ofInstantaneous radial force in time;
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle ofInstantaneous axial force in time;
and integrating the instantaneous tangential force, the instantaneous radial force and the instantaneous axial force on all the infinitesimal elements to obtain the cutting force of the whole cutter on the workpiece at the moment.
Further, TxyzIs a conversion matrix of a cutting edge coordinate system and a cutter coordinate system during forward milling:
wherein, KiFor the cutting edge at the infinitesimal position ziThe axial immersion angle is calculated as follows:
To achieve the above object, the present invention also provides a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the method as described in any of the preceding claims.
In order to achieve the above object, the present invention further provides an ACIS-based cutting force prediction apparatus for complex curved surface machining, comprising the computer-readable storage medium as described above and a processor for calling and processing a computer program stored in the computer-readable storage medium.
In general, compared with the prior art, the above technical solution contemplated by the present invention can obtain the following beneficial effects:
(1) geometric and position information of the workpiece and the cutter enveloping surface is restored by ACIS modeling in a computer, an actual cutting process is simulated, the engagement boundary of the workpiece and the cutter at any moment can be obtained by solid Boolean operation, and the cutting force of the complex curved surface is calculated based on a discretization method, so that the problem that the engagement boundary is difficult to accurately obtain for the complex curved surface in the prior art is solved.
(2) Through the conversion of the world coordinate system and the tool coordinate system, the relation between the meshing area and the tool can be accurately obtained, so that the radial immersion angle of the tool can be determined.
(3) By dispersing the meshing area, the upper and lower bounds of the radial immersion angle at each infinitesimal position in the cutting process can be accurately obtained, and the judgment of the instantaneous cutting position is facilitated.
(4) Through the decomposition of the instantaneous cutting force under the cutting edge coordinate system and the conversion of the instantaneous cutting force under the cutter coordinate system, the cutting force components of each infinitesimal under the cutter coordinate system are obtained through solving, finally, the total cutting force of the whole meshing area can be obtained through integration, and the prediction of the cutting force of the cutter at any moment is realized.
Drawings
FIG. 1 is a flow chart of the cutting force prediction method for ACIS-based complex curved surface machining according to the present invention;
FIG. 2 is a block diagram of a solid model of a workpiece and tool created using ACIS;
FIG. 3 is a diagram of the area of engagement of the workpiece and the tool using a Boolean algorithm;
FIG. 4 is a schematic representation of the rake and roll angles;
FIG. 5 is a schematic diagram of a coordinate system and engagement boundaries;
FIG. 6 is a schematic view of the engagement boundary abutting the tool surface;
FIG. 7 is a schematic diagram of a meshing boundary cut-in and cut-out position in a tool coordinate system;
FIG. 8 is a schematic view of a discrete method of cutting force modeling.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. In addition, the technical features involved in the embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
As shown in fig. 1, the method for predicting the cutting force of the present invention comprises the steps of:
1. constructing a solid model of the envelope surface of the workpiece and the cutter in a computer by utilizing an ACIS geometric modeling engine;
3Dmodeler (acis) is a three-dimensional geometric modeling engine with higher precision, supports a robust 3D modeling function, and is an ideal tool for constructing an application program with a hybrid modeling function. The work of modeling the envelope surfaces of the workpiece and tool can be performed using the SDK of the ACIS C + + version, the modeling portion of the example being shown in FIG. 2.
2. On the established entity model, obtaining a meshing area between the workpiece and the tool enveloping body by using a Boolean intersection algorithm, and extracting the boundary to obtain the meshing boundary of the tool workpiece under the global coordinate system;
in ACIS, all models are represented using the three-dimensional B-Rep method. The basic idea of the B-Rep method is to define an entity as a finite space enclosed by closed boundary surfaces, which can be represented by a subset of faces, which are defined by edges, edges by points, points and finally by three coordinate values. In the underlying data structure, the topological graph is used for representation, so the boolean intersection algorithm of the B-Rep method is to calculate the common part between the topological graphs of the two entity model underlayers, and the boolean operation result of the example is shown in fig. 3. And the boundary of the engagement region of the sheet is the engagement boundary.
3. Mapping the meshing boundary under the global coordinate system to the tool coordinate system
In total, three coordinate systems are established throughout the process, which are described below as the world coordinate system, the process coordinate system and the tool coordinate system.
The world coordinate system XYZ is a fixed coordinate system, and the whole modeling process is completed under the world coordinate system in step 1.
The process coordinate system FCN is an intermediate coordinate system established for the sake of conversion, wherein the origin of the process coordinate system is at the center of the sphere of the ball nose tool, and the three axes FCN are respectively parallel to the XYZ axes of the world coordinate system.
The tool coordinate system xyz origin is also at the tool position, and is obtained by rotating the process coordinate system around the C axis by the anteversion angle and then rotating the process coordinate system around the F axis by the laterial inclination angle.
A schematic of the process and tool coordinate systems is shown in fig. 4, and a schematic of the rake and roll angles is shown in fig. 5.
During the process of moving the ball end mill relative to a workpiece coordinate system (namely a world coordinate system), the ball end mill predicts the original point of the action point of milling force in a process coordinate system FCN, so that the action point of the predicted force makes relative translation motion in the world coordinate system. Meanwhile, because the force is a vector, the magnitude and the direction of the force only need to be considered, only the coordinate transformation from the process coordinate system to the tool coordinate system needs to be calculated, and because the original points are consistent, only two rotations need to be carried out. The coordinate system is transformed into the expression:
wherein T is a transformation matrix, l is a cutter rake angle, and T is a cutter roll angle.
Three coordinate systems are established in ACIS by the above method and the position of the engagement boundary in space is shown, as shown in fig. 6. The coordinates of the meshing boundary in space are converted into coordinates in the tool coordinate system using the above conversion formula, and thus can be attached to the tool surface, as shown in fig. 7.
4. The coordinate positions of the cutting-in position and the cutting-out position in the tool coordinate system can be obtained through the above process, as shown in fig. 8. Because of the general discrete cutting force model, the cutting force on each infinitesimal element is calculated firstly, and the cutting forces on all the infinitesimal elements are superposed, thereby obtaining the cutting force on the whole cutter. In this example, the selected discretization method discretizes the z-direction in the tool coordinate system as shown in FIG. 8. First, the upper and lower limits of the meshing boundary in the z-direction, such as zLow and zHigh in fig. 8, need to be obtained, and then the infinitesimal step length dz is set, so that the z-axis coordinate z of the ith infinitesimal isiThe expression of (a) is as follows:
zi=zLOW+i×dz (3)
in the formula,
i is an integer, zLow is the lower boundary of the meshing boundary, and zHigh is the upper boundary of the meshing boundary;
obtaining cut-in boundary and cut-out boundary on the meshing boundary contour through simulation, and obtaining coordinate z of ith infinitesimaliAfter, in ziMaking a plane parallel to the xy plane at the position to obtain the intersection point of the xy plane and the cut-in boundary and the cut-out boundary so as to obtain the radial immersion angle corresponding to the heightThe upper and lower bounds of (1) are startAngle and endAngle, respectively.
(5) Cutting at each instantIn the case of (2), the cutting force is generated only when the cutting edge is in the engagement region,
wherein, Ktc、Krc、Kac、Kte、Kre、KaeThe cutting coefficient can be obtained through a calibration test according to different cutter workpiece materials; the subscript j indicates the jth cutting edge on the tool; subscripts t, r, a are the t, r and a axes of a cutting edge coordinate system tra established in tangential, radial and axial directions of the cutting edge;is the cutting thickness h of the jth cutting edgejWith respect to z-coordinate ziAnd radial immersion angleA function of (a);
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle ofInstantaneous tangential force in time;
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle ofInstantaneous radial force in time;
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiRadial immersion angle ofIs composed ofInstantaneous axial force in time;
then, the following formula is used for conversion, and tangential, radial and axial forces of the cutting edge are converted into forces under a tool coordinate system;
wherein, TxyzThe matrix is a conversion matrix of a cutting edge coordinate system and a cutter coordinate system during forward milling or reverse milling;
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle ofInstantaneous tangential force in time;
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle ofInstantaneous radial force in time;
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle ofInstantaneous axial force in time;
and integrating the instantaneous tangential force, the instantaneous radial force and the instantaneous axial force on all the infinitesimal elements to obtain the cutting force of the whole cutter on the workpiece at the moment.
Taking the straight milling as an example, the conversion matrix of the cutting edge coordinate system and the tool coordinate system during the straight milling is as follows:
wherein, KiFor the cutting edge at the infinitesimal position ziThe axial immersion angle is calculated as follows:
The forward milling and the reverse milling have corresponding adjustment in the form of a conversion matrix due to different tool turning directions, but the overall calculation method is the same as that of the forward milling and is not repeated.
After the forces on all the infinitesimal elements are obtained, the forces on all the infinitesimal elements are integrated, and then the cutting force of the whole cutter on the workpiece at the moment can be obtained.
It will be understood by those skilled in the art that the foregoing is only a preferred embodiment of the present invention, and is not intended to limit the invention, and that any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (8)
1. The method for predicting the cutting force of the complex curved surface machining based on the ACIS is characterized by comprising the following steps of:
(1) adopting ACIS geometric modeling to construct a solid model of the envelope surface of the workpiece and the cutter;
(2) on the established entity model of the workpiece and the tool enveloping surface, obtaining a meshing area between the workpiece and the tool enveloping body by using a Boolean intersection algorithm, and extracting the boundary of the meshing area so as to obtain the meshing boundary of the tool and the workpiece under a global coordinate system;
(3) converting the meshing boundary under the global coordinate system into a meshing boundary under a cutter coordinate system, wherein the cutter coordinate system takes a cutter location point as an original point, a cutter axis vector as a z-axis, and a rotating shaft of a cutter forward inclination angle and a rotating shaft of a cutter side inclination angle as an x-axis and a y-axis respectively;
(4) discretizing the meshing boundary of the cutter and the workpiece in the z-axis direction of the cutter coordinate system, and obtaining the radial immersion angle range of each infinitesimal position obtained by discretization;
(5) the cutting force of each element is calculated according to the radial immersion angle, and then the elements in the meshing boundary are summed to obtain the cutting force of the whole cutter on the workpiece.
2. An ACIS-based cutting force prediction method for complex curved surface machining based on the claim 1, characterized in that, the world coordinate system XYZ is taken as the global coordinate system, and modeling is performed in the step (1) under the world coordinate system XYZ;
setting a process coordinate system FCN and a tool coordinate system xyz in the step (3), wherein:
the process coordinate system FCN takes the cutter location point as an origin, and an F axis, a C axis and an N axis are respectively parallel to an X axis, a Y axis and a Z axis of the world coordinate system;
the tool coordinate system xyz is obtained by rotating the process coordinate system FCN by the same angle around the axis C according to the size of the tool rake angle and then rotating by the same angle around the axis F according to the size of the tool side rake angle.
4. An ACIS-based cutting force prediction method for complex curved surface machining according to claim 2 or 3, wherein in step (4), the method for solving the radial immersion angle range is as follows:
setting a infinitesimal step length dz according to the meshing boundary outline under a tool coordinate system, decomposing a meshing region into a plurality of arc-shaped infinitesimals parallel to an xy plane along a z axis according to upper and lower boundaries zLow and zHigh of the meshing region in the z direction, and determining the z-axis coordinate z of the ith infinitesimaliThe expression of (a) is as follows:
zi=zLOW+i×dz (3)
in the formula,and i is an integer, zLow is the lower boundary of the meshing boundary, and zHigh is the upper boundary of the meshing boundary;
obtaining cut-in boundary and cut-out boundary on the meshing boundary contour through simulation, and obtaining coordinate z of ith infinitesimaliAfter, in ziMaking the intersection point of the cut-in boundary and the cut-out boundary at the position to obtain the corresponding radial immersion angle at the heightThe upper and lower bounds of (1) are startAngle and endAngle, respectively.
5. An ACIS-based cutting force prediction method for complex curved surface machining according to claim 4, wherein in step (5),
in the case of each instantaneous cut, the cutting force is only generated if the cutting edge is in the region of engagement, so that for each cutting edge at each infinitesimal position ziRadial immersion angle ofIf it is judged thatAt [ startAngle, endAngle ]]In the range, the following calculation formula of infinitesimal force is used for solving the tangential, radial and axial forces of the cutting edge;
wherein, Ktc、Krc、Kac、Kte、Kre、KaeIs the coefficient of cut; the subscript j indicates the jth cutting edge on the tool; subscripts t, r, a are the t, r and a axes of a cutting edge coordinate system tra established in tangential, radial and axial directions of the cutting edge;is the cutting thickness h of the jth cutting edgejWith respect to z-coordinate ziAnd radial immersion angleA function of (a);
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle ofInstantaneous tangential force in time;
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle ofInstantaneous radial force in time;
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle ofInstantaneous axial force in time;
then, the following formula is used for conversion, and tangential, radial and axial forces of the cutting edge are converted into forces under a tool coordinate system;
wherein, TxyzThe matrix is a conversion matrix of a cutting edge coordinate system and a cutter coordinate system during forward milling or reverse milling;
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle ofInstantaneous tangential force in time;
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle ofInstantaneous radial force in time;
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle ofInstantaneous axial force in time;
and integrating the instantaneous tangential force, the instantaneous radial force and the instantaneous axial force on all the infinitesimal elements to obtain the cutting force of the whole cutter on the workpiece at the moment.
6. The ACIS-based cutting force prediction method for complex curved surface machining according to claim 5, wherein T isxyzIs a conversion matrix of a cutting edge coordinate system and a cutter coordinate system during forward milling:
wherein, KiFor the cutting edge at the infinitesimal position ziThe axial immersion angle is calculated as follows:
7. A computer-readable storage medium, having stored thereon a computer program which, when executed by a processor, implements the method of any one of claims 1 to 6.
8. An ACIS-based complex machining cutting force prediction apparatus comprising the computer-readable storage medium of claim 7 and a processor for invoking and processing the computer program stored in the computer-readable storage medium.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911413493.0A CN111061218B (en) | 2019-12-31 | 2019-12-31 | Complex curved surface machining cutting force prediction method and device based on ACIS |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911413493.0A CN111061218B (en) | 2019-12-31 | 2019-12-31 | Complex curved surface machining cutting force prediction method and device based on ACIS |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111061218A true CN111061218A (en) | 2020-04-24 |
CN111061218B CN111061218B (en) | 2021-07-27 |
Family
ID=70306033
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911413493.0A Active CN111061218B (en) | 2019-12-31 | 2019-12-31 | Complex curved surface machining cutting force prediction method and device based on ACIS |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111061218B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112464399A (en) * | 2020-11-12 | 2021-03-09 | 华中科技大学 | Five-axis machining cutter envelope surface calculation method based on motion vector |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4138934A (en) * | 1977-05-31 | 1979-02-13 | Garden City Envelope Company | Envelope forming assemblies |
CN102566492A (en) * | 2012-01-13 | 2012-07-11 | 华中科技大学 | Method for forecasting maximum milling force for plunge milling of metal difficult-to-cut materials |
CN102581364A (en) * | 2012-03-05 | 2012-07-18 | 华中科技大学 | Method for calculating milling load of ball-end milling cutter |
CN102622489A (en) * | 2012-03-26 | 2012-08-01 | 上海交通大学 | Five-axis side milling cutting force predicting method based on ACIS platform |
CN103221587A (en) * | 2010-11-24 | 2013-07-24 | 康宁股份有限公司 | Methods of making an unsupported article of a semiconducting material using thermally active molds |
EP2919937A2 (en) * | 2012-11-16 | 2015-09-23 | Marcel Sobczyk | Method for determining the flank face contour of a hob peeling tool, hob peeling tool and use thereof |
CN105069257A (en) * | 2015-09-01 | 2015-11-18 | 华中科技大学 | Free-form surface micro-milling cutting force modeling method |
CN105242546A (en) * | 2015-11-04 | 2016-01-13 | 大连理工大学 | High-speed milling force modeling method based on material characteristics |
CN105701323A (en) * | 2016-03-23 | 2016-06-22 | 华中科技大学 | Plunge milling machining cutting force predicting and modeling method |
CN106808320A (en) * | 2016-12-12 | 2017-06-09 | 华中科技大学 | A kind of Prediction Method of Milling Forces for considering cutter angle between teeth, helical angle and bias |
CN107168245A (en) * | 2017-05-04 | 2017-09-15 | 武汉理工大学 | A kind of accurate Forecasting Methodology of chamfered edge circular bit cutting force for considering cutting edge effect |
CN108255134A (en) * | 2017-12-15 | 2018-07-06 | 武汉理工大学 | A kind of difficult-to-machine material high-speed turning prediction of Turning Force with Artificial method for considering chamfered edge geometry |
CN108563848A (en) * | 2018-03-27 | 2018-09-21 | 华中科技大学 | A kind of milling force modeling method of flat helical end mills |
CN109754332A (en) * | 2019-01-08 | 2019-05-14 | 西北工业大学 | The energy consumption model modeling method of lathe Milling Processes based on cutting force |
-
2019
- 2019-12-31 CN CN201911413493.0A patent/CN111061218B/en active Active
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4138934A (en) * | 1977-05-31 | 1979-02-13 | Garden City Envelope Company | Envelope forming assemblies |
CN103221587A (en) * | 2010-11-24 | 2013-07-24 | 康宁股份有限公司 | Methods of making an unsupported article of a semiconducting material using thermally active molds |
CN102566492A (en) * | 2012-01-13 | 2012-07-11 | 华中科技大学 | Method for forecasting maximum milling force for plunge milling of metal difficult-to-cut materials |
CN102581364A (en) * | 2012-03-05 | 2012-07-18 | 华中科技大学 | Method for calculating milling load of ball-end milling cutter |
CN102622489A (en) * | 2012-03-26 | 2012-08-01 | 上海交通大学 | Five-axis side milling cutting force predicting method based on ACIS platform |
CN102622489B (en) * | 2012-03-26 | 2014-01-15 | 上海交通大学 | Five-axis side milling cutting force predicting method based on ACIS platform |
US20190193176A1 (en) * | 2012-11-16 | 2019-06-27 | Profilator Gmbh & Co. Kg | Method for determining the flank face contour of a gear skiving tool, gear skiving tool and use thereof |
EP2919937A2 (en) * | 2012-11-16 | 2015-09-23 | Marcel Sobczyk | Method for determining the flank face contour of a hob peeling tool, hob peeling tool and use thereof |
CN105069257A (en) * | 2015-09-01 | 2015-11-18 | 华中科技大学 | Free-form surface micro-milling cutting force modeling method |
CN105242546A (en) * | 2015-11-04 | 2016-01-13 | 大连理工大学 | High-speed milling force modeling method based on material characteristics |
CN105701323A (en) * | 2016-03-23 | 2016-06-22 | 华中科技大学 | Plunge milling machining cutting force predicting and modeling method |
CN106808320A (en) * | 2016-12-12 | 2017-06-09 | 华中科技大学 | A kind of Prediction Method of Milling Forces for considering cutter angle between teeth, helical angle and bias |
CN107168245A (en) * | 2017-05-04 | 2017-09-15 | 武汉理工大学 | A kind of accurate Forecasting Methodology of chamfered edge circular bit cutting force for considering cutting edge effect |
CN108255134A (en) * | 2017-12-15 | 2018-07-06 | 武汉理工大学 | A kind of difficult-to-machine material high-speed turning prediction of Turning Force with Artificial method for considering chamfered edge geometry |
CN108563848A (en) * | 2018-03-27 | 2018-09-21 | 华中科技大学 | A kind of milling force modeling method of flat helical end mills |
CN109754332A (en) * | 2019-01-08 | 2019-05-14 | 西北工业大学 | The energy consumption model modeling method of lathe Milling Processes based on cutting force |
Non-Patent Citations (5)
Title |
---|
J.GRADISEK: ""Mechanistic identification of specific force coefficients for a general end mill"", 《INTERNATIONAL JOURNAL OF MACHINE TOOLS AND MANUFACTURE》 * |
LARUE A: ""Simulation of flank milfing processes"", 《INTERNATIONAL JOURNAL OF MACHINE TOOLS AND MANUFACTURE》 * |
付中涛: ""基于切削力预测模型的复杂曲面铣削进给速度优化研究"", 《中国博士学位论文全文数据库-工程科技I辑》 * |
周林: ""曲面微细铣削让刀误差预测、补偿与精度评价"", 《中国博士学位论文全文数据库工程科技I辑》 * |
张卫青: ""基于ACIS几何平台的复杂多曲面数控加工编程技术"", 《机械设计与制造》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112464399A (en) * | 2020-11-12 | 2021-03-09 | 华中科技大学 | Five-axis machining cutter envelope surface calculation method based on motion vector |
CN112464399B (en) * | 2020-11-12 | 2022-04-22 | 华中科技大学 | Five-axis machining cutter envelope surface calculation method based on motion vector |
Also Published As
Publication number | Publication date |
---|---|
CN111061218B (en) | 2021-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bedi et al. | Advanced interpolation techniques for NC machines | |
CN102528555B (en) | Geometry and mechanics integrated optimization information processing method of interference-free tool path on complex curved surface | |
EP2844427B1 (en) | High performance multi-axis milling | |
Imani et al. | Geometric simulation of ball-end milling operations | |
Li et al. | Interference-free inspection path generation for impeller blades using an on-machine probe | |
WO2005049284A1 (en) | Method for optimising the performance of a robot | |
Hao et al. | Tool path transplantation method for adaptive machining of large-sized and thin-walled free form surface parts based on error distribution | |
CN110989490B (en) | Method for acquiring optimal installation position of workpiece based on contour error | |
Ding et al. | Tool path generation and simulation of dynamic cutting process for five-axis NC machining | |
Sun et al. | Iso-planar feed vector-fields-based streamline tool path generation for five-axis compound surface machining with torus-end cutters | |
CN111061218B (en) | Complex curved surface machining cutting force prediction method and device based on ACIS | |
Leirmo et al. | Tolerancing from STL data: a legacy challenge | |
Lu et al. | Collision-free and smooth path generation method of five-axis ball-end milling considering effective tool diameter change | |
Si et al. | A solid-discrete-based method for extracting the cutter-workpiece engagement in five-axis flank milling | |
Gong et al. | Tool orientation optimization method based on ruled surface using genetic algorithm | |
Çakır et al. | Path planning for industrial robot milling applications | |
Rivière-Lorphèvre et al. | Cutting force prediction in robotic machining | |
EP4036520B1 (en) | Automatic generation of probe path for surface inspection and part alignment | |
CN114131605B (en) | Automatic registration adjustment device and method for ship part polishing track | |
Lu et al. | Geometric deviation evaluation for a five-axis flank milling tool path using the tool swept envelope | |
CN113276112B (en) | Mobile double-robot-based weak rigid member machining process planning method | |
Sarma et al. | An integrated system for NC machining of multi-patch surfaces | |
Wang et al. | Feasible spindle speed interval identification method for large aeronautical component robotic milling system | |
Inui et al. | Cutter engagement feature extraction by using dexel representation solid model | |
US20240176322A1 (en) | Machining Fixture Force Balancing |
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 |