CN103646141A - Cutting force modeling method for flat bottom spiral end mill orthogonal turning milling shaft parts - Google Patents

Cutting force modeling method for flat bottom spiral end mill orthogonal turning milling shaft parts Download PDF

Info

Publication number
CN103646141A
CN103646141A CN201310655052.8A CN201310655052A CN103646141A CN 103646141 A CN103646141 A CN 103646141A CN 201310655052 A CN201310655052 A CN 201310655052A CN 103646141 A CN103646141 A CN 103646141A
Authority
CN
China
Prior art keywords
msub
mtd
cutting
mrow
mtr
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
CN201310655052.8A
Other languages
Chinese (zh)
Other versions
CN103646141B (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
Original Assignee
Huazhong University of Science and Technology
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 filed Critical Huazhong University of Science and Technology
Priority to CN201310655052.8A priority Critical patent/CN103646141B/en
Publication of CN103646141A publication Critical patent/CN103646141A/en
Application granted granted Critical
Publication of CN103646141B publication Critical patent/CN103646141B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Numerical Control (AREA)

Abstract

The invention discloses a cutting force modeling method for flat bottom spiral end mill orthogonal turning milling shaft parts. The method includes the steps that a spiral end mill coordinate system and a shaft part coordinate system are respectively established, and coordinate transformation between the two coordinate systems is determined; in the coordinate systems, an orthogonal turning milling cutting motion trail is clear, and an axial feed, a radial feed and a total feed are determined; according to the orthogonal turning milling cutting motion trail obtained through calculation and the determined feeds, an orthogonal turning milling cutting geometry is determined; infinitesimal orthogonal turning milling cutting force calculation formulas are set up for discrete layers respectively, whether the discrete layers participate in cutting or not is judged according to the determined orthogonal turning milling cutting thicknesses, cutting depths and orthogonal turning milling cutting-in cutting-out angles, limits of integration of infinitesimal cutting force are determined, integrals participate in cutting infinitesimal cutting force, and the total cutting force of the current cutting state is obtained. According to the model, cutting force of the orthogonal turning milling shaft parts can be accurately obtained, and then technological parameter optimization on the basis of the cutting force of the orthogonal turning milling shaft parts is further achieved.

Description

Cutting force modeling method for flat-bottom spiral end mill orthogonal turning milling shaft part
Technical Field
The invention relates to the field of metal cutting machining, in particular to a cutting force modeling method for a flat-bottom spiral end mill orthogonal turn-milling shaft part, which is suitable for machining the flat-bottom spiral end mill orthogonal turn-milling shaft part.
Background
In recent years, with the development of the national aviation industry, the requirements on the functions and the performances of products are diversified by the research and the application of new models and new materials, and the integral forgings of various new materials are widely applied to aircraft landing gears, so that the processing difficulty of large-scale complex shaft parts is improved. The turning and milling machining center is a high-grade numerical control machine tool, can realize five-axis linkage machining, can finish machining of complex curved surfaces, and can obtain good surface quality, wherein orthogonal turning and milling machining is a machining mode for finishing forming and machining of complex shaft parts by combining workpiece rotary motion and cutter rotary motion, and has the advantage of meeting the requirement of one-time clamping and multi-process machining of large complex shaft parts. In order to orthogonally turn-mill the shaft parts with the maximum material removal rate under the normal working condition, the cutting force of the orthogonal turn-milling shaft parts needs to be modeled, the cutting force of the orthogonal turn-milling shaft parts under the given cutting parameters is predicted, and reasonable processing parameters are arranged.
Common models for cutting force prediction are: the model comprises an empirical model, a finite element numerical model, an artificial intelligence model and a mechanical model. The empirical model establishes a functional relationship between the cutting force and the cutting parameters by using a statistical principle according to the cutting parameters and the cutting force data obtained by the cutting force test experiment, and a large amount of data is needed for ensuring the accuracy of the model. The finite element model is suitable for simulating the cut-in and cut-out stages of cutting and discontinuous machining, but under the condition that a large amount of plastic deformation occurs in the machining, the grid deformation is too large, so that the simulation error is increased. The application potential of the artificial intelligence modeling technology in metal cutting research is huge, but more theoretical and experimental researches are needed for the problems of the precision, the application range and the like of the simulation calculation. The mechanical model takes the chip as a research entry point, and analyzes the rules of stress, strain, friction force and cutting force in the chip forming process through modeling and simulation of the chip forming process. The mechanical force model is a cutting force model which is mature and widely used at present, and one important content is the determination of cutting geometry. The cutting geometry of the orthogonal turn-milling has the characteristics of variable cutting thickness and variable cutting depth, the cutting thickness and the cutting depth are changed along with the change of the rotation angle of the cutter, and an accurate orthogonal turn-milling cutting geometry model needs to be established for modeling the orthogonal turn-milling cutting force.
Disclosure of Invention
The invention aims to provide a cutting force modeling method for a flat-bottom spiral end mill orthogonal turning and milling shaft part, which can accurately predict the cutting force of a steady-state orthogonal turning and milling shaft part through the analysis of the orthogonal turning and milling cutting geometry in the process of orthogonal turning and milling the shaft part, thereby optimizing processing parameters in the orthogonal turning and milling shaft part and improving the processing efficiency.
The specific technical scheme adopted for achieving the purpose of the invention is as follows:
a cutting force modeling method for a flat-bottom spiral end mill orthogonal turn-milling shaft part specifically comprises the following steps:
s1: respectively establishing a spiral end mill coordinate system and a shaft part coordinate system, and determining coordinate transformation between the two coordinate systems;
s2: in the established coordinate system, the orthogonal turn-milling cutting motion track is defined, and axial feed, radial feed and total feed are determined;
s3: in the established coordinate system, determining the orthogonal turn-milling cutting geometry according to the calculated orthogonal turn-milling motion track and the determined feed;
a: axially dispersing the spiral end mill into a plurality of layers, and determining the orthogonal turn-milling chip thickness of each discrete layer in an established coordinate system according to the obtained feeding and orthogonal turn-milling motion tracks;
b: determining the orthogonal turning and milling cutting depth of each corner position, and judging whether the discrete layer is within the cutting depth range or not;
c: determining a cutting-in and cutting-out angle of orthogonal turn-milling cutting according to the determined thickness and cutting depth of the orthogonal turn-milling cutting, and judging whether each cutter tooth is in cutting at each corner position;
s4: and establishing a infinitesimal orthogonal turning and milling cutting force calculation formula for each discrete layer, judging whether each discrete layer participates in cutting or not according to the determined orthogonal turning and milling chip thickness, the determined cutting depth and the orthogonal turning and milling cut-in and cut-out angle, determining an integral limit of the infinitesimal cutting force, integrating the infinitesimal cutting force participating in cutting, and obtaining the total cutting force of the current cutting state.
The coordinate transformation between the coordinate system of the spiral end mill and the coordinate system of the shaft part in the orthogonal turning, milling and cutting is obtained by the following formula:
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mi>w</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mi>w</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mi>w</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>sin</mi> <mi>&theta;</mi> </mtd> <mtd> <mi>cos</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mi>cos</mi> <mi>&theta;</mi> </mtd> <mtd> <mi>sin</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mi>F</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mi>F</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mi>F</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mrow> <mo>(</mo> <mi>R</mi> <mo>-</mo> <msub> <mi>a</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&theta;</mi> <mo>+</mo> <mi>e</mi> <mi>sin</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>(</mo> <mi>R</mi> <mo>-</mo> <msub> <mi>a</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mi>sin</mi> <mi>&theta;</mi> <mo>-</mo> <mi>e</mi> <mi>cos</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&lambda;N</mi> <msub> <mi>f</mi> <mi>a</mi> </msub> <mi>t</mi> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
wherein x isw,ywAnd zwIs the coordinate value, x, of the midpoint in the coordinate system of the shaft partF,yFAnd zFIs the coordinate value of the midpoint of the tool coordinate system, R is the workpiece radius, apFor a given depth of cut, theta is the workpiece angle, and e is z of the tool coordinate system with t equal to 0FAxial relative zwAxis in ywEccentricity in the negative direction of the axis, t cutting time, faFor axial feed, N is the number of teeth of the cutter, and lambda is the ratio of the rotational speed of the milling cutter to the rotational speed of the workpiece.
The orthogonal turn-milling total feed is calculated by the following formula:
f z = f a 2 + f t 2
wherein f istIn order to realize the tangential feeding,
Figure BDA0000430845480000033
n is the number of teeth of the cutter; f. ofaIn order to achieve the axial feed,
Figure BDA0000430845480000034
aelambda is the ratio of the rotational speed of the milling cutter to the rotational speed of the workpiece,nFfor rotational speed of milling cutter, nwIs the workpiece rotation speed, R is the workpiece radius, betawFor helical conveying of milling cutters relative to the workpieceThe lift angle of the moving track.
The orthogonal turning and milling chip thickness is as follows:
Figure BDA0000430845480000041
wherein, aeTo cut the width, ae=fa·λ·N·sinβw,βwIs the lead angle of the milling cutter relative to the spiral motion track of the workpiece, <math> <mrow> <msub> <mi>&beta;</mi> <mi>w</mi> </msub> <mo>=</mo> <mi>arctan</mi> <mfrac> <msub> <mi>f</mi> <mi>t</mi> </msub> <msub> <mi>f</mi> <mi>a</mi> </msub> </mfrac> <mo>,</mo> <mi>&psi;</mi> <mo>=</mo> <mfrac> <mi>&pi;</mi> <mn>2</mn> </mfrac> <mo>-</mo> <msub> <mi>&beta;</mi> <mi>w</mi> </msub> <mo>,</mo> </mrow> </math>
Figure BDA0000430845480000043
for the cut-in angle determined by the chip thickness,
Figure BDA0000430845480000044
for the cutting angle determined by the chip thickness,
Figure BDA0000430845480000045
<math> <mrow> <mi>x</mi> <mo>=</mo> <msqrt> <msubsup> <mi>f</mi> <mi>z</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msup> <mi>r</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>2</mn> <msub> <mi>f</mi> <mi>z</mi> </msub> <mo>&CenterDot;</mo> <mi>r</mi> <mo>&CenterDot;</mo> <mi>cos</mi> <mi>&theta;</mi> </msqrt> <mo>,</mo> <mi>&theta;</mi> <mo>=</mo> <mi>arcsin</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>r</mi> <mo>-</mo> <msub> <mi>a</mi> <mi>e</mi> </msub> </mrow> <mi>r</mi> </mfrac> <mo>)</mo> </mrow> <mo>.</mo> </mrow> </math>
the orthogonal turning and milling cutting depth is as follows:
wherein R is the radius of the tool, R is the radius of the workpiece,
Figure BDA0000430845480000048
is the contact angle of the ith infinitesimal layer on the blade j, apZ of the tool coordinate system for a given depth of cut with t equal to 0FAxial relative zwAxis in ywThe eccentricity existing in the negative direction of the shaft, and t is the cutting time.
The cutting-in and cutting-out angle of the orthogonal turn-milling cutting is as follows:
corner of cut
Figure BDA0000430845480000049
Cutting corner
Wherein,
Figure BDA00004308454800000411
for the cut angle determined by the depth of cut,
Figure BDA00004308454800000412
Figure BDA00004308454800000413
for the cutting angle determined by the cutting depth,
Figure BDA0000430845480000051
βwthe angle of the spiral line is increased, <math> <mrow> <msub> <mi>&beta;</mi> <mi>w</mi> </msub> <mo>=</mo> <mi>arctan</mi> <mfrac> <msub> <mi>f</mi> <mi>t</mi> </msub> <msub> <mi>f</mi> <mi>a</mi> </msub> </mfrac> <mo>,</mo> <mi>&psi;</mi> <mo>=</mo> <mfrac> <mi>&pi;</mi> <mn>2</mn> </mfrac> <mo>-</mo> <msub> <mi>&beta;</mi> <mi>w</mi> </msub> <mo>,</mo> </mrow> </math> zifor the ith infinitesimal layer in the tool axis zFThe coordinate values of (2).
The instantaneous orthogonal turning and milling cutting force is calculated by a mechanical model, an orthogonal turning and milling infinitesimal cutting force calculation formula is firstly obtained, then whether the discrete layer participates in cutting or not is judged according to the determined orthogonal turning and milling chip thickness, the determined cutting depth and the orthogonal turning and milling cut-in and cut-out angle, the integral limit of the infinitesimal cutting force is determined, the infinitesimal cutting force participating in cutting is integrated, and the total cutting force in the current cutting state is obtained.
Wherein, the algorithm of the total cutting force of the current cutting state is as follows:
according to the mechanical cutting force model of Altintas, the cutting force of the ith element of the jth cutter tooth can be expressed as:
Figure BDA0000430845480000053
in the formula, Ktc,KrcAnd KacAre respectively a shearing action pairCoefficient of action of tangential, radial and axial cutting forces, Kte,KreAnd KaeRespectively the action coefficients of the friction action of the cutting edge to tangential, radial and axial cutting forces,
Figure BDA0000430845480000054
is the chip thickness of the ith element of the jth cutter tooth.
Figure BDA0000430845480000055
Is the instantaneous contact angle of the ith element of the jth cutter tooth;
the infinitesimal tangential force, the infinitesimal radial force and the infinitesimal axial force in the formula (1) are converted into a milling cutter coordinate system xFyFzFThe infinitesimal cutting force component along the x-axis, the infinitesimal cutting force component along the y-axis, and the infinitesimal cutting force component along the z-axis in (b) to obtain:
Figure BDA0000430845480000056
integrating the infinitesimal cutting force along the axial direction of the milling cutter, and obtaining the cutting force according to the equations (1) and (2):
Figure BDA0000430845480000061
in the formula, apIn order to be the nominal depth of cut,
Figure BDA0000430845480000063
is the cutting-in angle of the milling cutter,
Figure BDA0000430845480000064
is the milling cutter cut angle.
The instantaneous cutting forces of all the helical cutting edges are summed to obtain the total instantaneous cutting force acting on the milling cutter:
<math> <mrow> <msub> <mi>F</mi> <mi>x</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>F</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>F</mi> <mi>y</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>F</mi> <mrow> <mi>y</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <msub> <mi>F</mi> <mi>z</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>F</mi> <mrow> <mi>z</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow> </math>
according to the cutting force modeling method for the orthogonal turning and milling shaft parts of the flat-bottom spiral end mill, the accurate cutting force model of the orthogonal turning and milling shaft parts can be obtained, the model fully considers the movement track, the cutting thickness, the cutting depth and the cutting-in and cutting-out angle of orthogonal turning and milling cutting in the process of orthogonal turning, milling and cutting the shaft parts, and the cutting force of the stable orthogonal turning and milling shaft parts can be accurately predicted. The method comprises the steps of dispersing a helical cutting edge into a plurality of tiny cutter tooth infinitesimal elements in the axial direction, establishing a calculation formula of instantaneous cutting force of orthogonal turning and milling shaft parts according to a mechanical force model, and predicting the orthogonal turning and milling cutting force under a given point cutting parameter by using a chip geometric model of the orthogonal turning and milling shaft parts.
According to the method, the cutting force of the orthogonal turning and milling shaft part is predicted according to the known material model, the geometric parameters of the cutter and the designed process parameters, and compared with the measured value of the cutting force in the experimental process of the orthogonal turning and milling shaft part adopting the same workpiece material, cutter and process parameters, the predicted cutting force can be well matched with the actually measured cutting force. The method can predict the cutting force of the steady-state orthogonal turning and milling shaft part, is beneficial to determining the processing parameters by technicians in the field and improving the processing efficiency under the normal working condition.
Drawings
Fig. 1 is a schematic diagram of orthogonal turn-milling coordinate transformation according to an embodiment of the present invention.
FIG. 2 is a schematic diagram of the orthogonal turn-milling feed motion of an embodiment of the present invention.
FIG. 3 is a discrete schematic view of a flat bottom helical end mill according to an embodiment of the present invention.
FIG. 4 is a schematic illustration of orthogonal turn-milling chip thicknesses according to an embodiment of the present invention.
FIG. 5 is a schematic diagram of the orthogonal turn-milling depth of cut according to an embodiment of the present invention.
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.
The cutting force modeling method for the orthogonal turn-milling shaft part of the flat-bottom spiral end mill comprises the following steps:
(one) coordinate transformation calculation
Orthogonal turn-milling is a combined motion of workpiece rotation and milling cutter rotation to complete the machining of the workpiece, and to obtain a better machined surface and increase the feed, the cutter may be offset so that the axis of the cutter is not coplanar with and orthogonal to the axis of the workpiece, and machining is performed in a downward-offset manner, as shown in fig. 1, for ease of analysis, assuming that the workpiece is stationary, the milling cutter is rotated at the workpiece rotation speed nwIn fig. 1(a), the workpiece to be machined by the milling cutter is viewed from the shank along the milling cutter axis at the time t equal to 0, and fig. 1(b) shows the milling cutter rotated by θ (that is, after the milling cutter is rotated by θ degrees)
Figure BDA0000430845480000071
Time) observing the milling cutter along the workpiece axis, establishing a workpiece coordinate system xwywzwWherein z iswThe axis being the workpiece axis, xwZ when t and t are equal to 0FWhen the direction of the axes is consistent and t is 0, z of the tool coordinate systemFAxial relative zwAxis in ywThe negative axis direction has a lower eccentricity e, and the direction in which the lower eccentricity increases is vertically downward as viewed in fig. 1(a), coinciding with the direction of rotation of the workpiece.
Any point (x) in the milling cutter coordinate systemF,yF,zF) Transforming into the coordinate system of the workpiece to obtain coordinates (x)w,yw,zw):
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mi>w</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mi>w</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mi>w</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>sin</mi> <mi>&theta;</mi> </mtd> <mtd> <mi>cos</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mi>cos</mi> <mi>&theta;</mi> </mtd> <mtd> <mi>sin</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mi>F</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mi>F</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mi>F</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mrow> <mo>(</mo> <mi>R</mi> <mo>-</mo> <msub> <mi>a</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&theta;</mi> <mo>+</mo> <mi>e</mi> <mi>sin</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>(</mo> <mi>R</mi> <mo>-</mo> <msub> <mi>a</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mi>sin</mi> <mi>&theta;</mi> <mo>-</mo> <mi>e</mi> <mi>cos</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&lambda;N</mi> <msub> <mi>f</mi> <mi>a</mi> </msub> <mi>t</mi> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math>
Wherein R is the workpiece radius, apTo a depth of cut, nFBy milling cutter speed, milling cutter-to-workpiece speed ratiofaThe axial feed amount of the milling cutter per tooth along the axial direction of the workpiece.
(II) determining the orthogonal turn-milling movement locus
In the process of orthogonal turning, milling and cutting shaft parts, the rotary motion of the shaft parts can generate feeding along the tangential direction of the shaft parts, and the feeding of a milling cutter along the axial direction of the shaft parts is combined to form the total feeding of the milling cutter, so that the motion track of the milling cutter relative to a workpiece is a continuous spiral line, and a helix angle exists. And respectively calculating tangential feed, axial feed and total feed of the orthogonal turn-milling machine in the determined coordinate system, and determining the helix angle of the motion track.
In orthogonal turn-milling, the milling cutter is fed f relative to each tooth of the workpiecezFed tangentially f by each toothtAnd axial feed per tooth faSynthesis, as shown in FIG. 2, of the computational expression
The path of movement of the milling cutter relative to the workpiece is a continuous helix, as shown in fig. 2, the helix angle βwFeeding tangentially to each tooth of the milling cutter ftAxial feed per tooth faCorrelation, expressed as follows <math> <mrow> <msub> <mi>&beta;</mi> <mi>w</mi> </msub> <mo>=</mo> <mi>arctan</mi> <mfrac> <msub> <mi>f</mi> <mi>t</mi> </msub> <msub> <mi>f</mi> <mi>a</mi> </msub> </mfrac> <mo>,</mo> </mrow> </math> Tangential feed per tooth of milling cutter <math> <mrow> <msub> <mi>f</mi> <mi>t</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> <mrow> <mo>(</mo> <mi>R</mi> <mo>-</mo> <msub> <mi>a</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&lambda;</mi> <mo>&CenterDot;</mo> <mi>N</mi> </mrow> </mfrac> <mo>,</mo> </mrow> </math> Axial feed per tooth of milling cutter f a = v fa n w , vfaIs the milling cutter axial feed speed.
Due to the helix angle, a rotary conversion seat is establishedMarker is x'Fy′Fz′FAnd coordinate system xFyFzFThe included angle therebetween is psi,
Figure BDA0000430845480000086
in a coordinate system of x'Fy′Fz′FAnd analyzing the thickness of the cutting chips of the shaft parts orthogonally turned and milled by the circumferential edge of the spiral end mill.
(III) determination of orthogonal turning and milling chip geometry
In an established coordinate system, axially dispersing a spiral end mill into a plurality of layers, and determining the orthogonal turn-milling chip thickness of each discrete layer at each corner position of the cutter according to the obtained feeding and orthogonal turn-milling motion tracks; then determining the orthogonal turning and milling cutting depth of each corner position, and judging whether the discrete layer is within the cutting depth range or not; and finally, determining a cutting-in and cutting-out angle of orthogonal turn-milling cutting according to the determined thickness and cutting depth of the orthogonal turn-milling cutting, and judging whether each cutter tooth is in cutting at each corner position.
As shown in fig. 3, the cutting edge of the helical end mill is a helical line, and it is very difficult to directly calculate the cutting force caused by orthogonal turning and milling of shaft parts, and the problem can be solved well by discretization. Therefore, as shown in fig. 3, the helical edge of the milling cutter is axially dispersed into several layers, that is, a plurality of plane cutting milling cutters perpendicular to the central axis of the milling cutter and having a distance Δ d, the sheet sandwiched between two adjacent planes is a discrete layer, the cutting of each discrete layer tooth infinitesimal (for example, when viewed along the axis of the milling cutter, four teeth are uniformly distributed on each discrete layer of a four-edge milling cutter, and each tooth is a tooth infinitesimal) can be regarded as the oblique cutting with the cutting edge inclination i equal to the milling cutter helix angle β, the orthogonal turning and milling cutting force calculation formula is derived from a mechanical force model for each infinitesimal cutting layer, and the integral is obtained in the cutting range to obtain the instantaneous orthogonal turning and milling cutting force. A in FIG. 3pFor cutting depth, Δ d is the discrete layer thickness,
Figure BDA0000430845480000091
n is a radicalAnd (5) the number of scattering layers.
(a) Calculating the orthogonal turn-milling chip thickness
To facilitate the analysis of the chip thickness of the peripheral edge of the helical end mill, a local coordinate system x 'is established in FIG. 4'Fy′Fz′FWherein x'FThe axis being perpendicular to the helix and to the coordinate system xFyFzFX ofFThe included angle between the axes is psi,
Figure BDA0000430845480000092
z′Faxis and zFAre superposed and y'FObtained by the right hand rule.
As shown in FIG. 4, in the milling cutter local coordinate system x'Fy′Fz′FIn the method, a calculation model of the orthogonal turning and milling chip thickness of the spiral end mill is established, and a coordinate system x ″, isFy′′Fz′′FThe milling cutter is rotated about the cutter axis, the blade continuously cuts into and cuts out the workpiece from y'FThe shaft calculates the cut-in and cut-out angle and the cut-in angle in the clockwise direction
Figure BDA0000430845480000093
And cutting out of cornersCan be given by:
Figure BDA0000430845480000095
Figure BDA0000430845480000096
in the formula, aeTo cut the width, ae=fa·λ·N·sinβwX 'is a coordinate system'Fy′Fz′FTo xFyFzFIn the coordinate system, x is obtainedFyFzFAngle of entry in a coordinate system
Figure BDA0000430845480000101
And cutting out of corners
Figure BDA0000430845480000102
Figure BDA0000430845480000103
Figure BDA0000430845480000104
When in use
Figure BDA0000430845480000105
When the ith element of the jth cutter tooth is in a cutting state, the chip thickness is calculated into two parts of CA and AB, and the following formula is shown as follows:
Figure BDA0000430845480000106
wherein, aeTo cut the width, ae=fa·λ·N·sinβw,βwIs the lead angle of the milling cutter relative to the spiral motion track of the workpiece, <math> <mrow> <msub> <mi>&beta;</mi> <mi>w</mi> </msub> <mo>=</mo> <mi>arctan</mi> <mfrac> <msub> <mi>f</mi> <mi>t</mi> </msub> <msub> <mi>f</mi> <mi>a</mi> </msub> </mfrac> <mo>,</mo> <mi>&psi;</mi> <mo>=</mo> <mfrac> <mi>&pi;</mi> <mn>2</mn> </mfrac> <mo>-</mo> <msub> <mi>&beta;</mi> <mi>w</mi> </msub> <mo>,</mo> </mrow> </math>
Figure BDA0000430845480000108
for the cut-in angle determined by the chip thickness,
Figure BDA0000430845480000109
for the cutting angle determined by the chip thickness,
Figure BDA00004308454800001011
<math> <mrow> <mi>x</mi> <mo>=</mo> <msqrt> <msubsup> <mi>f</mi> <mi>z</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msup> <mi>r</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>2</mn> <msub> <mi>f</mi> <mi>z</mi> </msub> <mo>&CenterDot;</mo> <mi>r</mi> <mo>&CenterDot;</mo> <mi>cos</mi> <mi>&theta;</mi> </msqrt> <mo>,</mo> <mi>&theta;</mi> <mo>=</mo> <mi>arcsin</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>r</mi> <mo>-</mo> <msub> <mi>a</mi> <mi>e</mi> </msub> </mrow> <mi>r</mi> </mfrac> <mo>)</mo> </mrow> <mo>.</mo> </mrow> </math>
(b) calculating orthogonal turn-milling cutting depth
Different from the common end milling, when the shaft parts are cut by orthogonal turn milling, the cutting depth a is constantly changed and is less than or equal to the given nominal cutting depth apIt is assumed that the contact angle of the ith bin of the instantaneous j teeth is, as shown in FIG. 5
Figure BDA00004308454800001013
When the corresponding milling cutter cutting depth is
Figure BDA00004308454800001014
It can be seen from the infinitesimal layering that the cutting is only participated when the infinitesimal is within the cutting depth range of the milling cutter, and the intersection point of the outer circle of the milling cutter and the outer circle of the workpiece is A at the moment, and the point is in the coordinate system x of the workpiecewywzwAnd milling cutter coordinate system xFyFzFRespectively is (x)A,w,yA,w,zA,w)、(xA,F,yA,F,zA,F) Wherein
Figure BDA0000430845480000111
Since point a is located on both the outer circle of the workpiece and the peripheral edge of the milling cutter, it can be deduced that the milling cutter cutting depth at this time is:
Figure BDA0000430845480000112
(c) calculating the cut-in and cut-out angle of orthogonal turning and milling shaft parts
And calculating the cut-in and cut-out angle of the orthogonal turning and milling shaft type part respectively according to the determined chip thickness and the determined cutting depth of the orthogonal turning and milling shaft type part, wherein the cut-in angle of the orthogonal turning and milling shaft type part is the larger value of the two calculated cut-in angles, and the cut-out angle is the smaller value of the two calculated cut-out angles.
Figure BDA0000430845480000113
Figure BDA0000430845480000114
Wherein,
Figure BDA0000430845480000115
for the cut angle determined by the depth of cut,
Figure BDA0000430845480000117
for the cutting angle determined by the cutting depth,
(IV) calculating the instantaneous cutting force of orthogonal turning and milling cutting shaft parts
According to the mechanical cutting force model of Altintas, the cutting force of the ith element of the jth cutter tooth can be expressed as:
Figure BDA0000430845480000119
in the formula, Ktc,KrcAnd KacCoefficient of action of the shearing action on the tangential, radial and axial cutting forces, Kte,KreAnd KaeThe coefficients of the tangential, radial and axial cutting forces are the blade friction.
Figure BDA00004308454800001110
Is the chip thickness of the ith element of the jth cutter tooth.
Figure BDA00004308454800001111
Is the instant contact angle of the ith element of the jth cutter tooth due to the milling cutter helix angle betaFThe instantaneous contact angle of different infinitesimal elements changes with the change of the cutting depth, and the instantaneous reference contact angle of a certain spiral blade bottom endpoint is assumed to be determined asThe instant contact angle of the ith element of the jth cutter tooth is
Figure BDA0000430845480000122
Wherein
Figure BDA0000430845480000123
The angle between the teeth is the angle between the teeth,
Figure BDA0000430845480000124
n is the number of teeth of the knife,
Figure BDA0000430845480000125
and r is the radius of the milling cutter.
The infinitesimal tangential force, the infinitesimal radial force and the infinitesimal axial force in the formula (1) are converted into a milling cutter coordinate system xFyFzFThe infinitesimal cutting force component along the x-axis, the infinitesimal cutting force component along the y-axis, and the infinitesimal cutting force component along the z-axis in (b) to obtain:
Figure BDA0000430845480000126
in order to obtain the total cutting force generated by the helical cutting edge j, the infinitesimal cutting force is integrated in the milling cutter axial direction in the range where the helical cutting edge participates in cutting, and the cutting force is obtained according to equations (10) and (11):
Figure BDA0000430845480000127
Figure BDA0000430845480000128
in the formula, apIs nominally cutDepth of cut (transient milling cutter depth of cut changes over time due to the kinematic nature of orthogonal turn-milling, which is less than or equal to the nominal depth of cut),
Figure BDA0000430845480000129
is the cutting-in angle of the milling cutter,is the milling cutter cut angle.
The instantaneous cutting forces of all the helical cutting edges are summed to obtain the total instantaneous cutting force acting on the milling cutter:
<math> <mrow> <msub> <mi>F</mi> <mi>x</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>F</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>F</mi> <mi>y</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>F</mi> <mrow> <mi>y</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <msub> <mi>F</mi> <mi>z</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>F</mi> <mrow> <mi>z</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow> </math>
the accurate instantaneous cutting force of the flat-bottom spiral end mill orthogonal turning-milling shaft type part can be obtained through the four steps.
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 (9)

1. A cutting force modeling method for a flat-bottom spiral end mill orthogonal turn-milling shaft part is characterized by comprising the following steps:
s1: respectively establishing a spiral end mill coordinate system and a shaft part coordinate system, and determining coordinate transformation between the two coordinate systems;
s2: in the established coordinate system, the orthogonal turn-milling cutting motion track is defined, and axial feed, radial feed and total feed are determined;
s3: in the established coordinate system, determining the orthogonal turn-milling chip geometry according to the calculated orthogonal turn-milling motion track and the determined feed:
a: axially dispersing the spiral end mill into a plurality of layers, and determining the orthogonal turn-milling chip thickness of each discrete layer in an established coordinate system according to the obtained feeding and orthogonal turn-milling motion tracks;
b: determining the orthogonal turning and milling cutting depth of each corner position, and judging whether the discrete layer is in the cutting depth range;
c: determining a cutting-in and cutting-out angle of orthogonal turn-milling cutting according to the determined thickness and cutting depth of the orthogonal turn-milling cutting, and judging whether each cutter tooth is in cutting at each corner position;
s4: and establishing a infinitesimal orthogonal turning and milling cutting force calculation formula for each discrete layer, judging whether each discrete layer participates in cutting or not according to the determined orthogonal turning and milling chip thickness, the determined cutting depth and the orthogonal turning and milling cut-in and cut-out angle, determining an integral limit of the infinitesimal cutting force, integrating the infinitesimal cutting force participating in cutting, and obtaining the total cutting force of the current cutting state.
2. The modeling method for the cutting force of the flat-bottom helical end mill orthogonal turning shaft part according to claim 1, wherein the coordinate transformation between the helical end mill coordinate system and the shaft part coordinate system in the orthogonal turning cutting is calculated by the following formula:
<math> <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mi>w</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mi>w</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mi>w</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>sin</mi> <mi>&theta;</mi> </mtd> <mtd> <mi>cos</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mi>cos</mi> <mi>&theta;</mi> </mtd> <mtd> <mi>sin</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mi>F</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mi>F</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mi>F</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mrow> <mo>(</mo> <mi>R</mi> <mo>-</mo> <msub> <mi>a</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&theta;</mi> <mo>+</mo> <mi>e</mi> <mi>sin</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>(</mo> <mi>R</mi> <mo>-</mo> <msub> <mi>a</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mi>sin</mi> <mi>&theta;</mi> <mo>-</mo> <mi>e</mi> <mi>cos</mi> <mi>&theta;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&lambda;N</mi> <msub> <mi>f</mi> <mi>a</mi> </msub> <mi>t</mi> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
wherein x isw,ywAnd zwIs the coordinate value, x, of the midpoint in the coordinate system of the shaft partF,yFAnd zFIs the coordinate value of the midpoint of the tool coordinate system, R is the workpiece radius, apFor a given depth of cut, theta is the workpiece angle, and e is z of the tool coordinate system with t equal to 0FAxial relative zwAxis in ywEccentricity in the negative direction of the axis, t cutting time, faFor axial feeding, N is the number of teeth of the cutter, lambda is the rotation speed ratio of the milling cutter to the workpiece,
Figure FDA0000430845470000021
nFfor rotational speed of milling cutter, nwThe workpiece rotation speed.
3. The modeling method for the cutting force of the flat-bottom helical end mill orthogonal turn-milling shaft part according to claim 2, wherein the total orthogonal turn-milling feed is calculated by the following formula:
f z = f a 2 + f t 2
wherein f istIn order to realize the tangential feeding,
Figure FDA0000430845470000023
fain order to achieve the axial feed,
Figure FDA0000430845470000024
vfais the axial feed speed of the milling cutter, R is the workpiece radius, nwIs the rotational speed of the workpiece, apFor a given depth of cut.
4. The method for modeling the cutting force of a flat-bottomed helical end mill orthogonal turning shaft part according to any one of claims 1 to 3, wherein the orthogonal turning chip geometry is determined by a variable chip thickness, a variable cutting depth and an incision angle of orthogonal turning cutting in a determined coordinate system and coordinate transformation.
5. The modeling method for the cutting force of the flat-bottom helical end mill orthogonal turning shaft part according to claim 3, wherein the chip thickness of the orthogonal turning is determined by the tool corner, the feed per tooth, the tool radius and the cutting width, and is obtained by the following calculation formula:
Figure FDA0000430845470000025
wherein, aeTo cut the width, ae=fa·λ·N·sinβw,βwIs the lead angle of the milling cutter relative to the spiral motion track of the workpiece, <math> <mrow> <msub> <mi>&beta;</mi> <mi>w</mi> </msub> <mo>=</mo> <mi>arctan</mi> <mfrac> <msub> <mi>f</mi> <mi>t</mi> </msub> <msub> <mi>f</mi> <mi>a</mi> </msub> </mfrac> <mo>,</mo> <mi>&psi;</mi> <mo>=</mo> <mfrac> <mi>&pi;</mi> <mn>2</mn> </mfrac> <mo>-</mo> <msub> <mi>&beta;</mi> <mi>w</mi> </msub> <mo>,</mo> </mrow> </math>
Figure FDA0000430845470000027
for the cut-in angle determined by the chip thickness,
Figure FDA0000430845470000031
for the cutting angle determined by the chip thickness, <math> <mrow> <mi>x</mi> <mo>=</mo> <msqrt> <msubsup> <mi>f</mi> <mi>z</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msup> <mi>r</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>2</mn> <msub> <mi>f</mi> <mi>z</mi> </msub> <mo>&CenterDot;</mo> <mi>r</mi> <mo>&CenterDot;</mo> <mi>cos</mi> <mi>&theta;</mi> </msqrt> <mo>,</mo> <mi>&theta;</mi> <mo>=</mo> <mi>arcsin</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>r</mi> <mo>-</mo> <msub> <mi>a</mi> <mi>e</mi> </msub> </mrow> <mi>r</mi> </mfrac> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </math> r is the radius of the cutter,the contact angle of the ith infinitesimal layer on blade j.
6. The method for modeling the cutting force of a flat-bottomed helical end mill orthogonal turn-milling shaft part according to any one of claims 2 to 5, wherein the orthogonal turn-milling cutting depth is as follows:
Figure FDA0000430845470000036
wherein R is the radius of the tool, R is the radius of the workpiece,
Figure FDA0000430845470000037
is the contact angle of the ith infinitesimal layer on the blade j, apZ of the tool coordinate system for a given depth of cut with t equal to 0FAxial relative zwAxis in ywThe eccentricity existing in the negative direction of the shaft, and t is the cutting time.
7. The method for modeling the cutting force of the flat-bottom helical end mill orthogonal turn-milling shaft part according to claim 6, wherein the cut-in and cut-out angles of the orthogonal turn-milling cutting are as follows:
corner of cut
Figure FDA0000430845470000038
Cutting corner
Figure FDA0000430845470000039
Wherein,for the cut angle determined by the depth of cut, <math> <mrow> <msub> <mi>&phi;</mi> <mrow> <mi>en</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>=</mo> <mi>arccos</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msqrt> <mn>2</mn> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mi>p</mi> </msub> <mo>-</mo> <msub> <mi>z</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mi>p</mi> </msub> <mo>-</mo> <msub> <mi>z</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <mo>-</mo> <mi>e</mi> </mrow> <mi>r</mi> </mfrac> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </math>
Figure FDA00004308454700000314
for the cutting angle determined by the cutting depth,
Figure FDA00004308454700000312
βwthe angle of the spiral line is increased, <math> <mrow> <msub> <mi>&beta;</mi> <mi>w</mi> </msub> <mo>=</mo> <mi>arctan</mi> <mfrac> <msub> <mi>f</mi> <mi>t</mi> </msub> <msub> <mi>f</mi> <mi>a</mi> </msub> </mfrac> <mo>,</mo> <mi>&psi;</mi> <mo>=</mo> <mfrac> <mi>&pi;</mi> <mn>2</mn> </mfrac> <mo>-</mo> <msub> <mi>&beta;</mi> <mi>w</mi> </msub> <mo>,</mo> </mrow> </math> zifor the ith infinitesimal layer in the tool axis zFThe coordinate values of (2).
8. The cutting force modeling method for the orthogonal turning and milling shaft part of the flat-bottom spiral end mill according to any one of claims 1 to 7, characterized in that the instantaneous orthogonal turning and milling cutting force is calculated by a mechanical model, an orthogonal turning and milling infinitesimal cutting force calculation formula is firstly obtained, then whether a discrete layer participates in cutting or not is judged according to the determined orthogonal turning and milling chip thickness, the determined cutting depth and the orthogonal turning and milling cutting-in and cutting-out angle, the integral limit of the infinitesimal cutting force is determined, and the infinitesimal cutting force participating in cutting is integrated to obtain the total cutting force in the current cutting state.
9. The modeling method of claim 8, wherein the current cutting state total cutting force algorithm is as follows:
according to the mechanical cutting force model of Altintas, the cutting force of the ith element of the jth cutter tooth can be expressed as:
in the formula, Ktc,KrcAnd KacCoefficient of action of the shearing action on the tangential, radial and axial cutting forces, Kte,KreAnd KaeRespectively the action coefficients of the friction action of the cutting edge to tangential, radial and axial cutting forces,
Figure FDA0000430845470000042
is the chip thickness of the ith infinitesimal of the jth cutter tooth, N is the cutter tooth number,
Figure FDA0000430845470000043
is the instantaneous contact angle of the ith element of the jth cutter tooth;
the infinitesimal tangential force, the infinitesimal radial force and the infinitesimal axial force in the formula (1) are converted into a milling cutter coordinate system xFyFzFThe infinitesimal cutting force component along the x-axis, the infinitesimal cutting force component along the y-axis, and the infinitesimal cutting force component along the z-axis in (b) to obtain:
Figure FDA0000430845470000044
integrating the infinitesimal cutting force along the axial direction of the milling cutter, and obtaining the cutting force according to the equations (1) and (2):
<math> <mrow> <msub> <mi>F</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mo>&Integral;</mo> <mn>0</mn> <msub> <mi>a</mi> <mi>p</mi> </msub> </msubsup> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mo>-</mo> <mi>cos</mi> <msub> <mi>&phi;</mi> <mi>ij</mi> </msub> </mtd> <mtd> <mo>-</mo> <mi>sin</mi> <msub> <mi>&phi;</mi> <mi>ij</mi> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <msub> <mi>&phi;</mi> <mi>ij</mi> </msub> </mtd> <mtd> <mo>-</mo> <mi>cos</mi> <msub> <mi>&phi;</mi> <mi>ij</mi> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>K</mi> <mi>tc</mi> </msub> </mtd> <mtd> <msub> <mi>K</mi> <mi>te</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>K</mi> <mi>rc</mi> </msub> </mtd> <mtd> <msub> <mi>K</mi> <mi>re</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>K</mi> <mi>ac</mi> </msub> </mtd> <mtd> <msub> <mi>K</mi> <mi>ae</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>&phi;</mi> <mi>ij</mi> </msub> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>&phi;</mi> <mi>ij</mi> </msub> <mo>)</mo> </mrow> <mi>dz</mi> <mo>,</mo> <mi>q</mi> <mo>=</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mi>z</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </math>
Figure FDA0000430845470000051
in the formula, apIn order to be the nominal depth of cut,
Figure FDA0000430845470000052
is the cutting-in angle of the milling cutter,
Figure FDA0000430845470000053
is the milling cutter cut-off angle;
the instantaneous cutting forces of all the helical cutting edges are summed to obtain the total instantaneous cutting force acting on the milling cutter:
<math> <mrow> <msub> <mi>F</mi> <mi>x</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>F</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>F</mi> <mi>y</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>F</mi> <mrow> <mi>y</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <msub> <mi>F</mi> <mi>z</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>F</mi> <mrow> <mi>z</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow> </math>
CN201310655052.8A 2013-12-05 2013-12-05 Cutting force modeling method for flat bottom spiral end mill orthogonal turning milling shaft parts Active CN103646141B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310655052.8A CN103646141B (en) 2013-12-05 2013-12-05 Cutting force modeling method for flat bottom spiral end mill orthogonal turning milling shaft parts

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310655052.8A CN103646141B (en) 2013-12-05 2013-12-05 Cutting force modeling method for flat bottom spiral end mill orthogonal turning milling shaft parts

Publications (2)

Publication Number Publication Date
CN103646141A true CN103646141A (en) 2014-03-19
CN103646141B CN103646141B (en) 2017-05-10

Family

ID=50251354

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310655052.8A Active CN103646141B (en) 2013-12-05 2013-12-05 Cutting force modeling method for flat bottom spiral end mill orthogonal turning milling shaft parts

Country Status (1)

Country Link
CN (1) CN103646141B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104239720A (en) * 2014-09-18 2014-12-24 北京理工大学 Complicated surface cutting force prediction method and parameter adjustment and tool path planning method
CN105868455A (en) * 2016-03-25 2016-08-17 哈尔滨理工大学 A layer cutting difference analysis method for left and right cutting edges of a cutter
CN107085412A (en) * 2017-03-30 2017-08-22 北京精密机电控制设备研究所 A kind of method of controlling curve turnery processing cutting force mutation
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
CN107971538A (en) * 2017-09-01 2018-05-01 徐州工程学院 A kind of emulation mode of orthogonal turn-milling finished surface microscopic appearance
CN108161573A (en) * 2017-11-28 2018-06-15 北京动力机械研究所 A kind of thin-walled workpiece bias orthogonal turn-milling machining deformation compensation method
CN108490881A (en) * 2018-04-15 2018-09-04 徐州工程学院 A method of orthogonal turn-milling Chip Shape and type are differentiated using emulation technology
CN108563848A (en) * 2018-03-27 2018-09-21 华中科技大学 A kind of milling force modeling method of flat helical end mills
CN108804749A (en) * 2018-04-19 2018-11-13 西北工业大学 Thin-wall part milling process Damping Modeling method
CN110262397A (en) * 2019-06-24 2019-09-20 北京理工大学 Turn-milling cutting spatially spiral trochoid motion profile and instantaneous Predictive Model of Cutting Force
CN112799299A (en) * 2020-12-24 2021-05-14 华中科技大学 Robot multi-axis milling stability model and construction method thereof

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3531786C2 (en) * 1985-09-06 1987-12-10 Rolf Klenk Hartmetallwerkzeugfabrik Gmbh & Co Kg, 7959 Balzheim, De
CN102059381A (en) * 2010-11-04 2011-05-18 西北工业大学 Machining error prediction method in peripheral milling process
CN102436527A (en) * 2011-11-10 2012-05-02 西北工业大学 Milling force modeling method of milling process of countersinking vertical milling cutter
CN102653013A (en) * 2012-04-16 2012-09-05 华中科技大学 Orthogonal turn-milling wide line processing method based on eccentric distance optimization
CN102681488A (en) * 2012-05-24 2012-09-19 南京航空航天大学 Modeling method for milling surface appearance of workpiece
CN103258095A (en) * 2013-05-14 2013-08-21 西北工业大学 Universal milling force modeling method for flat-bottom end mill

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3531786C2 (en) * 1985-09-06 1987-12-10 Rolf Klenk Hartmetallwerkzeugfabrik Gmbh & Co Kg, 7959 Balzheim, De
CN102059381A (en) * 2010-11-04 2011-05-18 西北工业大学 Machining error prediction method in peripheral milling process
CN102436527A (en) * 2011-11-10 2012-05-02 西北工业大学 Milling force modeling method of milling process of countersinking vertical milling cutter
CN102653013A (en) * 2012-04-16 2012-09-05 华中科技大学 Orthogonal turn-milling wide line processing method based on eccentric distance optimization
CN102681488A (en) * 2012-05-24 2012-09-19 南京航空航天大学 Modeling method for milling surface appearance of workpiece
CN103258095A (en) * 2013-05-14 2013-08-21 西北工业大学 Universal milling force modeling method for flat-bottom end mill

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
P. LEE ET AL: "PREDICTION OF BALL-END MILLING FORCES FROM ORTHOGONAL CUTTING DATA", 《TOOLS MANUFACT》, vol. 36, no. 9, 31 December 1996 (1996-12-31) *
闫兵 等: "球头铣刀铣削过程的物理建模与仿真", 《天津职业技术师范学院学报》, vol. 11, no. 4, 31 December 2001 (2001-12-31) *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104239720B (en) * 2014-09-18 2017-05-10 北京理工大学 Complicated surface cutting force prediction method and parameter adjustment and tool path planning method
CN104239720A (en) * 2014-09-18 2014-12-24 北京理工大学 Complicated surface cutting force prediction method and parameter adjustment and tool path planning method
CN105868455B (en) * 2016-03-25 2019-01-22 哈尔滨理工大学 Cutter cutting to left and right sword layered cutting difference analysis method
CN105868455A (en) * 2016-03-25 2016-08-17 哈尔滨理工大学 A layer cutting difference analysis method for left and right cutting edges of a cutter
CN107085412A (en) * 2017-03-30 2017-08-22 北京精密机电控制设备研究所 A kind of method of controlling curve turnery processing cutting force mutation
CN107168245B (en) * 2017-05-04 2019-08-23 武汉理工大学 A kind of accurate prediction technique of chamfered edge circular bit cutting force considering cutting edge effect
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
CN107971538A (en) * 2017-09-01 2018-05-01 徐州工程学院 A kind of emulation mode of orthogonal turn-milling finished surface microscopic appearance
CN108161573A (en) * 2017-11-28 2018-06-15 北京动力机械研究所 A kind of thin-walled workpiece bias orthogonal turn-milling machining deformation compensation method
CN108161573B (en) * 2017-11-28 2019-08-16 北京动力机械研究所 A kind of thin-walled workpiece bias orthogonal turn-milling machining deformation compensation method
CN108563848A (en) * 2018-03-27 2018-09-21 华中科技大学 A kind of milling force modeling method of flat helical end mills
CN108563848B (en) * 2018-03-27 2020-07-10 华中科技大学 Milling force modeling method of flat-bottom spiral end mill
CN108490881A (en) * 2018-04-15 2018-09-04 徐州工程学院 A method of orthogonal turn-milling Chip Shape and type are differentiated using emulation technology
CN108490881B (en) * 2018-04-15 2020-01-10 徐州工程学院 Method for judging shape and type of orthogonal turn-milling cutting by adopting simulation technology
CN108804749A (en) * 2018-04-19 2018-11-13 西北工业大学 Thin-wall part milling process Damping Modeling method
CN110262397A (en) * 2019-06-24 2019-09-20 北京理工大学 Turn-milling cutting spatially spiral trochoid motion profile and instantaneous Predictive Model of Cutting Force
CN110262397B (en) * 2019-06-24 2020-10-23 北京理工大学 Method for modeling spiral trochoid motion trajectory and instantaneous cutting force in turn-milling machining space
CN112799299A (en) * 2020-12-24 2021-05-14 华中科技大学 Robot multi-axis milling stability model and construction method thereof

Also Published As

Publication number Publication date
CN103646141B (en) 2017-05-10

Similar Documents

Publication Publication Date Title
CN103646141B (en) Cutting force modeling method for flat bottom spiral end mill orthogonal turning milling shaft parts
Wei et al. Cutting force prediction in ball end milling of sculptured surface with Z-level contouring tool path
CN102622489B (en) Five-axis side milling cutting force predicting method based on ACIS platform
CN110262397B (en) Method for modeling spiral trochoid motion trajectory and instantaneous cutting force in turn-milling machining space
Yücesan et al. Improved modelling of cutting force coefficients in peripheral milling
CN105701323A (en) Plunge milling machining cutting force predicting and modeling method
CN102129232A (en) Five-axis side milling machining process parameter design method
Han et al. Precise prediction of forces in milling circular corners
Pleta et al. Identification of optimal machining parameters in trochoidal milling of Inconel 718 for minimal force and tool wear and investigation of corresponding effects on machining affected zone depth
CN104794337B (en) A kind of orthogonal turn-milling judged based on boundary condition processes end face sword On Cutting Force Modeling
CN104462775B (en) A kind of Five-axis NC Machining Center Milling Force computational methods
CN106424969B (en) A kind of slotting accurate prediction technique of milling dynamic cutting force considering cutter deflection
CN103235556A (en) Feature-based numerical-control method for processing and manufacturing complicated parts
CN105414616B (en) Cutting force forecast and Convenient stable criterion during helical milling
CN106934170A (en) Chatter stability lobes flap figure modeling method based on rose cutter Yu workpiece contact zone
Desai et al. Process geometry modeling with cutter runout for milling of curved surfaces
CN106682281A (en) Method for predicting instantaneous cutting force of milling based on maximum cutting force
CN107168245A (en) A kind of accurate Forecasting Methodology of chamfered edge circular bit cutting force for considering cutting edge effect
Wang et al. Identification and analysis of cutting force coefficients in the helical milling process
CN106934172A (en) A kind of multiple-cutting-edge milling clearance computational methods of carbon fibre composite
CN111459096A (en) Method for calculating side milling force of end milling cutter with non-uniform wear arc head
CN105160128A (en) Prediction method for cutting force in curve end milling process
Shen et al. Grinding wheel parametric design for machining arbitrary grooves on the helical rake face of the tool
Chen et al. A study of design and manufacturing models for circular-arc ball-end milling cutters
Dikshit et al. Determining cutting force coefficients from instantaneous cutting forces in ball end milling

Legal Events

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