CN102681488B - Modeling method for milling surface appearance of workpiece - Google Patents

Modeling method for milling surface appearance of workpiece Download PDF

Info

Publication number
CN102681488B
CN102681488B CN201210162358.5A CN201210162358A CN102681488B CN 102681488 B CN102681488 B CN 102681488B CN 201210162358 A CN201210162358 A CN 201210162358A CN 102681488 B CN102681488 B CN 102681488B
Authority
CN
China
Prior art keywords
cutting
cutter
point
angle
trajectory
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.)
Expired - Fee Related
Application number
CN201210162358.5A
Other languages
Chinese (zh)
Other versions
CN102681488A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201210162358.5A priority Critical patent/CN102681488B/en
Publication of CN102681488A publication Critical patent/CN102681488A/en
Application granted granted Critical
Publication of CN102681488B publication Critical patent/CN102681488B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Numerical Control (AREA)
  • Milling Processes (AREA)

Abstract

The invention discloses a modeling method for milling the surface appearance of a workpiece, belonging to the field of digital milling. The modeling method comprises the steps of carrying out dispersing treatment on the cutting edge of a ball-end mill according to cutting movement locus of a milling cutter; establishing infinitesimal cutting locus equation of the cutting edge of the ball-end mill, wherein the surface appearance of the cutting locus equation is formed by the outermost side locus in the established locus equation; determining corresponding cutter turning angle range for milling the outermost side locus by judging the cutting infinitesimal position angle range corresponding to the outermost side locus; calculating corresponding milling time for milling the locus by combining with angular rotation speed obtained through a cutter spindle rotation speed; and obtaining the outermost side milling locus through a helical lag angle at the cutting infinitesimal position, thus obtaining the surface appearance of the milling workpiece. The modeling method can solve the generation problem on the surface appearance of the workpiece in milling process.

Description

