CN105608288B  One kind being based on overdamp effect milling parameter stability prediction method  Google Patents
One kind being based on overdamp effect milling parameter stability prediction method Download PDFInfo
 Publication number
 CN105608288B CN105608288B CN201610037067.1A CN201610037067A CN105608288B CN 105608288 B CN105608288 B CN 105608288B CN 201610037067 A CN201610037067 A CN 201610037067A CN 105608288 B CN105608288 B CN 105608288B
 Authority
 CN
 China
 Prior art keywords
 cutter
 milling
 damping
 directions
 kernel
 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
Classifications

 G—PHYSICS
 G06—COMPUTING; CALCULATING OR COUNTING
 G06F—ELECTRIC DIGITAL DATA PROCESSING
 G06F30/00—Computeraided design [CAD]
 G06F30/30—Circuit design
 G06F30/36—Circuit design at the analogue level
 G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

 G—PHYSICS
 G05—CONTROLLING; REGULATING
 G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
 G05B19/00—Programmecontrol systems
 G05B19/02—Programmecontrol systems electric
 G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or coordinated operations by means of programme data in numerical form
 G05B19/406—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or coordinated operations by means of programme data in numerical form characterised by monitoring or safety
Abstract
The present invention provides a kind of based on overdamp effect milling parameter stability prediction method,This method is in carrying out workpiece milling process,Obtain milling cutter geometrical structure parameter and milling process dynamic parameter,Determine the gross energy that plow power generates in workpiece milling process in a pirouette period,Utilize law of conservation of energy,Obtain equivalent linearity process damped coefficient,The damping of equivalent linearity process is changed into the equivalent processes damping of the directions x and the damping of the directions y equivalent processes,The directions x equivalent processes damped coefficient and the directions y equivalent processes damped coefficient input milling dynamics equation are obtained into the milling dynamics model of Kernelbased methods damping effect,The milling dynamics model of Kernelbased methods damping effect is solved using ZOA methods,Obtain the flutter stability model of Kernelbased methods damping effect,To obtain the relationship between marginal stability cuttingin and cutter rotating speed,And draw the flutter stability flap figure of Kernelbased methods damping effect.
Description
Technical field
The invention belongs to field of machining, and in particular to one kind being based on overdamp effect milling parameter stability prediction side
Method.
Background technology
Milling Process is widely used in the industries such as aerospace mold, and Regenerative Chatter is the main restriction for improving production efficiency
One of factor.The development of Highspeed machining theory keeps the incision frequency of cutter tooth synchronous with flutter frequency next by selecting the speed of mainshaft
Avoid the generation of flutter.In the highspeed milling stage of higher material removing rate, classical flutter flap figure provides accurate steady
Qualitative forecasting.However, in the low speed process segment, due to there are a large amount of complete vibration waves in a swing circle of main shaft
It is long, cause flap figure so dense so that classical flutter theory tends not to accurately predict stability.On the other hand, it tests
Observation shows that the stabilization cutting zone of system can dramatically increase when cutting speed is far below the intrinsic frequency of system of processing.This
Stablize dramatically increasing for cutting region when kind slow cutting, rubbing action institute between cutter flank and uneven workpiece surface can be attributed to
Caused cutting speed variation, i.e. process damping action.
The research of process damping effect is concentrated mainly in external experts and scholars' document, the domestic research phase in the field
To fewer.Professor Altintas thinks that process damping effect will be as most challenging research to the research of flutter stability
Project；Mentioning tool wear will cause process damped coefficient to increase, to increase flutter instability critical zone.Budak is mentioned
Process damping effect in low speed processing becomes one of research emphasis and difficult point with its complexity.The effect and cutting data are cut
It is closely related to cut the factors such as temperature, cutter material characteristic, shear surface variation；Heat height will be generated by cutting difficulttomachine material, easily be led
Cause tool wear, the finished work corrugated surface that the geometry for changing cutting edge is contacted with flank, to increase
Process damping effect.
Invention content
One kind being based on overdamp effect milling parameter stability prediction method, includes the following steps：
Step 1：In carrying out workpiece milling process, obtains milling cutter geometrical structure parameter and milling process dynamic is joined
Number；
The milling cutter geometrical structure parameter includes：Cutter number of teeth N_{f}, tool clearance λ and cutter diameter D；
The milling process dynamic parameter includes：Cutter flutter angular frequency_{c}, cutter amplitude A_{0}With tool angular velocity Ω.
Step 2：A pirouette period is determined according to milling cutter geometrical structure parameter and milling process dynamic parameter
The gross energy that plow power generates in interior workpiece milling process, using law of conservation of energy, the gross energy etc. which is generated
The energy that valence is generated in equivalent linearity process damping force in a cycle, obtains equivalent linearity process damped coefficient；
Step 2.1：The tangential displacement of cutter is determined according to milling cutter geometrical structure parameter and milling process dynamic parameter
With the radial displacement of cutter；
Step 2.2：The axial cuttingin a of cutter is divided into N_{z}A infinitesimal determines the tangential plow infinitesimal power of cutter cutter tooth
With the radial plow infinitesimal power of cutter tooth；
Step 2.3：Determine the energy e that plow power generates within a vibration period_{i}, i.e., cutter cutter tooth is in a vibration wave
Rimmer knife tooth footpath the sum of is done work to plow infinitesimal power and tangential plow infinitesimal power；
Step 2.4：According to cutter flutter angular frequency_{c}Determine that cutter pirouette stays in work in one week with tool angular velocity Ω
The oscillation mark number on part surface；
Step 2.5：Determine the gross energy of plow power, i.e., plow power is in the one week work(done of cutter pirouette；
Step 2.6：Using law of conservation of energy, the gross energy of plow power is equivalent to use linear viscous in a cycle
The energy that the process damping force of damping generates, obtains equivalent linearity process damped coefficient.
Step 3：The damping of equivalent linearity process is changed into the equivalent processes damping of the directions x and the damping of the directions y equivalent processes, is obtained
To the equivalent processes damped coefficient of the directions x equivalent processes damped coefficient and the directions y；
Step 4：The directions x equivalent processes damped coefficient and the directions y equivalent processes damped coefficient are inputted into milling dynamics side
Journey obtains the milling dynamics model of Kernelbased methods damping effect；
Step 5：The milling dynamics model that Kernelbased methods damping effect is solved using ZOA methods obtains Kernelbased methods damping
The flutter stability model of effect to obtain the relationship between marginal stability cuttingin and cutter rotating speed, and is drawn and was based on
The flutter stability flap figure of journey damping effect.
The flutter stability model of the Kernelbased methods damping effect is：
Wherein,For cutter and work piece interface product transfer function matrix,
The direct transmission function that the directions x are accumulated for cutter and work piece interface is that cutter accumulates the direct of the directions y with work piece interface
Transmission function, to intersect transmission function, a is the axial cuttingin of cutter, k_{t}For tangential cutting force constant,
N_{f}For cutter number of teeth, ω_{c}For cutter flutter angular frequency, T is cutter swing circle, and i is plural number.
Beneficial effects of the present invention：
The present invention proposes that one kind being based on overdamp effect milling parameter stability prediction method, in difficulttomachine material, complexity
Curved surface relatively low speed processing in, process damping effect play the role of in flutter stability it is very important, especially milling navigate
Empty material, have low heat conductivity titanium alloy when, will generate heat compare it is larger, so as to cause tool wear；Process damping system
Number will increase, proposed by the present invention steady based on overdamp effect milling parameter to change flutter instability region limit
Qualitative Forecast Methods are solved when low speed is processed, since flutter stability flap densification leads to not choose using stable region
The problem of cutting parameter, improves the accuracy of flutter instability regional prediction.Therefore in actual production, to joining in process
Several optimization improves production efficiency and has great importance.
Description of the drawings
Fig. 1 is the flow based on overdamp effect milling parameter stability prediction method in the specific embodiment of the invention
Figure；
Fig. 2 is that Kernelbased methods damp milling dynamics model in the specific embodiment of the invention；
Wherein, (a) is that cutter cuts part model, is (b) the vibration wavelength schematic diagram that cutter invades workpiece；
K_{x}For cutter x directional stiffness coefficients, C_{x}For the directions cutter x damped coefficient, K_{y}For cutter y directional stiffness coefficients, C_{y}For
The directions cutter y damped coefficient, F_{pt}For tangential plow power, F_{pr}For radial plow power, r_{j}For cutter cutter tooth j radial distances, u_{j}For knife
Have cutter tooth j tangential distances, λ is tool clearance, cutter_{j1}For 1 cutter tooth of cutter jth, cutter_{j}For jth of cutter tooth of cutter,
V is cutter linear velocity, A_{0}For cutter amplitude, L is the wavelength for the vibration that cutter invades workpiece, φ_{j}For cutter cutter tooth j and work
The contact angle of part, Ω are tool angular velocity；
Fig. 3 is that the gross energy that plow power generates is equivalent to equivalent linearity in a cycle in the specific embodiment of the invention
The energy that process damping force generates obtains the flow chart of equivalent linearity process damped coefficient；
Fig. 4 is not consider the flutter stability flap figure of process damping effect in the specific embodiment of the invention and be based on
The flutter stability flap figure comparison diagram of journey damping effect；
Fig. 5 is that the flutter of the Kernelbased methods damping effect in the case of different tool clearances in the specific embodiment of the invention is steady
Qualitative flap figure comparison diagram；
Fig. 6 is that the flutter of the Kernelbased methods damping effect in the case of different cutter amplitudes in the specific embodiment of the invention is steady
Qualitative flap figure comparison diagram；
Fig. 7 is that the flutter of the Kernelbased methods damping effect in the case of different tool stiffness in the specific embodiment of the invention is steady
Qualitative flap figure comparison diagram；
Fig. 8 quivers for the Kernelbased methods damping effect in the case of different cutter structures dampings in the specific embodiment of the invention
Shake the stability lobes diagram comparison diagram；
Fig. 9 is the flutter instability of the Kernelbased methods damping effect in the case of different number of teeth in the specific embodiment of the invention
Property threedimensional flap figure comparison diagram；
Figure 10 is the Kernelbased methods in the case of different radial the ratio between cuttingins and cutter diameter in the specific embodiment of the invention
The flutter stability threedimensional flap figure comparison diagram of damping effect.
Specific implementation mode
The specific embodiment of the invention is described in detail below in conjunction with the accompanying drawings.
The present invention provides one kind being based on overdamp effect milling parameter stability prediction method, the vibration decomposition of cutter
It for radial linear motion of reciprocating vibration and tangential, calculates within one vibration period of cutter tooth, since cutter flank invades
Then the invasiveness work done of workpiece machined surface is multiplied by one swing circle of cutter complete oscillation mark number, profit
With the linear damping of equal value of energy equivalence principle nonlinear dampling, flutter stability flap figure is solved using frequency domain method, low
In fast milling process, the flutter stability cuttingin limit can be efficiently and accurately predicted.Therefore to Optimizing Cutting Conditions, finished surface
Quality and processing efficiency provide theoretical direction.
One kind being based on overdamp effect milling parameter stability prediction method, as shown in Figure 1, including the following steps：
Step 1：In carrying out workpiece milling process, obtains milling cutter geometrical structure parameter and milling process dynamic is joined
Number.
In present embodiment, Kernelbased methods damp milling dynamics model as shown in Fig. 2, carbidetipped milling cutter geometry
Parameter includes：The cutter number of teeth N that helical angle is 0_{f}=4, tool clearance λ=10 degree, cutter diameter D=10mm；
Milling process dynamic parameter includes：Cutter flutter angular frequency_{c}, cutter amplitude A_{0}With tool angular velocity Ω.
Step 2：A pirouette period is determined according to milling cutter geometrical structure parameter and milling process dynamic parameter
The gross energy that plow power generates in interior workpiece milling process, using law of conservation of energy, the gross energy etc. which is generated
The energy that valence is generated in equivalent linearity process damping force in a cycle, obtains equivalent linearity process damped coefficient, such as Fig. 3 institutes
Show.
Step 2.1：The tangential displacement of cutter is determined according to milling cutter geometrical structure parameter and milling process dynamic parameter
With the radial displacement of cutter.
In present embodiment, according to cutter flutter angular frequency_{c}, cutter amplitude A_{0}, tool angular velocity Ω and cutter diameter D
The tangential displacement u (t) of determining cutter and the radial displacement r (t) of cutter are as shown in formula (1) and formula (2)：
R (t)=A_{0}sin(ω_{c}t) (2)
Wherein, v is cutter linear velocity, and t is the time.
Step 2.2：The axial cuttingin a of cutter is divided into N_{z}A infinitesimal determines the tangential plow infinitesimal power of cutter cutter tooth
With the radial plow infinitesimal power of cutter tooth.
In present embodiment, the axial cuttingin a of cutter is divided into N_{z}Shown in a infinitesimal such as formula (3)：
N_{z}=a/ Δs z (3)
Wherein, Δ z is infinitesimal.
The tangential plow infinitesimal power dF of determining cutter cutter tooth_{pr}[φ_{j}(t), k_{z}]With the radial plow infinitesimal of stage property cutter tooth
Power dF_{pt}[φ_{j}(t), k_{z}]As shown in formula (4) and formula (5)：
dF_{pr}[φ_{j}(t), k_{z}]=K_{sp}dV[φ_{j}(t), k_{z}]g[φ_{j}(t), k_{z}] (4)
dF_{pt}[φ_{j}(t), k_{z}]=μ dF_{pr}[φ_{j}(t), k_{z}] (5)
Wherein, φ_{j}(t)=Ω t+ (j1) 2 π/N_{f}, φ_{j}(t) it is the contact angle of cutter cutter tooth j and workpiece, by documents and materials
Know intrusion force constant K_{sp}=3000N/mm^{3}, V cutters intrusion workpiece volume, k_{z}∈N_{z}, by documents and materials know average friction coefficient μ=
0.3, g[φ_{j}(t), k_{z}]For unit jump function, for determining cutter tooth whether in cutting.
Unitstep function g[φ_{j}(t), k_{z}]As shown in formula (6)：
Wherein, φ_{st}For cutter entrance angle, φ_{ex}Angle is cut out for cutter.
Step 2.3：Determine the energy e that plow power generates within a vibration period_{i}, i.e., cutter cutter tooth is in a vibration wave
Rimmer knife tooth footpath the sum of is done work to plow infinitesimal power and tangential plow infinitesimal power.
In present embodiment, energy e that plow power generates within a vibration period_{i}Such as formula (7), formula (8) and formula (9) institute
Show：
e_{i}=e_{iu}+e_{ir} (7)
Wherein, e_{iu}For the tangential plow infinitesimal power work done of one vibration wavelength cutter tooth of cutter cutter tooth, e_{ir}For cutter cutter tooth
One vibration wavelength cutter tooth radial direction plow infinitesimal power work done, L=(2 π v)/ω_{c},
Step 2.4：According to cutter flutter angular frequency_{c}Determine that cutter pirouette stays in work in one week with tool angular velocity Ω
The oscillation mark number on part surface.
In present embodiment, according to cutter flutter angular frequency_{c}Cutter pirouette is determined with tool angular velocity Ω one week
Shown in the oscillation mark number such as formula (10) for staying in workpiece surface：
(ω_{c}/ Ω)=the π of η+ζ/2 (10)
Wherein, η is the integer number of one week oscillation mark for staying in workpiece surface of cutter pirouette, and the π of ζ/2 are to generate wave
The decimal of line number.
Step 2.5：Determine the gross energy of plow power, i.e., plow power is in the one week work(done of cutter pirouette.
In present embodiment, the gross energy E of plow power_{i}As shown in formula (11)：
E_{i}=E_{iu}+E_{ir}=η (e_{iu}+e_{ir}) (11)
Wherein, E_{iu}For one week tangential plow power infinitesimal power work done of pirouette, E_{ir}For one week radial plough of pirouette
Cut infinitesimal power work done.
Step 2.6：Using law of conservation of energy, the gross energy of plow power is equivalent to use linear viscous in a cycle
The energy that the process damping force of damping generates, obtains equivalent linearity process damped coefficient.
In present embodiment, the energy such as formula of the process damping force generation of linear viscous damping will be used in a cycle
(12) and shown in formula (13)：
Wherein, E^{p} _{ir}For equivalent radial process damping force within a pirouette period work done, E^{p} _{iu}It is equivalent
Tangential process damping force within a pirouette period work done.
Using law of conservation of energy, the gross energy of plow power is equivalent to use the mistake of linear viscous damping in a cycle
The energy that journey damping force generates, makes E^{p} _{ir}=E_{ir}, E^{p} _{iu}=E_{iu}, obtain radial equivalent linearity process damped coefficient c_{pr}Such as formula (14)
It is shown：
It can similarly obtain, tangential equivalent linearity process damped coefficient c_{pu}As shown in formula (15)：
Step 3：The damping of equivalent linearity process is changed into the equivalent processes damping of the directions x and the damping of the directions y equivalent processes, is obtained
To the equivalent processes damped coefficient of the directions x equivalent processes damped coefficient and the directions y.
In present embodiment, the directions x equivalent processes damped coefficient C is obtained_{px}With the equivalent processes damped coefficient C in the directions y_{py}
As shown in formula (16)：
Equivalent processes damped coefficient expression formula in the directions x and y is：
Step 4：The directions x equivalent processes damped coefficient and the directions y equivalent processes damped coefficient are inputted into milling dynamics side
Journey obtains the milling dynamics model of Kernelbased methods damping effect.
In present embodiment, the directions x equivalent processes damped coefficient and the input milling of the directions y equivalent processes damped coefficient are moved
Mechanical equation obtains shown in the milling dynamics model such as formula (17) of Kernelbased methods damping effect：
Wherein,For milling system quality coefficient, m_{xx}For cutter xaxis quality, m_{yy}For cutter yaxis matter
Amount,For tool stiffness, K_{x}For cutter x directional stiffness coefficients, K_{y}For cutter y directional stiffness coefficients,It is damped for milling system, C_{x}For the directions cutter x damped coefficient, C_{y}It is hindered for the directions cutter y
Buddhist nun's coefficient,For the dynamic displacement of cutter,Inside and outside cutter
The difference of cyclic shift, Δ x are the difference of the directions cutter x innerouter circulation displacement, and Δ y is the difference of the directions cutter y innerouter circulation displacement, x
(t) it is cutter x direction outer circulations in t moment displacement, y (t) is cutter y direction outer circulations in t moment displacement, and x (tT) is cutter
T moment displacement is circulated in the directions x, y (tT) is that t moment displacement is circulated in the directions cutter y, and T is cutter swing circle, k_{t}
=1596.3Mpa is tangential cutting force constant,For the direction coefficient of timevarying.
Shown in the direction coefficient A such as formulas (18) of timevarying：
Wherein, k_{r}=0.25 is radial cutting force constant, and φ is the instant contact angle of cutter and workpiece.
Step 5：The milling dynamics model that Kernelbased methods damping effect is solved using ZOA methods obtains Kernelbased methods damping
The flutter stability model of effect to obtain the relationship between marginal stability cuttingin and cutter rotating speed, and is drawn and was based on
The flutter stability flap figure of journey damping effect.
In present embodiment, the milling dynamics model of Kernelbased methods damping effect is solved using ZOA methods, was based on
Shown in the flutter stability model such as formula (19) of journey damping effect：
Wherein,For cutter and work piece interface product transfer function matrix,
The direct transmission function that the directions x are accumulated for cutter and work piece interface is that cutter accumulates the direct of the directions y with work piece interface
Transmission function, to intersect transmission function, i is plural number.
Obtain marginal stability cuttinginAs shown in formula (20)：
Wherein, Δ R^{p}It is the real part of consideration process damping milling dynamics formulation character value Δ, κ is the damping milling of consideration process
Cut the ratio of the imaginary part and real part of kinetics equation characteristic value Δ.
Wherein, consider shown in process damping milling dynamics formulation character value Δ such as formula (21)：
Wherein, a_{0}=φ^{p} _{xx}(iω_{c})φ^{p} _{vv}(iω_{c})(a_{xx}a_{yy}a_{xy}a_{yx}), a_{1}=φ^{p} _{xx}(iω_{c})φ^{p} _{yy}(iω_{c})
(a_{xx}a_{yy}a_{xy}a_{yx}),
ω_{nx}For the intrinsic frequency in the directions main shaft tooling system x, ω_{ny}For the intrinsic frequency in the directions main shaft tooling system y, ζ_{x}For main shaft knife
The structural damping ratio in the directions tool system x, ζ_{y}For the structural damping ratio in the directions main shaft tooling system y, ζ_{px}For the main shaft tooling system side x
To process damping ratio, ζ_{py}For the process damping ratio in the directions main shaft tooling system y.
It obtains shown in cutter range of speeds n such as formulas (22)：
Wherein, ψ=arctan κ are the phase shift that consideration process damps milling dynamics formulation character value Δ, l_{r}It is followed for vibration
Number of rings.
The modal parameter of main axle cutter system is as shown in table 1.
The modal parameter of 1 main axle cutter system of table
In present embodiment, go out not consider the flutter stability flap figure of process damping effect using Matlab Software on Drawing
Flutter stability flap figure comparison diagram with Kernelbased methods damping effect is as shown in figure 4, as shown in Figure 4, Kernelbased methods damping is imitated
The flutter stability flap figure answered can more accurately describe the relationship between marginal stability cuttingin and cutter rotating speed.
In present embodiment, the flutter stability flap figure pair of the Kernelbased methods damping effect in the case of different tool clearances
Than figure as shown in figure 5, when taking tool clearance λ=5 degree and λ=10 degree respectively, other parameters are constant, are programmed by MATLAB
Draw influence of the different relief angles to Kernelbased methods damping effect flutter stability.From figure 5 it can be seen that cuttingin stability limit with
The increase of tool clearance and reduce.
In present embodiment, the flutter stability flap figure pair of the Kernelbased methods damping effect in the case of different cutter amplitudes
Than scheming as shown in fig. 6, taking cutter amplitude A_{0}=20 μm and A_{0}At=40 μm, other parameters are constant, and picture is programmed by MATLAB
Go out influence of the cutter amplitude to the flutter stability of Kernelbased methods damping effect.From fig. 6 it can be seen that cuttingin stability limit with
It the increase of cutter amplitude and increases.
In present embodiment, the flutter stability flap figure pair of the Kernelbased methods damping effect in the case of different tool stiffness
Than scheming as shown in fig. 7, it is respectively 1 times of rigidity in table 1,1.5 times, 2 times to take tool stiffness, other parameters are constant, pass through
MATLAB is programmed the influence for drawing tool stiffness to the flutter stability of Kernelbased methods damping effect.It can from Fig. 7
Go out, cuttingin stability limit increases with the increase of tool stiffness.
In present embodiment, the flutter stability flap of the Kernelbased methods damping effect in the case of different cutter structure dampings
Figure comparison diagram is as shown in figure 8, it is respectively 1 times, 1.5 times, 2 times of structural damping in table 1 to take cutter structure damping, and other parameters are not
Become, the influence for drawing cutter structure damping to the flutter stability of Kernelbased methods damping effect is programmed by MATLAB.From
As can be seen that cuttingin stability limit increases with the increase that cutter structure damps in Fig. 8.
In present embodiment, the flutter stability threedimensional flap figure of the Kernelbased methods damping effect in the case of different number of teeth
For comparison diagram as shown in figure 9, when number of teeth being taken to be respectively 2,3,4,5, other parameters are constant, are programmed by MATLAB and draw knife
Influence of the number of teeth to the flutter stability of Kernelbased methods damping effect.It can be seen in figure 9 that cuttingin stability limit is with cutter
The increase of structural damping and reduce.
In present embodiment, the flutter of the Kernelbased methods damping effect in the case of different radial direction the ratio between cuttingins and cutter diameter
Stability threedimensional flap figure comparison diagram is as shown in Figure 10, take the ratio between radial cuttingin and cutter diameter be 0.5,0.7,0.9 when, other
Parameter constant is programmed the flutter for comparing Kernelbased methods damping effect for drawing radial cuttingin and cutter diameter by MATLAB
The influence of stability.It can be seen from fig. 10 that cuttingin stability limit with the increase of the ratio between radial cuttingin and cutter diameter and
Reduce.
Claims (3)
1. one kind being based on overdamp effect milling parameter stability prediction method, which is characterized in that include the following steps：
Step 1：In carrying out workpiece milling process, milling cutter geometrical structure parameter and milling process dynamic parameter are obtained；
Step 2：Work in a pirouette period is determined according to milling cutter geometrical structure parameter and milling process dynamic parameter
The gross energy that the plow power generates is equivalent to by the gross energy that plow power generates in part milling process using law of conservation of energy
The energy that equivalent linearity process damping force generates in a cycle, obtains equivalent linearity process damped coefficient, is as follows：
Step 2.1：The tangential displacement and knife of cutter are determined according to milling cutter geometrical structure parameter and milling process dynamic parameter
The radial displacement of tool；
Step 2.2：The axial cuttingin a of cutter is divided into N_{z}A infinitesimal determines the tangential plow infinitesimal power and cutter tooth of cutter cutter tooth
Radial plow infinitesimal power；
Step 2.3：Determine the energy e that plow power generates within a vibration period_{i}, i.e., cutter cutter tooth is in a vibration wave rimmer knife
Tooth footpath the sum of is done work to plow infinitesimal power and tangential plow infinitesimal power；
Step 2.4：According to cutter flutter angular frequency_{c}Determine that cutter pirouette stays in workpiece table in one week with tool angular velocity Ω
The oscillation mark number in face；
Step 2.5：Determine the gross energy of plow power, i.e., plow power is in the one week work(done of cutter pirouette；
Step 2.6：Using law of conservation of energy, the gross energy of plow power is equivalent to use linear viscous damping in a cycle
Process damping force generate energy, obtain equivalent linearity process damped coefficient；
Step 3：The damping of equivalent linearity process is changed into the equivalent processes damping of the directions x and the damping of the directions v equivalent processes, obtains x
The equivalent processes damped coefficient of direction equivalent processes damped coefficient and the directions y；
Step 4：The directions x equivalent processes damped coefficient and the directions y equivalent processes damped coefficient input milling dynamics equation are obtained
To the milling dynamics model of Kernelbased methods damping effect；
Step 5：The milling dynamics model that Kernelbased methods damping effect is solved using ZOA methods, obtains Kernelbased methods damping effect
Flutter stability model, to obtain the relationship between marginal stability cuttingin and cutter rotating speed, and draw Kernelbased methods resistance
The flutter stability flap figure of Buddhist nun's effect.
2. according to claim 1 be based on overdamp effect milling parameter stability prediction method, which is characterized in that described
Milling cutter geometrical structure parameter includes：Cutter number of teeth N_{f}, tool clearance λ and cutter diameter D；
The milling process dynamic parameter includes：Cutter flutter angular frequency_{c}, cutter amplitude A_{0}With tool angular velocity Ω.
3. according to claim 1 be based on overdamp effect milling parameter stability prediction method, which is characterized in that described
The flutter stability model of Kernelbased methods damping effect is：
Wherein,For cutter and work piece interface product transfer function matrix,For knife
The direct transmission function of tool and the directions work piece interface product x,The direct transmission in the directions y is accumulated for cutter and work piece interface
Function,To intersect transmission function, a is the axial cuttingin of cutter, k_{t}For tangential cutting force constant, N_{f}For
Cutter number of teeth, ω_{c}For cutter flutter angular frequency, T is cutter swing circle, and i is plural number.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201610037067.1A CN105608288B (en)  20160120  20160120  One kind being based on overdamp effect milling parameter stability prediction method 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201610037067.1A CN105608288B (en)  20160120  20160120  One kind being based on overdamp effect milling parameter stability prediction method 
Publications (2)
Publication Number  Publication Date 

CN105608288A CN105608288A (en)  20160525 
CN105608288B true CN105608288B (en)  20181023 
Family
ID=55988220
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201610037067.1A Active CN105608288B (en)  20160120  20160120  One kind being based on overdamp effect milling parameter stability prediction method 
Country Status (1)
Country  Link 

CN (1)  CN105608288B (en) 
Families Citing this family (11)
Publication number  Priority date  Publication date  Assignee  Title 

CN106843147B (en) *  20161215  20200317  太原科技大学  Method for predicting milling stability based on Hamming formula 
CN106647625B (en) *  20161215  20181221  太原科技大学  A method of based on Gear formula predictions milling stability 
CN106940746A (en) *  20170314  20170711  东北大学  The parallel time domain method of milling parameter stability prediction based on thinwall part 
CN107423489B (en) *  20170627  20190531  西北工业大学  Thinwall part milling process stability method for quick predicting 
CN107914183B (en) *  20171205  20200407  上海复合材料科技有限公司  Chatter stability prediction method for milling carbon fiber laminated plate 
CN108804749B (en) *  20180419  20190430  西北工业大学  Thinwall part milling process Damping Modeling method 
CN108647413B (en) *  20180427  20210706  北京理工大学  Comprehensive prediction method for position error and stability of fine surface 
CN108746795B (en) *  20180620  20200214  哈尔滨理工大学  Method for predicting flutter in numerical control milling of mold cavity 
CN109375578B (en) *  20181031  20200616  湖南工学院  Efficient machining control method for deep hole in engine oil pump shell 
CN111722586B (en) *  20200622  20230915  华中科技大学无锡研究院  Intermittent milling vibration sectional control method 
CN112621382B (en) *  20201217  20220524  北京工业大学  Online tool wear monitoring method based on cutting vibration damping characteristics 
Citations (4)
Publication number  Priority date  Publication date  Assignee  Title 

CN101493686A (en) *  20090226  20090729  上海交通大学  Cutting tool mode parameter uncertain curve fiveshaft numerical control process parameter optimizing method 
CN102554326A (en) *  20120217  20120711  上海交通大学  Milling finish machining method based on dynamic rigidity of impeller blade 
CN103823945A (en) *  20140313  20140528  大连理工大学  Flutter stability domain modeling approach for face cutting process 
CN104484506A (en) *  20141125  20150401  东北大学  Turning flutter forecasting method based on reliability lobe graph 

2016
 20160120 CN CN201610037067.1A patent/CN105608288B/en active Active
Patent Citations (4)
Publication number  Priority date  Publication date  Assignee  Title 

CN101493686A (en) *  20090226  20090729  上海交通大学  Cutting tool mode parameter uncertain curve fiveshaft numerical control process parameter optimizing method 
CN102554326A (en) *  20120217  20120711  上海交通大学  Milling finish machining method based on dynamic rigidity of impeller blade 
CN103823945A (en) *  20140313  20140528  大连理工大学  Flutter stability domain modeling approach for face cutting process 
CN104484506A (en) *  20141125  20150401  东北大学  Turning flutter forecasting method based on reliability lobe graph 
NonPatent Citations (2)
Title 

Altintas Y. Analytical prediction of chatter stability in milling;Budak E, Altintas Y;《Journal of Dynamic Systems, Measurement and Control》;19880131;第120卷(第1期);第2230页 * 
考虑过程阻尼的铣削稳定性建模及仿真分析;董亚峰;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20151215;第B022304页 * 
Also Published As
Publication number  Publication date 

CN105608288A (en)  20160525 
Similar Documents
Publication  Publication Date  Title 

CN105608288B (en)  One kind being based on overdamp effect milling parameter stability prediction method  
Seguy et al.  Suppression of period doubling chatter in highspeed milling by spindle speed variation  
Altintas et al.  Virtual high performance milling  
Ozlu et al.  Analytical modeling of chatter stability in turning and boring operations—part I: model development  
CN102184276B (en)  Method for optimizing feeding speed in cutting machining  
Kashyzadeh et al.  Study of chatter analysis in turning tool and control methods–a review  
Elhami et al.  Experimental study of surface roughness and tool flank wear during hybrid milling  
Zhai et al.  Influence of cutting parameters on force coefficients and stability in plunge milling  
Sun et al.  Modeling and simulation for micromilling mechanisms  
Rahman et al.  Effective corner machining via a constant feed rate looping tool path  
Dikshit et al.  Chatter and dynamic cutting force prediction in highspeed ball end milling  
Fromentin et al.  Effect of thread milling penetration strategies on the dimensional accuracy  
Baburaj et al.  Development and experimental validation of a mechanistic model of cutting forces in microball end milling of full slots  
Li et al.  Modeling of instantaneous cutting force for large pitch screw with vibration consideration of the machine tool  
Kiswanto et al.  Fundamental aspects in designing vibration assisted machining: a review  
Yemelyanov et al.  Improving machining accuracy of cnc machines with innovative design methods  
Zhu  Modeling of the machining process  
Hrytsay et al.  Load Parameters of the Gear Machining by Power Skiving and Their Influence on the Machining System  
Guo et al.  Study on the stability for nonuniform helix angle tools in the milling process  
Hadraba et al.  Virtual twin of the multispindle lathe for the chatter timedomain analysis  
Xu et al.  Research of gyroscopic effects on the stability of high speed milling  
Ozoegwu et al.  Time domain chatter stability comparison of turning and milling processes  
Cao et al.  CBRBased Cutting Parameter Selection System and Searching Algorithm for Metal Machining Operations  
Xiao et al.  Application of a new geometric modeling on machining cutting chatter analysis  
Nejman et al.  Relation between power and linear model of dynamic cutting coefficients 
Legal Events
Date  Code  Title  Description 

C06  Publication  
PB01  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
GR01  Patent grant  
GR01  Patent grant 