CN104680000A - Orthogonal polynomial-based milling stability prediction method - Google Patents

Orthogonal polynomial-based milling stability prediction method Download PDF

Info

Publication number
CN104680000A
CN104680000A CN201510068454.7A CN201510068454A CN104680000A CN 104680000 A CN104680000 A CN 104680000A CN 201510068454 A CN201510068454 A CN 201510068454A CN 104680000 A CN104680000 A CN 104680000A
Authority
CN
China
Prior art keywords
slideway
formula
represent
circuit
orthogonal polynomial
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
Application number
CN201510068454.7A
Other languages
Chinese (zh)
Other versions
CN104680000B (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201510068454.7A priority Critical patent/CN104680000B/en
Publication of CN104680000A publication Critical patent/CN104680000A/en
Application granted granted Critical
Publication of CN104680000B publication Critical patent/CN104680000B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention relates to the field of advanced manufacturing, in particular to an orthogonal polynomial-based milling stability prediction method. The orthogonal polynomial-based milling stability prediction method adopts orthogonal polynomials to approach a state term, a delay term and a periodical coefficient term in a kinetic equation, and adopts a plurality of known time points and responses thereof to fit the needed terms, so that the local errors of the calculation method are reduced, and thereby the precision of the prediction method is increased; meanwhile, in the process of obtaining a stability lobe diagram, an H matrix is introduced instead of directly substituting an F matrix for calculation, so that the number of iterations in the F matrix calculation process is reduced, consequently, the time of the calculation method is saved, and the efficiency of calculation is increased.

Description

A kind of milling stability Forecasting Methodology based on orthogonal polynomial
Technical field
The present invention relates to advanced manufacture field, be specifically related to a kind of milling stability Forecasting Methodology based on orthogonal polynomial.
Background technology
Along with China constantly incorporates globalizing market, the competitive pressure that manufacturing industry faces is increasing, High-speed Machining Technology is arisen at the historic moment, and be widely used in the complex parts manufacture in the fields such as Aeronautics and Astronautics, boats and ships, mould and automobile, flutter in working angles has a strong impact on the surface quality of workpiece, reduce lathe and cutting-tool's used life, and likely cause the destruction of lathe.The flutter stability of milling process is predicted, selects suitable processing conditions, the generation of flutter can be controlled, improve working (machining) efficiency, cut down finished cost.If do not carry out stability prediction, need to test many times in actual process, just can determine optimum machined parameters, cause a large amount of wastes of the wearing and tearing of cutter and workpiece material, greatly increase manufacturing cost.Therefore carry out prediction to the flutter stability in milling process to have great importance.
A kind of milling parameter stability prediction method based on semi-discrete method is disclosed in " Updated semi-discretization method for periodic delay-differential equations with discrete delay " literary composition on famous International Periodicals " International Journal For Numerical Methods In Engineering " magazine the 61st volume in 2004, use the method can obtain the flutter stability limit, but weak point is, the counting yield of the method is lower, obtain the stability lobes diagram and need long time, and the computational accuracy of this algorithm awaits promoting.
Summary of the invention
The milling stability Forecasting Methodology based on orthogonal polynomial that technical matters to be solved by this invention is to provide a kind of Computationally efficient, high computational accuracy is applicable to multiple parts.
For solving the problems of the technologies described above, the present invention adopts following technical scheme: the present invention includes following steps:
1. the kinetics equation in milling cutter single-degree-of-freedom milling process is set up:
(1)
Wherein, for constant coefficient matrix, for the matrix of coefficients of mechanical periodicity in time, represent that cutter exists the condition responsive in moment, represent the natural frequency of point of a knife point, represent relative damping, represent modal mass, represent axial cutting depth, represent time lag;
represent momentary-cutting-thick, its expression formula is:
(2) in formula (2), represent the cutter tooth number of milling cutter, with be respectively Cutting Force Coefficient that is tangential and normal direction, be the angular displacement of cutter tooth, expression formula is , window function definition is:
(3)
In formula (3), with be respectively the incision of cutter tooth and cut out angle, when adopting climb cutting, ; When adopting upmilling, , for radial direction immerses ratio, the i.e. ratio of radial cutting-in/tool diameter;
2. by the time lag item of the milling process kinetics equation (1) of single-degree-of-freedom be equally divided into individual minizone, then time step is , wherein any one time minizone is expressed as
By equation (1) in time minizone on carry out integration, obtain
(4)
3. build by structure the status items that orthogonal polynomial carrys out fit procedure 2. Chinese style (4) , hangover state item several with time dependent periodic system , detailed process is as follows:
Orthogonal polynomial race coefficient by (6) determine;
Formula (6) be with for power is about point set orthogonal polynomial, wherein
(7)
(8)
By step 2. in the status items of equation (4) be expressed as by orthonormal polynomial approximation:
(9)
Wherein it is polynomial expression coefficient, pass through (10) determine, wherein be the condition responsive in moment;
Get in (6) formula pass through time point with its response come fit procedure 2. in status items,
(7) formula of utilization calculates , and substitute into can obtain:
(11)
So
(12)
Obtained by (10) formula
(13)
(14)
(12), (13), (14) formula are substituted into (10) formula and obtain
(15)
Arrange
(16)
In like manner time lag item can be expressed as:
(17)
(18);
4. build Floquet transition matrix, (16), (17), (18) formula substituted into (4) formula, can obtain:
(19)
Wherein
(20a)
(20b)
(20c)
(21a)
(21b)
(21c)
(21d)
Equation (19) can be written as
(22)
Wherein
(23)
By equation (22), can obtain current cutter tooth and a upper cutter tooth respond between mapping relations, as follows by matrix representation:
(24)
Wherein
(25a)
(25b)
(25c)
The Floquet transition matrix of system can be expressed as
(26)
Wherein
(27)
5. Floquet transition matrix is calculated eigenwert, by the stability of the mould decision-making system of eigenwert, concrete decision criteria is as follows: .
Step of the present invention 3. in also can adopt status items based on building Quadratic Orthogonal polynomial method and come fit procedure 2. Chinese style (4) , hangover state item several with time dependent periodic system ; Detailed process is as follows:
In modus ponens (6) with , pass through time point with the response of its correspondence carry out matching status items , wherein
(29a)
(29b)
(29c)
By (7) Shi Ke get
(30)
So
(31)
(32)
(33)
Order , and substitute into can obtain (30), (31), (32), (33) formula obtain
(34)
(35)
(36)
(37)
(34), (35), (36), (37) formula are substituted into (29c) and obtain
(38)
(29a), (36) and (38) formula are substituted into (10) formula and obtain
So
Arrange
The present invention in order to accurately obtain step 1. in modal parameter, need bracing or strutting arrangement be adopted;
Described bracing or strutting arrangement comprises brace table, the first circuit, the second circuit, X to slideway, the cooperation button being provided with through hole, the first slideway, the second slideway, the 3rd slideway, the 4th slideway and the supporting mechanism for fixed support thin-wall part;
Described first circuit and the second circuit are with one heart and be successively set on brace table, described X also passes through the center of circle of the first circuit and the second circuit on brace table to installation of slide track, described first slideway and the second slideway are symmetricly set on X to slideway both sides, described 3rd slideway and the 4th slideway are symmetricly set on X to slideway both sides, and described cooperation button is located at the intersection of the first circuit, the second circuit and the first slideway, the second slideway, the 3rd slideway, the 4th slideway;
Described supporting mechanism comprises train wheel bridge, lower plate, spring and pole, and described lower plate is arranged on pole, and described train wheel bridge is connected with train wheel bridge by spring, and described pole lower end is provided with shoulder, and described pole is connect with coordinating to link by its lower end shoulder.
Good effect of the present invention is as follows: the present invention adopts orthogonal polynomial several to status items, time lag item and the periodic system of approaching in kinetics equation, multiple known time point and response thereof is adopted to carry out matching the required entries, reduce the local error of computing method, thus improve the precision of Forecasting Methodology; Simultaneously in the process obtaining the stability lobes diagram, introduce H matrix, instead of directly substitution F matrix calculates, and reduces the iterations in F matrix computation process, thus save the time of computing method, improve counting yield.
Accompanying drawing explanation
Fig. 1 is that the present invention's radial direction immerses than stability diagram when being 0.05;
Fig. 2 is that the present invention's radial direction immerses than stability diagram when being 0.5;
Fig. 3 is that the present invention's radial direction immerses than stability diagram when being 1;
Fig. 4 is support table structure schematic diagram of the present invention;
Fig. 5 is lower plate structural representation of the present invention;
In the drawings: 1 first circuit, 2 second circuits, 3 brace tables, 4 coordinate button, 5X to slideway, 6 through holes, 7-1 first slideway, 7-2 second slideway, 8-1 the 3rd slideway, 8-2 the 4th slideway, 9-1 train wheel bridge, 9-2 lower plate, 10 springs, 11 poles.
Embodiment
The present invention includes following steps:
1. the kinetics equation in milling cutter single-degree-of-freedom milling process is set up:
(1)
Wherein, for constant coefficient matrix, for the matrix of coefficients of mechanical periodicity in time, represent that cutter exists the condition responsive in moment, represent the natural frequency of point of a knife point, represent relative damping, represent modal mass, represent axial cutting depth, represent time lag;
represent momentary-cutting-thick, its expression formula is:
(2) in formula (2), represent the cutter tooth number of milling cutter, with be respectively Cutting Force Coefficient that is tangential and normal direction, be the angular displacement of cutter tooth, expression formula is , window function definition is:
(3)
In formula (3), with be respectively the incision of cutter tooth and cut out angle, when adopting climb cutting, ; When adopting upmilling, , for radial direction immerses ratio, the i.e. ratio of radial cutting-in/tool diameter;
2. by the time lag item of the milling process kinetics equation (1) of single-degree-of-freedom be equally divided into individual minizone, then time step is , wherein any one time minizone is expressed as
By equation (1) in time minizone on carry out integration, obtain
(4)
3. build by structure the status items that orthogonal polynomial carrys out fit procedure 2. Chinese style (4) , hangover state item several with time dependent periodic system , detailed process is as follows:
Orthogonal polynomial race coefficient by (6) determine;
Formula (6) be with for power is about point set orthogonal polynomial, wherein
(7)
(8)
By step 2. in the status items of equation (4) be expressed as by orthonormal polynomial approximation:
(9)
Wherein it is polynomial expression coefficient, pass through (10) determine, wherein be the condition responsive in moment;
Get in (6) formula pass through time point with its response come fit procedure 2. in status items,
(7) formula of utilization calculates , and substitute into can obtain:
(11)
So
(12)
Obtained by (10) formula
(13)
(14)
(12), (13), (14) formula are substituted into (10) formula and obtain
(15)
Arrange
(16)
In like manner time lag item can be expressed as:
(17)
(18);
4. build Floquet transition matrix, (16), (17), (18) formula substituted into (4) formula, can obtain:
(19)
Wherein
(20a)
(20b)
(20c)
(21a)
(21b)
(21c)
(21d)
Equation (19) can be written as
(22)
Wherein
(23)
By equation (22), can obtain current cutter tooth and a upper cutter tooth respond between mapping relations, as follows by matrix representation:
(24)
Wherein
(25a)
(25b)
(25c)
The Floquet transition matrix of system can be expressed as
(26)
Wherein
(27)
5. Floquet transition matrix is calculated eigenwert, by the stability of the mould decision-making system of eigenwert, concrete decision criteria is as follows: .
the present inventiondescribed step 3. in also can adopt status items based on building Quadratic Orthogonal polynomial method and come fit procedure 2. Chinese style (4) , hangover state item several with time dependent periodic system ; Detailed process is as follows:
In modus ponens (6) with , pass through time point with the response of its correspondence carry out matching status items , wherein
(29a)
(29b)
(29c)
By (7) Shi Ke get
(30)
So
(31)
(32)
(33)
Order , and substitute into can obtain (30), (31), (32), (33) formula obtain
(34)
(35)
(36)
(37)
(34), (35), (36), (37) formula are substituted into (29c) and obtain
(38)
(29a), (36) and (38) formula are substituted into (10) formula and obtain
So
Arrange
The present invention in order to accurately obtain step 1. in modal parameter, bracing or strutting arrangement need be adopted, as shown in Figure 4,5, described bracing or strutting arrangement comprises brace table 3, first circuit 1, second circuit 2, X to slideway 5, cooperation button 4, the first slideway 7-1, the second slideway 7-2 that are provided with through hole 6, the 3rd slideway 8-1, the 4th slideway 8-2 and the supporting mechanism for fixed support thin-wall part;
Described first circuit 1 and the second circuit 2 are with one heart and be successively set on brace table 3, described X to be arranged on brace table 3 and by the center of circle of the first circuit 1 and the second circuit 2 to slideway 5, described first slideway 7-1 and the second slideway 7-2 is symmetricly set on X to slideway 5 both sides, described 3rd slideway 8-1 and the 4th slideway 8-2 is symmetricly set on X to slideway 5 both sides, and described cooperation button 4 is located at the intersection of the first circuit 1, second circuit 2 and the first slideway 7-1, the second slideway 7-2, the 3rd slideway 8-1, the 4th slideway 8-2;
Described supporting mechanism comprises train wheel bridge 9-1, lower plate 9-2, spring 10 and pole 11, described lower plate 9-2 is arranged on pole 11, described train wheel bridge 9-1 is connected with train wheel bridge 9-1 by spring 10, described pole 11 lower end is provided with shoulder, and described pole 11 is connected with cooperation button 4 by its lower end shoulder.
The present invention coordinates the support attitude that can complete different spaces shape on brace table 3 by coordinating button 4 with circuit, slideway, effectively can complete and the support of thin-wall part is fixed, utilize bracing or strutting arrangement to carry out clamping and can obtain more accurate modal parameter thus make by utilizing the limit of stability of orthogonal polynomial prediction more accurate.
When given relevant parameter: milling cutter tooth number is 2, and tangential force coefficient and normal force coefficient are respectively with , the first natural frequency of point of a knife point is , modal damping is 0.011, and modal mass is 0.03993kg, climb cutting.By time lag be divided into 40 minizones, the plane be made up of the speed of mainshaft and the radial cutting degree of depth is divided into grid.
Above-mentioned steps and parameter are carried out programming by Matlab software and draws stability Lobe figure, the stability in milling process is predicted by stability diagram, choosing different radial directions and immerse ratio, obtaining stability diagram for not getting 0.05,0.5,1, as shown in Figure 1, 2, 3.
The above embodiment is only the preferred embodiments of the present invention, and and the feasible enforcement of non-invention exhaustive.For persons skilled in the art, to any apparent change done by it under the prerequisite not deviating from the principle of the invention and spirit, all should be contemplated as falling with within claims of the present invention.

Claims (3)

1., based on a milling stability Forecasting Methodology for orthogonal polynomial, it is characterized in that it comprises the following steps:
1. the kinetics equation in milling cutter single-degree-of-freedom milling process is set up:
(1)
Wherein, for constant coefficient matrix, for the matrix of coefficients of mechanical periodicity in time, represent that cutter exists the condition responsive in moment, represent the natural frequency of point of a knife point, represent relative damping, represent modal mass, represent axial cutting depth, represent time lag;
represent momentary-cutting-thick, its expression formula is:
(2) in formula (2), represent the cutter tooth number of milling cutter, with be respectively Cutting Force Coefficient that is tangential and normal direction, be the angular displacement of cutter tooth, expression formula is , window function definition is:
(3)
In formula (3), with be respectively the incision of cutter tooth and cut out angle, when adopting climb cutting, ; When adopting upmilling, , for radial direction immerses ratio, the i.e. ratio of radial cutting-in/tool diameter;
2. by the time lag item of the milling process kinetics equation (1) of single-degree-of-freedom be equally divided into individual minizone, then time step is , wherein any one time minizone is expressed as
By equation (1) in time minizone on carry out integration, obtain
(4)
3. build by structure the status items that orthogonal polynomial carrys out fit procedure 2. Chinese style (4) , hangover state item several with time dependent periodic system , detailed process is as follows:
Orthogonal polynomial race coefficient by (6) determine;
Formula (6) be with for power is about point set orthogonal polynomial, wherein
(7)
(8)
By step 2. in the status items of equation (4) be expressed as by orthonormal polynomial approximation:
(9)
Wherein it is polynomial expression coefficient, pass through (10) determine, wherein be the condition responsive in moment;
Get in (6) formula pass through time point with its response come fit procedure 2. in status items,
(7) formula of utilization calculates , and substitute into can obtain:
(11)
So
(12)
Obtained by (10) formula
(13)
(14)
(12), (13), (14) formula are substituted into (10) formula and obtain
(15)
Arrange
(16)
In like manner time lag item can be expressed as:
(17)
(18);
4. build Floquet transition matrix, (16), (17), (18) formula substituted into (4) formula, can obtain:
(19)
Wherein
(20a)
(20b)
(20c)
(21a)
(21b)
(21c)
(21d)
Equation (19) can be written as
(22)
Wherein
(23)
By equation (22), can obtain current cutter tooth and a upper cutter tooth respond between mapping relations, as follows by matrix representation:
(24)
Wherein
(25a)
(25b)
(25c)
The Floquet transition matrix of system can be expressed as
(26)
Wherein
(27)
5. Floquet transition matrix is calculated eigenwert, by the stability of the mould decision-making system of eigenwert, concrete decision criteria is as follows: .
2. a kind of milling stability Forecasting Methodology based on orthogonal polynomial according to claim 1, is characterized in that: described step 3. in also can adopt status items based on building Quadratic Orthogonal polynomial method and come fit procedure 2. Chinese style (4) , hangover state item several with time dependent periodic system ; Detailed process is as follows:
In modus ponens (6) with , pass through time point with the response of its correspondence carry out matching status items , wherein
(29a)
(29b)
(29c)
By (7) Shi Ke get
(30)
So
(31)
(32)
(33)
Order , and substitute into can obtain (30), (31), (32), (33) formula obtain
(34)
(35)
(36)
(37)
(34), (35), (36), (37) formula are substituted into (29c) and obtain
(38)
(29a), (36) and (38) formula are substituted into (10) formula and obtain
So
Arrange
3. a kind of milling stability Forecasting Methodology based on orthogonal polynomial according to claim 1, is characterized in that: in order to accurately obtain step 1. in modal parameter, need bracing or strutting arrangement be adopted;
Described bracing or strutting arrangement comprises brace table (3), the first circuit (1), the second circuit (2), X to slideway (5), the cooperation button (4) being provided with through hole (6), the first slideway (7-1), the second slideway (7-2), the 3rd slideway (8-1), the 4th slideway (8-2) and the supporting mechanism for fixed support thin-wall part;
Described first circuit (1) and the second circuit (2) are with one heart and be successively set on brace table (3), described X is arranged on brace table (3) to slideway (5) and goes up and the center of circle passing through the first circuit (1) and the second circuit (2), described first slideway (7-1) and the second slideway (7-2) are symmetricly set on X to slideway (5) both sides, described 3rd slideway (8-1) and the 4th slideway (8-2) are symmetricly set on X to slideway (5) both sides, described cooperation button (4) is located at the first circuit (1), second circuit (2) and the first slideway (7-1), second slideway (7-2), 3rd slideway (8-1), the intersection of the 4th slideway (8-2),
Described supporting mechanism comprises train wheel bridge (9-1), lower plate (9-2), spring (10) and pole (11), described lower plate (9-2) is arranged on pole (11), described train wheel bridge (9-1) is connected with train wheel bridge (9-1) by spring (10), described pole (11) lower end is provided with shoulder, and described pole (11) is detained (4) by its lower end shoulder with cooperation and is connected.
CN201510068454.7A 2015-02-10 2015-02-10 A kind of milling stability Forecasting Methodology based on orthogonal polynomial Expired - Fee Related CN104680000B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510068454.7A CN104680000B (en) 2015-02-10 2015-02-10 A kind of milling stability Forecasting Methodology based on orthogonal polynomial

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510068454.7A CN104680000B (en) 2015-02-10 2015-02-10 A kind of milling stability Forecasting Methodology based on orthogonal polynomial

Publications (2)

Publication Number Publication Date
CN104680000A true CN104680000A (en) 2015-06-03
CN104680000B CN104680000B (en) 2018-03-09

Family

ID=53315033

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510068454.7A Expired - Fee Related CN104680000B (en) 2015-02-10 2015-02-10 A kind of milling stability Forecasting Methodology based on orthogonal polynomial

Country Status (1)

Country Link
CN (1) CN104680000B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105184077A (en) * 2015-09-06 2015-12-23 河南师范大学 Excessively-close-range particle-swarm exponential method for optimizing efficiency of resonant electric energy transmitting system
CN106774148A (en) * 2017-01-12 2017-05-31 太原科技大学 A kind of milling stability Forecasting Methodology based on Bull formula
CN110162733A (en) * 2018-02-13 2019-08-23 中国商用飞机有限责任公司 Milling stability analysis method based on whole discrete strategies
CN110750891A (en) * 2019-10-12 2020-02-04 北京理工大学 Parallel synchronous orthogonal turning and milling flutter stability lobe graph prediction method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060271231A1 (en) * 2005-05-26 2006-11-30 Nejat Olgac System and method for chatter stability prediction and control in simultaneous machining applications
CN102873381A (en) * 2012-09-29 2013-01-16 西安交通大学 High-speed milling process parameter optimizing method based on dynamic model
CN103559550A (en) * 2013-09-09 2014-02-05 西北工业大学 Milling stability domain prediction method under multi-modal coupling
CN103846690A (en) * 2012-11-29 2014-06-11 大连凯泓科技有限公司 Movable positioning device
CN203726149U (en) * 2014-03-04 2014-07-23 栗生锐 Flexible auxiliary supporting and clamping device used for large thin-wall annular members

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060271231A1 (en) * 2005-05-26 2006-11-30 Nejat Olgac System and method for chatter stability prediction and control in simultaneous machining applications
CN102873381A (en) * 2012-09-29 2013-01-16 西安交通大学 High-speed milling process parameter optimizing method based on dynamic model
CN103846690A (en) * 2012-11-29 2014-06-11 大连凯泓科技有限公司 Movable positioning device
CN103559550A (en) * 2013-09-09 2014-02-05 西北工业大学 Milling stability domain prediction method under multi-modal coupling
CN203726149U (en) * 2014-03-04 2014-07-23 栗生锐 Flexible auxiliary supporting and clamping device used for large thin-wall annular members

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
QICHANG ZHANG 等: "Stability Prediction of Milling Process with Variable Pitch Cutter", 《MATHEMATICAL PROBLEMS IN ENGINEERING》 *
YE DING 等: "A full-discretization method for prediction of milling stability", 《INTERNATIONAL JOURNAL OF MACHINE TOOLS AND MANUFACTURE》 *
YILONG LIU 等: "An efficient full-discretization method for prediction of milling stability", 《INTERNATIONAL JOURNAL OF MACHINE TOOLS AND MANUFACTURE》 *
丁汉 等: "铣削过程稳定性分析的时域法研究进展", 《科学通报》 *
夏省祥 等: "《常用数值算法及其MATLAB实现》", 30 April 2014, 北京:清华大学出版社 *
李忠群: "复杂切削条件高速铣削加工动力学建模、仿真与切削参数优化研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105184077A (en) * 2015-09-06 2015-12-23 河南师范大学 Excessively-close-range particle-swarm exponential method for optimizing efficiency of resonant electric energy transmitting system
CN105184077B (en) * 2015-09-06 2018-07-31 河南师范大学 Cross short distance low-resonance electric energy transmission system improving efficiency population index method
CN106774148A (en) * 2017-01-12 2017-05-31 太原科技大学 A kind of milling stability Forecasting Methodology based on Bull formula
CN106774148B (en) * 2017-01-12 2018-11-13 太原科技大学 A kind of milling stability prediction technique based on Bull formula
CN110162733A (en) * 2018-02-13 2019-08-23 中国商用飞机有限责任公司 Milling stability analysis method based on whole discrete strategies
CN110162733B (en) * 2018-02-13 2023-06-02 中国商用飞机有限责任公司 Milling stability analysis method based on integral discrete strategy
CN110750891A (en) * 2019-10-12 2020-02-04 北京理工大学 Parallel synchronous orthogonal turning and milling flutter stability lobe graph prediction method
CN110750891B (en) * 2019-10-12 2021-09-21 北京理工大学 Parallel synchronous orthogonal turning and milling flutter stability lobe graph prediction method

Also Published As

Publication number Publication date
CN104680000B (en) 2018-03-09

Similar Documents

Publication Publication Date Title
CN104680000A (en) Orthogonal polynomial-based milling stability prediction method
CN105488282A (en) Cutting parameter segmentation and variable cutting depth optimizing method based on dynamic machining feature
CN105160059A (en) BP and GA based blade machining cutting quantity optimization selection method
CN104657606B (en) A kind of milling stability Forecasting Methodology based on cubic polynomial
CN105425727A (en) Five-axis side milling machining cutter path smoothing method
CN104007697B (en) Five axle multirow Flank machining tool path scheduling methods
CN104898554A (en) Composite cutter path generation method based on discrete cutter location point
CN108628247A (en) Curved surface subregion Toolpath Generation method based on the residual high constraint in boundary
CN107451382B (en) Control method for surface appearance of high-speed cutting workpiece
CN104318022A (en) Method for predicting workpiece surface roughness and increasing cutting efficiency
CN105955195A (en) Milling force prediction-based curved surface processing trajectory generation method
CN104708458A (en) Supporting device and machining precision improving method for thin-wall part
CN104794337A (en) Modeling method for judging cutting force of orthogonal turn-milling machining end face on basis of boundary conditions
CN105414616A (en) Method for cutting force forecasting and stability judging in spiral hole milling process
Sun et al. Chatter stability of orthogonal turn-milling analyzed by complete discretization method
CN105739436A (en) Adaptive rapid interpolation method for numerical control machine tool spline curve
CN103971019A (en) Method for predicting workpiece machining energy consumption based on geometrical characteristics
Scippa et al. Milled surface generation model for chip thickness detection in peripheral milling
CN106126778A (en) Thin-wall part week milling stability prediction method with curved surface
CN104657607A (en) Thin-wall part supporting device and milling stability forecasting method
Wang et al. A model of deformation of thin-wall surface parts during milling machining process
CN106980720A (en) Thin-wall part milling distortion inaccuracy Forecasting Methodology
CN110516373A (en) A kind of method for milling of circular saw roughing special-shaped level
CN112883505B (en) Ultra-precise end face turning surface modeling method considering relative vibration of cutter workpiece
CN104014633B (en) A kind of analogy method of removal bar core defect based on limited element analysis technique and the hole-punching method based on this analogy method

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

Granted publication date: 20180309

Termination date: 20200210

CF01 Termination of patent right due to non-payment of annual fee