A kind of modeling method of Milling Process workpiece surface appearance
Technical field
The invention belongs to numerical control milling manufacture field, the workpiece surface appearance modeling field that particularly in numerical control milling processing, tool wear causes.
Background technology
Along with the development of modern mechanical manufacturing industry, more and more higher to the requirement of part processing precision.And in actual Milling Processes, owing to being subject to the impact of various factors inevitably can produce mismachining tolerance.Mismachining tolerance is to affect the very important factor of workpiece crudy, can make Workpiece Machining Accuracy obviously reduce, and excessive mismachining tolerance even can cause part to scrap, and has a strong impact on working (machining) efficiency and benefit.The microscopic appearance of surface of the work and workpiece surface roughness have close contact, the wearing quality of workpiece and assembly precision are had to material impact, also be the important indicator of reflection surface of the work crudy simultaneously, by the pattern of prediction surface of the work, can obtain the coordinate figure of any location point on workpiece, the theoretical coordinate value of itself and this point is contrasted, can obtain the mismachining tolerance value at this some place.
In order to improve the crudy of surface of the work, reduce mismachining tolerance, reduce roughness, domestic and international many scholars are studied workpiece surface appearance from microcosmic angle, have obtained some achievements.Successively propose some and set up the modeling method of workpiece surface appearance model.
The research that the factors such as the cutter distortion mostly modeling method of existing workpiece surface appearance is to cause with regard to cutting parameter selection, cutter positioning error, cutting force are carried out the impact of workpiece surface appearance, also seldom has the research of the impact of workpiece surface appearance being correlated with for the cutter Dynamic wear in rose cutter milling process at present.
For this reason, the present invention has overcome the deficiency of the modeling method existence of above-mentioned workpiece surface appearance, according to the feature of rose cutter cutting edge, sets up Cutting trajectory equation, and further tool wear is taken into account the impact of Cutting trajectory, the workpiece surface appearance after Cutter wear is studied.Analyze the change of shape of bulb cutting edge outline sword line after tool wear, can further analyze accordingly the Cutting trajectory of the relative workpiece of cutting edge after wearing and tearing, set up the Cutting trajectory equation of considering wearing and tearing, and then the workpiece surface appearance after being processed.
Summary of the invention
Object of the present invention, is to provide a kind of modeling method of Milling Process workpiece surface appearance, and it can solve modeling and the visual simulating problem of the workpiece surface appearance that in Milling Processes, tool wear causes.
In order to reach above-mentioned purpose, solution of the present invention is:
A modeling method for Milling Process workpiece surface appearance, is characterized in that comprising the steps:
(1) according to Milling Process Tool in Cutting movement locus, dividing cutting edge discrete ball head is a series of cutting infinitesimals, because difference cutting infinitesimal linear velocity when cutting is different, and in process, there is straight-line feed motion and self rotatablely moving around cutter spindle simultaneously, therefore, the movement locus that removes other point on bulb point of a knife ectosphere head cutting edge will form a series of trochoid, the pattern forming on workpiece in order to analyze Tool in Cutting track, the cutting infinitesimal P point of usining after discrete is as research object, motion by analysis site P in feed process, obtain putting the Cutting trajectory that P finally stays on workpiece after process finishing, remove the track that the cut trajectory of falling of work in-process obtains finally forming workpiece pattern, consider straight-line trajectory and rotational trajectory, when cutter weares and teares, the desired cut equation of locus that P is ordered as the formula (1):
Figure GDA0000389611900000021
Wherein, X ', Y ' for P point through cut the coordinate after the time, X 0, Y 0for P point origin coordinates, the tool feeding speed that v is the unit interval, t 1for cut time, t 2for the impact on this process time of the spiral drag angle of P point present position,
Figure GDA0000389611900000022
t=t 1-t 2, for cutter spiral drag angle corresponding to P point,
Figure GDA0000389611900000024
for the residing position angle of P point the cutter radial radius at place,
Figure GDA0000389611900000026
the rose cutter radius of R when not wearing and tearing, ω is cutter spindle angular velocity of rotation,
Figure GDA0000389611900000027
for the line of P point and bulb center O and the angle of cutter spindle Z, the angle that is description cutting infinitesimal position is called position angle;
When the impact of tool wear is taken into account, along with cutter constantly weares and teares, the corresponding radially radius of point on cutting edge on sustained height and corresponding position angle
Figure GDA0000389611900000029
all will change, therefore, during tool abrasion in considering process, the Cutting trajectory that P is ordered can be represented by formula (2):
Figure GDA00003896119000000210
Wherein, X, Y for consider tool wear P point through cut the coordinate after the time, X 0, Y 0for P point origin coordinates, the tool feeding speed that v is the unit interval, t 1for cut time, t 2for the impact on this process time of the spiral drag angle of P point present position, and
Figure GDA00003896119000000211
for cutter spiral drag angle corresponding to P point, ω is cutter spindle angular velocity of rotation;
for the bulb of P point At The Height after tool wear radius radially, by formula (3), solved;
Figure GDA0000389611900000031
In formula, R is the radius of rose cutter while not wearing and tearing,
Figure GDA0000389611900000032
for the cutting plane height at cutting infinitesimal place, VB is
Figure GDA0000389611900000033
the tool flank wear of At The Height Tool in Cutting sword;
(2) the milling track obtaining according to milling equation of locus is not all to form workpiece surface appearance, part trajectory work in-process can be cut, only have outermost trajectory to form workpiece surface appearance, in order to determine outermost trajectory, need to solve the position of Tool in Cutting sword corresponding to outermost trajectory
Figure GDA0000389611900000034
angle span
Figure GDA0000389611900000035
obtain in interval each
Figure GDA0000389611900000037
the Tool in Cutting trajectory that angle is corresponding, the set of these Cutting trajectory lines has just formed the workpiece surface appearance after processing;
Cutter is along the feeding of x axle, and the trajectory that y axle both sides are finally left to surface of the work is respectively P 1p 1' and P 2p 2', wanting to guarantee the existence of outermost track, the coordinate figure that cutting edge track must meet trajectory and feed shaft intersection point is greater than the straight-line feed amount f that cutter rotates a circle p, obtain thus formula (5):
Figure GDA0000389611900000038
Above formula obtains
Figure GDA0000389611900000039
for the minimum value satisfying condition;
The maximal value of position angle is the position angle of the cutting point that the residual high highest point of workpiece surface appearance is corresponding, namely section angle of cut of twice feed in vertical feed direction axial cross section before and after cutter;
For direction of feed right side P 1p 1' this section of track,
Figure GDA00003896119000000310
size to discuss in two kinds of situation;
1) feed for the first time, can be directly according to cutting depth a pcalculate formula (6) with the geometric relationship of Z value maximum point in cutter under this cutting depth and workpiece contact point:
Figure GDA00003896119000000312
2) from feed for the second time,
Figure GDA00003896119000000313
with feed line-spacing f psize have direct relation, can draw shown in the relational expression (7) between them:
Figure GDA00003896119000000314
P for direction of feed left side 2p 2' track,
Figure GDA00003896119000000315
size and P 1p 1' there is relative relation, for last feed,
Figure GDA00003896119000000316
size is shown in formula (8):
Figure GDA00003896119000000317
Feed is for the first time in the process of second from the bottom feed, and the outmost trajectory in left side is corresponding
Figure GDA0000389611900000047
value be all shown in formula (9):
Figure GDA0000389611900000041
(3) in order to obtain each position
Figure GDA0000389611900000048
the Cutting trajectory line that angle is corresponding, need to obtain respectively outermost trajectory P 1p 1' section and P 2p 2the tool motion time in ' section is interval, by obtaining the corner that P is ordered in this is interval, then try to achieve corresponding process time in conjunction with the speed of mainshaft of cutter, consider the spiral drag angle at place, cutting infinitesimal position, can obtain outermost pattern track, thereby obtain the surface topography of workpiece;
In order to set up P 1and P 2corresponding corner
Figure GDA0000389611900000049
with
Figure GDA00003896119000000410
formula, analyze P 1and P 2the corner of point in milling process;
For P 1point, from the negative y axle of starting point, x=0 place forwards P to 1time, P 1point horizontal ordinate meets formula (10) below:
x P 1 = 1 2 f p - - - ( 10 )
That is:
Figure GDA0000389611900000043
Wherein, v is the amount of feeding of bulb cutter unit interval, and ω is cutter shaft rotational speed, and R is tool radius,
Figure GDA00003896119000000411
for P 1the position angle of point,
Figure GDA00003896119000000412
for position angle
Figure GDA00003896119000000413
locate corresponding cutter corner, and meet
Figure GDA00003896119000000414
For P 2point, from the negative y axle of starting point, x=0 place forwards P to 2time, P 2point horizontal ordinate meets formula (11) below:
x p 2 = f p - - - ( 11 )
That is:
Wherein, v is the amount of feeding of bulb cutter unit interval, and ω is cutter shaft rotational speed, and R is tool radius,
Figure GDA00003896119000000415
for P 2the position angle of point,
Figure GDA00003896119000000416
for position angle
Figure GDA00003896119000000417
locate corresponding cutter corner, and meet
Figure GDA00003896119000000418
Use Newton iteration method to solve nonlinear equation (10) and (11), obtain correspondence position the cutter corner value at angle, in conjunction with the angular velocity of rotation of cutter, can be in the hope of the cutting time t in this segment of cutting 1=φ/ω, φ is the cutter anglec of rotation; ;
Because the spiral cutting sword of rose cutter will cause spiral hysteresis phenomenon, the time that will participate in cutting to place, diverse location angle cutting point exerts an influence, and then affects Cutting trajectory, each position angle corresponding cutter spiral drag angle can be expressed as formula (12):
Wherein, β 0for the maximum spiral drag angle of cutter,
Figure GDA0000389611900000055
for the position angle of cutting point on Tool in Cutting sword, the time effects that spiral lags behind and produces
Figure GDA0000389611900000056
By the position angle of trying to achieve
Figure GDA0000389611900000057
scope, cutter tried to achieve
Figure GDA0000389611900000058
interval in corner φ, time that cutter turns over this corner be updated in the Cutting trajectory equation (2) of considering tool wear, can obtain the workpiece surface appearance of rose cutter milling, be expressed as formula (13);
Figure GDA0000389611900000052
Adopt after such scheme, the present invention carries out discrete processes according to Milling Process Tool in Cutting movement locus by rose cutter cutting edge, set up rose cutter Tool in Cutting sword infinitesimal milling equation of locus, Cutting trajectory equation surface topography outermost track in the equation of locus of above-mentioned foundation forms, by judging the positional angle range of the corresponding cutting of these outermost tracks infinitesimal, determine the cutter angle range that these outermost tracks of milling are corresponding, in conjunction with process time corresponding to angular velocity of rotation calculating these tracks of milling being obtained by speed of cutter spindle, consider the spiral drag angle at place, cutting infinitesimal position, obtain these outermost milling tracks, thereby obtain milling workpiece surface appearance, solve the workpiece surface appearance problem being caused by tool wear in Milling Processes.
Accompanying drawing explanation
Fig. 1 is the discrete schematic diagram of bulb cutting edge in the present invention;
Fig. 2 is rose cutter cutting edge Cutting trajectory in the present invention;
Fig. 3 is the tracing point that forms surface topography in the present invention;
Fig. 4 cuts infinitesimal position view before and after tool wear in the present invention;
Fig. 5 is rose cutter monodentate Cutting trajectory schematic diagram in the present invention;
Fig. 6 is that in the present invention, adjacent twice feed cuts angle of cut schematic diagram;
Fig. 7 is the cutter corner schematic diagram in tool cutting process of the present invention;
Number in the figure title: 1-the first cutting edge; 2-the second cutting edge, 3-Cutting trajectory, the point that 4-remains, the point that 5-cuts away, 6-workpiece is surface finally, the point that 7-remains.
Embodiment
Below with reference to accompanying drawing, technical scheme of the present invention is elaborated.
A modeling method for Milling Process workpiece surface appearance, comprises the steps:
(1) rose cutter cutting edge is discrete, sets up rose cutter Tool in Cutting sword infinitesimal milling equation of locus.
According to Milling Process Tool in Cutting movement locus, dividing cutting edge discrete ball head is a series of cutting infinitesimals, be illustrated in figure 1 the discrete schematic diagram of bulb cutting edge in the present invention, in figure, dotted line is depicted as a spiral cutting sword in bulb part, P is a cutting infinitesimal on it, O is ball head branch center
Figure GDA0000389611900000062
for the line of P point and bulb center O and the angle of cutter spindle Z, the angle that is description cutting infinitesimal position is called position angle.Because difference cutting infinitesimal linear velocity when cutting is different, and in process, there is straight-line feed motion and self rotatablely moving around cutter spindle simultaneously, therefore, except the movement locus of other point on bulb point of a knife ectosphere head cutting edge is actual, it is a series of trochoid, the Cutting trajectory that trochoid forms as shown in Figure 2, solid line represents certain some Cutting trajectory on rose cutter the first cutting edge, dotted line represents the Cutting trajectory of sustained height point on rose cutter the second cutting edge, the outermost tracing point in track above and below has formed workpiece pattern, by the tracing point in this region amplify show after as shown in Figure 3.In Fig. 3, the point that hollow dots representative is cut in cutting process, solid dot is the point remaining, and namely forms the point of workpiece surface appearance.
The pattern forming on workpiece in order to analyze Tool in Cutting track, the cutting infinitesimal P point of usining after discrete is as research object, motion by analysis site P in feed process, obtain putting the Cutting trajectory that P finally stays on workpiece after process finishing, remove the track that the cut trajectory of falling of work in-process obtains finally forming workpiece pattern.On cutter, P space of points movement locus consists of rectilinear motion and rotatablely moving, therefore when cutter weares and teares, the desired cut equation of locus that P is ordered as the formula (1):
Wherein, X ', Y ' for P point through cut the coordinate after the time, X 0, Y 0for P point origin coordinates, the tool feeding speed that v is the unit interval, t 1for cut time, t 2for the impact on this process time of the spiral drag angle of P point present position,
Figure GDA0000389611900000063
t=t 1-t 2,
Figure GDA0000389611900000064
for cutter spiral drag angle corresponding to P point,
Figure GDA0000389611900000065
for the residing position angle of P point
Figure GDA0000389611900000066
the cutter radial radius at place,
Figure GDA0000389611900000067
the rose cutter radius of R when not wearing and tearing, ω is cutter spindle angular velocity of rotation.
When the impact of tool wear is taken into account, along with cutter constantly weares and teares, the corresponding radially radius of point on cutting edge on sustained height and corresponding position angle
Figure GDA0000389611900000075
all will change, therefore, during tool abrasion in considering process, the Cutting trajectory that P is ordered can be represented by formula (2):
Figure GDA0000389611900000071
Wherein, X, Y for consider tool wear P point through cut the coordinate after the time, X 0, Y 0for P point origin coordinates, the tool feeding speed that v is the unit interval, t 1for cut time, t 2for the impact on this process time of the spiral drag angle of P point present position, and
Figure GDA0000389611900000077
for cutter spiral drag angle corresponding to P point,
Figure GDA0000389611900000078
for the bulb of P point At The Height after tool wear radius radially, ω is cutter spindle angular velocity of rotation.
Fig. 4 has shown position angle
Figure GDA0000389611900000079
the variation of the cutting infinitesimal P that chooses of place present position before and after cutting edge wearing and tearing, while there is not tool wear, P point is positioned on the profile of bulb outermost, and at cutter, after processing after a while, ball head divides and weares and teares, identical with P point place height
Figure GDA00003896119000000710
section, the physical location after the wearing and tearing of P point becomes P ', the part of PP ' for being worn in process.Therefore, work as position angle
Figure GDA00003896119000000711
while getting different numerical value, just can solve the position coordinates of the cutting infinitesimal after a series of cutting edge wearing and tearing, thereby set up the geometric relationship model of the rear cutter outline sword line of wearing and tearing.
In Fig. 4, while not there are not wearing and tearing, P point place level
Figure GDA00003896119000000712
locate radially radius of corresponding rose cutter can be expressed as:
Figure GDA0000389611900000072
In formula, R 0radius while not wearing and tearing for rose cutter,
Figure GDA00003896119000000714
position angle corresponding to P point while not wearing and tearing for cutter.
After processing a period of time, cutter bulb partly will wear and tear, rose cutter sustained height
Figure GDA00003896119000000715
having there is variation in the actual Working position of the cutting point on place's cutting edge, is changed to P ' point, as shown in Figure 4 by P point.Now, rose cutter cross section diameter corresponding to this plane highly of living in to radius by theoretical value
Figure GDA00003896119000000716
be changed to
Figure GDA00003896119000000717
therefore, after wearing and tearing
Figure GDA00003896119000000718
the radially radius of the cutter that At The Height is corresponding can be expressed as:
Figure GDA0000389611900000073
In formula, R is the radius of rose cutter while not wearing and tearing,
Figure GDA0000389611900000082
for the cutting plane height at cutting infinitesimal place, VB is
Figure GDA0000389611900000083
the tool flank wear of At The Height Tool in Cutting sword.
(2) solve the position of Tool in Cutting sword corresponding to outermost trajectory
Figure GDA0000389611900000084
angle span
Figure GDA0000389611900000085
The milling track obtaining according to milling equation of locus is not all to form workpiece surface appearance, part trajectory work in-process can be cut, only have outermost trajectory to form workpiece surface appearance, in order to determine outermost trajectory, need to solve the position of Tool in Cutting sword corresponding to outermost trajectory
Figure GDA0000389611900000086
angle span
Figure GDA0000389611900000087
Figure GDA0000389611900000088
obtain
Figure GDA0000389611900000089
in interval each
Figure GDA00003896119000000810
the Tool in Cutting trajectory that angle is corresponding, the set of these Cutting trajectory lines has just formed the workpiece surface appearance after processing;
Analysis chart 5, cutter is along the feeding of x axle, and the trajectory that y axle both sides are finally left to surface of the work is respectively P 1p 1' and P 2p 2' these outermost tracks.In order to try to achieve this part final track that forms workpiece surface appearance, remove the cut track of all the other work in-processes, need to obtain respectively Tool in Milling P 1p 1' and P 2p 2the cutting-edge angle that ' time is corresponding
Figure GDA00003896119000000811
scope, illustrate below and how to try to achieve P 1p 1' and P 2p 2' corresponding
Figure GDA00003896119000000812
angle value.
Want to guarantee P 2p 2' exist, the track of cutting edge is this situation of prolate cycloid, must meet Fig. 5 mid point P 3x coordinate figure be greater than the straight-line feed amount f that cutter rotates a circle p.Shown in Fig. 5, solid line represents that position angle is
Figure GDA00003896119000000813
the Cutting trajectory of certain point, the Cutting trajectory of this point finally can form workpiece surface appearance; Dotted line represents that position comparatively
Figure GDA00003896119000000814
the Cutting trajectory of certain point, the Cutting trajectory of this point finally can not form workpiece surface appearance.While meeting following the relation, Cutting trajectory finally can form workpiece surface appearance, and namely corner is that 90 ° of corresponding x coordinates are greater than the x coordinate that corner is (360+270) °, that is:
Figure GDA0000389611900000081
Formula (5) obtains
Figure GDA00003896119000000815
for the minimum value satisfying condition.
The maximal value of position angle is the position angle of the cutting point that the residual high highest point of workpiece surface appearance is corresponding, namely section angle of cut of twice feed in vertical feed direction axial cross section before and after cutter.Direction of feed in Fig. 6 (X-direction) is pointed in paper by O point.
Can find out as shown in Figure 6: for direction of feed right side P 1p 1' this section of track,
Figure GDA00003896119000000816
size to discuss in two kinds of situation.
1) feed for the first time, can be directly according to cutting depth a pcalculate formula (6) with the geometric relationship of Z value maximum point in cutter under this cutting depth and workpiece contact point:
Figure GDA0000389611900000091
2) from feed for the second time, with feed line-spacing f psize have direct relation, can draw shown in the relational expression (7) between them:
Figure GDA0000389611900000092
P for direction of feed left side 2p 2' track, size and P 1p 1' there is relative relation, for last feed,
Figure GDA00003896119000000918
size is shown in formula (8):
Figure GDA0000389611900000093
Feed is for the first time in the process of second from the bottom feed, and the outmost trajectory in left side is corresponding
Figure GDA00003896119000000910
value be all shown in formula (9):
Figure GDA0000389611900000094
(3) outermost trajectory P 1p 1' section and P 2p 2the calculating of tool motion time in ' section is interval
In order to obtain each position
Figure GDA00003896119000000911
the Cutting trajectory line that angle is corresponding, need to obtain respectively outermost trajectory P 1p 1' section and P 2p 2the tool motion time in ' section is interval, by obtaining the corner that P is ordered in this is interval, then try to achieve corresponding process time in conjunction with the speed of mainshaft of cutter, consider the spiral drag angle at place, cutting infinitesimal position, can obtain outermost pattern track, thereby obtain the surface topography of workpiece.
In order to set up P in Fig. 7 1and P 2corresponding corner
Figure GDA00003896119000000912
with
Figure GDA00003896119000000913
formula, analyze P 1and P 2the corner of point in milling process.
For P 1point, from the negative y axle of starting point, x=0 place forwards P to 1time, P 1point horizontal ordinate meets formula (10) below:
x P 1 = 1 2 f p - - - ( 10 )
That is:
Wherein, v is the amount of feeding of bulb cutter unit interval, and ω is cutter shaft rotational speed, and R is tool radius,
Figure GDA00003896119000000914
for P 1the position angle of point,
Figure GDA00003896119000000915
for position angle
Figure GDA00003896119000000916
locate corresponding cutter corner, and meet
Figure GDA00003896119000000917
For P 2point, from the negative y axle of starting point, x=0 place forwards P to 2time, P 2point horizontal ordinate meets formula (11) below:
x p 2 = f p - - - ( 11 )
That is:
Figure GDA0000389611900000102
Wherein, v is the amount of feeding of bulb cutter unit interval, and ω is cutter shaft rotational speed, and R is tool radius,
Figure GDA0000389611900000105
for P 2the position angle of point, for position angle
Figure GDA0000389611900000107
locate corresponding cutter corner, and meet
Figure GDA0000389611900000108
Use Newton iteration method to solve nonlinear equation (10) and (11), obtain correspondence position
Figure GDA0000389611900000109
the cutter corner value at angle, in conjunction with the angular velocity of rotation of cutter, can be in the hope of the cutting time t in this segment of cutting 1=φ/ω, φ is the cutter anglec of rotation; .
Because the spiral cutting sword of rose cutter will cause spiral hysteresis phenomenon, the time that will participate in cutting to place, diverse location angle cutting point exerts an influence, and then affects Cutting trajectory, each position angle
Figure GDA00003896119000001010
corresponding cutter spiral drag angle
Figure GDA00003896119000001011
can be expressed as formula (12):
Wherein, β 0for the maximum spiral drag angle of cutter,
Figure GDA00003896119000001012
for the position angle of cutting point on Tool in Cutting sword, the time effects that spiral lags behind and produces
Figure GDA00003896119000001013
By the position angle of trying to achieve scope, cutter tried to achieve
Figure GDA00003896119000001015
interval in corner φ, time that cutter turns over this corner be updated in the Cutting trajectory equation (2) of considering tool wear, can obtain the workpiece surface appearance of rose cutter milling, be expressed as formula (13);
Figure GDA0000389611900000104
Above example only, for explanation technological thought of the present invention, can not limit protection scope of the present invention with this, every technological thought proposing according to the present invention, and any change of doing on technical scheme basis, within all falling into protection domain of the present invention.

