CN102581364B - Method for calculating milling load of ball-end milling cutter - Google Patents

Method for calculating milling load of ball-end milling cutter Download PDF

Info

Publication number
CN102581364B
CN102581364B CN201210055859.3A CN201210055859A CN102581364B CN 102581364 B CN102581364 B CN 102581364B CN 201210055859 A CN201210055859 A CN 201210055859A CN 102581364 B CN102581364 B CN 102581364B
Authority
CN
China
Prior art keywords
cutter
milling
infinitesimal
cutting
workpiece
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201210055859.3A
Other languages
Chinese (zh)
Other versions
CN102581364A (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.)
HUST Wuxi Research Institute
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 CN201210055859.3A priority Critical patent/CN102581364B/en
Publication of CN102581364A publication Critical patent/CN102581364A/en
Application granted granted Critical
Publication of CN102581364B publication Critical patent/CN102581364B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Numerical Control (AREA)

Abstract

The invention discloses a method for calculating the milling load of a ball-end milling cutter, and the method comprises the following steps of: obtaining a relative position parameter model of the cutter and a workpiece, and determining all cutter micro-units participating in cutting according to the relative position parameter model of the cutter and the workpiece, calculating the cut-in angle, the cut-out angle and the instantaneous cutting thickness of each side-tilting milling micro-unit participating in cutting, calculating the cut-in angle, the cut-out angle and the instantaneous cutting thickness of each forward-tilting milling micro-unit participating in cutting each involved cutting, superimposing the instantaneous cutting thickness of each side-tilting milling micro-unit and the instantaneous cutting thickness of each forward-tilting milling micro-unit participating in cutting so as to obtain the instantaneous cutting thickness of a cutter micro-unit, and summing the instantaneous cutting thickness of all cutter micro-units participating in cutting so as to obtain the instantaneous milling load of the ball-end milling cutter. With the method of the invention, the instantaneous cutting state of a cutting edge and the workpiece and the cutting thickness can be obtained during side-tilting and forward-tilting milling of the cutter, so that predication about a cutting force during five-axis milling processing can be realized.

Description

