CN104647132B - A kind of milling parameter Active Control Method based on magnetic suspension bearing electric chief axis - Google Patents

A kind of milling parameter Active Control Method based on magnetic suspension bearing electric chief axis Download PDF

Info

Publication number
CN104647132B
CN104647132B CN201410805910.7A CN201410805910A CN104647132B CN 104647132 B CN104647132 B CN 104647132B CN 201410805910 A CN201410805910 A CN 201410805910A CN 104647132 B CN104647132 B CN 104647132B
Authority
CN
China
Prior art keywords
amp
displacement
described
matrix
milling
Prior art date
Application number
CN201410805910.7A
Other languages
Chinese (zh)
Other versions
CN104647132A (en
Inventor
陈智勇
吴越
赵欢
张海涛
赵杰
丁汉
赵鑫
张国强
Original Assignee
华中科技大学
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 华中科技大学 filed Critical 华中科技大学
Priority to CN201410805910.7A priority Critical patent/CN104647132B/en
Publication of CN104647132A publication Critical patent/CN104647132A/en
Application granted granted Critical
Publication of CN104647132B publication Critical patent/CN104647132B/en

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23QDETAILS, 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
    • B23Q15/00Automatic control or regulation of feed movement, cutting velocity or position of tool or work
    • B23Q15/007Automatic control or regulation of feed movement, cutting velocity or position of tool or work while the tool acts upon the workpiece
    • B23Q15/12Adaptive control, i.e. adjusting itself to have a performance which is optimum according to a preassigned criterion

Abstract

The invention discloses a kind of milling parameter Active Control Method based on magnetic suspension bearing electric chief axis, comprise the steps: current displacement signal when gathering electro spindle generation milling parameter; Obtain the state vector Q (t) be made up of electro spindle displacement and speed; According to electro spindle displacement [F xf y] twith trigonometric function column vector obtain harmonic displacement obtain the first self adaptation weight coefficient with the second self adaptation weight coefficient according to harmonic displacement h (t), the first self adaptation weight coefficient λ 1with the second self adaptation weight coefficient λ 2obtain the first adaptive rate with the second adaptive rate according to the first adaptive rate second adaptive rate and harmonic displacement h (t) obtains Self Adaptive Control electric current by Self Adaptive Control electric current Q ct () acts on radial magnetic bearing, and produce corresponding magnetic field force, thus realizes suppressing the flutter of main shaft.The present invention eliminates the flutter produced in milling process, and then ensure that crudy and improve working (machining) efficiency.

Description

A kind of milling parameter Active Control Method based on magnetic suspension bearing electric chief axis

Technical field

The invention belongs to Computerized Numerical Control processing technology field, more specifically, relate to a kind of milling parameter Active Control Method based on magnetic suspension bearing electric chief axis.

Background technology

In complex-curved class part five-shaft numerical control processing Flutter Suppression, conventional method is based on Passive Control, these class methods are by regulating the technological parameters such as the speed of mainshaft, tool feeding cycle, maximum cutting-in, cutter length, angle between teeth, make system works in the below of salient angle effect point, thus obtain higher rotating speed and cutting-in under the prerequisite ensureing milling stability.Although passive control methods design and installation is simple, the stable region that it can adjust is very little, limited to the improvement of closed loop power.Therefore the control device changing the intrinsic dynamic characteristic of main axle unit must be adopted to suppress flutter, and this control device is namely based on the active control technology of magnetic suspension bearing electric chief axis.

The thought of magnetic suspension bearing electric chief axis active control technology is by configuring magnetic suspension bearing brake and sensor in spindle motor unit proper site, utilize the effect of FEEDBACK CONTROL, the overall dynamics performance of closed-loop system is changed, thus the closed-loop stabilization region of effective extension system, therefrom find the point that performance is more excellent, improve maximum metal milling rate.Active control technology has the ability of radical change system dynamics behavior, and therefore it suppresses the ability of flutter to be better than Passive Control technology far away.

Through the retrieval to existing document, find that there is following prior art: application number is 201310475810.8, publication No. is CN103769945A, name is called " vibration suppressing method and lathe ", this technology is by changing the speed of mainshaft to reach the effect suppressing flutter, belong to Passive Control, also limited for the ability improving working (machining) efficiency.

Summary of the invention

For the defect of prior art, the object of the present invention is to provide a kind of milling parameter Active Control Method based on magnetic suspension bearing electric chief axis, be intended to solve the workpiece surface quality caused because of self-excitation chatter phenomenon in existing Digit Control Machine Tool Milling Process decline, the problem that milling efficiency is low.

The invention provides a kind of milling parameter Active Control Method based on magnetic suspension bearing electric chief axis, comprise the steps:

(1) when electro spindle generation milling parameter, the current displacement signal of electro spindle is gathered;

(2), after filter and amplification and differential process being carried out to described current displacement signal, the state vector Q (t) be made up of electro spindle displacement and speed is obtained; Described state vector Q (t) is the column vector of 4 dimensions;

(3) electro spindle displacement [F is obtained according to described state vector Q (t) xf y] t=C q[Q (t-T/N)-Q (t)];

Wherein C q = 1 0 0 0 0 1 0 0 , The electro spindle displacement that Q (t-T/N) is the t-T/N moment and the state vector that speed is formed, T represents main shaft swing circle, and N represents the number of teeth of milling cutter, the electro spindle displacement F in x direction xwith the electro spindle displacement F in y direction yare all vectors of 1 dimension, the subscript T of vector represents the transposition to vector;

(4) Fourier transformation is carried out to each element in cutting force dynamic matrix K (t), obtain cutting force harmonic wave dynamic matrix K α(t), and according to cutting force harmonic wave dynamic matrix K αt () obtains harmonic wave

Wherein cutting force harmonic wave dynamic matrix K α(t) for each element in cutting force dynamic matrix K (t) carry out Fourier expansion front 2 α+1 item parts and, ω=2 π N/T is harmonic frequency, and α is the suitable positive integer chosen;

(5) according to described electro spindle displacement [F xf y] twith described trigonometric function column vector obtain harmonic displacement wherein harmonic displacement h (t) is 4 α+2 dimensional vectors;

(6) according to state displacement term coefficient matrix A and Lyapunov Equation Χ A+A tΧ=-I obtains symmetric positive definite matrix Χ, obtains state weight coefficient [Λ according to symmetric positive definite matrix Χ 1Λ 2]=XB, according to state weight coefficient [Λ 1Λ 2] and described state vector Q (t) obtain the first self adaptation weight coefficient with the second self adaptation weight coefficient

Wherein, I is the unit matrix that 4 row 4 arrange, Λ ifour dimensional vectors, λ ithe vector of one dimension, i=1,2;

(7) according to described harmonic displacement h (t), described first self adaptation weight coefficient λ 1with described second self adaptation weight coefficient λ 2obtain the first adaptive rate with the second adaptive rate

Wherein, the first adaptive change rate second adaptive change rate Φ 1and Φ 2be respectively the symmetric positive definite matrix on two 4 α+2 rank of getting arbitrarily, the first adaptive rate with the second adaptive rate be 4 α+2 dimensional vectors;

(8) according to described first adaptive rate described second adaptive rate and described harmonic displacement h (t) obtains Self Adaptive Control electric current wherein Q ct () is 2 dimensional vectors;

(9) by described Self Adaptive Control electric current Q ct () acts on radial magnetic bearing, and produce corresponding magnetic field force, thus realizes suppressing the flutter of main shaft.

Further, described state displacement term coefficient matrix A is: A = 0 0 1 0 0 0 0 1 - ω n 2 0 - 2 ζ ω n 0 0 - ω n 2 0 - 2 ζ ω n , ζ is the damping ratio of milling cutter system architecture, ω nfor the intrinsic frequency of milling cutter system architecture.

Further, described Cutting Force Coefficient matrix B is: B = 0 0 0 0 1 / m 0 0 1 / m , M is the modal mass of milling cutter system architecture.

Further, described cutting force dynamic matrix is: K ( t ) = Σ j = 1 N g j ( φ j ( t ) ) S ( t ) [ K t K r ] T [ sin φ j ( t ) cos φ j ( t ) ] ; Wherein φ jt () represents the anglec of rotation of a jth cutter tooth of milling cutter, K tand K rrepresent Cutting Force Coefficient, g jj(t)) and S (t) represent screening function and trigonometric function battle array respectively, expression formula is: S ( t ) = cos φ j ( t ) sin φ j ( t ) - sin φ j ( t ) cos φ j ( t ) ; Wherein φ sand φ erepresent the incision of cutter tooth opposite piece respectively and cut out angle.