Claims (1)

1. a modeling method for Milling Process workpiece surface appearance, is characterized in that comprising the steps:
(1) according to Milling Process Tool in Cutting movement locus, dividing cutting edge discrete ball head is a series of cutting infinitesimals, because difference cutting infinitesimal linear velocity when cutting is different, and in process, there is straight-line feed motion and self rotatablely moving around cutter spindle simultaneously, therefore, the movement locus that removes other point on bulb point of a knife ectosphere head cutting edge will form a series of trochoid, the pattern forming on workpiece in order to analyze Tool in Cutting track, the cutting infinitesimal P point of usining after discrete is as research object, motion by analysis site P in feed process, obtain putting the Cutting trajectory that P finally stays on workpiece after process finishing, remove the track that the cut trajectory of falling of work in-process obtains finally forming workpiece pattern, consider straight-line trajectory and rotational trajectory, when cutter weares and teares, the desired cut equation of locus that P is ordered as the formula (1):
Figure FDA0000389611890000011
Wherein, X ', Y ' for P point through cut the coordinate after the time, X 0, Y 0for P point origin coordinates, the tool feeding speed that v is the unit interval, t 1for cut time, t 2for the impact on this process time of the spiral drag angle of P point present position,
Figure FDA0000389611890000012
t=t 1-t 2,
Figure FDA0000389611890000013
for cutter spiral drag angle corresponding to P point,
Figure FDA0000389611890000014
for the residing position angle of P point
Figure FDA0000389611890000015
the cutter radial radius at place,
Figure FDA0000389611890000016
the rose cutter radius of R when not wearing and tearing, ω is cutter spindle angular velocity of rotation,
Figure FDA0000389611890000017
for the line of P point and bulb center O and the angle of cutter spindle Z, the angle that is description cutting infinitesimal position is called position angle;
When the impact of tool wear is taken into account, along with cutter constantly weares and teares, the corresponding radially radius of point on cutting edge on sustained height
Figure FDA0000389611890000018
and corresponding position angle
Figure FDA0000389611890000019
all will change, therefore, during tool abrasion in considering process, the Cutting trajectory that P is ordered can be represented by formula (2):
Figure FDA00003896118900000110
Wherein, X, Y for consider tool wear P point through cut the coordinate after the time, X 0, Y 0for P point origin coordinates, the tool feeding speed that v is the unit interval, t 1for cut time, t 2for the impact on this process time of the spiral drag angle of P point present position, and
Figure FDA00003896118900000111
for cutter spiral drag angle corresponding to P point, ω is cutter spindle angular velocity of rotation;
Figure FDA0000389611890000021
for the bulb of P point At The Height after tool wear radius radially, by formula (3), solved;
In formula, R is the radius of rose cutter while not wearing and tearing,
Figure FDA0000389611890000023
for the cutting plane height at cutting infinitesimal place, VB is
Figure FDA0000389611890000024
the tool flank wear of At The Height Tool in Cutting sword;
(2) the milling track obtaining according to milling equation of locus is not all to form workpiece surface appearance, part trajectory work in-process can be cut, only have outermost trajectory to form workpiece surface appearance, in order to determine outermost trajectory, need to solve the position of Tool in Cutting sword corresponding to outermost trajectory
Figure FDA0000389611890000025
angle span
Figure FDA0000389611890000026
obtain
Figure FDA0000389611890000027
in interval each
Figure FDA0000389611890000028
the Tool in Cutting trajectory that angle is corresponding, the set of these Cutting trajectory lines has just formed the workpiece surface appearance after processing;
Cutter is along the feeding of x axle, and the trajectory that y axle both sides are finally left to surface of the work is respectively P 1p 1' and P 2p 2', wanting to guarantee the existence of outermost track, the coordinate figure that cutting edge track must meet trajectory and feed shaft intersection point is greater than the straight-line feed amount f that cutter rotates a circle p, obtain thus formula (5):
Figure FDA0000389611890000029
Above formula obtains
Figure FDA00003896118900000210
for the minimum value satisfying condition;
The maximal value of position angle is the position angle of the cutting point that the residual high highest point of workpiece surface appearance is corresponding, namely section angle of cut of twice feed in vertical feed direction axial cross section before and after cutter;
For direction of feed right side P 1p 1' this section of track,
Figure FDA00003896118900000211
size to discuss in two kinds of situation;
1) feed for the first time,
Figure FDA00003896118900000212
can be directly according to cutting depth a pcalculate formula (6) with the geometric relationship of Z value maximum point in cutter under this cutting depth and workpiece contact point:
Figure FDA00003896118900000213
2) from feed for the second time,
Figure FDA00003896118900000214
with feed line-spacing f psize have direct relation, can draw shown in the relational expression (7) between them:
Figure FDA00003896118900000215
P for direction of feed left side 2p 2' track, size and P 1p 1' there is relative relation, for last feed, size is shown in formula (8):
Figure FDA0000389611890000031
Feed is for the first time in the process of second from the bottom feed, and the outmost trajectory in left side is corresponding
Figure FDA0000389611890000032
value be all shown in formula (9):
Figure FDA0000389611890000033
(3) in order to obtain each position
Figure FDA0000389611890000034
the Cutting trajectory line that angle is corresponding, need to obtain respectively outermost trajectory P 1p 1' section and P 2p 2the tool motion time in ' section is interval, by obtaining the corner that P is ordered in this is interval, then try to achieve corresponding process time in conjunction with the speed of mainshaft of cutter, consider the spiral drag angle at place, cutting infinitesimal position, can obtain outermost pattern track, thereby obtain the surface topography of workpiece;
In order to set up P 1and P 2corresponding corner
Figure FDA0000389611890000035
with
Figure FDA0000389611890000036
formula, analyze P 1and P 2the corner of point in milling process;
For P 1point, from the negative y axle of starting point, x=0 place forwards P to 1time, P 1point horizontal ordinate meets formula (10) below:
Figure FDA0000389611890000037
That is:
Figure FDA0000389611890000038
Wherein, v is the amount of feeding of bulb cutter unit interval, and ω is cutter shaft rotational speed, and R is tool radius,
Figure FDA0000389611890000039
for P 1the position angle of point,
Figure FDA00003896118900000310
for position angle
Figure FDA00003896118900000311
locate corresponding cutter corner, and meet
Figure FDA00003896118900000319
For P 2point, from the negative y axle of starting point, x=0 place forwards P to 2time, P 2point horizontal ordinate meets formula (11) below:
Figure FDA00003896118900000312
That is:
Figure FDA00003896118900000313
Wherein, v is the amount of feeding of bulb cutter unit interval, and ω is cutter shaft rotational speed, and R is tool radius,
Figure FDA00003896118900000314
for P 2the position angle of point,
Figure FDA00003896118900000315
for position angle
Figure FDA00003896118900000316
locate corresponding cutter corner, and meet
Figure FDA00003896118900000320
Use Newton iteration method to solve nonlinear equation (10) and (11), obtain correspondence position
Figure FDA00003896118900000317
the cutter corner value at angle, in conjunction with the angular velocity of rotation of cutter, can be in the hope of the cutting time in this segment of cutting
Figure FDA00003896118900000318
φ is the cutter anglec of rotation;
Because the spiral cutting sword of rose cutter will cause spiral hysteresis phenomenon, the time that will participate in cutting to place, diverse location angle cutting point exerts an influence, and then affects Cutting trajectory, each position angle corresponding cutter spiral drag angle
Figure FDA0000389611890000042
can be expressed as formula (12):
Figure FDA0000389611890000043
Wherein, β 0for the maximum spiral drag angle of cutter, for the position angle of cutting point on Tool in Cutting sword, the time effects that spiral lags behind and produces
Figure FDA0000389611890000048
By the position angle of trying to achieve
Figure FDA0000389611890000045
scope, cutter tried to achieve
Figure FDA0000389611890000046
interval in corner φ, time that cutter turns over this corner be updated in the Cutting trajectory equation (2) of considering tool wear, can obtain the workpiece surface appearance of rose cutter milling, be expressed as formula (13);
Figure FDA0000389611890000047
CN201210162358.5A 2012-05-24 2012-05-24 Modeling method for milling surface appearance of workpiece Expired - Fee Related CN102681488B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210162358.5A CN102681488B (en) 2012-05-24 2012-05-24 Modeling method for milling surface appearance of workpiece

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210162358.5A CN102681488B (en) 2012-05-24 2012-05-24 Modeling method for milling surface appearance of workpiece

Publications (2)

Publication Number Publication Date
CN102681488A CN102681488A (en) 2012-09-19
CN102681488B true CN102681488B (en) 2014-01-29

Family

ID=46813563

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210162358.5A Expired - Fee Related CN102681488B (en) 2012-05-24 2012-05-24 Modeling method for milling surface appearance of workpiece

Country Status (1)

Country Link
CN (1) CN102681488B (en)

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103042434B (en) * 2012-12-08 2014-12-10 华中科技大学 Ultraprecise milling surface topography texture control method
CN103286636B (en) * 2013-06-20 2015-04-08 广东省南方彩色制版有限公司 Printing carving roller turning and grinding integrated machine cutter angle laser measuring and controlling device and adjusting method
CN103646141B (en) * 2013-12-05 2017-05-10 华中科技大学 Cutting force modeling method for flat bottom spiral end mill orthogonal turning milling shaft parts
CN103713576B (en) * 2013-12-31 2016-08-17 南京航空航天大学 Multi-axis milling processing workpiece surface appearance modeling method
CN104090528B (en) * 2014-06-30 2017-02-15 华南理工大学 Method suitable for cycloid high-speed milling machining path
CN106033025B (en) * 2015-03-11 2019-02-12 香港理工大学 A kind of Tool Wear Monitoring method and system
CN108256174B (en) * 2017-12-27 2021-05-11 西安建筑科技大学 Shaving engagement modeling method for analyzing influence of shaving process on tooth profile quality
CN108241784A (en) * 2018-01-08 2018-07-03 华南理工大学 The cavity multicutter method for milling that a kind of cycloidal path is combined with ring cutting track
CN108515217B (en) * 2018-04-09 2019-05-31 吉林大学 A kind of ball-end milling free form surface surface topography emulation mode
CN108705692B (en) * 2018-05-25 2020-06-30 哈尔滨工业大学 Micro-milling repair process method for laser damage on surface of large-caliber KDP crystal element
CN109765841A (en) * 2019-01-09 2019-05-17 西北工业大学 The space-time mapping method of online monitoring data and part Working position
WO2020217498A1 (en) * 2019-04-26 2020-10-29 三菱電機株式会社 Machining condition determination assistance device and machining condition determination device
CN110161963B (en) * 2019-06-14 2022-01-07 哈尔滨理工大学 Simulation model and verification method for milling cutter cutting error forming process
CN110262397B (en) * 2019-06-24 2020-10-23 北京理工大学 Method for modeling spiral trochoid motion trajectory and instantaneous cutting force in turn-milling machining space
CN110727246B (en) * 2019-10-29 2021-03-19 江苏长虹智能装备股份有限公司 Tool and workpiece instantaneous contact contour extraction method based on tool position file
CN111633468B (en) * 2020-05-27 2021-07-16 武汉理工大学 Method and device for determining contact condition of round-edge cutter based on cutting force
CN113070524B (en) * 2021-04-06 2022-04-29 湖南大学 Cutting processing method of special-shaped section bar
CN113210688B (en) * 2021-05-18 2023-09-05 西安航天华阳机电装备有限公司 Milling groove track determination method for throat flow channel groove of drill shank of spray pipe part
CN113204852B (en) * 2021-06-16 2023-03-24 山东大学 Method and system for predicting milling surface appearance of ball-end milling cutter
CN113704928B (en) * 2021-09-15 2023-11-14 上海交通大学 Inclined milling surface morphology generation method based on machining physical process
CN114444297A (en) * 2022-01-21 2022-05-06 安徽理工大学 Method for predicting surface topography of threaded workpiece in cyclone milling

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102004818B (en) * 2010-11-04 2012-08-29 西北工业大学 Instantaneous milling force modeling method in free curve shape part peripheral milling process

