CN106965032A - Thin-wall part milling parameter suppressing method - Google Patents
Thin-wall part milling parameter suppressing method Download PDFInfo
- Publication number
- CN106965032A CN106965032A CN201710171598.4A CN201710171598A CN106965032A CN 106965032 A CN106965032 A CN 106965032A CN 201710171598 A CN201710171598 A CN 201710171598A CN 106965032 A CN106965032 A CN 106965032A
- Authority
- CN
- China
- Prior art keywords
- milling
- workpiece
- matrix
- thin
- mass
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000003801 milling Methods 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000012545 processing Methods 0.000 claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims description 57
- 239000000463 material Substances 0.000 claims description 46
- 238000005520 cutting process Methods 0.000 claims description 33
- 238000013016 damping Methods 0.000 claims description 24
- 238000012360 testing method Methods 0.000 claims description 19
- 238000002474 experimental method Methods 0.000 claims description 12
- 238000010606 normalization Methods 0.000 claims description 9
- 230000002068 genetic effect Effects 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000005316 response function Methods 0.000 claims description 5
- 238000011438 discrete method Methods 0.000 claims description 4
- 230000035772 mutation Effects 0.000 claims description 4
- 230000001629 suppression Effects 0.000 claims description 4
- 238000010998 test method Methods 0.000 claims description 4
- 239000000654 additive Substances 0.000 claims description 2
- 230000000996 additive effect Effects 0.000 claims description 2
- 238000010079 rubber tapping Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 description 5
- 238000003754 machining Methods 0.000 description 4
- 229910000831 Steel Inorganic materials 0.000 description 3
- 229910045601 alloy Inorganic materials 0.000 description 3
- 239000000956 alloy Substances 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 229910000838 Al alloy Inorganic materials 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000003631 expected effect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 229910001008 7075 aluminium alloy Inorganic materials 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q11/00—Accessories fitted to machine tools for keeping tools or parts of the machine in good working condition or for cooling work; Safety devices specially combined with or arranged in, or specially adapted for use in connection with, machine tools
- B23Q11/0032—Arrangements for preventing or isolating vibrations in parts of the machine
- B23Q11/0035—Arrangements for preventing or isolating vibrations in parts of the machine by adding or adjusting a mass, e.g. counterweights
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Numerical Control (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
The invention discloses a kind of thin-wall part milling parameter suppressing method, the technical problem for solving existing milling stability Forecasting Methodology poor practicability.Technical scheme is that the part of thin-walled parts kinetic parameter is changed by additional mass, sets up a kind of efficient processing technology to improve the stable region of Milling Process, reliable parameter range of choice is provided for thin-wall part high-rate wireless LAN；Final chosen using optimized algorithm can realize no flutter, efficient machined parameters, realize the high speed of thin-wall part without flutter Milling Process.The present invention is changed by the part to thin-walled parts kinetic parameter, a kind of efficient processing technology is set up to improve the stable region of Milling Process, preferably resolve workpiece starting and final position two ends poor rigidity, stable region scope is small, seriously the problem of the selection of restriction milling process machined parameters；Reliable parameter range of choice is provided for thin-wall part high-rate wireless LAN, the high speed of thin-wall part is realized without flutter Milling Process, practicality is good.
Description
Technical Field
The invention belongs to the field of thin-wall part manufacturing, and particularly relates to a milling chatter suppression method for a thin-wall part.
Background
Document 1 "Song Q, Liu Z, Wan Y, et al application of Sherman-Morrison-Woodbury for a glass in an alternating dynamic of a peripheral milling for a thin-walled component [ J ]. International Journal of Mechanical Sciences,2015,96-97: 79-90" discloses a milling stability prediction method using the Sherman-Morrison-Woodbury formula to take into account the effect of material removal on the kinetic parameters of thin-walled parts during milling. The method comprises the steps of dispersing a milling process, obtaining a variation rule of a kinetic parameter of a dispersed thin-wall part along with material removal in the milling process through a Sherman-Morrison-Woodbury formula so as to obtain a corresponding kinetic parameter, then obtaining a relation between the axial rotation speed and the axial cutting depth in each dispersing process by utilizing a stability solving equation, and finally obtaining a prediction three-dimensional graph of influence of material removal on stability in the milling process.
Document 2 "Yang Y, Zhang WH, Ma YC, et al, character prediction for using a thin-walled workpiece with a curved surface [ J ]. International journal of Machine Tools and manufacturing, 2016,109: 36-48" discloses a milling stability prediction method that considers the variation of the dynamic parameters of the workpiece at different tool milling positions and axial heights at the same time. The method comprises the steps of firstly dispersing a cutter and a workpiece along the axial direction to obtain a kinetic equation of each dispersed part. And then dispersing the milling process to obtain a kinetic equation of each milling process. And finally, solving the kinetic equation to obtain a milling stability prediction method considering the milling position and the axial height kinetic parameter change.
The above documents all consider the influence of material removal on the dynamic parameters of the workpiece in the milling process, and predict the stable regions at different tool positions; however, the problems of poor workpiece rigidity and low milling stability region at the initial position and the final position of the milling process are not effectively solved, so that the selectable range of parameters of the milling process is small, and the processing efficiency cannot be improved.
Disclosure of Invention
The invention provides a thin-wall part milling chatter suppression method, aiming at overcoming the defect that the existing milling stability prediction method is poor in practicability. The method establishes an efficient processing technique method to improve the stable domain of milling processing through local modification of the dynamic parameters of the thin-wall part by additional mass, and provides a reliable parameter selection range for high-speed milling processing of the thin-wall part; finally, the optimization algorithm is used for selecting machining parameters capable of achieving flutter-free and high efficiency, and high-speed flutter-free milling machining of the thin-wall part is achieved. According to the invention, an efficient processing method is established to improve the stable region of milling by locally modifying the dynamic parameters of the thin-wall part, so that the problems that the rigidity of two ends of the starting position and the ending position of a workpiece is poor, the stable region range is small, and the selection of the processing parameters in the milling process is seriously restricted are solved; the method provides a reliable parameter selection range for the high-speed milling of the thin-wall part, and realizes the high-speed non-flutter milling of the thin-wall part.
The technical scheme adopted by the invention for solving the technical problems is as follows: a thin-wall part milling chatter suppression method is characterized by comprising the following steps:
firstly, establishing a milling dynamic model which is in multipoint contact and considers the deformation of a cutter and a workpiece; the motion equation of the milling system is as follows:
wherein,t(t) is a vector representing the modal displacement of the tool;wp(t) is a vector representing modal displacement of the workpiece; zetatDiagonal matrix, ζ, representing the damping ratio of the toolwpA diagonal matrix representing a damping ratio of the workpiece; omegatA diagonal matrix representing the natural frequency of the tool; omegawpA diagonal matrix representing the natural frequency of the workpiece; u shapetRepresenting the modal shape, U, of the tool after mass normalizationwpRepresenting the modal shape of the workpiece after mass normalization; q represents the number of nodes for differentiating the tool from the workpiece in the axial direction; f (t) represents a matrix of milling forces at each contact point;
step two, after the selected milling cutter is installed on a main shaft of a machine tool, measuring the feeding direction of the cutter and modal parameters vertical to the feeding direction by adopting a multi-point tapping test method and a linear interpolation method; zetam,t,x，ζm,t,yRepresents the damping ratio of the mth order; omegam,t,x，ωm,t,yA system natural frequency representing the mth order; u shapem,t,x，Um,t,yIs a matrix of dimension q × 1, representing the mode shape of the mth order of the system mass normalization;
step three, establishing a finite element model of the thin-wall workpiece; in the process of establishing the model, a set of removed materials, final workpieces and additional mass is respectively established; giving the model corresponding material properties andadding constraints and loads which are consistent with actual processing to the workpiece, then carrying out finite element analysis to obtain the natural frequency of the whole workpiece, extracting the mass matrix and the rigidity matrix of each unit, finally assembling to obtain a whole workpiece model and the mass matrix and the rigidity matrix of the additional mass, and further obtaining the natural frequency omega of each order of the whole workpiece modelwpSum mode vibration type Uwp；
Step four, assuming that the damping ratio is unchanged, extracting the damping characteristic of the workpiece by using a hammering mode experiment; determining a damping ratio matrix zeta of a workpiece by measuring the frequency response function of the workpiece at different points and fitting itwp；
Step five, using the modal parameters obtained by the test in the step two and the natural frequency omega of each order of the whole workpiece model obtained in the step three and the step fourwpNormal mode vibration type UwpAnd damping ratio matrix ζ of workpiecewpSubstituting the obtained result into the step I, respectively solving the state equations of the tool at the milling initial position and the milling final position by using a semi-discrete method to obtain the axial cutting depth apAnd a stability lobe graph of the milling initial position and the milling final position with the main shaft rotating speed n as a variable;
step six, the maximum material removal rate MRR is an objective function, and the MRR is ap×ae×n×N×f；apIndicating axial cutting depth, aeThe radial cutting depth is represented, N represents the rotating speed of a main shaft, N represents the number of cutter teeth, and f represents the feed amount of each tooth; radial cutting depth a during millingeThe number N of cutter teeth and the feed amount f of each tooth are determined in advance; by axial cutting-in of apAnd a main shaft rotating speed n variable, setting population size, random seed generation probability, mutation probability, cross probability and genetic algebra parameters by respectively taking the stability lobe graphs of the milling initial position and the milling end position obtained in the step five as constraints, and respectively obtaining optimized axial cutting depths a of the milling initial position and the milling end position by utilizing a genetic algorithmpAnd a spindle speed n; selecting a corresponding to the initial position and the end positionpThe smallest of which is the selected processing parameter;
step seven, modifyingThe material density rho and the Young modulus E in the material attribute change the mass matrix and the rigidity matrix of the additional mass, the material density rho and the Young modulus E in the material attribute are modified by using a Taguchi method orthogonal test to change the mass matrix and the rigidity matrix of the additional mass, and then the dynamic parameters of the modified workpiece are obtained through assembly; substituting the new kinetic parameters into the fifth step to obtain a stability lobe graph corresponding to the orthogonal test; through the sixth step, the MRR corresponding to the Taguchi method orthogonal experiment and the axial cutting depth a corresponding to the MRR are obtainedpAnd a spindle speed n; finally, analyzing the data obtained by the Taguchi test, and selecting the optimal additional mass combination which can maximize the MRR and the corresponding axial cutting depth apAnd a spindle speed n.
The invention has the beneficial effects that: the method establishes an efficient processing technique method to improve the stable domain of milling processing through local modification of the dynamic parameters of the thin-wall part by additional mass, and provides a reliable parameter selection range for high-speed milling processing of the thin-wall part; finally, the optimization algorithm is used for selecting machining parameters capable of achieving flutter-free and high efficiency, and high-speed flutter-free milling machining of the thin-wall part is achieved. According to the invention, an efficient processing method is established to improve the stable region of milling by locally modifying the dynamic parameters of the thin-wall part, so that the problems that the rigidity of two ends of the starting position and the ending position of a workpiece is poor, the stable region range is small, and the selection of the processing parameters in the milling process is seriously restricted are solved; the method provides a reliable parameter selection range for the high-speed milling of the thin-wall part, and realizes the high-speed non-flutter milling of the thin-wall part.
The present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.
Drawings
FIG. 1 is a schematic diagram of a thin-wall part milling dynamic model considering deformation of a cutter and a workpiece in the method, wherein the cutter and the workpiece are respectively differentiated into q micro elements along the axial direction;
FIG. 2 is a schematic diagram of a multi-tap test tool frequency response function, PX1, PX2, PX3, PX4 being four measurement points in the tool feed direction, PY1, PY2, PY3, PY4 being four measurement points perpendicular to the tool feed direction;
FIG. 3 is a sheet and additive mass model validated in an embodiment of the method of the invention;
FIG. 4 is a graph of stability lobes for initial and final milling positions in an embodiment of the method of the present invention.
Detailed Description
The following examples refer to fig. 1-4.
Example 1 sheet size 115mm × 36mm × 3.5.5 mm, material 7075 aluminum alloy, modulus of elasticity 71GPa, density 2810kg/m3The Poisson's ratio is 0.33; the cutter is a hard alloy milling cutter with the number of cutting edges of 2, the diameter of 15.875mm and the helix angle of 30 degrees, and the extension length of the cutter is 78 mm.
Firstly, respectively differentiating a cutter and a workpiece into 40 infinitesimals along the axial direction by utilizing a thin-wall part milling dynamic model of the deformation of the cutter and the workpiece, namely q is 41; respectively establishing a milling kinetic equation at each infinitesimal position, and solving a dynamic milling force:
secondly, after the cutter is installed on a main shaft of the machine tool, measuring the feeding direction of the cutter and modal parameters vertical to the feeding direction by adopting a multi-point knocking test method and a linear interpolation method; zetam,t,x，ζm,t,yRepresents the damping ratio of the mth order; omegam,t,x，ωm,t,yA system natural frequency representing the mth order; u shapem,t,x，Um,t,yIs a matrix of q × 1 dimension, representing the m-th order modal vibration of system quality normalizationMolding; in a modal knock experiment, because of the influence of the size of an accelerometer, 4 points are selected and are measured by the accelerometer to obtain a modal vibration mode, and modal displacements at the other infinitesimal points are obtained by a linear interpolation method; four measuring points PX1, PX2, PX3, PX4 in the tool feed direction are located at distances of 5, 23, 50, 72mm from the tool tip point, respectively, and four measuring points PY1, PY2, PY3, PY4 perpendicular to the tool feed direction are located at distances of 0, 29, 43, 72mm from the tool tip point, respectively; selecting the first three modes (m is 1,2 and 3) which are most easily deformed, and providing modal analysis and identification modal parameters below; from U'm,t,x，U'm,t,y(m-th order modal shape of tool system quality normalization obtained by four-point test) is subjected to linear interpolation to obtain Um,t,x，Um,t,yFinally obtain Ut；
Ut＝[U1,t…U1,m,t]
Establishing a finite element model of the thin-wall workpiece, respectively establishing a set of removed materials (the size is 115 × 36 × 0.5.5 mm), a final workpiece (the size is 115 × 36 × 3mm) and three pieces of additional mass (the size is 18 × 18 × 3mm) in the process of establishing the model, giving corresponding material attributes to the model, adding constraints and loads which are consistent with actual processing to the workpiece, performing finite element analysis to obtain the inherent frequency of the whole workpiece, extracting the mass matrix and the rigidity matrix of each unit, finally assembling to obtain the mass matrix and the rigidity matrix of the whole workpiece model and the additional mass, and further obtaining the inherent frequency omega of each order of the whole workpiece modelwpSum mode vibration type Uwp；
Fourthly, assuming that the damping ratio is unchanged, extracting the damping characteristic of the workpiece by using a hammering mode experiment; determining a damping ratio matrix zeta of a workpiece by measuring the frequency response function of the workpiece at different points and fitting itwp；
Fifthly, using the modal parameters obtained by the test in the second step and the natural frequency omega of each order of the whole workpiece model obtained in the third and fourth stepswpNormal mode vibration type UwpAnd damping ratio matrix ζ of workpiecewpSubstituting the obtained result into the step I, respectively solving the state equations of the tool at the milling initial position and the milling final position by using a semi-discrete method to obtain the axial cutting depth ap(mm) stability lobe plot of milling initial position (solid line 1) and end position (dashed line 2) with spindle speed n (rpm) as variables;
sixthly, taking the maximum material removal rate MRR as an objective function, wherein MRR is ap×ae×n×N×f，apIndicating axial cutting depth, aeThe radial cutting depth is represented, N represents the rotating speed of a main shaft, N represents the number of cutter teeth, and f represents the feed amount of each tooth; during the milling processMedium radial cutting depth aeThe cutter tooth number N is 2, and the feed amount f of each tooth is 0.1 mm/revolution/tooth; these parameters have been determined in advance; by axial cutting-in of apAnd a variable of the rotation speed n of the main shaft, respectively taking the stability lobe graphs of the milling initial position and the milling final position obtained in the step five as constraints, and setting the range of the axial cutting depth to be 7000 a or more according to actual processing parameterspLess than or equal to 12000; setting the population size to be 20, the random seed generation probability to be 0.12221, the mutation probability to be 0.7, the cross probability to be 0.8 and the genetic algebra parameter to be 30, and respectively obtaining the optimized axial cutting depths a of the initial milling position and the final milling position by utilizing a genetic algorithmpAnd a spindle speed n; selecting a corresponding to the initial position and the end positionpThe smallest of which is the selected processing parameter;
starting position: a isp1.852 n 12000 end position: a isp＝1.675 n＝12000
The parameters finally selected are as follows: a isp＝1.675 n＝12000 MRR＝2010
Seventhly, adding mass blocks with different masses on the thin plate to locally modify the mass matrix and the rigidity matrix of the workpiece; according to the material parameters of several common materials given in the table I, the method L of the Taguo method is utilized16(44) Modifying the material density rho and the Young modulus E in the material attribute through an orthogonal test to change the mass matrix and the rigidity matrix of three additional masses, and then assembling to obtain the dynamic parameters of the modified workpiece; substituting the new kinetic parameters into the fifth step to obtain a stability lobe graph corresponding to the orthogonal test; through the sixth step, the product is obtained according to the 'Tiankou method' L16(44) Table three data corresponding to the orthogonal experiment; finally, three mass matrixes and stiffness matrixes of the additional mass under different material attributes, corresponding MRRs and corresponding axial cutting depths a are obtainedpAnd a spindle speed n; k1,K2,K3,K4Shows the influence of various factors on the indexes of the horizontal test, and obtains the optimal additional mass block pasting mode according to the influence, namely 343 the combination, namely the material 45 for the mass block 1 at the starting position#Steel, middle massCopper for the mass 2 and 45 for the terminating mass 3#Steel; the processing parameters obtained were: axial cutting depth ap3.247 and the main shaft speed n 12000, the material removal rate MRR 3896.4; compared with the processing parameters and the material removal rate obtained by not sticking the mass block, the material removal efficiency is improved by 93.8 percent; the method is proved to achieve good expected effect and have good practicability.
Table one: material parameters of several common materials
Table two: use of Tiankou L16(44) Data from orthogonal experiments
Example 2 sheet size 100mm × 40mm × 4.5.5 mm, material aluminum alloy 7075, elastic modulus 71GPa, density 2810kg/m3The Poisson's ratio is 0.33; the cutter is a hard alloy milling cutter with the number of cutting edges of 2, the diameter of 15.875mm and the helix angle of 30 degrees, and the extension length of the cutter is 78 mm.
Firstly, respectively differentiating a cutter and a workpiece into 30 infinitesimals along the axial direction by using a thin-wall part milling dynamic model of the deformation of the cutter and the workpiece, namely q is 31; respectively establishing a milling kinetic equation at each infinitesimal position, and solving a dynamic milling force:
secondly, after the tool is arranged on a main shaft of the machine tool, the feeding direction of the tool is measured and the feeding direction of the tool is perpendicular to the feeding direction of the tool by adopting a multipoint knocking test method and a linear interpolation methodModal parameters of the feed direction; zetam,t,x，ζm,t,yRepresents the damping ratio of the mth order; omegam,t,x，ωm,t,yA system natural frequency representing the mth order; u shapem,t,x，Um,t,yThe matrix is q × 1 dimension matrix, representing the m-th order modal shape of system mass normalization, in a modal knock experiment, selecting 4 points to obtain the modal shape by accelerometer measurement, and obtaining modal displacement at the rest infinitesimal points by a linear interpolation method, wherein four measurement points PX1, PX2, PX3 and PX4 in the cutter feeding direction are respectively located at the positions 5, 23, 50 and 72mm away from the cutter point, four measurement points PY1, PY2, PY3 and PY4 perpendicular to the cutter feeding direction are respectively located at the positions 0, 29, 43 and 72mm away from the cutter point, the most deformable first three orders (m is 1,2 and 3) are selected, modal analysis and identification modal parameters are given below, and U's modal shape is used'm,t,x，U'm,t,y(m-th order modal shape of tool system quality normalization obtained by four-point test) is subjected to linear interpolation to obtain Um,t,x，Um,t,yFinally obtain Ut；
Ut＝[U1,t…U1,m,t]
Establishing a finite element model of the thin-wall workpiece, respectively establishing a set of three pieces of additional mass (the size is 20 × 20 × 4mm) of the removed material (100mm × 40mm × 0.5.5 mm) and the final workpiece (100mm × 40mm × 4mm 354 mm) in the process of establishing the model, giving corresponding material attributes to the model, adding constraints and loads which are consistent with actual processing to the workpiece, performing finite element analysis to obtain the inherent frequency of the whole workpiece, extracting the mass matrix and the rigidity matrix of each unit, and finally assembling to obtain the mass matrix and the rigidity matrix of the whole workpiece model and the additional mass, thereby obtaining the inherent frequency omega of each order of the whole workpiece modelwpSum mode vibration type Uwp；
Fourthly, assuming that the damping ratio is unchanged, extracting the damping characteristic of the workpiece by using a hammering mode experiment; determining a damping ratio matrix zeta of a workpiece by measuring the frequency response function of the workpiece at different points and fitting itwp；
Fifthly, using the modal parameters obtained by the test in the second step and the natural frequency omega of each order of the whole workpiece model obtained in the third and fourth stepswpNormal mode vibration type UwpAnd damping ratio matrix ζ of workpiecewpSubstituting the obtained result into the step I, respectively solving the state equations of the tool at the milling initial position and the milling final position by using a semi-discrete method to obtain the axial cutting depth ap(mm) stability lobe plot of milling initial position (solid line 1) and end position (dashed line 2) with spindle speed n (rpm) as variables;
sixthly, taking the maximum material removal rate MRR as an objective function, wherein MRR is ap×ae×n×N×f；apIndicating axial cutting depth, aeThe radial cutting depth is represented, N represents the rotating speed of a main shaft, N represents the number of cutter teeth, and f represents the feed amount of each tooth; radial cutting depth a during millingeThe cutter tooth number N is 2, and the feed amount f of each tooth is 0.1 mm/revolution/tooth; these parameters have been determined beforehand, with axial cutting depth apAnd a variable of the rotation speed n of the main shaft, respectively taking the stability lobe graphs of the milling initial position and the milling final position obtained in the step five as constraints, and setting the range of the axial cutting depth to be 7000 a or more according to actual processing parameterspLess than or equal to 12000; setting the population size to be 20, the random seed generation probability to be 0.12221, the mutation probability to be 0.7, the cross probability to be 0.8 and the genetic algebra parameter to be 30, and respectively obtaining the optimized axial cutting depths a of the initial milling position and the final milling position by utilizing a genetic algorithmpAnd a spindle speed n; selecting a corresponding to the initial position and the end positionpThe smallest of which is the selected processing parameter;
starting position: a isp2.833 n-11200 end position: a isp＝2.763 n＝11400
The parameters finally selected are as follows: a isp＝2.763 n＝11400 MRR＝3149.82
Seventhly, adding mass blocks with different masses on the thin plate to locally modify the mass matrix and the rigidity matrix of the workpiece; according to the material parameters of several common materials given in the table I, the method L of the Taguo method is utilized16(44) Modifying the material density rho and the Young modulus E in the material attribute through an orthogonal test to change the mass matrix and the rigidity matrix of three additional masses, and then assembling to obtain the dynamic parameters of the modified workpiece; substituting the new kinetic parameters into the fifth step to obtain a stability lobe graph corresponding to the orthogonal test; through the sixth step, the product is obtained according to the 'Tiankou method' L16(44) Table three data corresponding to the orthogonal experiment; finally, a mass matrix and a rigidity matrix of the three additional masses under different material properties, corresponding MRRs and corresponding axial directions thereof are obtainedCutting depth apAnd a spindle speed n; k1,K2,K3,K4The influence of various factors on the indexes of the horizontal test is shown, and the optimal additional mass block pasting mode can be obtained by 423 combination, namely the mass block 1 at the starting position is made of copper, the mass block 2 at the middle position is made of aluminum alloy, and the mass block 3 at the ending position is made of 45#Steel; the processing parameters obtained were: axial cutting depth ap5.098 and the main shaft speed n is 10950, the material removal rate MRR is 5582.31; compared with the processing parameters and the material removal rate obtained by not sticking the mass block, the material removal efficiency is improved by 77.2 percent; the method is proved to achieve good expected effect and have good practicability.
Table three: use of Tiankou L16(44) Data from orthogonal experiments
Claims (1)
1. A milling chatter suppression method for a thin-wall part is characterized by comprising the following steps:
firstly, establishing a milling dynamic model which is in multipoint contact and considers the deformation of a cutter and a workpiece; the motion equation of the milling system is as follows:
wherein,t(t) is a vector representing the modal displacement of the tool;wp(t) is a vector representing modal displacement of the workpiece; zetatDiagonal matrix, ζ, representing the damping ratio of the toolwpA diagonal matrix representing a damping ratio of the workpiece; omegatA diagonal matrix representing the natural frequency of the tool; omegawpA diagonal matrix representing the natural frequency of the workpiece; u shapetRepresenting the modal shape, U, of the tool after mass normalizationwpRepresenting the modal shape of the workpiece after mass normalization; q represents the number of nodes for differentiating the tool from the workpiece in the axial direction; f (t) represents a matrix of milling forces at each contact point;
step two, after the selected milling cutter is installed on a main shaft of a machine tool, measuring the feeding direction of the cutter and modal parameters vertical to the feeding direction by adopting a multi-point tapping test method and a linear interpolation method; zetam,t,x，ζm,t,yRepresents the damping ratio of the mth order; omegam,t,x，ωm,t,yA system natural frequency representing the mth order; u shapem,t,x，Um,t,yIs a matrix of dimension q × 1, representing the mode shape of the mth order of the system mass normalization;
step three, establishing a finite element model of the thin-wall workpiece; in the process of establishing the model, the removing material, the final workpiece and the attaching material are respectively establishedA set of additive weights; giving corresponding material attributes to the model, adding constraints and loads which are consistent with actual processing to the workpiece, performing finite element analysis to obtain the natural frequency of the whole workpiece, extracting the mass matrix and the rigidity matrix of each unit, finally assembling to obtain the whole workpiece model and the mass matrix and the rigidity matrix of the additional mass, and further obtaining the natural frequency omega of each order of the whole workpiece modelwpSum mode vibration type Uwp；
Step four, assuming that the damping ratio is unchanged, extracting the damping characteristic of the workpiece by using a hammering mode experiment; determining a damping ratio matrix zeta of a workpiece by measuring the frequency response function of the workpiece at different points and fitting itwp；
Step five, using the modal parameters obtained by the test in the step two and the natural frequency omega of each order of the whole workpiece model obtained in the step three and the step fourwpNormal mode vibration type UwpAnd damping ratio matrix ζ of workpiecewpSubstituting the obtained result into the step I, respectively solving the state equations of the tool at the milling initial position and the milling final position by using a semi-discrete method to obtain the axial cutting depth apAnd a stability lobe graph of the milling initial position and the milling final position with the main shaft rotating speed n as a variable;
step six, the maximum material removal rate MRR is an objective function, and the MRR is ap×ae×n×N×f；apIndicating axial cutting depth, aeThe radial cutting depth is represented, N represents the rotating speed of a main shaft, N represents the number of cutter teeth, and f represents the feed amount of each tooth; radial cutting depth a during millingeThe number N of cutter teeth and the feed amount f of each tooth are determined in advance; by axial cutting-in of apAnd a main shaft rotating speed n variable, setting population size, random seed generation probability, mutation probability, cross probability and genetic algebra parameters by respectively taking the stability lobe graphs of the milling initial position and the milling end position obtained in the step five as constraints, and respectively obtaining optimized axial cutting depths a of the milling initial position and the milling end position by utilizing a genetic algorithmpAnd a spindle speed n; selecting a corresponding to the initial position and the end positionpThe smallest of which is the selected processing parameter;
step seven, modifying the material density rho and the Young's modulus in the material propertyChanging the mass matrix and the rigidity matrix of the additional mass by the quantity E, modifying the material density rho and the Young modulus E in the material attribute by utilizing a Taguchi method orthogonal test to change the mass matrix and the rigidity matrix of the additional mass, and then assembling to obtain the modified dynamic parameters of the workpiece; substituting the new kinetic parameters into the fifth step to obtain a stability lobe graph corresponding to the orthogonal test; through the sixth step, the MRR corresponding to the Taguchi method orthogonal experiment and the axial cutting depth a corresponding to the MRR are obtainedpAnd a spindle speed n; finally, analyzing the data obtained by the Taguchi test, and selecting the optimal additional mass combination which can maximize the MRR and the corresponding axial cutting depth apAnd a spindle speed n.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710171598.4A CN106965032B (en) | 2017-03-22 | 2017-03-22 | Thin-wall part milling parameter suppressing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710171598.4A CN106965032B (en) | 2017-03-22 | 2017-03-22 | Thin-wall part milling parameter suppressing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106965032A true CN106965032A (en) | 2017-07-21 |
| CN106965032B CN106965032B (en) | 2018-01-16 |
Family
ID=59329526
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710171598.4A Active CN106965032B (en) | 2017-03-22 | 2017-03-22 | Thin-wall part milling parameter suppressing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106965032B (en) |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107479499A (en) * | 2017-09-29 | 2017-12-15 | 西北工业大学 | Thin-wall part cuts interface mismachining tolerance compensating Modeling and penalty coefficient learning control method |
| CN107807526A (en) * | 2017-10-31 | 2018-03-16 | 上海交通大学 | A kind of method for intelligently suppressing processing flutter based on Simulation of stability |
| CN107914183A (en) * | 2017-12-05 | 2018-04-17 | 上海复合材料科技有限公司 | The flutter stability Forecasting Methodology of milling carbon fiber layer plywood |
| CN108595874A (en) * | 2018-05-07 | 2018-09-28 | 西南石油大学 | A kind of large thin-wall slewing parts Milling Process vibration mechanical model construction method |
| CN108732995A (en) * | 2018-04-25 | 2018-11-02 | 西北工业大学 | The fast acquiring method of milling process workpiece kinetic parameter |
| CN108846242A (en) * | 2018-07-09 | 2018-11-20 | 西北工业大学 | The thin-wall part milling parameter suppressing method applied based on pre-tensile stress |
| CN109093130A (en) * | 2018-08-08 | 2018-12-28 | 西北工业大学 | Annular thin-wall parts roughing technique optimization method based on subregion processing |
| CN110153781A (en) * | 2019-05-24 | 2019-08-23 | 西北工业大学 | Thin-wall part processing equipment for inhibiting of vibration and method based on flexure type actuator |
| CN110245402A (en) * | 2019-05-31 | 2019-09-17 | 西安交通大学 | A kind of structural member machined parameters calculation method based on Combined machining feature |
| CN110516340A (en) * | 2019-08-21 | 2019-11-29 | 西北工业大学 | Process based on milling system stiffness variation damps unified model modeling method |
| CN113059221A (en) * | 2021-03-10 | 2021-07-02 | 华中科技大学 | Method and device for improving machining stability of thin-wall part based on mirror image magnetic adsorption of tool nose |
| CN113894333A (en) * | 2021-09-26 | 2022-01-07 | 西北工业大学 | Titanium alloy thin-wall structure precision milling surface state robustness process control method |
| CN114895566A (en) * | 2022-05-24 | 2022-08-12 | 西北工业大学 | Chatter prediction method for milling process by adopting transfer matrix reduction technology |
| CN117103280A (en) * | 2023-10-19 | 2023-11-24 | 中国长江电力股份有限公司 | Material reduction processing method and system for large-sized water turbine top cover on-site robot |
| CN118123097A (en) * | 2024-05-08 | 2024-06-04 | 南京航空航天大学 | Thin-wall part milling flutter suppression method based on magneto-rheological damper |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS57121410A (en) * | 1981-01-10 | 1982-07-28 | Mitsubishi Heavy Ind Ltd | Processing of thin wall cylinder |
| CN1515382A (en) * | 2001-12-19 | 2004-07-28 | 北京工业大学 | Machine cutting flutter on-line intelligent control system |
| EP1716969A1 (en) * | 2005-04-28 | 2006-11-02 | Okuma Corporation | A machine tool with an active vibration suppression means |
| JP2011067922A (en) * | 2009-09-28 | 2011-04-07 | Okuma Corp | Vibration control tool and turning method for thin wall cylindrical workpiece in which vibration control tool is inserted |
| CN105935795A (en) * | 2016-06-12 | 2016-09-14 | 哈尔滨理工大学 | Adjustable damping type vibration-attenuation and chatter-suppression device in milling process of thin-walled workpieces and vibration-attenuation method |
-
2017
- 2017-03-22 CN CN201710171598.4A patent/CN106965032B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS57121410A (en) * | 1981-01-10 | 1982-07-28 | Mitsubishi Heavy Ind Ltd | Processing of thin wall cylinder |
| CN1515382A (en) * | 2001-12-19 | 2004-07-28 | 北京工业大学 | Machine cutting flutter on-line intelligent control system |
| EP1716969A1 (en) * | 2005-04-28 | 2006-11-02 | Okuma Corporation | A machine tool with an active vibration suppression means |
| JP2011067922A (en) * | 2009-09-28 | 2011-04-07 | Okuma Corp | Vibration control tool and turning method for thin wall cylindrical workpiece in which vibration control tool is inserted |
| CN105935795A (en) * | 2016-06-12 | 2016-09-14 | 哈尔滨理工大学 | Adjustable damping type vibration-attenuation and chatter-suppression device in milling process of thin-walled workpieces and vibration-attenuation method |
Cited By (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107479499A (en) * | 2017-09-29 | 2017-12-15 | 西北工业大学 | Thin-wall part cuts interface mismachining tolerance compensating Modeling and penalty coefficient learning control method |
| CN107479499B (en) * | 2017-09-29 | 2019-07-19 | 西北工业大学 | Thin-wall part cuts interface mismachining tolerance compensating Modeling and penalty coefficient learning control method |
| CN107807526A (en) * | 2017-10-31 | 2018-03-16 | 上海交通大学 | A kind of method for intelligently suppressing processing flutter based on Simulation of stability |
| CN107914183A (en) * | 2017-12-05 | 2018-04-17 | 上海复合材料科技有限公司 | The flutter stability Forecasting Methodology of milling carbon fiber layer plywood |
| CN107914183B (en) * | 2017-12-05 | 2020-04-07 | 上海复合材料科技有限公司 | Chatter stability prediction method for milling carbon fiber laminated plate |
| CN108732995B (en) * | 2018-04-25 | 2019-05-31 | 西北工业大学 | The fast acquiring method of milling process workpiece kinetic parameter |
| CN108732995A (en) * | 2018-04-25 | 2018-11-02 | 西北工业大学 | The fast acquiring method of milling process workpiece kinetic parameter |
| CN108595874A (en) * | 2018-05-07 | 2018-09-28 | 西南石油大学 | A kind of large thin-wall slewing parts Milling Process vibration mechanical model construction method |
| CN108846242B (en) * | 2018-07-09 | 2019-05-07 | 西北工业大学 | The thin-wall part milling parameter suppressing method applied based on pre-tensile stress |
| CN108846242A (en) * | 2018-07-09 | 2018-11-20 | 西北工业大学 | The thin-wall part milling parameter suppressing method applied based on pre-tensile stress |
| CN109093130A (en) * | 2018-08-08 | 2018-12-28 | 西北工业大学 | Annular thin-wall parts roughing technique optimization method based on subregion processing |
| CN110153781A (en) * | 2019-05-24 | 2019-08-23 | 西北工业大学 | Thin-wall part processing equipment for inhibiting of vibration and method based on flexure type actuator |
| CN110153781B (en) * | 2019-05-24 | 2020-11-20 | 西北工业大学 | Thin-wall part machining vibration suppression device and method based on bending actuator |
| CN110245402A (en) * | 2019-05-31 | 2019-09-17 | 西安交通大学 | A kind of structural member machined parameters calculation method based on Combined machining feature |
| CN110516340A (en) * | 2019-08-21 | 2019-11-29 | 西北工业大学 | Process based on milling system stiffness variation damps unified model modeling method |
| CN110516340B (en) * | 2019-08-21 | 2022-04-12 | 西北工业大学 | Process damping unified model modeling method based on milling system rigidity change |
| CN113059221A (en) * | 2021-03-10 | 2021-07-02 | 华中科技大学 | Method and device for improving machining stability of thin-wall part based on mirror image magnetic adsorption of tool nose |
| CN113894333A (en) * | 2021-09-26 | 2022-01-07 | 西北工业大学 | Titanium alloy thin-wall structure precision milling surface state robustness process control method |
| CN113894333B (en) * | 2021-09-26 | 2022-07-22 | 西北工业大学 | Titanium alloy thin-wall structure precision milling surface state robustness process control method |
| CN114895566A (en) * | 2022-05-24 | 2022-08-12 | 西北工业大学 | Chatter prediction method for milling process by adopting transfer matrix reduction technology |
| CN114895566B (en) * | 2022-05-24 | 2023-11-10 | 西北工业大学 | Chatter prediction method for milling process by adopting transmission matrix reduction technology |
| CN117103280A (en) * | 2023-10-19 | 2023-11-24 | 中国长江电力股份有限公司 | Material reduction processing method and system for large-sized water turbine top cover on-site robot |
| CN117103280B (en) * | 2023-10-19 | 2023-12-22 | 中国长江电力股份有限公司 | Material reduction processing method and system for large-sized water turbine top cover on-site robot |
| CN118123097A (en) * | 2024-05-08 | 2024-06-04 | 南京航空航天大学 | Thin-wall part milling flutter suppression method based on magneto-rheological damper |
| CN118123097B (en) * | 2024-05-08 | 2024-09-20 | 南京航空航天大学 | Thin-wall part milling flutter suppression method based on magneto-rheological damper |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106965032B (en) | 2018-01-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106965032A (en) | Thin-wall part milling parameter suppressing method | |
| CN107914183B (en) | Chatter stability prediction method for milling carbon fiber laminated plate | |
| Budak | Analytical models for high performance milling. Part II: process dynamics and stability | |
| Budak et al. | Maximizing chatter free material removal rate in milling through optimal selection of axial and radial depth of cut pairs | |
| Huang et al. | Milling force vibration analysis in high-speed-milling titanium alloy using variable pitch angle mill | |
| CN106802630B (en) | The Parallel frequency domain method of milling parameter stability prediction based on thin-wall part | |
| Kishawy et al. | An energy based analysis of broaching operation: Cutting forces and resultant surface integrity | |
| CN104484506A (en) | Turning flutter forecasting method based on reliability lobe graph | |
| CN113168491B (en) | Method for simulating surface appearance of flutter-free milling | |
| CN107423502B (en) | Milling optimization method and optimization system based on milling cutter's helix angle and curvature effect | |
| CN106940746A (en) | The parallel time domain method of milling parameter stability prediction based on thin-wall part | |
| Li et al. | Development of a hybrid cutting force model for micromilling of brass | |
| CN106126778A (en) | Thin-wall part week milling stability prediction method with curved surface | |
| CN110102787B (en) | Amplitude modulation-based variable spindle rotating speed turning chatter suppression method | |
| CN104298799B (en) | Slotting cutter key geometric parameter design method based on processing vibration | |
| Kukiełka | Ecological aspects of the implementation of new technologies processing for machinery parts | |
| Markopoulos et al. | Molecular dynamics modeling of a single diamond abrasive grain in grinding | |
| CN114509991A (en) | Numerical control machine tool cutting stability prediction and optimization method considering parameter uncertainty | |
| CN108846242B (en) | The thin-wall part milling parameter suppressing method applied based on pre-tensile stress | |
| CN112643059B (en) | Design method of turning large-pitch thread vibration reduction tool | |
| No et al. | Scanning and modeling for non-standard edge geometry endmills | |
| CN112016203B (en) | Method for predicting milling stability based on segmented Hermite interpolation polynomial and integral discrete strategy | |
| CN103823787A (en) | Multi-turning-tool parallel turning stability judgment method based on differential quadrature method | |
| CN108520117A (en) | A method of obtaining the stability lobes diagram using approximate shceme method | |
| CN109048466B (en) | Milling flutter suppression method based on multi-frequency variable rotation speed |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |