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 PDF

Info

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
Application number
CN201911413493.0A
Other languages
Chinese (zh)
Other versions
CN111061218B (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.)
Huazhong University of Science and Technology
AECC Guizhou Liyang Aviation Power Co Ltd
Original Assignee
Huazhong University of Science and Technology
AECC Guizhou Liyang Aviation Power Co Ltd
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 Huazhong University of Science and Technology, AECC Guizhou Liyang Aviation Power Co Ltd filed Critical Huazhong University of Science and Technology
Priority to CN201911413493.0A priority Critical patent/CN111061218B/en
Publication of CN111061218A publication Critical patent/CN111061218A/en
Application granted granted Critical
Publication of CN111061218B publication Critical patent/CN111061218B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical 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/19Numerical 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
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/35Nc in input of data, input till input file format
    • G05B2219/35349Display 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

Complex curved surface machining cutting force prediction method and device based on ACIS
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:
Figure BDA0002350589080000021
Figure BDA0002350589080000031
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,
Figure BDA0002350589080000032
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 height
Figure BDA0002350589080000035
The 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 of
Figure BDA0002350589080000036
If it is judged that
Figure BDA0002350589080000037
At [ 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;
Figure BDA0002350589080000033
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;
Figure BDA0002350589080000034
is the cutting thickness h of the jth cutting edgejWith respect to z-coordinate ziAnd radial immersion angle
Figure BDA0002350589080000041
A function of (a);
Figure BDA0002350589080000042
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle of
Figure BDA0002350589080000043
Instantaneous tangential force in time;
Figure BDA0002350589080000044
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle of
Figure BDA0002350589080000045
Instantaneous radial force in time;
Figure BDA0002350589080000046
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle of
Figure BDA0002350589080000047
Instantaneous 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;
Figure BDA0002350589080000048
wherein, TxyzThe matrix is a conversion matrix of a cutting edge coordinate system and a cutter coordinate system during forward milling or reverse milling;
Figure BDA0002350589080000049
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle of
Figure BDA00023505890800000410
Instantaneous tangential force in time;
Figure BDA00023505890800000411
is the j-th cutting edge under the coordinate system of the toolAt infinitesimal position ziHas a radial immersion angle of
Figure BDA00023505890800000412
Instantaneous radial force in time;
Figure BDA00023505890800000413
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle of
Figure BDA00023505890800000414
Instantaneous 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:
Figure BDA0002350589080000051
wherein, KiFor the cutting edge at the infinitesimal position ziThe axial immersion angle is calculated as follows:
Figure BDA0002350589080000052
Figure BDA0002350589080000053
wherein,
Figure BDA0002350589080000054
is the position z in the engagement areaiRadius of the arc-shaped infinitesimal at (a).
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;
3D
Figure BDA0002350589080000061
modeler (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:
Figure BDA0002350589080000081
Figure BDA0002350589080000082
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,
Figure BDA0002350589080000083
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 height
Figure BDA0002350589080000084
The 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,
Figure BDA0002350589080000091
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;
Figure BDA0002350589080000092
is the cutting thickness h of the jth cutting edgejWith respect to z-coordinate ziAnd radial immersion angle
Figure BDA0002350589080000093
A function of (a);
Figure BDA0002350589080000094
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle of
Figure BDA0002350589080000095
Instantaneous tangential force in time;
Figure BDA0002350589080000096
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle of
Figure BDA0002350589080000097
Instantaneous radial force in time;
Figure BDA0002350589080000098
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiRadial immersion angle ofIs composed of
Figure BDA0002350589080000099
Instantaneous 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;
Figure BDA00023505890800000910
wherein, TxyzThe matrix is a conversion matrix of a cutting edge coordinate system and a cutter coordinate system during forward milling or reverse milling;
Figure BDA00023505890800000911
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle of
Figure BDA00023505890800000912
Instantaneous tangential force in time;
Figure BDA00023505890800000913
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle of
Figure BDA00023505890800000914
Instantaneous radial force in time;
Figure BDA0002350589080000101
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle of
Figure BDA0002350589080000102
Instantaneous 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:
Figure BDA0002350589080000103
wherein, KiFor the cutting edge at the infinitesimal position ziThe axial immersion angle is calculated as follows:
Figure BDA0002350589080000104
Figure BDA0002350589080000105
wherein,
Figure BDA0002350589080000106
is the position z in the engagement areaiRadius of the arc-shaped infinitesimal at (a).
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.
3. An ACIS-based cutting force prediction method for complex curved surface machining according to claim 2, wherein the coordinate system is transformed into the expression:
Figure FDA0002350589070000021
Figure FDA0002350589070000022
wherein T is a transformation matrix, l is a cutter rake angle, and T is a cutter roll 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,
Figure FDA0002350589070000023
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 height
Figure FDA0002350589070000024
The 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 of
Figure FDA0002350589070000025
If it is judged that
Figure FDA0002350589070000026
At [ 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;
Figure FDA0002350589070000031
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;
Figure FDA0002350589070000032
is the cutting thickness h of the jth cutting edgejWith respect to z-coordinate ziAnd radial immersion angle
Figure FDA0002350589070000033
A function of (a);
Figure FDA0002350589070000034
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle of
Figure FDA0002350589070000035
Instantaneous tangential force in time;
Figure FDA0002350589070000036
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle of
Figure FDA0002350589070000037
Instantaneous radial force in time;
Figure FDA0002350589070000038
is the position z of the j-th cutting edge at infinitesimal position under the cutting edge coordinate systemiHas a radial immersion angle of
Figure FDA0002350589070000039
Instantaneous 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;
Figure FDA00023505890700000310
wherein, TxyzThe matrix is a conversion matrix of a cutting edge coordinate system and a cutter coordinate system during forward milling or reverse milling;
Figure FDA00023505890700000311
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle of
Figure FDA00023505890700000312
Instantaneous tangential force in time;
Figure FDA00023505890700000313
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle of
Figure FDA00023505890700000314
Instantaneous radial force in time;
Figure FDA00023505890700000315
is the position z of the j-th cutting edge at infinitesimal position under the tool coordinate systemiHas a radial immersion angle of
Figure FDA0002350589070000041
Instantaneous 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:
Figure FDA0002350589070000042
wherein, KiFor the cutting edge at the infinitesimal position ziThe axial immersion angle is calculated as follows:
Figure FDA0002350589070000043
Figure FDA0002350589070000044
wherein,
Figure FDA0002350589070000045
is the position z in the engagement areaiRadius of the arc-shaped infinitesimal at (a).
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.
CN201911413493.0A 2019-12-31 2019-12-31 Complex curved surface machining cutting force prediction method and device based on ACIS Active CN111061218B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (16)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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