Further, described cutting force dynamic matrix K (t) take T/N as the periodic function matrix in cycle, i.e. K (t+T/N)=K (t), and then Fourier transformation can be carried out to each element in cutting force dynamic matrix K (t).

Further, it is larger that described positive integer α chooses, and approximate is more accurate, but computation complexity is larger.

Further, described 4 α+2 rank symmetric positive definite matrix Φ 1and Φ 2, can be taken as when α=1: Φ 1 = Φ 2 = 1 2 2 1 2 2 * 10 11 .

By the above technical scheme that the present invention conceives, compared with prior art, this technical scheme utilizes this feature of cutting force dynamic matrix cyclically-varying, employing Fourier space is similar to, design the self-adjusting Self Adaptive Control electric current along with main shaft change in displacement, and then suppress flutter by the intrinsic dynamic characteristic of magnetic suspension bearing change main axle unit, the beneficial effect carrying out Milling Process when axial cutting-in is larger and whole milling process still can be kept stable can be obtained, thus while ensure that workpiece surface quality, also improve milling efficiency.

Accompanying drawing explanation

Fig. 1 is milling Stability diagram, and transverse axis represents rotating speed size, and unit is rev/min; The longitudinal axis represents axial cutting-in size, and unit is millimeter;

Fig. 2 is closed loop control framework schematic diagram;

Fig. 3 is the agent structure schematic diagram of Milling Process Active Flutter Suppression System;

In figure, the implication of each label is as follows: 1 is milling cutter; 2 is handle of a knife; 3 is magnetic suspension bearing rotor displacement Autonomous test parts; 4 is radial magnetic bearing; 5 is magnetic suspension bearing rotor displacement Autonomous test parts; 6 is radial magnetic bearing; 7 is power amplifier; 8 is power amplifier; 9 is magnetically suspended bearing; 10 is PC; 11 is wave filter; 12 is workbench; 13 is workpiece.

Detailed description of the invention

In order to make object of the present invention, technical scheme and advantage clearly understand, below in conjunction with drawings and Examples, the present invention is further elaborated.Should be appreciated that specific embodiment described herein only in order to explain the present invention, be not intended to limit the present invention.

Below in conjunction with accompanying drawing, the specific embodiment of the present invention is described further.It should be noted that at this, the explanation for these embodiments understands the present invention for helping, but does not form limitation of the invention.In addition, if below in described each embodiment of the present invention involved technical characteristic do not form conflict each other and just can mutually combine.

In actual Milling Process, to the speed of mainshaft Ω specified arbitrarily, when axial cutting-in b is very little, dynamic cutting force F tt () is very little on the impact of whole milling process, namely whole milling process is all stable, but now milling efficiency is very low, does not meet our requirement; So wish that larger axial cutting-in is to obtain higher milling efficiency, but along with the increase of axial cutting-in b, whole system can be subject to dynamic cutting force F tthe impact of (t), the instability therefore become, produces chatter phenomenon, have impact on the crudy of workpiece; As shown in Figure 1, transverse axis represents the size of speed of mainshaft Ω, and the longitudinal axis represents the size of axial cutting-in b, with solid-line curve in scheming for boundary, represent stability region (flutter not occurring) below curve, above curve, represent unstable region (generation flutter).

The present invention, in order to improve milling efficiency and ensure workpiece processing quality, the basis of Dynamic Model of Milling Process proposes a kind of milling parameter Active Control Method based on magnetic suspension bearing electric chief axis; Specifically comprise the following steps:

(1) set up milling dynamics equation, milling dynamics equation arranged and carries out state space transformation, obtaining state space equation;

(2) according to the milling dynamics equation set up, design Self Adaptive Control electric current Q c(t).

More specifically, in step (), when main shaft is static, obtain the modal mass m of milling cutter system architecture, the dampingratioζ of milling cutter system architecture and the intrinsic pi ω of milling cutter system architecture by identification n, the milling cutter system architecture modal parameter obtained when recycling main shaft is static replaces modal parameter during main shaft rotation milling, and then set up corresponding 2DOF system model to describe the vibration of main shaft, the spindle vibration model of acquisition is:

m q · · 1 ( t ) m q · · 2 ( t ) + 2 ζ ω n 0 0 2 ζ ω n q · 1 ( t ) q · 2 ( t ) + ω 2 0 0 ω 2 q 1 ( t ) q 2 ( t ) = 0

Wherein q 1(t) and q 2t () is the displacement in main shaft two frees degree.

Further, go out Cutting Force Coefficient by milling Experimental Identification, carry out dynamic cutting force F t(t) modeling; In conjunction with the cutting force COEFFICIENT K that identification obtains twith normal direction Cutting Force Coefficient K r, the dynamic cutting force model of foundation is F t ( t ) = bf z K ( t ) 1 0 + bK ( t ) q 1 ( t - T / N ) - q 1 ( t ) q 2 ( t - T / N ) - q 2 ( t ) , Wherein b represents axial cutting-in, f zrepresent every cutter tooth amount of feeding, T represents main shaft swing circle, and N represents the number of teeth of milling cutter, K ( t ) = Σ j = 1 N g j ( φ j ( t ) ) S ( t ) [ K t K r ] T [ sin φ j ( t ) cos φ j ( t ) ] Be cutting force dynamic matrix, cutting force dynamic matrix K (t) take T/N as the periodic function matrix in cycle, i.e. K (t+T/N)=K (t), φ jt () represents the anglec of rotation of a jth cutter tooth of milling cutter, the subscript T of vector represents the transposition (lower same) to vector, g jj(t)) and S (t) represent screening function and trigonometric function battle array respectively, expression formula is:

S ( t ) = cos φ j ( t ) sin φ j ( t ) - sin φ j ( t ) cos φ j ( t )

Wherein φ sand φ erepresent the incision of cutter tooth opposite piece respectively and cut out angle.

And then can described in establishment step () milling dynamics equation be:

m q · · 1 ( t ) m q · · 2 ( t ) + 2 ζ ω n 0 0 2 ζ ω n q · 1 ( t ) q · 2 ( t ) + ω n 2 0 0 ω n 2 q 1 ( t ) q 2 ( t ) = bf z K ( t ) 1 0 + bK ( t ) q 1 ( t - T / N ) - q 1 ( t ) q 2 ( t - T / N ) - q 2 ( t )

And state space equation is as follows described in step ():

q · 1 ( t ) q · 2 ( t ) q · · 1 ( t ) q · · 2 ( t ) = 0 0 1 0 0 0 0 1 - ω n 2 0 - 2 ζ ω n 0 0 - ω n 2 0 - 2 ζ ω n q 1 ( t ) q 2 ( t ) q · 1 ( t ) q · 2 ( t ) +

b 0 0 0 0 1 / m 0 0 1 / m f z K ( t ) { 1 0 + K ( t ) 1 0 0 0 0 1 0 0 q 1 ( t - T / N - q 1 ( t ) q 2 ( t - T / N ) - q 2 ( t ) q · 1 ( t - T / N ) - q · 1 ( t ) q · 2 ( t - T / N ) - q · 2 ( t ) }

Order,

Q ( t ) q 1 ( t ) q 2 ( t ) q · 1 ( t ) q · 2 ( t ) , A = 0 0 1 0 0 0 0 1 - ω n 2 0 - 2 ζ ω n 0 0 - ω n 2 0 - 2 ζ ω n , B = 0 0 0 0 1 / m 0 0 1 / m

Then above-mentioned formula state space equation becomes:

Q · ( t ) = AQ ( t ) + bB { f z K ( t ) 1 0 T + K ( t ) C q [ Q ( t - T / N ) - Q ( t ) ] }

Wherein, represent state velocity, A represents state displacement term coefficient matrix, and B represents Cutting Force Coefficient matrix, C qrepresent hangover state displacement term coefficient matrix.

In step (two), Self Adaptive Control electric current Q ct the design process of () is as follows:

(1) Fourier transformation is carried out to each element in cutting force dynamic matrix K (t), obtain cutting force harmonic wave dynamic matrix K α(t), and according to cutting force harmonic wave dynamic matrix K αt () obtains harmonic wave

Wherein cutting force harmonic wave dynamic matrix K α(t) for each element in cutting force dynamic matrix K (t) carry out Fourier expansion front 2 α+1 item parts and, expression formula is: the constant vector that 2 unknown α+1 tie up can be set to; ω=2 π N/T is harmonic frequency, and α is the suitable positive integer chosen, and α is larger, and approximate is more accurate, but computation complexity is larger.At this, we choose α=2, can set further harmonic displacement as wherein electro spindle displacement [F xf y] t=C q[Q (t-T/N)-Q (t)], the electro spindle displacement F in x direction xwith the electro spindle displacement F in y direction yare all vectors of 1 dimension, x direction and y direction are perpendicular in the plane of main shaft, and x direction is mutually vertical with y direction, and h (t) is 4 α+2 dimensional vectors.

(2) according to state displacement term coefficient matrix A and Lyapunov Equation Χ A+A tΧ=-I obtains symmetric positive definite matrix Χ, obtains state weight coefficient [Λ according to symmetric positive definite matrix Χ 1Λ 2]=XB, according to state weight coefficient [Λ 1Λ 2] and described state vector Q (t) obtain the first self adaptation weight coefficient with the second self adaptation weight coefficient wherein I is the unit matrix that 4 row 4 arrange, Λ 1and Λ 2all four dimensional vectors, λ ithe vector of one dimension, i=1,2.

(3) the first adaptive rate is designed with the second adaptive rate get arbitrarily two 4 α+2 rank symmetric positive definite matrixs, use Φ respectively at this 1and Φ 2represent, then make the first adaptive change rate be second adaptive change rate is again by separating the adaptive change rate of differential equation form, the first adaptive rate can be obtained with the second adaptive rate wherein the first adaptive rate with the second adaptive rate be 4 α+2 dimensional vectors.

(4) according to the first adaptive rate second adaptive rate and harmonic displacement h (t) obtains Self Adaptive Control electric current Q c(t), its expression formula is: wherein Q ct () is 2 dimensional vectors.

Further, the Self Adaptive Control electric current Q will calculated ct () is multiplied by corresponding magnetic suspension bearing electric chief axis COEFFICIENT K amb, the size obtaining the power that magnetic suspension bearing produces is: F a(t)=K ambq c(t), thus the closed-loop control system (as shown in Figure 2) that formation one is stable, namely Q · ( t ) = AQ ( t ) + bB { f z K ( t ) 1 0 T + K ( t ) C q [ Q ( t - T / N ) - Q ( t ) ] } + F a ( t ) , Wherein K amb4 row 2 column matrix.

Secondly, the online ACTIVE CONTROL based on Milling Process active flutter surppression device is carried out according to the Self Adaptive Control electric current of design; As shown in Figure 3, the critical piece needed for online ACTIVE CONTROL is made up of magnetic suspension bearing milling electro spindle, magnetic suspension bearing rotor displacement Autonomous test parts, magnetically suspended bearing, signal APU four major part; Wherein, magnetic suspension bearing milling electro spindle comprises the spindle units such as two radial magnetic bearings 4 and 6, cutter 1, handle of a knife 2; Magnetic suspension bearing rotor displacement Autonomous test parts comprise band-pass circuit, one-level amplifying circuit, differential detection circuit, demodulator circuit, low-pass filter circuit etc., mainly integrated in the drawings 3 and 5 in; Magnetically suspended bearing is an integrated manipulator 9 based on DSP; Signal APU comprises wave filter 11, power amplifier 7 and 8, A/D and D/A signaling conversion circuit.Above-mentioned parts are adopted to realize the online active flutter surppression of Digit Control Machine Tool Milling Process, its process is: spindle motor carry its tools 1 rotates and carries out Milling Process to workpiece 13, in the process when axial cutting-in b is larger, very easily there is flutter, our chatter phenomenon by producing in magnetic suspension bearing rotor displacement Autonomous test parts 3 and 5 Real-Time Monitoring process; When detecting that chatter phenomenon occurs, magnetic suspension bearing rotor displacement Autonomous test parts 3 and 5 amplify through power amplifier 8 after vibration signal being passed to wave filter 11, then magnetically suspended bearing 9 is passed to by A/D conversion, and in PC10, show this vibration signal, vibratory output size is brought into the Self Adaptive Control electric current Q of offline design above simultaneously ct () calculates Self Adaptive Control electric current Q in PC10 cthe size of (t), the Self Adaptive Control electric current Q finally will calculated ct output signal that () obtains after magnetically suspended bearing 9 processes is processed through power amplifier 7 after being changed by D/A, be converted into current signal, act on radial magnetic bearing 4 and 6, produce the flutter of corresponding magnetic field force to main shaft to suppress, ensure that the stability of Milling Process accordingly.Point curve wherein in the visible Fig. 1 of closed-loop control simulated effect, the stability region that the stability region that this closed loop point curve is formed is formed much larger than open loop solid-line curve, namely when axial cutting-in is larger, closed-loop system is also stable, and then improves working (machining) efficiency.

Those skilled in the art will readily understand; the foregoing is only preferred embodiment of the present invention; not in order to limit the present invention, all any amendments done within the spirit and principles in the present invention, equivalent replacement and improvement etc., all should be included within protection scope of the present invention.

Claims (6)

1., based on a milling parameter Active Control Method for magnetic suspension bearing electric chief axis, it is characterized in that, comprise the steps:
(1) current displacement signal during electro spindle generation milling parameter is gathered;
(2), after filter and amplification and differential process being carried out to described current displacement signal, the state vector Q (t) be made up of electro spindle displacement and speed is obtained;
(3) electro spindle displacement [F is obtained according to described state vector Q (t) xf y] t=C q[Q (t-T/N)-Q (t)];
Wherein C q = 1 0 0 0 0 1 0 0 , The electro spindle displacement that Q (t-T/N) is the t-T/N moment and the state vector that speed is formed, T represents main shaft swing circle, and N represents the number of teeth of milling cutter, F xfor the electro spindle displacement in x direction, F yfor the electro spindle displacement in y direction, subscript T represents the transposition to vector;
(4) Fourier transformation is carried out to each element in cutting force dynamic matrix K (t), obtain cutting force harmonic wave dynamic matrix K α(t), and according to cutting force harmonic wave dynamic matrix K αt () obtains trigonometric function column vector
Wherein cutting force harmonic wave dynamic matrix K α(t) for each element in cutting force dynamic matrix K (t) carry out Fourier expansion front 2 α+1 item parts and, ω=2 π N/T is harmonic frequency, and α is positive integer;
(5) according to described electro spindle displacement [F xf y] twith described trigonometric function column vector obtain harmonic displacement
(6) according to state displacement term coefficient matrix A and Lyapunov Equation Χ A+A tΧ=-I obtains symmetric positive definite matrix Χ, obtains state weight coefficient [Λ according to symmetric positive definite matrix Χ 1Λ 2]=XB, according to state weight coefficient [Λ 1Λ 2] and described state vector Q (t) obtain the first self adaptation weight coefficient with the second self adaptation weight coefficient b is Cutting Force Coefficient matrix;
(7) according to described harmonic displacement h (t), described first self adaptation weight coefficient λ 1with described second self adaptation weight coefficient λ 2obtain the first adaptive rate with the second adaptive rate
Wherein, the first adaptive change rate second adaptive change rate Φ 1and Φ 2be respectively the symmetric positive definite matrix on two 4 α+2 rank of getting arbitrarily, the first adaptive rate with the second adaptive rate be 4 α+2 dimensional vectors;
(8) according to described first adaptive rate described second adaptive rate and described harmonic displacement h (t) obtains Self Adaptive Control electric current (9) by described Self Adaptive Control electric current Q ct () acts on radial magnetic bearing, and produce corresponding magnetic field force, thus realizes suppressing the flutter of main shaft.
2. milling parameter Active Control Method as claimed in claim 1, it is characterized in that, described state displacement term coefficient matrix A is:
A = 0 0 1 0 0 0 0 1 - ω n 2 0 - 2 ζω n 0 0 - ω n 2 0 - 2 ζω n , ζ is the damping ratio of milling cutter system architecture, ω nfor the intrinsic frequency of milling cutter system architecture.
3. milling parameter Active Control Method as claimed in claim 1, it is characterized in that, described Cutting Force Coefficient matrix B is: B = 0 0 0 0 1 / m 0 0 1 / m , M is the modal mass of milling cutter system architecture.
4. milling parameter Active Control Method as claimed in claim 1, it is characterized in that, described cutting force dynamic matrix is: K ( t ) = Σ j = 1 N g j ( φ j ( t ) ) S ( t ) K t K r T sinφ j ( t ) cosφ j ( t ) ;
Wherein φ jt () represents the anglec of rotation of a jth cutter tooth of milling cutter, K tand K rrepresent Cutting Force Coefficient, g jj(t)) and S (t) represent screening function and trigonometric function battle array respectively; S ( t ) = cosφ j ( t ) sinφ j ( t ) - sinφ j ( t ) cosφ j ( t ) ; φ sand φ erepresent the incision of cutter tooth opposite piece respectively and cut out angle.
5. milling parameter Active Control Method as claimed in claim 1, is characterized in that, described cutting force dynamic matrix K (t) take T/N as the periodic function matrix in cycle, K (t+T/N)=K (t).
6. milling parameter Active Control Method as claimed in claim 1, is characterized in that, described 4 α+2 rank symmetric positive definite matrix Φ 1and Φ 2, when α=1, Φ 1 = Φ 2 = 1 2 2 1 2 2 * 10 11 .
CN201410805910.7A 2014-12-22 2014-12-22 A kind of milling parameter Active Control Method based on magnetic suspension bearing electric chief axis CN104647132B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410805910.7A CN104647132B (en) 2014-12-22 2014-12-22 A kind of milling parameter Active Control Method based on magnetic suspension bearing electric chief axis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410805910.7A CN104647132B (en) 2014-12-22 2014-12-22 A kind of milling parameter Active Control Method based on magnetic suspension bearing electric chief axis

Publications (2)

Publication Number Publication Date
CN104647132A CN104647132A (en) 2015-05-27
CN104647132B true CN104647132B (en) 2016-01-20

Family

ID=53239042

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410805910.7A CN104647132B (en) 2014-12-22 2014-12-22 A kind of milling parameter Active Control Method based on magnetic suspension bearing electric chief axis

Country Status (1)

Country Link
CN (1) CN104647132B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105750570B (en) * 2016-04-13 2017-10-20 西安交通大学 A kind of electro spindle milling parameter time delay Active Control Method and its system
CN105700476B (en) * 2016-04-28 2018-05-01 上海交通大学 A kind of flutter Active Control Method under model-free parameter and driver saturation
CN106392753B (en) * 2016-05-11 2018-06-26 西安交通大学 A kind of intelligence milling electro spindle active cutting-vibration suppression system and method online
CN106363450B (en) * 2016-09-07 2018-10-09 北京理工大学 A kind of online suppressing method of milling parameter
CN108490789A (en) * 2018-05-08 2018-09-04 西安交通大学 A kind of electro spindle milling parameter robust Active Control Method and its active control system
CN109396948A (en) * 2018-11-16 2019-03-01 华中科技大学 A kind of Active Control Method and system improving five-axis milling system processing stability

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010023162A (en) * 2008-07-17 2010-02-04 Okuma Corp Chatter vibration suppression method of machine tool and device used for the same
CN102387892A (en) * 2009-04-10 2012-03-21 日本Nt工程技术株式会社 Method and device for suppressing chattering of work machine
CN102554685A (en) * 2010-10-13 2012-07-11 大隈株式会社 Vibration suppressing method and vibration suppressing device for use in machine tool
EP2540442A1 (en) * 2010-03-30 2013-01-02 Okuma Corporation Tool trajectory generation device, tool trajectory computation method, and tool trajectory generation program
CN103252680A (en) * 2013-04-24 2013-08-21 华中科技大学 Milling machining chatter active control system and method thereof
CN103769945A (en) * 2012-10-18 2014-05-07 大隈株式会社 Chatter vibration suppressing method and machine tool

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010023162A (en) * 2008-07-17 2010-02-04 Okuma Corp Chatter vibration suppression method of machine tool and device used for the same
CN102387892A (en) * 2009-04-10 2012-03-21 日本Nt工程技术株式会社 Method and device for suppressing chattering of work machine
EP2540442A1 (en) * 2010-03-30 2013-01-02 Okuma Corporation Tool trajectory generation device, tool trajectory computation method, and tool trajectory generation program
CN102554685A (en) * 2010-10-13 2012-07-11 大隈株式会社 Vibration suppressing method and vibration suppressing device for use in machine tool
CN103769945A (en) * 2012-10-18 2014-05-07 大隈株式会社 Chatter vibration suppressing method and machine tool
CN103252680A (en) * 2013-04-24 2013-08-21 华中科技大学 Milling machining chatter active control system and method thereof

