CN112372372B - Efficient milling cutter accumulated friction wear boundary identification and verification method - Google Patents
Efficient milling cutter accumulated friction wear boundary identification and verification method Download PDFInfo
- Publication number
- CN112372372B CN112372372B CN202011439558.1A CN202011439558A CN112372372B CN 112372372 B CN112372372 B CN 112372372B CN 202011439558 A CN202011439558 A CN 202011439558A CN 112372372 B CN112372372 B CN 112372372B
- Authority
- CN
- China
- Prior art keywords
- cutter
- milling
- face
- tooth
- boundary
- 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
Links
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q17/00—Arrangements for observing, indicating or measuring on machine tools
- B23Q17/09—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool
- B23Q17/0904—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool before or after machining
- B23Q17/0914—Arrangements for measuring or adjusting cutting-tool geometry machine tools
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Abstract
The invention discloses a method for identifying and verifying an accumulated friction and wear boundary of a high-efficiency milling cutter, belongs to the technical field of milling cutters, and aims to research the defects that the maximum wear width of a rear cutter face of a cutter tooth changes along with a milling stroke and the formation process of the friction and wear boundary of the rear cutter face of the cutter tooth under the action of milling vibration and cutter tooth error cannot be disclosed by the conventional milling experiment method. The invention provides a method for measuring and characterizing the frictional wear boundary of the rear cutter face of a cutter tooth of an efficient milling cutter, a method for simulating and identifying the instantaneous frictional wear boundary of the rear cutter face of the cutter tooth, a method for calculating the cumulative frictional wear boundary of the rear cutter face of the cutter tooth and an effective model for revealing the forming process of the cumulative frictional wear boundary of the rear cutter face of the cutter tooth.
Description
The technical field is as follows:
the invention belongs to the technical field of milling cutters, and particularly relates to a high-efficiency milling cutter accumulated friction wear boundary identification and verification method.
The background art comprises the following steps:
the efficient milling cutter is widely applied to the processing of large titanium alloy structural members in the aerospace field due to the efficient cutting performance, is influenced by milling vibration and cutter tooth errors in the milling process, and uncertainty exists in the formation and evolution of accumulated friction and wear boundaries of the rear cutter face of the cutter tooth of the milling cutter, so that the whole service life of the efficient milling cutter is difficult to accurately evaluate and predict, the processing quality of workpieces is unstable, and the production cost is increased.
The uncertainty of the forming process of the frictional wear boundary of the rear cutter face of the cutter tooth of the milling cutter under the influence of milling vibration and cutter tooth errors is the key for predicting the service life of the cutter. At present, although research on frictional wear of a rear cutter face of a milling cutter tooth is available, the research mainly focuses on revealing the change characteristic of the frictional wear width of the rear cutter face of the milling cutter tooth along with the milling stroke through a milling experiment, but the research on the frictional wear width of the rear cutter face of the milling cutter tooth cannot reveal the contact relation between the rear cutter face of the milling cutter tooth and a workpiece along with the change of the milling, and due to experiment limitations, an instantaneous frictional wear boundary of the rear cutter face of the milling cutter tooth cannot be obtained through the milling experiment, so that the relation between the instantaneous and accumulated frictional wear boundary curves of the rear cutter face of the milling cutter tooth and the forming process of the accumulated frictional wear boundary curves cannot be revealed through the milling experiment.
The invention content is as follows:
the invention provides a method for identifying and verifying the accumulated frictional wear boundary of a high-efficiency milling cutter, aiming at overcoming the defects that the forming and evolution process of the accumulated frictional wear boundary of the rear cutter face of the cutter tooth of the traditional milling cutter is difficult to accurately evaluate and predict the service life of the whole high-efficiency milling cutter.
The technical scheme adopted by the invention is as follows: a method for identifying and verifying accumulated friction and wear boundaries of a high-efficiency milling cutter specifically comprises the following steps:
step A: providing a method for measuring and characterizing the frictional wear boundary of the rear cutter face of the cutter tooth of the efficient milling cutter, constructing a measuring coordinate system of the frictional wear boundary of the rear cutter face of the cutter tooth of the milling cutter, and characterizing the curve shape and the position of the frictional wear boundary of the rear cutter face of the cutter tooth of the milling cutter in the measuring coordinate system;
and B: considering a milling mode, milling vibration, cutter tooth errors and differences of surfaces to be machined of different layers as simulation boundary conditions, and providing an axial layered milling simulation model of the efficient milling cutter and a boundary condition construction method;
and C: carrying out axial layered milling thermal coupling field simulation on the efficient milling cutter by using an axial layered milling simulation model of the efficient milling cutter, extracting an instantaneous temperature, an equivalent stress and an equivalent strain distribution field of a rear cutter face of a cutter tooth of the milling cutter, and providing a method for identifying an instantaneous frictional wear simulation boundary of the rear cutter face of the cutter tooth according to the process of forming the instantaneous frictional wear boundary of the rear cutter face of the cutter tooth of the milling cutter;
step D: and C, acquiring instantaneous frictional wear boundaries of the rear cutter face of the cutter tooth at different moments, analyzing the formation process of the accumulated frictional wear boundaries of the rear cutter face of the cutter tooth, providing a simulation boundary calculation method for the accumulated frictional wear boundaries of the rear cutter face of the cutter tooth of the milling cutter, calculating the gray correlation degree and average relative error of the curve of the accumulated frictional wear boundaries of the rear cutter face of the cutter tooth of the milling cutter obtained through simulation and experiments, and verifying the correctness of the calculation method for the accumulated frictional wear boundaries of the rear cutter face of the cutter tooth of the milling cutter.
Preferably, the step a further comprises the following steps:
step A1: the milling mode of the efficient milling cutter under the vibration action is as follows: in a workpiece coordinate system, reflecting the influence of milling vibration on the cutting process of the efficient milling cutter by using the cutting state of the milling cutter under the biasing and vibration effects of the milling cutter coordinate system;
step A2: determining a milling cutter tooth error and a method for measuring the frictional wear of the rear cutter face of the milling cutter according to the structural characteristics of the efficient milling cutter;
step A3: measuring a relation matrix between a coordinate system by utilizing a workpiece coordinate system, a milling cutter coordinate system, a cutter tooth coordinate system and a rear cutter face friction and wear boundary under the influence of milling vibration and cutter tooth error, and reflecting the influence of the milling vibration and the cutter tooth error on the instantaneous contact relation between the rear cutter face of the cutter tooth and the workpiece;
step A4: and sequentially measuring the frictional wear boundary end points and inflection points of the rear cutter face of the cutter tooth in a cutter tooth rear cutter face wear boundary measurement coordinate system, constructing an upper boundary curve and a lower boundary curve of the frictional wear of the rear cutter face of the cutter tooth of the milling cutter, and determining the frictional wear boundary of the rear cutter face of the cutter tooth.
Preferably, the step B further comprises the following steps:
step B1: considering a milling mode, a cutter tooth error, milling vibration and differences of surfaces to be machined in different layers as simulation boundary conditions, designing an integral structure of an axial layered milling simulation model of the efficient milling cutter, and determining and resolving design variables of the simulation model;
step B2: constructing a simulation model milling unit, solving a milling cutter motion track and an instantaneous cutting attitude angle under the action of vibration, determining a milling unit cross section model, and constructing the milling unit by considering milling vibration and cutter tooth errors;
step B3: in the simulation process, according to the milling cutter motion trail under the vibration action in the axial layered milling experiment, determining the discrete time interval of the milling cutter motion trail points, and resolving the milling cutter motion trail points under the vibration action in the simulation process.
Preferably, the step C further comprises the following steps:
step C1: carrying out thermal coupling field simulation through an efficient milling cutter axial layered milling simulation model, and extracting the distribution conditions of an equivalent stress field, an equivalent strain field and a temperature field of the free instantaneous cutting of the rear cutter face of the cutter tooth of the milling cutter;
step C2: extracting node equivalent stress, equivalent strain and temperature values at different positions of a rear cutter face of the cutter tooth of the milling cutter, constructing a cutter face equivalent stress, equivalent strain and temperature distribution curved surface of the cutter tooth, sequentially making section planes which are perpendicular to the distribution curved surface through nodes on cutting edges in the direction of a midpoint normal vector of the cutting edges on a projection plane of the rear cutter face of the cutter tooth in parallel, and constructing section plane equivalent stress, equivalent strain and temperature distribution curves of the nodes on the section planes at different positions of the rear cutter face of the cutter tooth;
step C3: sequentially fitting equivalent stress, equivalent strain and temperature distribution curves at each section plane of the rear cutter face of the cutter tooth by using a unitary high-order polynomial, and obtaining a node change rate distribution curve by derivation;
step C4: identifying critical characteristic nodes with equivalent stress larger than and closest to yield on all section planes by taking the yield strength of the cutter tooth material as a criterion, wherein characteristic points on the section planes at different positions jointly form an instantaneous frictional wear simulation upper boundary of the cutter tooth rear cutter face in a cutter tooth rear cutter face frictional wear measurement coordinate system;
step C5: and sequentially identifying nodes with obvious abrupt changes of temperature, equivalent stress and equivalent strain change rate at each section plane of the rear cutter face of the cutter tooth, constructing a characteristic curve of the temperature, the equivalent stress and the equivalent strain abrupt changes of the nodes on the rear cutter face of the cutter tooth by taking the abrupt changes as characteristic points in a measurement coordinate system of the frictional wear boundary of the rear cutter face of the cutter tooth, and forming an outermost profile jointly formed by the characteristic curves of the temperature, the equivalent strain and the equivalent strain abrupt changes of the nodes to form an instantaneous frictional wear simulation lower boundary of the rear cutter face of the cutter tooth.
Preferably, the step D further comprises the steps of:
step D1: by means of the instant frictional wear boundary curve of the rear cutter face of the cutter tooth at different milling moments, taking the accumulated frictional wear lower boundary curve as an example, a method for calculating the accumulated frictional wear boundary of the rear cutter face of the cutter tooth of the milling cutter is provided;
step D2: designing a simulation and experiment comparison scheme of upper and lower boundaries of accumulated frictional wear of the rear cutter face of the cutter tooth, constructing a simulation upper and lower boundary curve of the accumulated frictional wear of the rear cutter face of the cutter tooth of the milling cutter by adopting a calculation method of the boundary curve of the accumulated frictional wear of the rear cutter face of the cutter tooth of the milling cutter, and comparing the experiment upper and lower boundary curve of the accumulated frictional wear with the simulation upper and lower boundary curve;
step D3: and calculating a grey correlation value and a mean relative error between simulation and experiment accumulated friction and wear boundary curves, analyzing shape similarity and position contact ratio between the simulation and experiment accumulated friction and wear boundary curves, and verifying the correctness of the accumulated friction and wear boundary calculation method.
The invention has the beneficial effects that:
1. the invention fully considers the characteristic information of the whole abrasion boundary of the rear cutter face of the cutter tooth of the milling cutter, constructs a characteristic point curve of the upper and lower boundaries of the frictional abrasion of the rear cutter face of the milling cutter by establishing a frictional abrasion boundary measuring coordinate system of the rear cutter face of the cutter tooth and taking the end points and inflection points on the left and right sides of the frictional abrasion boundary as characteristic points, can represent the shape and position information of the frictional abrasion boundary curve of the rear cutter face of the cutter tooth more completely, and is beneficial to disclosing the contact relation of a cutter worker and the formation process of the frictional abrasion boundary in the milling process by a relation matrix among a workpiece coordinate system, a milling cutter coordinate system, a cutter tooth coordinate system and a rear cutter face frictional abrasion boundary measuring coordinate system under the influence of milling vibration and cutter tooth error.
2. The invention fully considers the influence of the milling mode, the milling vibration, the cutter tooth error and the difference of different layers of surfaces to be processed on the structure of the simulation model, considers the milling mode, the milling vibration, the cutter tooth error and the different layers of surfaces to be processed as simulation boundary conditions, and provides a novel design method of the axial layered milling simulation model of the efficient milling cutter.
3. The invention provides a method for identifying upper and lower boundaries of the instantaneous frictional wear simulation of the rear cutter face of the cutter tooth of a milling cutter based on the forming process of the instantaneous frictional wear simulation boundary of the rear cutter face of the cutter tooth of the milling cutter and according to the distribution characteristics of the temperature field, the equivalent stress and the equivalent strain of the instantaneous cutting node of the rear cutter face of the cutter tooth, a method for calculating the cumulative frictional wear simulation boundary of the rear cutter face of the cutter tooth of the milling cutter, a forming process of the cumulative frictional wear boundary of the rear cutter face of the cutter tooth of the milling cutter and a method for verifying the method for calculating the cumulative frictional wear boundary of the rear cutter face of the cutter tooth of the milling cutter.
Description of the drawings:
FIG. 1 is a flow chart of a method for calculating the accumulated frictional wear boundary of the rear cutter face of a high-efficiency milling cutter tooth;
FIG. 2 is a state diagram of efficient milling by the milling cutter under the action of vibration;
FIG. 3 is a diagram showing the instantaneous cutting of the milling cutter under vibration;
fig. 4 is a structure of a high-efficiency milling cutter and a tooth error measuring method thereof, wherein fig. 4(a) is a schematic diagram of the overall structure of the high-efficiency milling cutter, and fig. 4(b) is a schematic diagram of the tooth structure of the milling cutter;
FIG. 5 is a schematic diagram of a method for measuring a frictional wear boundary and characterizing a wear depth of a flank of a cutter tooth, wherein FIG. 5(a) is a method for measuring a frictional wear boundary of a flank of a cutter tooth, and FIG. 5(b) is a wear depth of a flank of a cutter tooth;
FIG. 6 is a graph of the upper and lower boundaries of the frictional wear of the flank surface of a milling cutter tooth;
FIG. 7 is a schematic diagram of a vibration acceleration time domain signal of an effective milling stroke of a layer 4 in a layered milling mode;
fig. 8 is a wear state diagram of a bottom edge flank of a high-efficiency milling cutter tooth with a milling stroke of 5m, wherein fig. 8(a) is a wear state diagram of a first flank of the tooth, fig. 8(b) is a wear state diagram of a second flank of the tooth, and fig. 8(c) is a wear state diagram of a third flank of the tooth;
fig. 9 is a distribution graph of upper and lower boundary characteristic points of the frictional wear of the rear face of the milling cutter tooth, in which fig. 9(a) is a distribution graph of upper boundary characteristic points of the frictional wear of the rear face of the milling cutter tooth, and fig. 9(b) is a distribution graph of lower boundary characteristic points of the frictional wear of the rear face of the milling cutter tooth;
FIG. 10 is a diagram of an overall structure and a unit cross-section model of an axial layered milling simulation workpiece of the high-efficiency milling cutter, wherein FIG. 10(a) is a diagram of an overall structure of an axial layered milling simulation workpiece model of the high-efficiency milling cutter, and FIG. 10(b) is a diagram of an i (i is not less than 2) th layer of cross-section model of a milling unit;
FIG. 11 is an identification view of the cut-in of the ith (i.gtoreq.2) layer of the milling cutter;
FIG. 12 is a schematic diagram of the path of the milling cutter movement under vibration;
fig. 13 is a signal diagram of an instantaneous cutting attitude angle of the milling cutter of the i-th layer, in which fig. 13(a) is a deflection angle of the milling cutter in a feeding direction and fig. 13(b) is a deflection angle of the milling cutter in a cutting width direction;
FIG. 14 is a solid view of the milling unit of the i-th layer, wherein FIG. 14(a) is a schematic surface view of the milling unit of the i-th layer, FIG. 14(B) is a cross-sectional view taken along line A-A in FIG. 14(a), and FIG. 14(c) is a cross-sectional view taken along line B-B in FIG. 14 (a);
fig. 15 is a simulation diagram of the axial layered milling process of the high-efficiency milling cutter, wherein fig. 15(a) is a simulation diagram of a cutting-in stage of the milling cutter, fig. 15(b) is a simulation diagram of a cutting middle section of the milling cutter, and fig. 15(c) is a simulation diagram of a cutting-out stage of the milling cutter;
FIG. 16 is a graph of instantaneous equivalent stress, equivalent strain field and temperature field distribution for the flank of a milling cutter tooth, wherein FIG. 16(a) is a graph of equivalent stress distribution, FIG. 16(b) is a graph of equivalent strain field distribution, and FIG. 16(c) is a graph of temperature field distribution;
FIG. 17 is a sectional plane parallel to the normal vector direction of the midpoint of the cutting edge of the cutter tooth and perpendicular to the curved surface of the cutter tooth rear face for equivalent stress, equivalent strain and temperature distribution, wherein fig. 17(a1) is a distribution diagram of equivalent stress of the rear face of the tooth, fig. 17(a2) is a sectional plane perpendicular to the curved surface of fig. 17(a1) sequentially made by crossing the upper node of the cutting edge in a direction parallel to the midpoint normal vector of the cutting edge on the projection plane of the rear face of the tooth, fig. 17(b1) is a distribution diagram of equivalent strain of the rear face of the tooth, fig. 17(b2) is a sectional plane perpendicular to the curved surface of fig. 17(b1) sequentially made by crossing the upper node of the cutting edge in a direction parallel to the middle normal vector of the cutting edge on the projection plane of the rear face of the tooth, fig. 17(c1) is a temperature distribution diagram of the rear face of the tooth, and fig. 17(c2) is a sectional plane perpendicular to the curved surface of fig. 17(c1) sequentially made by crossing the upper node of the cutting edge in a direction parallel to the middle normal vector of the projection plane of the cutting edge on the rear face of the tooth;
fig. 18 is a graph of equivalent stress distribution at nodes at different positions of the flank of the tooth, in which fig. 18(a) is a graph of stress distribution at a node at a section a, fig. 18(B) is a graph of stress distribution at a node at a section B, and fig. 18(C) is a graph of stress distribution at a node at a section C;
FIG. 19 is a graph showing the equivalent strain and temperature distribution at the node of the section B of the flank of the tooth, wherein FIG. 19(a) is a graph showing the equivalent strain distribution at the section B and FIG. 19(B) is a graph showing the temperature distribution at the section B;
FIG. 20 is a graph showing the distribution of equivalent stress, equivalent strain and temperature change rate at section B of the flank of a tooth, in which FIG. 20(a) is a graph showing the distribution of equivalent stress change rate, FIG. 20(B) is a graph showing the distribution of equivalent strain change rate, and FIG. 20(c) is a graph showing the distribution of temperature change rate;
FIG. 21 is a graph of an upper boundary of an instantaneous frictional wear simulation of a flank face of a milling cutter tooth;
FIG. 22 is a graph of simulated lower boundary recognition of instantaneous frictional wear of the flank face of a milling cutter tooth;
FIG. 23 is a graph of the evolution process of the cumulative frictional wear boundary formation on the flank face of a milling cutter tooth;
FIG. 24 is a graph showing the upper and lower boundary of instantaneous frictional wear of a first flank of a milling cutter tooth, wherein FIG. 24(a) is a graph showing the upper boundary of instantaneous frictional wear of the first flank of the milling cutter tooth, and FIG. 24(b) is a graph showing the lower boundary of instantaneous frictional wear of the second flank of the milling cutter tooth;
FIG. 25 is a graph comparing upper and lower boundary curves of a simulated and experimental cumulative frictional wear of a milling cutter tooth I.
The specific implementation mode is as follows:
as shown in fig. 1, the method for identifying and verifying the accumulated friction and wear boundary of the milling cutter of the present invention is characterized in that a relation matrix among a workpiece coordinate system, a milling cutter coordinate system, a cutter tooth coordinate system and a rear cutter face friction and wear boundary measurement coordinate system under the influence of milling vibration and cutter tooth error is used for reflecting the influence of milling vibration and cutter tooth error on the instantaneous contact relation between the rear cutter face of the cutter tooth and the workpiece, a cutter tooth rear cutter face friction and wear boundary measurement coordinate system is established, a method for measuring and characterizing the cutter tooth rear cutter face friction and wear boundary is provided by using end points and inflection points on the left side and the right side of the friction and wear boundary as characteristic points, and the curve shape and the position of the cutter tooth rear cutter face friction and wear boundary are completely characterized in the measurement coordinate system; considering milling vibration, cutter tooth error, milling mode and difference of surfaces to be machined of different layers as simulation boundary conditions, and providing an axial layered milling simulation model of the efficient milling cutter and a boundary condition construction method; providing an upper boundary simulation criterion and a lower boundary simulation criterion of the instantaneous frictional wear of the rear cutter face of the milling cutter tooth based on the formation process of the instantaneous frictional wear boundary of the rear cutter face of the cutter tooth and the distribution characteristics of the instantaneous node temperature field, the equivalent stress and the equivalent strain of the rear cutter face of the cutter tooth; the method for calculating the accumulated friction and wear boundary of the rear cutter face of the milling cutter tooth is provided based on an instantaneous friction and wear simulation criterion, the similarity degree of the shape and the position between simulation and experiment accumulated friction and wear boundary curves is analyzed by calculating the gray correlation value and the average relative error between the simulation and experiment accumulated friction and wear boundary curves, and the correctness of the method for calculating the accumulated friction and wear boundary of the rear cutter face of the milling cutter tooth is verified.
The method for identifying and verifying the accumulated friction and wear boundary of the efficient milling cutter specifically comprises the following steps:
step A: method for measuring and characterizing frictional wear boundary of rear cutter face of cutter tooth of efficient milling cutter
Milling mode of efficient milling cutter under vibration action of step A1
In a workpiece coordinate system, the influence of milling vibration on the cutting process of the high-efficiency milling cutter is reflected by using the offset and the cutting state of the milling cutter under the action of the milling cutter coordinate system, as shown in figures 2 and 3, and variable parameters in figures 2 and 3 are explained as shown in a table 1-1.
TABLE 1-1 efficient milling cutter milling mode and instantaneous cutting state variable interpretation under vibration
Step A2 cutter tooth error and flank friction wear measurement
In order to research the characteristic of influence of the cutter tooth error of the milling cutter on the frictional wear boundary of the rear cutter face of the cutter tooth, the radial error of the cutter tooth is measured by taking the maximum excircle radius of the cutter tooth of the milling cutter as a reference in a milling cutter coordinate system, and the axial error of the cutter tooth is measured by taking the lowest point of the cutter tooth of the milling cutter as a reference. As shown in fig. 4, the explanation of the variable parameters in fig. 4 is shown in tables 1-2.
TABLE 1-2 advanced milling cutter structure and variable interpretation of tooth error measurement method thereof
Making the mounting reference surface of the cutter tooth perpendicular to the measuring horizontal plane, and on the projection of the rear cutter surface of the cutter tooth on the measuring horizontal plane, making the normal vector passing through the cutter point outside the cutter tooth and parallel to the middle point of the cutting edge and the reverse intersection point O of the tangent vector along the middle point of the cutting edge j Establishing a frictional wear boundary measurement coordinate system O of the rear cutter face of the cutter tooth of the milling cutter as an original point j UVS, axis U parallel to the vector tangential to the midpoint of the cutting edge of the tooth on the projection plane, axis V parallel to the vector normal to the midpoint of the cutting edge on the projection plane, and axis S perpendicular to the measurement plane. And the flank face frictional wear upper and lower boundaries are measured on a parallel plane of the measurement coordinate system parallel to the plane of the midpoint tangent of the tooth cutting edge, as shown in fig. 5 (a). The method for characterizing the frictional wear depth of the rear face of the cutter tooth is shown in fig. 5 (b).
The explanation of the variable parameters in fig. 5 is shown in tables 1-3.
TABLE 1-3 knife teeth flank friction wear boundary measurement method variable interpretation
Step A3 influence of milling vibration and cutter tooth error on instantaneous contact relation between cutter tooth flank and workpiece
And reflecting the influence of the milling vibration and the cutter tooth error on the instantaneous contact relation between the rear cutter face of the cutter tooth and the workpiece by utilizing a relation matrix among a workpiece coordinate system, a milling cutter coordinate system, the cutter tooth coordinate system and the rear cutter face frictional wear measurement coordinate system under the influence of the milling vibration and the cutter tooth error. The relation matrix between the milling cutter coordinate system and the workpiece coordinate system under the vibration-free action is shown as a formula (1-1), the relation matrix between the milling cutter coordinate system under the vibration-free action is shown as a formula (1-3), the relation matrix between the cutter tooth coordinate system and the milling cutter coordinate system is shown as a formula (1-5), and the relation matrix between the cutter tooth coordinate system and the rear cutter face friction wear measurement coordinate system is shown as a formula (1-7).
[x g y g z g 1] T ＝M 1 [a b c 1] T (1-1)
[a b c 1] T ＝M 2 T 1 [a v b v c v 1] T (1-3)
[a b c 1] T ＝M 3 T 2 [a q b q c q 1] T (1-5)
[a q b q c q 1] T ＝M 4 T 4 T 3 [U V S 1] T (1-7)
Wherein M is 1 、M 2 、M 3 、M 4 For translation matrices, T 1 、T 2 、T 3 、T 4 For a rotation matrix, M 1 Wherein x (T), y (T), z (T) are the positions of the milling cutter center in the workpiece coordinates at time T, T 4 In the step (1), the first step,
is an included angle between the cutter tooth positioning surface and the milling cutter groove mounting surface,
is pi/2.
Step A4 cutter tooth flank frictional wear boundary determination
In order to completely represent the shape and position information of the upper and lower boundaries of the frictional wear of the rear cutter face of the milling cutter tooth, a measuring coordinate system O of the wear boundary of the rear cutter face of the milling cutter tooth measured under a microscope with ultra-depth of field j And (3) sequentially identifying left and right end points and inflection points of the frictional wear upper and lower boundaries of the rear cutter face of the cutter tooth along the cutting edge direction on a projection plane on the UV plane, sequentially measuring coordinate values of the end points and the inflection points along the V-axis direction in a cutter tooth rear cutter face wear boundary measurement coordinate system, and constructing a milling cutter tooth rear cutter face frictional wear upper and lower boundary curve by taking the left and right end points and the inflection points as characteristic points. As shown in fig. 6, wherein V ═ V 0 (U) is the profile curve V of the original cutting edge of the rear face of the milling cutter tooth rq (U, L) (q is more than or equal to 1 and less than or equal to 3) is a characteristic point curve of the upper boundary of the frictional wear of the rear face of the milling cutter tooth q with the milling stroke of L, and V is V hq And (U, L) is a lower boundary characteristic point curve of the frictional wear of the rear cutter face of the milling cutter tooth q with the milling stroke L.
The method comprises the following implementation steps of carrying out an experiment of milling titanium alloy by axial layering of the high-efficiency milling cutter to obtain the friction and wear result of the rear cutter face of the cutter tooth of the high-efficiency milling cutter, wherein the milling scheme is shown in tables 1-4.
TABLE 1-4 Experimental scheme for milling titanium alloy by using efficient milling cutter
In tables 1-4, L is the total effective milling stroke of 5m, Δ c z i Is the axial error, Δ c, of tooth i j i Is the radial error of the cutter tooth i. And calculating the axial and radial errors of each cutter tooth of the milling cutter by taking the reference of the lowest point of the cutter tooth of the milling cutter and the maximum radius of the cutter point of the cutter tooth of the milling cutter as the reference.
After the cutter tooth error measurement is finished, a titanium alloy milling experiment is carried out on a triaxial milling machining center, wherein the length and width of a workpiece are 250 multiplied by 100 multiplied by 10mm, the axial layered milling experiment is finished in 20 layers, in the milling experiment process, an acceleration sensor is adopted to detect a vibration signal generated by the excitation of cutting force on the surface of the workpiece close to a cutting area, and a DH5922 transient signal test analysis system is adopted to carry out data analysis on the collected vibration signal. As shown in fig. 7.
After the experiment is finished, the ultra-depth-of-field microscope is used for sequentially detecting the frictional wear boundaries of the rear tool faces of the cutter teeth of the milling cutter, as shown in fig. 8.
The characteristic points of the upper and lower boundary curves of the wear of the rear cutter face of the cutter tooth of the milling cutter in fig. 8 are sequentially extracted by adopting the method for identifying the characteristic points of the upper and lower boundary curves of the friction wear, as shown in fig. 9.
As can be seen from FIG. 9, the measurement method can completely characterize the shape and the position of the upper and lower boundary of the wear of the rear cutter face of the milling cutter tooth; under the influence of cutter tooth errors and milling vibration, the distribution of the upper and lower boundaries of the wear of the rear cutter faces of different cutter teeth of the milling cutter presents obvious difference.
And B: efficient milling cutter axial layered milling simulation model and boundary condition construction method
Step B1 high-efficiency milling cutter axial layered milling simulation model overall structure design
Considering the milling mode, the cutter tooth error, the milling vibration and the difference of the surfaces to be machined of different layers as simulation boundary conditions, designing the overall structure of the axial layered milling simulation model of the efficient milling cutter, and determining and calculating the design variables of the simulation model, as shown in fig. 10.
In fig. 10(a), the simulated workpiece is composed of a plurality of milling units, each layer of milling unit corresponds to a workpiece form in an axial layered milling process, rectangular groove structures are formed among the milling units, the profiles of the surfaces to be processed of different layers have differences but are positioned on the same horizontal height, and L is a height i (i+1) And H' is the height of the rectangular groove, and is the distance between the (i-1) th layer milling unit and the ith layer milling unit.
In FIG. 10(b), O i The center of a circle where a bottom edge cutting edge of a milling cutter tooth is positioned, r 0 Radius of the circle on which the cutting edge of the milling cutter tooth is located, d m As axis C of milling cutter 1 C 2 Distance from work piece in cutting width direction, J d For the lowest point of participation in cutting on the cutting edge of the milling cutter tooth, J h Is the highest contact point, delta r, of the side edge of the cutter tooth and the surface to be machined z max And Δ r j max The maximum axial error and the maximum radial error of the cutter teeth of the milling cutter are respectively; h i The height of the milling unit of the ith (i ≧ 2) layer is represented, wherein the cross section of the milling unit of the first layer is the same size as the initial cross section of the workpiece.
The control parameters of the simulation workpiece are as follows: l is a radical of an alcohol i (i+1) 、H'、L i-1 i 、H i 、W、L 0 、d m 、a p 、Δr z max 、Δr j max 、r 0 、r、α、J d 、J h 。L i (i+1) The solving method of (2) is as follows:
wherein，T i i+1 For milling cutters in actual machining processes from P 3i To P 2(i+1) The time of (c).
H i ＝H+(i-1)a p (i≥2) (2-3)
Step B2 simulation model milling unit construction
(1) According to the steps A shown in the figures 2 and 3, the instantaneous deflection angles delta of the milling cutter along the feed speed direction and the cutting width direction at the time t in the milling process are solved 1 (t)、δ 2 (t) is represented by the following formulae (2-4) and (2-5).
Solving the track equation f of the central point of the milling cutter in the workpiece coordinate system under the action of vibration v ＝(x gv ,y gv ,z gv And t) is shown as formula (2-6).
In the formula (2-6), A x (t)、A y (t)、A z (t) the vibration acceleration signal is obtained by performing secondary integration on the vibration acceleration signal measured in the experimental process and then using 8 times of Sun of sin Function fitting in Matlab. t is t e For milling cutters from P 4i To P 1(i+1) The time used is as shown in formula (2-7), T i For milling cutters from P 1(i+1) Move to P 1(i+1) The time period of (2) is as shown in the formula (2-8).
T i ＝[T si ,T ei ]＝[(i-1)(T c +t e ),iT c +(i-1)t e ] (2-8)
Wherein, T si For milling cutter along x g Moment of starting point of tool path for ith axial feed, T ei Is the milling cutter edge x g The tool lifting time of the ith feeding in the axial direction; t is a unit of c Milling cutter from P 1i Move to P 4i The time used is as in formula (2-9).
(2) Identifying the cutting-in and cutting-out time in the ith (i is more than or equal to 2) layer milling vibration signal, and solving the ith layer cutting section T according to the formula (2-6) as shown in FIG. 11 m Equation f of motion track of milling cutter v ＝(x gv ,y gv ,z gv T), as shown in fig. 12.
(3) The i-layer milling cutter instantaneous cutting attitude angle was solved using matlab according to equations (2-4) and (2-5), and the result is shown in fig. 13.
(4) As shown in fig. 14, according to the milling unit cross section model in fig. 10(b), the i-th layer milling unit cross section is constructed, and NJ is added h The section is swept along a milling cutter track curve to form a surface to be machined of the milling unit, an i-th layer of milling cutter instantaneous cutting attitude angle is introduced in the sweeping process, and finally other surfaces of the milling unit are constructed in sequence to form an entity, wherein parameter explanation of variables in the graph 14 is shown in a table 2-1.
TABLE 2-1 milling Unit variable resolution for ith layer
B3 axial layered milling simulation motion trail calculation of efficient milling cutter under vibration action
In the simulation process, in order to introduce the milling cutter vibration displacement in the milling process, according to the milling cutter motion trail under the vibration action in the axial layered milling experiment, on the premise of not changing the feeding speed, the discrete time interval of the milling cutter motion trail points is determined, and the milling cutter motion trail points under the vibration action in the simulation process are solved, wherein the method comprises the following steps of:
(1) determining discrete time interval of trace points, and determining minimum unit size L of minimum grid unit of workpiece in order to make simulation result more accurate min 1/3 is the simulation step size, t d As shown in formulas (2-10).
(2) And (4) re-calculating the coordinate value of the cutting-in to cutting-out of any ith (i is more than or equal to 2) layer of the milling cutter, as shown in a formula (2-11).
(3) And (3) the coordinates of milling cutter track points between any ith layer of cutting points and the (i + 1) th layer of cutting points are supplemented, the time step number from the ith layer of cutting points to the (i + 1) th layer of cutting points is calculated as formula (2-12), the step distance between adjacent track points from the ith layer of cutting points to the (i + 1) th layer of cutting points is calculated as formula (2-13), and the coordinate calculation method between the adjacent two layers of track points is calculated as formula (2-14).
And C: method for identifying instantaneous frictional wear boundary of rear cutter face of cutter tooth
C1 simulation of axial layered milling thermal coupling field of efficient milling cutter
And C, in order to obtain the instantaneous frictional wear simulation boundary of the rear cutter face of the cutter tooth of the milling cutter, utilizing the milling experimental parameters in the step A, milling a simulation model in an axial layered mode according to the efficient milling cutter designed in the step B, and carrying out thermal coupling field simulation, wherein a Johnson-Cook constitutive model adopted by the simulation is shown as a formula (3-1), and the titanium alloy milling experimental processing parameters of the efficient milling cutter in the table 1-4 are adopted as the simulation milling parameters.
In the formula: σ — equivalent flow stress; a, B, n, c, m-yield stress strength, strain strengthening constant, strain strengthening index, strain rate strengthening parameter and temperature strain rate sensitivity; t, T r ，T m -real time temperature, room temperature, 25 ℃ and material melting temperature/° c; the sum of the distances of epsilon,
material equivalent strain, equivalent strain rate/s -1 The reference strain rate takes 1s -1 (ii) a Wherein the high-efficiency milling cutter base material is hard alloy WC, the cutter coating material is TiN, the thickness is 5 mu m, and the yield strength sigma is s The intrinsic parameters of the TC4 titanium alloy material J-C are shown in the table 3-1 when the thickness is 5.9 GPa. The milling simulation process is shown in fig. 15.
TABLE 3-1 TC4 constitutive parameter table of titanium alloy J-C
After the simulation is finished, an equivalent stress field, an equivalent strain field and a temperature field of any instantaneous cutting of the rear cutter face of the cutter tooth of the milling cutter are extracted, as shown in fig. 16.
C2 curve of equivalent stress, equivalent strain and temperature change of instantaneous node at different positions of cutter tooth rear cutter face
Extracting node equivalent stress, equivalent strain and temperature values at different positions of the rear cutter face of the cutter tooth of the milling cutter, constructing a cutter tooth rear cutter face equivalent stress, equivalent strain and temperature distribution curved surface, and sequentially making a plurality of nodes perpendicular to the middle point normal vector direction of the upper cutting edge on the projection plane of the rear cutter face of the cutter tooth to pass through the upper cutting edge in parallelThe section plane of the distribution curved surface is, as shown in fig. 17, constructed by the equivalent stress distribution curve of the node on the section plane at different positions of the flank face and identified that the equivalent stress on the section plane is greater than and closest to the yield strength σ s Fig. 18 shows that a curve of equivalent strain and temperature distribution of the nodes on the section plane at different positions of the flank face is constructed, and fig. 19 shows an example of the section plane B.
Step C3 curve of stress, strain and temperature change rate at different positions of cutter tooth rear cutter face
In the cutting process of the high-efficiency milling cutter, shear and friction exist between the rear cutter face of the cutter tooth and the transition surface of a workpiece, so that stress concentration is formed inside the cutter tooth in a contact area and is diffused and distributed to a non-contact area, and the stress value is sharply reduced at a contact and non-contact boundary node; meanwhile, the temperature of the contact surface of the rear cutter face of the cutter tooth is sharply increased by heat generated in the shearing and rubbing process, the generated heat is transferred to a non-contact area at the rear cutter face of the cutter tooth, and the temperature value is sharply reduced at a contact and non-contact boundary node, so that a very large temperature gradient is generated on the rear cutter face of the cutter tooth; in the milling process, the interaction between the rear cutter face of the cutter tooth and the transition surface of the workpiece generates extrusion deformation, so that the contact boundary node of the rear cutter face of the cutter tooth and the transition surface of the workpiece generates strain concentration, and the non-contact boundary generates micro strain; based on the analysis, the numerical values of stress, strain and temperature at the nodes of the instantaneous frictional wear boundary of the cutter tooth rear cutter face are mutated, so that the node temperature, equivalent stress and equivalent strain distribution curves at all the sections are sequentially fitted through a unitary high-order polynomial, the node change rate distribution curve is obtained by derivation through the formulas (3-2), (3-3) and (3-4), and the node inflection points at which the temperature, stress and strain change rate are mutated at all the sections of the rear cutter face are respectively identified by taking the section B as an example, as shown in fig. 20.
Step C4 cutter tooth flank instantaneous friction wear simulation upper boundary identification
Based on the mechanical property of the material, when the equivalent stress on the cutting edge of the cutter tooth is greater than the yield strength, the material on the cutting edge is peeled off and taken away by the cuttings, so that the shape of the edge is changed to form an instantaneous frictional wear upper boundary of the rear cutter face of the cutter tooth, and based on the analysis, the yield strength sigma of the cutter tooth material is used s The simulation upper boundary identification judgment for the instantaneous frictional wear of the rear cutter face of the milling cutter tooth is shown as a formula (3-5), a critical characteristic node with equivalent stress larger than and closest to yield strength on all the section planes is identified, and in a rear cutter face frictional wear measurement coordinate system, a boundary surrounded by the characteristic nodes on the section planes is the upper boundary of the instantaneous frictional wear of the rear cutter face of the milling cutter tooth, as shown in fig. 21.
σ≥σ s (3-5)
Step C5 cutter tooth flank instantaneous frictional wear simulation lower boundary recognition
In the instantaneous cutting process, the rear cutter face of the milling cutter tooth and the transition surface of the workpiece are sheared, extruded and rubbed to form an instantaneous frictional wear lower boundary of the rear cutter face of the milling cutter tooth, the temperature, equivalent stress and equivalent strain between contact and non-contact areas of the rear cutter face of the milling cutter tooth and the transition surface of the workpiece are obviously different, based on the analysis, a simulation lower boundary identification method for the instantaneous frictional wear of the rear cutter face of the milling cutter tooth is provided, the temperature, the equivalent stress and the nodes where the equivalent strain is obviously mutated at each section plane of the rear cutter face of the milling cutter tooth are sequentially identified by adopting the method, the temperature, the equivalent stress and the equivalent strain mutation curves of the nodes on the rear cutter face of the milling cutter tooth are shown in a graph 22, and the calculation method for the instantaneous frictional wear simulation lower boundary curves of the rear cutter face of the milling cutter tooth is shown in a formula (3-6) and a formula (3-7).
V m (U i ,t)＝max{V T (U i ,t),V σ (U i ,t),V ε (U i ,t)}(0≤i≤10.7) (3-6)
V m (U,t)＝{V m (U 1 ,t),V m (U 2 ,t)...V m (U i ,t)} (3-7)
As shown in fig. 22, V ═ V 0 (U) is the boundary of the original profile of the milling cutter, V ═ V T (U,t)、V＝V σ (U,t)、V ε (U, t) are respectively the curves of the temperature, stress and strain sudden change node of the rear cutter face of the milling cutter tooth at the time t, and in the formula (3-6), V T (U i ,t)、V σ (U i ,t)、V ε (U i T) is the temperature, stress and strain abrupt change node curve U of the rear cutter face of the cutter tooth at the time t i The V-direction coordinate value at the position. V in the formula (3-7) m (U i T) is a lower boundary curve U of instantaneous frictional wear of the rear cutter face of the milling cutter tooth at the time t i Coordinate value of V direction at position, V m And (U, t) is a lower boundary curve of instantaneous frictional wear of the rear cutter face of the milling cutter tooth at the time t.
Step D: method for calculating and verifying accumulated friction and wear boundary of cutter tooth rear cutter face
Step D1 calculation of accumulated friction and wear boundary of cutter tooth rear face of milling cutter
In the milling process, under the influence of milling vibration and cutter tooth errors, instantaneous contact boundaries between rear cutter faces of cutter teeth of the milling cutter and transition surfaces of a workpiece dynamically change along with time, accumulated friction and wear boundaries of the rear cutter faces of the cutter teeth of the milling cutter are formed by accumulating a plurality of instantaneous friction and wear boundaries, based on the analysis, by taking the lower boundary simulated by the accumulated friction and wear of the rear cutter faces of the cutter teeth of the milling cutter as an example, instantaneous friction and wear lower boundary curves of the rear cutter faces of the cutter teeth at different milling moments are extracted, and as shown in fig. 23, a calculation method for the accumulated friction and wear boundaries of the rear cutter faces of the cutter teeth of the milling cutter is provided, and is shown in a formula (4-1) and a formula (4-2).
V a (U i ,Δt)＝max{V m (U i ,t 1 ),V m (U i ,t 2 )...V m (U i ,t n )}Δt＝t n -t 1 (4-1)
V a (U,Δt)＝{V a (U 1 ,Δt),V a (U 2 ,Δt)...V a (U i ,Δt)} (4-2)
In the formula (4-1), V a (U i And delta t) is the accumulated friction wear boundary curve of the rear cutter face of the milling cutter tooth in the time period of delta t at U i Coordinate value of V direction at position, in formula (4-2), V a (U, delta t) is a boundary curve of accumulated friction wear of the rear face of the milling cutter tooth in a delta t time period.
D2 cutter tooth flank accumulated friction wear lower boundary simulation and experiment comparison scheme
In order to verify the correctness of the method for calculating the accumulated frictional wear boundary of the efficient milling cutter, firstly, the simulation recognition criteria of the upper and lower boundary of the instantaneous frictional wear of the rear cutter face of the cutter tooth of the milling cutter in the steps C4 and C5 are adopted, and the simulation upper and lower boundaries of the instantaneous frictional wear of the cutter face of each cutter tooth of the milling cutter at the time of 0.5m, 1m, 1.5m, 2m, 2.5m, 3m, 3.5m, 4m, 4.5m and 5m in the simulation of the thermal-force coupling field in the step C1 are extracted, taking the cutter tooth I as an example, as shown in FIG. 24.
And D1, constructing a simulation upper and lower boundary curve of the accumulated frictional wear of the rear cutter face of the milling cutter tooth by adopting the calculation method of the boundary curve of the accumulated frictional wear of the rear cutter face of the milling cutter tooth in the step D1, and comparing the simulation upper and lower boundary curve of the accumulated frictional wear of the rear cutter face of the milling cutter tooth with the experimental frictional wear upper and lower boundary curve of the rear cutter face of the milling cutter tooth with the milling stroke of 5m in the step A4, taking the cutter tooth I as an example, as shown in FIG. 25.
Step D3 simulation and experiment cumulative friction wear boundary shape and position similarity analysis
Calculating a gray correlation value and a mean relative error between simulation and experimental accumulated friction and wear boundary curves, analyzing shape similarity and position contact degree between the simulation and experimental wear boundary curves, and verifying the correctness of the accumulated friction and wear boundary calculation method, wherein the gray correlation analysis method comprises the following steps:
firstly, respectively taking coordinate values of V axes of characteristic points of wear boundaries of rear cutter faces of experimental and simulated cutter teeth as a reference sequence and a comparison sequence, wherein the coordinate values are expressed as formulas (4-3) and (4-4).
V q ＝(V q (1),V q (2),...,V q (m),...,V q (n)) (4-3)
V q *＝(V q (1)*,V q (2)*,...,V q (m)*,...,V q (n)*) (4-4)
V in the formulae (4-3) and (4-4) q (m) and V q (m) is the coordinate value V of the mth characteristic point V along the cutting edge direction on the boundary curve of the wear of the flank of the q (q is more than or equal to 0 and less than or equal to 3) th cutter tooth after the experiment and simulation respectively m 。
The method for calculating the average relative error of the wear curves of the rear cutter face of the simulated and experimental cutter teeth is shown as the formula (4-5).
In the formula (4-5), δ q The relative error of the wear boundary curve of the q-th cutter tooth rear cutter surface is simulated and tested.
By adopting the correlation and average relative error calculation method, the correlation and relative error of the upper and lower boundary curves of the wear of the rear cutter face of the simulated and experimental cutter tooth are shown in the table 4-1, wherein gamma (V) q ,V q (1 is more than or equal to q is less than or equal to 3) represents the grey correlation degree between the simulation of the cutter tooth q and the experimental accumulated friction wear boundary curve, delta q The average relative error between the tooth q simulation and the experimental cumulative wear boundary curve is shown.
TABLE 4-1 simulation and experiment cutter tooth back cutter face abrasion upper and lower boundary curve correlation and relative error
From table 4-1, the correlation values of the upper and lower boundary curves of the simulation and the experiment accumulated frictional wear of each cutter tooth rear cutter face of the milling cutter are all above 0.75, and the average relative errors are all less than 20%, and the results show that the curves and the positions of the simulation and the experiment accumulated frictional wear boundary curves of each cutter tooth rear cutter face of the milling cutter have high similarity, the effectiveness of the design method of the high-efficiency milling cutter axial layered milling simulation model is achieved, and meanwhile, the correctness of the calculation method of the accumulated frictional wear boundary of the cutter tooth rear cutter face of the milling cutter is verified.
The above description is only a preferred embodiment of the present invention, and these embodiments are based on different implementations of the present invention, and the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (4)
1. A high-efficiency milling cutter accumulated friction and wear boundary identification and verification method is characterized by comprising the following steps:
step A: providing a method for measuring and characterizing the frictional wear boundary of the rear cutter face of the cutter tooth of the efficient milling cutter, constructing a measuring coordinate system of the frictional wear boundary of the rear cutter face of the cutter tooth of the milling cutter, and characterizing the shape and the position of the frictional wear boundary of the rear cutter face of the cutter tooth on a plane parallel to a measuring coordinate system parallel to a mid-point tangent plane of the cutting edge of the cutter tooth;
the step A further comprises the following steps:
step A1: the milling mode of the efficient milling cutter under the vibration action is as follows: in a workpiece coordinate system, reflecting the influence of milling vibration on the cutting process of the efficient milling cutter by using the cutting state of the milling cutter under the biasing and vibration effects of the milling cutter coordinate system;
step A2: determining a milling cutter tooth error and a method for measuring a frictional wear boundary of a rear cutter face of the milling cutter according to the structural characteristics of the efficient milling cutter;
step A3: measuring a relation matrix between a coordinate system by utilizing a workpiece coordinate system, a milling cutter coordinate system, a cutter tooth coordinate system and a rear cutter face friction and wear boundary under the influence of milling vibration and cutter tooth error, and reflecting the influence of the milling vibration and the cutter tooth error on the instantaneous contact relation between the rear cutter face of the cutter tooth and the workpiece;
step A4: sequentially measuring frictional wear boundary end points and inflection points of the rear cutter face of the cutter tooth in a cutter tooth rear cutter face wear boundary measurement coordinate system, constructing an upper boundary curve and a lower boundary curve of frictional wear of the rear cutter face of the cutter tooth of the milling cutter, and determining a frictional wear boundary of the rear cutter face of the cutter tooth;
and B: considering a milling mode, milling vibration, cutter tooth errors and differences of surfaces to be machined of different layers as simulation boundary conditions, and providing an axial layered milling simulation model of the efficient milling cutter and a boundary condition construction method;
and C: carrying out axial layered milling thermal coupling field simulation on the efficient milling cutter by using an axial layered milling simulation model of the efficient milling cutter, extracting an instantaneous temperature, an equivalent stress and an equivalent strain distribution field of a rear cutter face of a cutter tooth of the milling cutter, and providing a method for identifying an instantaneous frictional wear simulation boundary of the rear cutter face of the cutter tooth according to the process of forming the instantaneous frictional wear boundary of the rear cutter face of the cutter tooth of the milling cutter;
step D: and C, acquiring instantaneous frictional wear boundaries of the rear cutter face of the cutter tooth at different moments, analyzing the formation process of the accumulated frictional wear boundaries of the rear cutter face of the cutter tooth, providing a simulation boundary calculation method for the accumulated frictional wear boundaries of the rear cutter face of the cutter tooth of the milling cutter, calculating the gray correlation degree and average relative error of the curve of the accumulated frictional wear boundaries of the rear cutter face of the cutter tooth of the milling cutter obtained through simulation and experiments, and verifying the correctness of the calculation method for the accumulated frictional wear boundaries of the rear cutter face of the cutter tooth of the milling cutter.
2. A high efficiency milling cutter cumulative frictional wear boundary identification and verification method as claimed in claim 1, wherein: the step B further comprises the following steps:
step B1: considering a milling mode, a cutter tooth error, milling vibration and differences of surfaces to be machined of different layers as simulation boundary conditions, designing an integral structure of an axial layered milling simulation model of the efficient milling cutter, and determining and resolving design variables of the simulation model;
step B2: constructing a simulation model milling unit, solving a milling cutter motion track and an instantaneous cutting attitude angle under the action of vibration, determining a milling unit cross section model, and constructing the milling unit by considering milling vibration and cutter tooth errors;
step B3: in the simulation process, according to the milling cutter motion trail under the vibration action in the axial layered milling experiment, determining the discrete time interval of the milling cutter motion trail points, and resolving the milling cutter motion trail points under the vibration action in the simulation process.
3. A high efficiency milling cutter cumulative frictional wear boundary identification and verification method as claimed in claim 2, wherein: the step C further comprises the following steps:
step C1: carrying out thermal coupling field simulation through an efficient milling cutter axial layered milling simulation model, and extracting the distribution conditions of an equivalent stress field, an equivalent strain field and a temperature field of the free instantaneous cutting of the rear cutter face of the cutter tooth of the milling cutter;
step C2: extracting node equivalent stress, equivalent strain and temperature values at different positions of a rear cutter face of the cutter tooth of the milling cutter, constructing a cutter face equivalent stress, equivalent strain and temperature distribution curved surface of the cutter tooth, sequentially making section planes which are perpendicular to the distribution curved surface through nodes on cutting edges in the direction of a midpoint normal vector of the cutting edges on a projection plane of the rear cutter face of the cutter tooth in parallel, and constructing section plane equivalent stress, equivalent strain and temperature distribution curves of the nodes on the section planes at different positions of the rear cutter face of the cutter tooth;
step C3: sequentially fitting equivalent stress, equivalent strain and temperature distribution curves of each section plane of the rear cutter face of the cutter tooth by adopting a unary high-order polynomial, and obtaining a node change rate distribution curve by derivation;
step C4: identifying critical characteristic nodes with equivalent stress larger than and closest to yield on all section planes by taking the yield strength of the cutter tooth material as a criterion, wherein characteristic points on the section planes at different positions in a friction and wear measurement coordinate system of the cutter tooth rear cutter face jointly form an instantaneous friction and wear simulation upper boundary of the cutter tooth rear cutter face;
step C5: and sequentially identifying nodes with obvious abrupt changes in temperature, equivalent stress and equivalent strain change rate at each section plane of the rear cutter face of the cutter tooth, constructing a node temperature, equivalent stress and equivalent strain abrupt change characteristic curve on the rear cutter face of the cutter tooth in a frictional wear boundary measurement coordinate system by taking the abrupt change nodes as characteristic points, and forming an outermost profile jointly formed by the node temperature, the equivalent strain and the equivalent strain abrupt change characteristic curve to form an instantaneous frictional wear simulation lower boundary of the rear cutter face of the cutter tooth.
4. A high efficiency milling cutter cumulative frictional wear boundary identification and verification method as claimed in claim 1, wherein: the step D further comprises the following steps:
step D1: by means of the instant frictional wear boundary curve of the rear cutter face of the cutter tooth at different milling moments, taking the accumulated frictional wear lower boundary curve as an example, a method for calculating the accumulated frictional wear boundary of the rear cutter face of the cutter tooth of the milling cutter is provided;
step D2: designing a simulation and experiment comparison scheme of upper and lower boundaries of accumulated friction wear of the rear cutter face of the cutter tooth, constructing a simulation upper and lower boundary curve of the accumulated friction wear of the rear cutter face of the cutter tooth of the milling cutter by adopting a calculation method of the curve of the accumulated friction wear boundary of the rear cutter face of the cutter tooth of the milling cutter, and comparing the experiment upper and lower boundary curve of the accumulated friction wear with the simulation upper and lower boundary curve;
step D3: and calculating a grey correlation value and a mean relative error between simulation and experiment accumulated friction and wear boundary curves, analyzing shape similarity and position contact ratio between the simulation and experiment accumulated friction and wear boundary curves, and verifying the correctness of the accumulated friction and wear boundary calculation method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011439558.1A CN112372372B (en) | 2020-12-07 | 2020-12-07 | Efficient milling cutter accumulated friction wear boundary identification and verification method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011439558.1A CN112372372B (en) | 2020-12-07 | 2020-12-07 | Efficient milling cutter accumulated friction wear boundary identification and verification method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112372372A CN112372372A (en) | 2021-02-19 |
| CN112372372B true CN112372372B (en) | 2022-09-02 |
Family
ID=74590798
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011439558.1A Active CN112372372B (en) | 2020-12-07 | 2020-12-07 | Efficient milling cutter accumulated friction wear boundary identification and verification method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112372372B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114769692B (en) * | 2022-04-19 | 2023-01-03 | 哈尔滨理工大学 | Efficient cross-scale identification method for friction wear boundary of rear cutter face of milling cutter |
| CN114952420B (en) * | 2022-05-23 | 2023-04-11 | 大连理工大学 | Wear prediction calculation method for milling CFRP (carbon fiber reinforced plastics) of ball end mill based on path dispersion |
| CN116061004B (en) * | 2022-10-14 | 2023-08-08 | 哈尔滨理工大学 | Identification method for change characteristics of transient cutting energy efficiency distribution of milling cutter under vibration effect |
| CN117521374B (en) * | 2023-06-17 | 2024-04-09 | 哈尔滨理工大学 | Method for calculating propagation and attenuation characteristics of friction stress wave |
| CN116690313B (en) * | 2023-08-01 | 2023-12-12 | 成都飞机工业（集团）有限责任公司 | Failure monitoring method for machining cutter of web plate of aircraft structural member |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5251144A (en) * | 1991-04-18 | 1993-10-05 | Texas Instruments Incorporated | System and method utilizing a real time expert system for tool life prediction and tool wear diagnosis |
| CN109434564A (en) * | 2018-12-21 | 2019-03-08 | 哈尔滨理工大学 | A kind of cutter wear state monitoring method based on deep neural network |
| CN109940460A (en) * | 2019-04-10 | 2019-06-28 | 哈尔滨理工大学 | The detection method of milling machined surface geometric error distribution character under a kind of effect of vibration |
| CN109940461A (en) * | 2019-04-10 | 2019-06-28 | 哈尔滨理工大学 | A kind of detection method of high feed milling cutter cutter tooth wear of the tool flank characteristic |
| CN110103079A (en) * | 2019-06-17 | 2019-08-09 | 中国科学院合肥物质科学研究院 | The on-line monitoring method of tool wear in a kind of micro- milling process |
| CN110704969A (en) * | 2019-09-24 | 2020-01-17 | 哈尔滨理工大学 | Analytical model construction and verification method for revealing abrasion difference of cutter teeth of high-feed milling cutter |
| CN111002104A (en) * | 2019-12-16 | 2020-04-14 | 哈尔滨理工大学 | Method for detecting and calculating frictional wear boundary of rear cutter face of cutter tooth of high-feed milling cutter |
| CN111644900A (en) * | 2020-05-21 | 2020-09-11 | 西安交通大学 | Tool damage real-time monitoring method based on spindle vibration characteristic fusion |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE50001492D1 (en) * | 2000-05-15 | 2003-04-24 | Prometec Gmbh | Method and device for monitoring the state of wear of a tool |
| KR100517880B1 (en) * | 2003-11-20 | 2005-09-28 | 학교법인 포항공과대학교 | Off-line feed rate scheduling for reduction of machining time and enhancement of machining accuracy in cnc machining |
| EP3074748B1 (en) * | 2013-11-29 | 2021-10-13 | Safran Aircraft Engines | Method and device for determining the wear of a cutting tool flank |
| US11630435B2 (en) * | 2018-10-12 | 2023-04-18 | Tata Consultancy Services Limited | Method and system for monitoring tool wear to estimate RUL of tool in machining |
-
2020
- 2020-12-07 CN CN202011439558.1A patent/CN112372372B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5251144A (en) * | 1991-04-18 | 1993-10-05 | Texas Instruments Incorporated | System and method utilizing a real time expert system for tool life prediction and tool wear diagnosis |
| CN109434564A (en) * | 2018-12-21 | 2019-03-08 | 哈尔滨理工大学 | A kind of cutter wear state monitoring method based on deep neural network |
| CN109940460A (en) * | 2019-04-10 | 2019-06-28 | 哈尔滨理工大学 | The detection method of milling machined surface geometric error distribution character under a kind of effect of vibration |
| CN109940461A (en) * | 2019-04-10 | 2019-06-28 | 哈尔滨理工大学 | A kind of detection method of high feed milling cutter cutter tooth wear of the tool flank characteristic |
| CN110103079A (en) * | 2019-06-17 | 2019-08-09 | 中国科学院合肥物质科学研究院 | The on-line monitoring method of tool wear in a kind of micro- milling process |
| CN110704969A (en) * | 2019-09-24 | 2020-01-17 | 哈尔滨理工大学 | Analytical model construction and verification method for revealing abrasion difference of cutter teeth of high-feed milling cutter |
| CN111002104A (en) * | 2019-12-16 | 2020-04-14 | 哈尔滨理工大学 | Method for detecting and calculating frictional wear boundary of rear cutter face of cutter tooth of high-feed milling cutter |
| CN111644900A (en) * | 2020-05-21 | 2020-09-11 | 西安交通大学 | Tool damage real-time monitoring method based on spindle vibration characteristic fusion |
Non-Patent Citations (1)
| Title |
|---|
| 高进给铣刀刀齿后刀面磨损差异性研究;左林晗;《中国优秀硕士学位论文全文数据库 工程科技I辑（月刊）》;20190815(第08期);7-65 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112372372A (en) | 2021-02-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112372372B (en) | Efficient milling cutter accumulated friction wear boundary identification and verification method | |
| US11761930B2 (en) | Prediction method of part surface roughness and tool wear based on multi-task learning | |
| CN104392090B (en) | The construction method of aluminum alloy materials end mill cutting force and machining distorted pattern | |
| Desai et al. | On cutter deflection surface errors in peripheral milling | |
| CN102566492B (en) | Method for forecasting maximum milling force for plunge milling of metal difficult-to-cut materials | |
| CN111002104B (en) | Method for detecting and calculating frictional wear boundary of rear cutter face of cutter tooth of high-feed milling cutter | |
| CN114429064A (en) | Identification method for fractal characteristics of friction boundary of rear cutter face of cutter tooth of high-energy-efficiency milling cutter | |
| Ortiz-de-Zarate et al. | Experimental and FEM analysis of surface integrity when broaching Ti64 | |
| Denkena et al. | Kinematic and stochastic surface topography of machined TiAl6V4-parts by means of ball nose end milling | |
| CN105269402A (en) | Method for predicating surface roughness of titanium alloy material based on milling | |
| Shandilya et al. | Prediction of surface roughness during wire electrical discharge machining of SiCp/6061 Al metal matrix composite | |
| CN107168245A (en) | A kind of accurate Forecasting Methodology of chamfered edge circular bit cutting force for considering cutting edge effect | |
| Kurra et al. | Multi-objective optimisation of single point incremental sheet forming using Taguchi-based grey relational analysis | |
| CN105242546B (en) | A kind of High Speed Milling Force modeling method based on material property | |
| Jia et al. | A new cutting force prediction method in ball-end milling based on material properties for difficult-to-machine materials | |
| CN112528535B (en) | Tongue-and-groove broaching process simulation analysis method based on thermal-force-flow multi-field coupling | |
| Liu et al. | Surface roughness prediction in ball screw whirlwind milling considering elastic-plastic deformation caused by cutting force: Modelling and verification | |
| CN109857061A (en) | A kind of workpiece surface residual stress regulation method based on thermal influence zone | |
| CN109648397B (en) | Broach life prediction method based on correlation between width of cutting edge strip and broaching load | |
| Hilligardt et al. | A new approach for local cutting force modeling enabling the transfer between different milling conditions and tool geometries | |
| Guo et al. | Studies on residue stress and deformation behavior of GCr15 subjected to whirlwind milling | |
| CN104476326A (en) | Ceramic tool groove wear prediction method | |
| Essa | Finite element prediction of deformation mechanics in incremental forming processes | |
| Jiang | Intelligent prediction of surface roughness of milling aluminium alloy based on least square support vector machine | |
| Schulze et al. | FE analysis on the influence of sequential cuts on component conditions for different machining strategies |
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 |