Calculate the method for milling load of ball-end milling cutter
Technical field
The present invention relates to computer numerical control processing technique field, be specifically related to calculate in a kind of five-shaft numerical control processing the method for milling load of ball-end milling cutter.
Background technology
Rose cutter five axle millings are the complex-curved accurately machined technical ways of aero-engine impeller, blade space-like curved surface and mould class.Due to the complexity of part, thereby the axis direction that conventionally need to change cutter is avoided interfering.For predicting the cutting force of five axle milling process, reduce distortion and vibration in process, first need to calculate the momentary cutting thick of milling cutter cutting workpiece.
Modeling of Milling Forces for Ball-end Milling Cutter is generally divided into some discoid cutter infinitesimals by cutter along axis direction, by calculating cutting state and the infinitesimal cutting force of infinitesimal, then the cutting force of infinitesimal is superposeed and obtains whole Milling Force.Yet existing model is all generally in the situation that the relative position of cutter and workpiece more special (only have and roll or only lean forward) carries out Milling Force prediction, and adopts vector calculation method to obtain cutting load.Yet this method cannot reflect that infinitesimal cuts workpiece and cut out the true process of workpiece in one-period, also cannot calculate the milling load that rose cutter not only rolls but also leans forward.
Summary of the invention
The object of the invention is to the defect existing in existing milling force modeling, a kind of method of calculating milling load of ball-end milling cutter is provided, it can obtain instantaneous cutting state and the thickness of cutting of cutting edge and the workpiece of cutter when cutter rolls and leans forward milling, thereby realizes the prediction of five axle Milling Forces.
The present invention is achieved by the following technical solutions:
A method of calculating milling load of ball-end milling cutter, comprises the steps:
(1) obtain the relative position parameter model of cutter and workpiece;
(2) according to the relative position parameter model of cutter and workpiece, determine the whole cutter infinitesimals that participate in cutting, cutter infinitesimal comprises inclination milling infinitesimal and leans forward milling infinitesimal;
(3) according to the relative position parameter model of cutter and workpiece, calculate each entrance angle that participates in the inclination milling infinitesimal of cutting, cut out angle and momentary cutting thick;
(4) according to the relative position parameter model of cutter and workpiece, calculate each entrance angle that participates in the milling infinitesimal that leans forward of cutting, cut out angle and momentary cutting thick;
(5) momentary cutting thick of the milling infinitesimal that leans forward of the momentary cutting thick of each inclination milling infinitesimal and each participation cutting is superposeed, to obtain the momentary cutting thick of cutter infinitesimal;
(6) momentary cutting thick that participates in all cutter infinitesimals of cutting is sued for peace, to obtain the Instantaneous Milling load of rose cutter.
The present invention has following advantage and technique effect: by rose cutter working angles is decomposed, disclose the incision of cutter infinitesimal and the process that cuts out workpiece of the different axial locations of cutter, and the thickness of the instantaneous cutting of infinitesimal, obtain the cutting point of instantaneous participation cutting, thereby calculate the thickness of cutting in cutter inclination situation, to predict that five axles process instantaneous cutting force.
Accompanying drawing explanation
Fig. 1 is the flow chart that the present invention calculates the method for milling load of ball-end milling cutter
Fig. 2 is the relative tertiary location schematic diagram of cutter and workpiece.
Fig. 3 (a) is that to Fig. 3 (c) cutter sphere participates in the part of cutting and the relation of direction of feed.
Fig. 4 is the decomposing schematic representation of cutter infinitesimal direction of feed.
Fig. 5 is for rolling milling schematic diagram.
Fig. 6 (a) is the milling schematic diagram that leans forward.
Fig. 6 (b) is to the different milling view of Fig. 6 (e) for the milling infinitesimal that leans forward.
Fig. 7 (a) and (b) the milling load result that adopts the inventive method to calculate is shown.
The specific embodiment
Below in conjunction with drawings and Examples, the scheme of accepting of the present invention is described in further detail, but the present embodiment is not limited to the present invention.
The parameter of the rose cutter that the present invention adopts is as follows: 10 millimeters of diameters, and sword number is 4, helical angle 30 degree, surface to be machined is plane, and cutting way adopts groove milling mode, and cutting depth is 1 millimeter, and per tooth feeding is 0.1 millimeter.
As shown in Figure 1, the method that the present invention calculates milling load of ball-end milling cutter comprises the following steps:
1, obtain the relative position parameter model of cutter and workpiece;
As shown in Figure 2, T is tool axis direction, and cutting direction of feed is identical with X-direction, and the incline direction of cutter and angle are shown in right figure, O is the centre of sphere point of cutter, OA represents tool axis, and α represents tool axis and finished surface normal direction angulation, and φ represents the corner of cutter infinitesimal, be that tool axis is at the angle of finished surface projection and direction of feed X, A ' expression A point, in the projection (as shown in Fig. 3 (a) to (c)) of finished surface, when φ=0, is the milling of leaning forward; When φ=pi/2, for rolling milling.
2, according to the relative position parameter model of cutter and workpiece, determine the whole cutter infinitesimals that participate in cutting, cutter infinitesimal comprises inclination milling infinitesimal and the milling infinitesimal that leans forward, and the feed motion of cutter infinitesimal can be decomposed into the feed motion of inclination milling infinitesimal and lean forward the feed motion of milling infinitesimal;
As shown in Figure 4, the feeding of cutter infinitesimal can be divided into inclination milling feed f in the feed motion of cutter infinitesimal cwith the milling feed f that leans forward q, and calculated by following formula (1):
f c = f sin φ f q = f cos φ - - - ( 1 )
The feeding that in formula, f is rose cutter.
As shown in Figure 5, θ represents the position angle of cutter infinitesimal, and establishing the region minimum bit angle setting that participates in cutting is θ s, that maximum is θ e,
θ s = α - arccos | OH | R θ e = α + arccos | OH | R - - - ( 2 )
Work as θ s< θ < θ etime, infinitesimal participates in cutting.The radius that wherein R is rose cutter, | OH| represents that centre of sphere O is to the line segment length of intersection point H.
3, according to the relative position parameter model of cutter and workpiece, calculate each entrance angle that participates in the inclination milling infinitesimal of cutting, cut out angle and momentary cutting thick;
Roll milling cross section as shown in Figure 5, direction of feed is outside perpendicular to paper.For the different infinitesimals that participate in cutting, it is not identical that its incision cuts out angle.Position angle is the cutter infinitesimal at θ place, and its radius of clean-up is | AD|, cuts out angle Φ ex, entrance angle Φ stexpression formula is:
&Phi; st = arccos ( - | AC | | AD | ) &Phi; ex = 2 &pi; - &Phi; st - - - ( 3 )
In formula, | AC| with | the value of AD| can be passed through similar Δ BAC and Δ BHO:
| AC | = | OH | &CenterDot; | AB | | HB | = | OH | &CenterDot; | OB | - | OA | | OH &CenterDot; tan &alpha; | = | OH &CenterDot; sec &alpha; | - | R &CenterDot; cos &theta; | tan &alpha; | AD | = R &CenterDot; sin &theta; - - - ( 4 )
In rolling milling, because tool axis is perpendicular to direction of feed, when the corner of cutter infinitesimal meet time, its thickness of cutting expression formula is:
Due to in variation desirable on the occasion of and negative value, so h cmay be negative value, now represent that this blade is away from workpiece cutting surface.
4, according to the relative position parameter model of cutter and workpiece, calculate each entrance angle that participates in the milling infinitesimal that leans forward of cutting, cut out angle and momentary cutting thick;
As Fig. 6 (a), in the milling of leaning forward, the position angle of establishing the cutter infinitesimal of current participation cutting is θ, and the cutter infinitesimal position angle at this place of front cutting at one time is θ ', perpendicular to the feed engagement f of cutter axis orientation θcan be expressed as:
±f θ=Rsinθ+f qcosα-Rsinθ′(6)
Rcosθ=Rcosθ′+f qsinα (7)
Be f θmeet equation with θ:
( R sin &theta; + f q cos &alpha; + &OverBar; f &theta; ) 2 + ( R cos &theta; - f q sin &alpha; ) 2 = R 2 - - - ( 8 )
F in formula qfor the feed engagement of the milling of leaning forward, in formula (6), in darker regions, get positive sign, in light areas, get negative sign, in formula (8), in contrast.
In the milling mode leaning forward, the instantaneous thickness of cutting of not being out of shape of infinitesimal is as Fig. 6.In Fig. 6 (a), dotted outline is the tooth feeding position of cutter before, the outline position that solid line is current cutting tool.While leaning forward milling, there are four exemplary position, four positions as shown in Fig. 6 (a) 1,2,3,4, Fig. 6 (b), (c), (d), (e) in cutter infinitesimal be A to plane.Position 1, before the corresponding infinitesimal cutting of current cutter, this place's workpiece was not cut by last position knives, and therefore the cutting state of corresponding infinitesimal is as shown in Fig. 6 (b), and hatching is expressed as orthogonal cutting; Position 2, as Fig. 6 (c), dotted line represents to be cut by last position knives, and before the corresponding infinitesimal cutting of current cutter, this place's workpiece was cut by last position knives, but current infinitesimal still can switch to surface of the work, and hatching is expressed as orthogonal cutting; Position 3, the cutting state of corresponding infinitesimal is as shown in Fig. 6 (d), before current corresponding infinitesimal cutting, this place's workpiece was cut by last position knives, but current infinitesimal can not switch to surface of the work, hatching is orthogonal cutting region, and cutter tooth is cut and cut out is negative cutting zone constantly; During position 4, as Fig. 6 (e), current cutter infinitesimal can not switch to workpiece, and cutter tooth is in negative cutting zone.
As shown in Fig. 6 (b), 1 place in position, the angular position theta of cutter infinitesimal meets following formula:
arccos R cos &theta; e + f q sin &alpha; R &le; &theta; < &theta; e
In formula, θ eby formula (2), provided.
The O2 that sets up an office is d to horizontal distance θ, be equal in formula (4) | AC|:
d &theta; = | OH | &CenterDot; sec &alpha; - R cos &theta; tan &alpha;
Cutting out angle can be expressed as:
&Phi; st = &pi; - arccos d &theta; R sin &theta; &Phi; ex = &pi; + arccos d &theta; R sin &theta;
When cutter infinitesimal corner meet time, its thickness of cutting expression formula is:
Because be less than zero, so above formula positive sign.
The incision that in like manner can obtain position 2, position 3 and position 4 cutters cuts out angle and momentary cutting thick.
5, the momentary cutting thick of the milling infinitesimal that leans forward of the momentary cutting thick of each inclination milling infinitesimal and each participation cutting is superposeed, to obtain the momentary cutting thick of cutter infinitesimal;
Make α=45 °, φ=45 °, result of calculation is as Fig. 7 (a), and the curve in figure is the rotate a circle state of diverse location cutting workpiece of a blade of cutter, being less than or equal to null representation does not have cutting workpiece, is greater than the momentary cutting thick of null representation cutting workpiece.Visible in figure, the cutting maximum ga(u)ge of infinitesimal occurs in the instantaneous of incision workpiece, and thickness of cutting reduces in time gradually afterwards, until cut out workpiece.Make α=45 °, φ=10 °, result of calculation is as Fig. 7 (b), and the cutting maximum ga(u)ge of infinitesimal occurs in a certain moment after incision workpiece, can calculate the thickness of cutting of arbitrfary point on any time blade simultaneously.
6, the momentary cutting thick that participates in all cutter infinitesimals of cutting is sued for peace, to obtain the Instantaneous Milling load of rose cutter.

Claims (1)

1. a method of calculating rose cutter five axle milling loads, is characterized in that, comprises the steps:
(1) obtain the relative position parameter model of cutter and workpiece, in this step, the cutter centre of sphere point of rose cutter is defined as to O, tool axis is defined as OA, the projection of tool axis OA on finished surface is defined as OA ', between tool axis and finished surface normal direction, angulation is expressed as α, and the described projection OA ' of tool axis on finished surface and the angle between direction of feed are expressed as φ, build in this way the relative position parameter model between cutter and workpiece;
(2) according to the above-mentioned relative position parameter model of cutter and workpiece, determine the whole cutter infinitesimals that participate in cutting, described cutter infinitesimal comprises inclination milling infinitesimal simultaneously and leans forward milling infinitesimal, thus the feed motion of all cutter infinitesimals is decomposed into separately to the inclination milling infinitesimal feed motion and the milling infinitesimal feed motion of leaning forward that by following formula (), are represented:
(1)
Wherein, f crepresent to roll the feed motion amount of milling infinitesimal, f qrepresent the to lean forward feed motion amount of milling infinitesimal, the corner that represents each cutter infinitesimal, is also formed angle respectively between the projection of its corresponding tool axis on finished surface and direction of feed, and f represents the total feed of rose cutter;
(3) according to the above-mentioned relative position parameter model of cutter and workpiece calculate described each participate in cutting inclination milling infinitesimal entrance angle and cut out angle, based on following formula (two), calculate the momentary cutting thick h that each rolls milling infinitesimal simultaneously c:
(2)
Wherein, f crepresent that each rolls the feed motion amount of milling infinitesimal, its numerical value calculates in step (2);
(4) according to the relative position parameter model of described cutter and workpiece calculate described each participate in cutting the milling infinitesimal that leans forward entrance angle and cut out angle, based on following formula (three) and (four), calculate the momentary cutting thick h of each milling infinitesimal that leans forward simultaneously q:
(3)
d &theta; = | OH | &CenterDot; sec &alpha; - R cos &theta; tan &alpha; (4)
Wherein, R represents the radius of rose cutter, and θ represents the position angle of each cutter infinitesimal, | OH| represents that the cutter centre of sphere of milling cutter is to the line segment length between its intersection point;
(5) momentary cutting thick of the milling infinitesimal that leans forward of the momentary cutting thick of described each inclination milling infinitesimal and described each participation cutting is superposeed, to obtain the momentary cutting thick of described cutter infinitesimal;
(6) momentary cutting thick that participates in all cutter infinitesimals of cutting is sued for peace, to obtain the Instantaneous Milling load of rose cutter.
CN201210055859.3A 2012-03-05 2012-03-05 Method for calculating milling load of ball-end milling cutter Active CN102581364B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210055859.3A CN102581364B (en) 2012-03-05 2012-03-05 Method for calculating milling load of ball-end milling cutter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210055859.3A CN102581364B (en) 2012-03-05 2012-03-05 Method for calculating milling load of ball-end milling cutter