Also Published As

Publication number Publication date
CN102681488A (en) 2012-09-19

Similar Documents

Publication Publication Date Title
CN102681488B (en) Modeling method for milling surface appearance of workpiece
Zhu et al. Research on rotary surface topography by orthogonal turn-milling
CN103537743B (en) The method of multiaxis Surface NC Machining complex curved surface parts
CN101518843B (en) Method and apparatus for working a screw rotor, end mill for working , and method of manufacturing a screw compressor
CN111644909B (en) Method for solving grinding track of rear cutter face of woodworking forming milling cutter
CN103713576A (en) Modeling method for workpiece surface appearance machined through multi-axis milling
CN106934170B (en) Chatter stability lobes flap figure modeling method based on rose cutter and workpiece contact zone
CN106216747A (en) A kind of integral wheel 5-shaft linkage numerical control cut track path processing method
CN110262397B (en) Method for modeling spiral trochoid motion trajectory and instantaneous cutting force in turn-milling machining space
CN103586519B (en) Ladder-type trough layered milling rough machining method
US9459166B2 (en) Cutting resistance analysis device, cutting and machining device equipped with same, and cutting resistance analysis program
CN102004818A (en) Instantaneous milling force modeling method in free curve shape part peripheral milling process
CN103008986A (en) Numerical control boring and milling machining method of internal spherical surface
CN103645674A (en) A method for generating a mixed path of rough-semifine-fine milling of an integrated impeller blade
CN102581705A (en) Division compensation method for machining error during milling of workpiece
CN107544433B (en) Numerical-controlled machine tool machining process rose cutter and workpiece contact zone semi analytic modeling
CN102905847A (en) Adaptive control of machining process
CN103862346A (en) Non-instantaneous-pole envelope grinding method of spiral curved surface of superfine milling cutter
CN102393647B (en) Method for controlling cutting quantity of spiral bevel gear
CN106270680A (en) The processing method of profile of tooth seamed edge rounded corner inside a kind of duplex rack rails wheel
CN103286360B (en) A kind of tangential offset orthogonal turn-milling cutting method
CN100424601C (en) Method for processing shoe tree
CN101298104A (en) Method for milling turbo wheel
CN110516373A (en) A kind of method for milling of circular saw roughing special-shaped level
CN101763048A (en) Method for selecting roughing tool for aircraft complex component based on Voronoi Mountain variation

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140129

Termination date: 20150524

EXPY Termination of patent right or utility model