Also Published As

Publication number Publication date
CN104647132A (en) 2015-05-27

Similar Documents

Publication Publication Date Title
Sebastian et al. Characterization of the unsteady aerodynamics of offshore floating wind turbines
CN104242769B (en) Permanent magnet synchronous motor speed composite control method based on continuous terminal slip form technology
Bossanyi Wind turbine control for load reduction
Zorzi et al. Finite element simulation of rotor-bearing systems with internal damping
Tarng et al. Chatter suppression in turning operations with a tuned vibration absorber
Yang et al. Design and implementation of two-degree-of-freedom tuned mass damper in milling vibration mitigation
Munoa et al. Chatter suppression in ram type travelling column milling machines using a biaxial inertial actuator
CN102873381B (en) High-speed milling process parameter optimizing method based on dynamic model
Strah et al. Speed and active power control of hydro turbine unit
Faassen et al. An improved tool path model including periodic delay for chatter prediction in milling
CN103823945A (en) Flutter stability domain modeling approach for face cutting process
Noshadi et al. System identification and robust control of multi-input multi-output active magnetic bearing systems
CN106126778B (en) Thin-wall part week with curved surface mills stability prediction method
CN103267085A (en) Hydraulic multifrequency dynamic vibration absorber for rotary mechanical rotor
Bekakra et al. Optimal tuning of PI controller using PSO optimization for indirect power control for DFIG based wind turbine with MPPT
CN102554326A (en) Milling finish machining method based on dynamic rigidity of impeller blade
CN105312679A (en) Three-dimensional ultrasound elliptical vibration assisted cutting device and elliptical orbit generation method
JP2013000850A (en) Controller and control method of machine tool
Subrahmanyam et al. Vibration and buckling of rotating, pretwisted, preconed beams including Coriolis effects
Choi et al. Control of linear motor machine tool feed drives for end milling: Robust MIMO approach
Zheng et al. Active balancing control of AMB-rotor systems using a phase-shift notch filter connected in parallel mode
Sivrioglu et al. Sliding mode control with time-varying hyperplane for AMB systems
CN104458128A (en) Turbocharger rotor unbalance amount control method based on dynamic characteristics
CN107457609B (en) Milling parameter suppressing method and milling parameter optimization system based on stiffness variation
Ullah et al. Adaptive fractional order terminal sliding mode control of a doubly fed induction generator-based wind energy system

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model