CN102581364B - Method for calculating milling load of ball-end milling cutter - Google Patents
Method for calculating milling load of ball-end milling cutter Download PDFInfo
- 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
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
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):
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,
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:
In formula, | AC| with | the value of AD| can be passed through similar Δ BAC and Δ BHO:
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 θ:
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:
In formula, θ
eby formula (2), provided.
The O2 that sets up an office is d to horizontal distance
θ, be equal in formula (4) | AC|:
Cutting out angle can be expressed as:
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)
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.
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)
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)
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)
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 |
-
2012
- 2012-03-05 CN CN201210055859.3A patent/CN102581364B/en active Active
Patent Citations (5)
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)
Title |
---|
富宏亚,张翔,韩振宇,孙雅洲.微径球头铣刀铣削力建模与仿真.《计算机集成制造系统》.2011,第17卷(第7期),1448-1453. * |
数控铣削加工物理建模与仿真;苑辉;《哈尔滨理工大学硕士学位论文》;20090331;8-16 * |
苑辉.数控铣削加工物理建模与仿真.《哈尔滨理工大学硕士学位论文》.2009,8-16. |
Cited By (1)
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 |