Publications (2)

Publication Number Publication Date
CN102581364A CN102581364A (en) 2012-07-18
CN102581364B true CN102581364B (en) 2014-08-27

Family

ID=46470904

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210055859.3A Active CN102581364B (en) 2012-03-05 2012-03-05 Method for calculating milling load of ball-end milling cutter

Country Status (1)

Country Link
CN (1) CN102581364B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI562870B (en) * 2014-12-30 2016-12-21 Tongtai Machine & Tool Co Ltd Method for determining processing parameters of ball end mill in cutting process

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462775B (en) * 2014-11-12 2018-07-20 天津大学 A kind of Five-axis NC Machining Center Milling Force computational methods
CN104794337B (en) * 2015-04-17 2017-10-20 北京航空航天大学 A kind of orthogonal turn-milling judged based on boundary condition processes end face sword On Cutting Force Modeling
CN109656192B (en) * 2018-12-18 2020-02-21 华中科技大学 Machining optimization method based on milling force of ball-end milling cutter
CN111061218B (en) * 2019-12-31 2021-07-27 华中科技大学 Complex curved surface machining cutting force prediction method and device based on ACIS

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5383131A (en) * 1992-07-14 1995-01-17 Mazda Motor Corporation Method of determining cutting trajectory in N/C machining system
US5827020A (en) * 1995-07-11 1998-10-27 Toshiba Kikai Kabushiki Kaisha Numerical control apparatus
US6824336B2 (en) * 2000-12-28 2004-11-30 Mazda Motor Corporation Method for controlling cutting machine
US6942436B2 (en) * 2002-03-11 2005-09-13 Yoshiaki Kakino NC program generating method, NC apparatus, computer memory product, and computer program product
CN101269423A (en) * 2008-05-14 2008-09-24 北京航空航天大学 Multi-coordinate end milling process tool position optimization method using longitude line division tool bit

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5383131A (en) * 1992-07-14 1995-01-17 Mazda Motor Corporation Method of determining cutting trajectory in N/C machining system
US5827020A (en) * 1995-07-11 1998-10-27 Toshiba Kikai Kabushiki Kaisha Numerical control apparatus
US6824336B2 (en) * 2000-12-28 2004-11-30 Mazda Motor Corporation Method for controlling cutting machine
US6942436B2 (en) * 2002-03-11 2005-09-13 Yoshiaki Kakino NC program generating method, NC apparatus, computer memory product, and computer program product
CN101269423A (en) * 2008-05-14 2008-09-24 北京航空航天大学 Multi-coordinate end milling process tool position optimization method using longitude line division tool bit

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
富宏亚,张翔,韩振宇,孙雅洲.微径球头铣刀铣削力建模与仿真.《计算机集成制造系统》.2011,第17卷(第7期),1448-1453. *
数控铣削加工物理建模与仿真;苑辉;《哈尔滨理工大学硕士学位论文》;20090331;8-16 *
苑辉.数控铣削加工物理建模与仿真.《哈尔滨理工大学硕士学位论文》.2009,8-16.

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI562870B (en) * 2014-12-30 2016-12-21 Tongtai Machine & Tool Co Ltd Method for determining processing parameters of ball end mill in cutting process

Also Published As

Publication number Publication date
CN102581364A (en) 2012-07-18

Similar Documents

Publication Publication Date Title
CN102581364B (en) Method for calculating milling load of ball-end milling cutter
CN102004818A (en) Instantaneous milling force modeling method in free curve shape part peripheral milling process
EP1868054B1 (en) A method for milling grooves in bored workpieces
CN102637216B (en) Method for generating numerical-control side milling machining tool path for complicated curved surfaces
CN102681488B (en) Modeling method for milling surface appearance of workpiece
CN104102171B (en) A kind of ball head knife Multi-axis Machining optimal tool orientation method
CN110355614B (en) Grinding wheel track solving method for integral end mill cutter body gap grinding process
US20140277686A1 (en) Method of controlling tool orientation and step-over distance in face milling of curvilinear surfaces
CN103042262B (en) Integrated finish machining method for inner profile and corner of pocket
CN113927378B (en) Numerical control grinding track calculation method for peripheral tooth chip dividing groove of rough milling cutter
CN105701323A (en) Plunge milling machining cutting force predicting and modeling method
CN102825308A (en) Aircraft structural part dwell angle back gouging method
Lauwers et al. Development of a five-axis milling tool path generation algorithm based on faceted models
CN104493305A (en) Car gear machining method based on regulation of gear cutting machine
CN107544433B (en) Numerical-controlled machine tool machining process rose cutter and workpiece contact zone semi analytic modeling
CN105160128A (en) Prediction method for cutting force in curve end milling process
CN105171077B (en) A kind of composite bevel milling method of high-precision control surface roughness
Gdula Adaptive method of 5-axis milling of sculptured surfaces elements with a curved line contour
CN105911957A (en) Control method for tooth surface scallop height of gear processed by adopting disc-like milling cutter generating method
CN103042434A (en) Ultraprecise milling surface topography texture control method
RU2429949C1 (en) Procedure for processing mono-wheel
CN106599389B (en) A kind of blade contact method of discrimination of square end mill wide-line processing
CN105759717B (en) A kind of method that anti-mistake of cutter path for five-shaft numerical control processing is cut
Kuo et al. The analysis of NC machining efficiency for marine propellers
CN106077775B (en) A kind of processing method of three axis machine tooling back-off class workpiece

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200615

Address after: No.329, Yanxin Road, Huishan Economic Development Zone, Huishan District, Wuxi City, Jiangsu Province

Patentee after: HUST-WUXI Research Institute

Address before: 430074 Hubei Province, Wuhan city Hongshan District Luoyu Road No. 1037

Patentee before: HUAZHONG University OF SCIENCE AND TECHNOLOGY