CN108107734A - A kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method - Google Patents

A kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method Download PDF

Info

Publication number
CN108107734A
CN108107734A CN201711411778.1A CN201711411778A CN108107734A CN 108107734 A CN108107734 A CN 108107734A CN 201711411778 A CN201711411778 A CN 201711411778A CN 108107734 A CN108107734 A CN 108107734A
Authority
CN
China
Prior art keywords
msub
mrow
mechanical
permanent magnetic
msubsup
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
CN201711411778.1A
Other languages
Chinese (zh)
Other versions
CN108107734B (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201711411778.1A priority Critical patent/CN108107734B/en
Publication of CN108107734A publication Critical patent/CN108107734A/en
Application granted granted Critical
Publication of CN108107734B publication Critical patent/CN108107734B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Linear Motors (AREA)

Abstract

The present invention provides a kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method, comprises the following steps:1. utilize spectrum signature of the electromagnetic energy method analysis with calculating permanent magnetic linear synchronous motor servo thrust output;2. establish permanent magnetic linear synchronous motor feed system mechanical kinetics model using Lagrange's equation;3. analyze the interactively that intercouples between permanent magnetic linear synchronous motor feed system servo-drive and mechanical system;4. according to the interactively that intercouples between servo-drive and mechanical system, permanent magnetic linear synchronous motor feed system various factors coupling electromechanical Coupling Model is established.The present invention can quickly and effectively in analysis system all kinds of coupled problems essential laws.Influence size of all kinds of mechanical-electric couplings of evaluation fast and reliable simultaneously to system motion stationarity.It further supplements due to new change caused by mechanical-electric coupling problem, the perfect integrated modelling approach of complex electromechanical systems.

Description

A kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method
Technical field
The present invention relates to High Speed NC Machine Tools Dynamic Performance Analysis fields, are specially a kind of permanent magnetic linear synchronous motor feeding System mechanical-electric coupling modeling method.
Background technology
In recent years, as the efficient height of digital control processing is refined the increasingly raising of demand, permanent-magnetism linear motor feed system is with it The advantages that outstanding high thrust, high speed, high acceleration and high-precision, before there is wide application in high-grade, digitally controlled machine tools Scape.The Zero-drive Chain structure of Linear motor feeding system uniqueness so that motor thrust directly acts on mechanical system, and external interference Etc. mechanical oscillation caused by factors also directly react on motor.Relation is intercoupled between servo-drive and mechanical system more Step up it is close, have become influence Linear motor feeding system robust motion an important factor for.Especially it is applied to numerical control The Linear motor feeding system of lathe, mechanical system and nonideal single inertia system, are influenced by each dynamic engaging portion characteristic, are had There is multistage Oscillatory mode shape, the coupled relation between system can be aggravated, deteriorate system motion precision.
It researchs and analyses at present in work, due to ignoring the dynamic characteristic of mechanical system, permanent magnetic linear synchronous motor feeding system Mechanical-electric coupling problem in system obtains the deep analysis and research of system not yet.And influenced by many factors, servo is driven It is dynamic between mechanical system there is not a kind of mechanical-electric coupling phenomenon, also exist between each coupling circuit and influence each other, increase The difficulty researched and analysed.How the modeling of system is carried out to mechanical-electric coupling problem with analyzing for further controlling and changing The influence of kind coupling has great importance.
The content of the invention
For problems of the prior art, the present invention provides a kind of permanent magnetic linear synchronous motor feed system electromechanics coupling Modeling method is closed, can parsing characterization quickly and effectively be carried out to a few quasi-representative mechanical-electric coupling phenomenons present in system, is disclosed The essential laws of coupling, at the same can quantitative description its influence to system motion stationarity, for commenting for mechanical-electric coupling problem Valency and control, which improve, to have great importance.
The present invention is to be achieved through the following technical solutions:
A kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method, comprises the following steps:
Step 1 analyzes the frequency spectrum spy with calculating permanent magnetic linear synchronous motor servo thrust output using electromagnetic energy method Sign;
Step 2 establishes permanent magnetic linear synchronous motor feed system mechanical kinetics model using Lagrange's equation;
Step 3, the feed system mechanical kinetics mould of servo thrust output spectrum signature and step 2 based on step 1 Type, from mechanical oscillation caused by the variation of linear motor air gap, thrust harmonic wave caused by mechanical torsional oscillation and mechanical twisting vibration Three aspects of influence to grating scale feedback signal, analyze permanent magnetic linear synchronous motor feed system servo-drive and mechanical system Between the interactively that intercouples;
Step 4 according to the interactively that intercouples between the servo-drive and mechanical system obtained in step 3, is built Vertical permanent magnetic linear synchronous motor feed system various factors coupling electromechanical Coupling Model is as follows,
Wherein, F is permanent magnetic linear synchronous motor thrust, and k=a, b, c represent a, b, c three-phase of electrical-coil;ElkAnd Emk Respectively the armature back-emf of three-phase coil and unloaded back-emf, ikFor servo three-phase current, v is motor movement speed;F0' be Consider the motor name thrust of air gap variation, Fr' original thrust the harmonic wave of motor changed for consideration air gap, FgTo consider that air gap becomes The new thrust harmonic components changed, xiIt is command signal, xoFor output signal of displacement, Gp(s) letter is passed for position ring control, s is micro- Divide operator, Gv(s) letter, K are passed for speed ring controlAFor electric current loop equieffective ratio gain, KFFor motor thrust constant, FmFor motor The nominal thrust of output, FrFor motor output harmonic wave power, GmF (s) be mechanical system direction of feed model, Gmi(s) for machinery other Directional dynamics model;△ δ are that mechanical twisting vibrates the encoder errors generated, Ge(s) it is encoder errors and mechanical twisting Biography letter between vibration;xδ0To consider the feedback signal of encoder errors.
Preferably, in step 1, according to the structure composition of driving circuit, considering driving circuit, to include dead zone humorous All kinds of non-linear factors and electric machine structure of ripple, modulation harmonic wave and counter electromotive force harmonic wave include slot effect, end effect Should be with all kinds of non-linear factors of magnetic linkage harmonic wave, using electromagnetic energy method, the frequency spectrum for calculating permanent magnetic linear synchronous motor thrust F is special Sign is as follows,
F=∑s (Elk+Emk)·ik/v。
Preferably, in step 2, ignore the flexibility of workbench, take workbench along three axis and the fortune around three axis It moves, i.e. { x, y, z, θx, θy, θz}TFor generalized coordinates;Ignore the effect of intercoupling between each vibration shape, utilize Lagrange's equation The mechanical kinetics equation of Linear motor feeding system is established,
Wherein, X=[x y z θx θy θz]TFor worktable displacement,The respectively single order of X, second dervative, M, C, K are respectively mass matrix, damping matrix and stiffness matrix, and F is external force matrix.
Preferably, when considering that linear motor air gap changes caused by being vibrated due to mechanical twisting in step 3, introduce opposite Permeance function describes influence of the mechanical twisting vibration to motor gas-gap magnetic field, and utilizes electromagnetic energy method, and consideration gas is calculated The motor thrust output of gap variation,
F=F '0+F′r+Fg
Preferably, caused by considering thrust harmonic wave in step 3 during mechanical oscillation, establish based on PID three close-loop control models, Wherein, position ring is controlled using P, and speed ring is controlled using PI, and electric current loop is equivalent to proportional gain, the motor that step 1 is obtained Thrust output is introduced in the form of interference;By the dynamic characteristic in other directions of the mechanical system in addition to direction of feed respectively etc. It imitates and introduces model for relatively independent two-order oscillation system, the parameters in model are obtained by Experimental Identification, it is final to establish The electromechanical integrated model of system under the effect of thrust harmonic wave,
[(xi-x0)·Gp(s)-x0·s]·Gv(s)·KA·KF=F '0
(F′0+F′r+Fg)·Gmf(s)+∑(F′0+F′r+Fg)·Gmi(s)=x0
Preferably, when influence of the mechanical twisting vibration to grating scale feedback signal is considered in step 3, based on grating scale work Make principle, the encoder errors that mechanical twisting vibration generates be calculated,
Δ δ=∑ (F0′+Fr′+Fg)·Gmi(s)·Ge(s),
And obtained encoder feedback error is introduced into feedback control loop, obtain considering the feedback letter of encoder errors Number,
xδ0=x0+Δδ。
Compared with prior art, the present invention has technique effect beneficial below:
The present invention is for the typical mechanical-electric coupling phenomenon of three classes present in system, it is proposed that permanent magnetic linear synchronous motor into Give system mechanical-electric coupling modeling method, can quickly and effectively in analysis system all kinds of coupled problems essential laws.Profit simultaneously System displacement fluctuation caused by can accurately calculating all kinds of mechanical-electric coupling problems with modeling method proposed by the present invention, quickly may be used Influence size of all kinds of mechanical-electric couplings of evaluation leaned on to system motion stationarity.The present invention is on the basis of traditional integrated modelling approach On, it further supplements due to new change caused by mechanical-electric coupling problem, the perfect integrated modelling approach of complex electromechanical systems, It traces to the source for the kinematic error of Linear motor feeding system, the analysis and optimization of mechanical-electric coupling problem, which improves, has important meaning Justice.
Description of the drawings
Fig. 1 is the schematic diagram of permanent magnetic linear synchronous motor feeding system structure described in present example.
Fig. 2 is permanent magnetic linear synchronous motor feed system electromechanical Coupling Model of the present invention.
Fig. 3 is the thrust result of spectrum analysis of experimental bench motor described in present example.
Fig. 4 is the result of spectrum analysis of experimental bench displacement fluctuation described in present example.
In figure:1 is workbench;2 be line slideway and sliding block;3 be grating scale scale;4 be grating ruler reading head;5 be electricity Motor-driven son;6 be motor stator;7 be lathe bed.
Specific embodiment
With reference to specific embodiment, the present invention is described in further detail, it is described be explanation of the invention and It is not to limit.
A kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method of the present invention considers servo-drive spy The interactively that intercouples of property, mechanical dynamics property and the two, establishes electromechanical Coupling Model.Specifically include following step Suddenly:
Step 1 considers all kinds of non-linear factor (such as dead zone harmonic wave, modulation harmonic wave and anti-electricity in driving circuit Kinetic potential harmonic wave etc.) and electric machine structure non-linear factor (such as slot effect, end effect, magnetic linkage harmonic wave), utilize electromagnetic energy Amount method calculates permanent magnetic linear synchronous motor thrust spectrum signature, i.e.,
F=∑s (Elk+Emk)·ik/v
Wherein, k=a, b, c represent a, b, c three-phase of electrical-coil;ElkAnd EmkThe armature of respectively three-phase coil is anti-electric Gesture and unloaded back-emf;ikFor servo three-phase current;V is motor movement speed.
Step 2 ignores the flexibility of workbench, takes workbench along three axis and the movement around three axis, i.e., and x, y, Z, θx, θy, θz}TFor generalized coordinates.Ignore the effect of intercoupling between each vibration shape, straight-line electric is established using Lagrange's equation The mechanical kinetics equation of machine feed system, i.e.,
Wherein, X=[x y z θx θy θz]TFor worktable displacement,The respectively single order of X, second dervative, M, C, K is respectively mass matrix, damping matrix and stiffness matrix, and F is external force matrix.
Step 3 considers that linear motor air gap changes caused by being vibrated due to mechanical twisting, introduces relative magnetic permeability letter first Number description mechanical twisting vibrates the influence to motor gas-gap magnetic field, and is calculated using electromagnetic energy method and considers air gap variation Motor thrust output, i.e.,
F=F '0+F′r+Fg
Wherein, F '0To consider the motor name thrust of air gap variation, F 'rTo consider that the original thrust of motor of air gap variation is humorous Ripple, FgTo consider the new thrust harmonic components of air gap variation.
Secondly, mechanical oscillation caused by considering thrust harmonic wave is established based on PID three close-loop control models, wherein, position ring is adopted It is controlled with P, speed ring is controlled using PI, and electric current loop is equivalent to proportional gain, the motor thrust output that step 1 is obtained, with dry The form disturbed introduces;Mechanical system is equivalent to relatively independent two-order oscillation system and introduces model, is obtained by Experimental Identification Parameters in model finally establish the electromechanical integrated model of the system under the effect of thrust harmonic wave, i.e.,
[(xi-x0)·Gp(s)-x0·s]·Gv(s)·KA·KF=F '0
(F′0+F′r+Fg)·Gmf(s)+∑(F′0+F′r+Fg)·Gmi(s)=x0
Wherein, xiIt is command signal, xoFor output signal of displacement, Gp(s) letter being passed for position ring control, s is differential operator, Gv(s) letter, K are passed for speed ring controlAFor electric current loop equieffective ratio gain, KFFor motor thrust constant, FmName is exported for motor Thrust, FrFor motor output harmonic wave power, GmF (s) be mechanical system direction of feed model, Gmi(s) it is other mechanical direction power Learn model.
Additionally, it is contemplated that mechanical twisting vibrates the influence to grating scale feedback signal, based on grating scale operation principle, calculate The encoder errors generated to mechanical twisting vibration, i.e.,
Δ δ=∑ (F '0+F′r+Fg)·Gmi(s)·Ge(s)
Wherein, Ge(s) be encoder errors and mechanical twisting vibration between biography letter.
Obtained encoder feedback error is introduced into feedback control loop, obtains considering the feedback letter of encoder errors Number, i.e.,
xδ0=x0+Δδ
Step 4, the interaction relationship between the servo-drive and mechanical system that are obtained based on step 3, establishes permanent magnetism Linear synchronous motor feed system various factors coupling electromechanical Coupling Model, i.e.,
Specifically, permanent magnetic linear synchronous motor feed system is as shown in Figure 1, the uniaxial linear motor feeding experiment of selection one Platform is test cases, and displacement signal is acquired by laser interferometer, sample frequency 10KHz.Motor thrust is carried by servo Monitoring of software is acquired, sample frequency 1KHz.In experiment test, feed speed 20m/min.It is as follows:
1) all kinds of non-linear factor (such as dead zone harmonic wave, modulation harmonic wave and counter electromotive force in driving circuit are considered Harmonic wave etc.) and electric machine structure non-linear factor (such as slot effect, end effect, magnetic linkage harmonic wave), using electromagnetic energy method, Permanent magnetic linear synchronous motor thrust spectrum signature is calculated, i.e.,
F=∑s (Elk+Emk)·ik/ v=F0+Fr+Fc+Fe+Fq+FL (1)
Wherein F0It is motor name thrust, FrIt is motor ripple thrust, FcIt is Slot force, FeIt is End Force, FqIt is ripple tooth socket coupling With joint efforts, FLIt is inductance asymmetry HARMONIC FORCE.
2) ignore the flexibility of workbench, take workbench along three axis and the movement around three axis, i.e. { x, y, z, θx, θy, θz}TFor generalized coordinates.Ignore the effect of intercoupling between each vibration shape, using Lagrange's equation establish linear motor into To the mechanical kinetics equation of system, i.e.,
Wherein, wherein:X=[x y z θx θy θz]TFor worktable displacement,Respectively the single order of X, second order are led Number, M, C, K are respectively mass matrix, damping matrix and stiffness matrix, and F is external force matrix.
In the test cases, mechanical structure is simple, and the coupling between each mode of machinery can be ignored.In addition, Since workbench y is to very high to intrinsic frequency with z, and the damping in the two directions is larger, and oscillatory response amplitude is relatively low, point It can be not considered in analysis feed system Self excitation Response.Therefore in subsequent analysis Main Analysis workbench in the vibration of direction of feed And the twisting vibration around each axis.
3) first, consider that linear motor air gap changes caused by being vibrated due to mechanical twisting, introduces relative magnetic permeability function and retouches Influence of the mechanical twisting vibration to motor gas-gap magnetic field is stated, and the motor for considering air gap variation is calculated using electromagnetic energy method Thrust output, i.e.,
F=F '0+F′r+Fg (3)
Wherein, F '0To consider the motor name thrust of air gap variation, F 'rTo consider the original thrust harmonic wave of motor of air gap variation, FgTo consider the new thrust harmonic components of air gap variation.
Due in present case, mechanical twisting oscillation amplitude very little, therefore motor name thrust caused by it and original thrust Harmonic wave varies less, and when analyzing influence of the air gap fluctuation to feed system robust motion, only considers what air gap variation generated New motor thrust harmonic wave.
Secondly, mechanical oscillation caused by considering thrust harmonic wave is established based on PID three close-loop control models, and wherein position ring is adopted It is controlled with P, speed ring is controlled using PI, and electric current loop is equivalent to proportional gain, the motor thrust output that step 1 is obtained, with dry The form disturbed introduces;Mechanical system direction of feed is equivalent to single inertia system, and the torsional oscillation in three directions of mechanical system is equivalent Model is introduced for relatively independent two-order oscillation system, the parameters in model are obtained by Experimental Identification, final establish pushes away The electromechanical integrated model of system under the effect of power harmonic wave, i.e.,
[(xi-x0)·Gp(s)-x0·s]·Gv(s)·KA·KF=F '0 (4)
(F′0+F′r+Fg)·Gmf(s)+∑(F′0+F′r+Fg)·Gmi(s)=x0 (5)
Wherein, xiIt is command signal, xoFor output signal of displacement, Gp(s) letter is passed for position ring control, s is differential operator, Gv(s) It is controlled for speed ring and passes letter, KAFor electric current loop equieffective ratio gain, KFFor motor thrust constant, FmNominal thrust is exported for motor, FrFor motor output harmonic wave power, GmF (s) be mechanical system direction of feed model, Gmi(s) it is other mechanical directional dynamics moulds Type.
Additionally, it is contemplated that mechanical twisting vibrates the influence to grating scale feedback signal, based on grating scale operation principle, calculate The encoder errors generated to mechanical twisting vibration, i.e.,
Δ δ=∑ (F '0+F′r+Fg)·Gmi(s)·Ge(s) (6)
Wherein, Ge(s) be encoder errors and mechanical twisting vibration between biography letter.
Obtained encoder feedback error is introduced into feedback control loop, obtains considering the feedback letter of encoder errors Number, i.e.,
xδ0=x0+Δδ (7)
4) the comprehensive model being previously obtained, establishes permanent magnetic linear synchronous motor feed system various factors coupling mechanical-electric coupling mould Type is as shown in Fig. 2, i.e.
In Fig. 2, xiFor command signal, xoIt exports and responds for system, KpFor position loop gain, KvFor speed ring gain, TvFor The speed ring time of integration, KFFor thrust constant, FrFor thrust harmonic wave caused by encoder errors, m is driving load, Jx,Jy,JzFor The rotary inertia in three directions of mechanical system, Cθx,Cθy,CθzFor the damping of three torsional oscillations of mechanical system, Kθx,Kθy,KθzFor The torsion stiffness of mechanical three torsional oscillations;My,Mp,MrFor equieffective ratio coefficient of the motor thrust on three torsional directions; Sy,Sp,SrFor the ratio conversion factor between three torsional oscillations and feeding displacement fluctuation;Gp(s),Gy(s) and Gr(s) it is respectively Transmission function between mechanical three direction torsional oscillations and encoder errors.
5) choose the obtained thrust of experiment test and displacement fluctuation data steady-state process same section, carries out spectrum analysis, As a result as shown in Figure 3 and Figure 4.In Fig. 3, the deviation of digitized representation theoretical calculation and measured result in thrust frequency spectrum bracket; In Fig. 4, bracket inner digital represents the frequency and amplitude of corresponding points.It, will for main thrust harmonic components according to Fig. 3 and Fig. 4 The theoretical calculation displacement fluctuation result that integrated model obtains is compared with experiment test displacement fluctuation result, as a result such as 1 institute of table Show.
Theoretical calculation and the experimental results comparison of 1 displacement fluctuation of table
It can be obtained by Fig. 3, Fig. 4 and table 1, a variety of mechanical-electric coupling phenomenons are implicitly present in Linear motor feeding system, are made The stable state displacement fluctuation that can not ignore into system.The wherein maximum deviation of displacement fluctuation the calculated results and measured result is 5.3%, it was confirmed that the accuracy and reliability of electromechanical Coupling Model result of calculation.

Claims (6)

1. a kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method, which is characterized in that comprise the following steps:
Step 1 analyzes the spectrum signature with calculating permanent magnetic linear synchronous motor servo thrust output using electromagnetic energy method;
Step 2 establishes permanent magnetic linear synchronous motor feed system mechanical kinetics model using Lagrange's equation;
Step 3, the feed system mechanical kinetics model of servo thrust output spectrum signature and step 2 based on step 1, From mechanical oscillation caused by the variation of linear motor air gap, thrust harmonic wave caused by mechanical torsional oscillation and mechanical twisting vibration to light Three aspects of influence of grid ruler feedback signal, are analyzed between permanent magnetic linear synchronous motor feed system servo-drive and mechanical system The interactively that intercouples;
Step 4 according to the interactively that intercouples between the servo-drive and mechanical system obtained in step 3, is established forever Magnetic-synchro Linear motor feeding system various factors coupling electromechanical Coupling Model is as follows,
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mi>F</mi> <mo>=</mo> <mi>&amp;Sigma;</mi> <mo>(</mo> <msub> <mi>E</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>E</mi> <mrow> <mi>m</mi> <mi>k</mi> </mrow> </msub> <mo>)</mo> <mo>&amp;CenterDot;</mo> <msub> <mi>i</mi> <mi>k</mi> </msub> <mo>/</mo> <mi>v</mi> <mo>=</mo> <msubsup> <mi>F</mi> <mn>0</mn> <mo>&amp;prime;</mo> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mi>r</mi> <mo>&amp;prime;</mo> </msubsup> <mo>+</mo> <msub> <mi>F</mi> <mi>g</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>x</mi> <mrow> <mi>&amp;delta;</mi> <mn>0</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>G</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>x</mi> <mrow> <mi>&amp;delta;</mi> <mn>0</mn> </mrow> </msub> <mo>&amp;CenterDot;</mo> <mi>s</mi> <mo>&amp;rsqb;</mo> <mo>&amp;CenterDot;</mo> <msub> <mi>G</mi> <mi>v</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>K</mi> <mi>A</mi> </msub> <mo>&amp;CenterDot;</mo> <msub> <mi>K</mi> <mi>F</mi> </msub> <mo>=</mo> <msubsup> <mi>F</mi> <mn>0</mn> <mo>&amp;prime;</mo> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>(</mo> <msubsup> <mi>F</mi> <mn>0</mn> <mo>&amp;prime;</mo> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mi>r</mi> <mo>&amp;prime;</mo> </msubsup> <mo>+</mo> <msub> <mi>F</mi> <mi>g</mi> </msub> <mo>)</mo> <mo>&amp;CenterDot;</mo> <msub> <mi>G</mi> <mrow> <mi>m</mi> <mi>f</mi> </mrow> </msub> <mo>(</mo> <mi>s</mi> <mo>)</mo> <mo>+</mo> <mi>&amp;Sigma;</mi> <mo>(</mo> <msubsup> <mi>F</mi> <mn>0</mn> <mo>&amp;prime;</mo> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mi>r</mi> <mo>&amp;prime;</mo> </msubsup> <mo>+</mo> <msub> <mi>F</mi> <mi>g</mi> </msub> <mo>)</mo> <mo>&amp;CenterDot;</mo> <msub> <mi>G</mi> <mrow> <mi>m</mi> <mi>i</mi> </mrow> </msub> <mo>(</mo> <mi>s</mi> <mo>)</mo> <mo>=</mo> <msub> <mi>x</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;Delta;</mi> <mi>&amp;delta;</mi> <mo>=</mo> <mi>&amp;Sigma;</mi> <mrow> <mo>(</mo> <msubsup> <mi>F</mi> <mn>0</mn> <mo>&amp;prime;</mo> </msubsup> <mo>+</mo> <msubsup> <mi>F</mi> <mi>r</mi> <mo>&amp;prime;</mo> </msubsup> <mo>+</mo> <msub> <mi>F</mi> <mi>g</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>G</mi> <mrow> <mi>m</mi> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>G</mi> <mi>e</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mrow> <mi>&amp;delta;</mi> <mn>0</mn> </mrow> </msub> <mo>=</mo> <msub> <mi>x</mi> <mn>0</mn> </msub> <mo>+</mo> <mi>&amp;Delta;</mi> <mi>&amp;delta;</mi> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
Wherein, F is permanent magnetic linear synchronous motor thrust, and k=a, b, c represent a, b, c three-phase of electrical-coil;ElkAnd EmkRespectively For the armature back-emf of three-phase coil and unloaded back-emf, ikFor servo three-phase current, v is motor movement speed;F0' it is to consider The motor name thrust of air gap variation, Fr' original thrust the harmonic wave of motor changed for consideration air gap, FgTo consider air gap variation New thrust harmonic components, xiIt is command signal, xoFor output signal of displacement, Gp(s) letter is passed for position ring control, s calculates for differential Son, Gv(s) letter, K are passed for speed ring controlAFor electric current loop equieffective ratio gain, KFFor motor thrust constant, FmIt is exported for motor Nominal thrust, FrFor motor output harmonic wave power, GmF (s) be mechanical system direction of feed model, Gmi(s) it is other mechanical directions Kinetic model;△ δ are that mechanical twisting vibrates the encoder errors generated, Ge(s) it is that encoder errors are vibrated with mechanical twisting Between biography letter;xδ0To consider the feedback signal of encoder errors.
2. a kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method according to claim 1, feature It is, in step 1, according to the structure composition of driving circuit, considering driving circuit includes dead zone harmonic wave, modulation harmonic wave And all kinds of non-linear factors and electric machine structure of counter electromotive force harmonic wave include slot effect, end effect and magnetic linkage harmonic wave All kinds of non-linear factors, using electromagnetic energy method, it is as follows to calculate the spectrum signature of permanent magnetic linear synchronous motor thrust F,
F=∑s (Elk+Emk)·ik/v。
3. a kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method according to claim 1, feature Be, in step 2, ignore the flexibility of workbench, take workbench along three axis and the movement around three axis, i.e., x, y, z, θx, θy, θz}TFor generalized coordinates;Ignore the effect of intercoupling between each vibration shape, linear motor is established using Lagrange's equation The mechanical kinetics equation of feed system,
<mrow> <mi>M</mi> <mover> <mi>X</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>C</mi> <mover> <mi>X</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>K</mi> <mi>X</mi> <mo>=</mo> <mi>F</mi> </mrow>
Wherein, X=[x y z θx θy θz]TFor worktable displacement,The respectively single order of X, second dervative, M, C, K points Not Wei mass matrix, damping matrix and stiffness matrix, F be external force matrix.
4. a kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method according to claim 1, feature It is, when considering that linear motor air gap changes caused by being vibrated due to mechanical twisting in step 3, introduces relative magnetic permeability function and retouch Influence of the mechanical twisting vibration to motor gas-gap magnetic field is stated, and utilizes electromagnetic energy method, the electricity for considering air gap variation is calculated Machine thrust output,
F=F0′+Fr′+Fg
5. a kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method according to claim 1, feature It is, caused by considering thrust harmonic wave in step 3 during mechanical oscillation, establishes based on PID three close-loop control models, wherein, position ring It being controlled using P, speed ring is controlled using PI, and electric current loop is equivalent to proportional gain, the motor thrust output that step 1 is obtained, with The form of interference introduces;The dynamic characteristic in other directions of the mechanical system in addition to direction of feed is equivalent to respectively relatively independent Two-order oscillation system introduce model, the parameters in model are obtained by Experimental Identification, finally establish thrust harmonic wave effect Under system electromechanical integrated model,
[(xi-x0)·Gp(s)-x0·s]·Gv(s)·KA·KF=F0
(F0′+Fr′+Fg)·Gmf(s)+∑(F0′+Fr′+Fg)·Gmi(s)=x0
6. a kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method according to claim 1, feature It is, when influence of the mechanical twisting vibration to grating scale feedback signal is considered in step 3, based on grating scale operation principle, calculates The encoder errors that mechanical twisting vibration generates are obtained,
Δ δ=∑ (F0′+Fr′+Fg)·Gmi(s)·Ge(s),
And obtained encoder feedback error is introduced into feedback control loop, obtain considering the feedback signal of encoder errors,
xδ0=x0+Δδ。
CN201711411778.1A 2017-12-23 2017-12-23 Electromechanical coupling modeling method for permanent magnet synchronous linear motor feeding system Expired - Fee Related CN108107734B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711411778.1A CN108107734B (en) 2017-12-23 2017-12-23 Electromechanical coupling modeling method for permanent magnet synchronous linear motor feeding system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711411778.1A CN108107734B (en) 2017-12-23 2017-12-23 Electromechanical coupling modeling method for permanent magnet synchronous linear motor feeding system

Publications (2)

Publication Number Publication Date
CN108107734A true CN108107734A (en) 2018-06-01
CN108107734B CN108107734B (en) 2020-06-19

Family

ID=62212531

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711411778.1A Expired - Fee Related CN108107734B (en) 2017-12-23 2017-12-23 Electromechanical coupling modeling method for permanent magnet synchronous linear motor feeding system

Country Status (1)

Country Link
CN (1) CN108107734B (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109581967A (en) * 2018-12-29 2019-04-05 西安交通大学 A kind of high-speed feeding system speed anti-vibration planing method based on mechanical electric coupling
CN110083944A (en) * 2019-04-28 2019-08-02 重庆理工大学 Power drive system mechanical-electric coupling hybrid modeling method
CN110231775A (en) * 2019-07-02 2019-09-13 西安交通大学 A kind of accurate prediction technique of direct-drive high-speed feed system kinematic accuracy considering mechanical twisting oscillation
CN111086004A (en) * 2020-01-08 2020-05-01 山东理工大学 Human-simulated flexible joint arm electromechanical coupling modeling method
CN111123834A (en) * 2019-12-26 2020-05-08 西安交通大学 Method for evaluating electromechanical coupling strength of linear motor feeding system
CN111752213A (en) * 2019-12-26 2020-10-09 西安交通大学 Electromechanical coupling analysis method applied to high-speed high-precision direct-drive feeding system
CN112446105A (en) * 2020-10-15 2021-03-05 昆明理工大学 Method for establishing Lagrange equation of diesel engine
CN114619451A (en) * 2022-03-31 2022-06-14 西安交通大学 Electromechanical control inductive coupling modeling method of serial industrial robot
WO2023179060A1 (en) * 2022-03-24 2023-09-28 亚龙智能装备集团股份有限公司 Linear electric motor module

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102022201150A1 (en) * 2022-02-03 2023-08-03 Volkswagen Aktiengesellschaft Electronically commutated electrical machine with several sub-machines and method for operating such an electrical machine

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101340173A (en) * 2008-04-30 2009-01-07 沈阳工业大学 Method for suppressing pulse of linear motor pushing force system
RU2430397C2 (en) * 2009-11-17 2011-09-27 Открытое акционерное общество "Научно-исследовальский институт "Субмикрон" Electrohydraulic servo drive adaptive control system
CN202713223U (en) * 2012-06-11 2013-01-30 桂林电子科技大学 Non-linear interference control system for permanent magnetism synchronization linear motor
CN106130432A (en) * 2016-07-07 2016-11-16 中国矿业大学 A kind of permanent magnetic linear synchronous motor force oscillation suppression technology based on complex controll
CN106849791A (en) * 2017-01-05 2017-06-13 江苏大学 It is a kind of to suppress the sliding mode speed control method based on internal model that permanent-magnetism linear motor mismatches disturbance
CN107203184A (en) * 2017-06-20 2017-09-26 南京理工大学 The dynamic control method of straight line steering wheel Electric Loading System
CN108021039A (en) * 2017-12-23 2018-05-11 西安交通大学 A kind of Linear motor feeding system electromechanics integrated modelling approach

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101340173A (en) * 2008-04-30 2009-01-07 沈阳工业大学 Method for suppressing pulse of linear motor pushing force system
RU2430397C2 (en) * 2009-11-17 2011-09-27 Открытое акционерное общество "Научно-исследовальский институт "Субмикрон" Electrohydraulic servo drive adaptive control system
CN202713223U (en) * 2012-06-11 2013-01-30 桂林电子科技大学 Non-linear interference control system for permanent magnetism synchronization linear motor
CN106130432A (en) * 2016-07-07 2016-11-16 中国矿业大学 A kind of permanent magnetic linear synchronous motor force oscillation suppression technology based on complex controll
CN106849791A (en) * 2017-01-05 2017-06-13 江苏大学 It is a kind of to suppress the sliding mode speed control method based on internal model that permanent-magnetism linear motor mismatches disturbance
CN107203184A (en) * 2017-06-20 2017-09-26 南京理工大学 The dynamic control method of straight line steering wheel Electric Loading System
CN108021039A (en) * 2017-12-23 2018-05-11 西安交通大学 A kind of Linear motor feeding system electromechanics integrated modelling approach

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
XIAOJUN YANG,DUN LU,WANHUA ZHAO,: "Decoupling and effects of the mechanical vibration on the dynamic precision for the direct-driven machine tool", 《THE INTERNATION JOURNAL OF ADVANCED MANUFACTURING TECHNOLOGY》 *
XIAOJUN YANG,DUNLU,CHENGFANGMA, JUNZHANG, WANHUAZHAO: "Analysis on the multi-dimensional spectrum of the thrust force for the linear motor feed drive system in machine tools", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》 *
XIAOJUN YANG,DUNLU,JUNZHANG,WANHUAZHAO: "Dynamic electromechanical coupling resulting from the air-gap fluctuation of the linear motor in machine tools", 《INTERNATIONAL JOURNAL OF MACHINE TOOLS & MANUFACTURE》 *
孙立宁: "高加速精密直线电机运动平台键合图建模与分析", 《机器人》 *
宁建荣: "高精数控机床永磁直线进给系统机电耦合法向振动研究", 《中国博士学位论文全文数据库》 *
曾理湛: "永磁同步直线电机复杂机电耦合系统分析与控制方法", 《中国博士学位论文全文数据库》 *
杨晓君,赵万华,刘辉,张会杰: "直线电机进给系统机械系统动态特性研究", 《西安交通大学学报》 *
武轲: "永磁直线同步电机的建模及控制方法研究", 《中国优秀硕士学位论文全文数据库》 *
马成方,杨晓君,张俊,吕盾,赵万华: "面向运动平稳性的直线电机推力频谱特征分析", 《西安交通大学学报》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109581967A (en) * 2018-12-29 2019-04-05 西安交通大学 A kind of high-speed feeding system speed anti-vibration planing method based on mechanical electric coupling
CN110083944A (en) * 2019-04-28 2019-08-02 重庆理工大学 Power drive system mechanical-electric coupling hybrid modeling method
CN110083944B (en) * 2019-04-28 2023-04-07 重庆理工大学 Electromechanical coupling hybrid modeling method for electric drive system
CN110231775A (en) * 2019-07-02 2019-09-13 西安交通大学 A kind of accurate prediction technique of direct-drive high-speed feed system kinematic accuracy considering mechanical twisting oscillation
CN111123834A (en) * 2019-12-26 2020-05-08 西安交通大学 Method for evaluating electromechanical coupling strength of linear motor feeding system
CN111752213A (en) * 2019-12-26 2020-10-09 西安交通大学 Electromechanical coupling analysis method applied to high-speed high-precision direct-drive feeding system
CN111086004A (en) * 2020-01-08 2020-05-01 山东理工大学 Human-simulated flexible joint arm electromechanical coupling modeling method
CN111086004B (en) * 2020-01-08 2022-09-13 山东理工大学 Human-simulated flexible joint arm electromechanical coupling modeling method
CN112446105A (en) * 2020-10-15 2021-03-05 昆明理工大学 Method for establishing Lagrange equation of diesel engine
WO2023179060A1 (en) * 2022-03-24 2023-09-28 亚龙智能装备集团股份有限公司 Linear electric motor module
CN114619451A (en) * 2022-03-31 2022-06-14 西安交通大学 Electromechanical control inductive coupling modeling method of serial industrial robot
CN114619451B (en) * 2022-03-31 2023-12-05 西安交通大学 Electromechanical control inductive coupling modeling method for serial industrial robots

Also Published As

Publication number Publication date
CN108107734B (en) 2020-06-19

Similar Documents

Publication Publication Date Title
CN108107734A (en) A kind of permanent magnetic linear synchronous motor feed system mechanical-electric coupling modeling method
CN108021039B (en) Electromechanical integration modeling method for linear motor feeding system
Zhang et al. Tracking control of ball screw drives using ADRC and equivalent-error-model-based feedforward control
Bai et al. Time-optimal freeform S-curve profile under positioning error and robustness constraints
Chow et al. Disturbance and response time improvement of submicrometer precision linear motion system by using modified disturbance compensator and internal model reference control
CN101639667A (en) Method for designing servo system
Inoue et al. Spatial disturbance suppression of a flexible system based on wave model
Ulu et al. Design and analysis of a modular learning based cross-coupled control algorithm for multi-axis precision positioning systems
Ghaffari et al. Calculating the frequency of oscillation for servo axes distressed by clearance or preloading
Takeuchi et al. Motion reproduction with time-adaptation control for dealing with variations of environmental location
CN112207331B (en) Dynamic integrated modeling method for milling process of single-shaft feeding system
Mekid et al. High precision linear slide. Part II: control and measurements
Peng et al. The design and control of a rigid-flexible coupling positioning stage for enhanced settling performance
CN111610753B (en) Differential micro-feeding servo system of double linear motors and control method
Jastrzebski et al. Analysis of a segmented axial active magnetic bearing for multi-MW compressor applications
Zhang et al. Novel dynamic microactuation method for tracking-error reduction of permanent magnet linear synchronous motors
Nakamura et al. Torque control of a series elastic actuator using an ultrasonic motor with angular‐velocity saturation
Liu et al. The immersion and invariance adaptive control for a class of linear motor systems with its application
Xu et al. Robust μ parameterization with low tuning complexity of cascaded control for feed drives
Kwon et al. Simultaneous reduction of normal-and tangential-direction force ripples in permanent-magnet linear synchronous motors
Ke et al. Tangential velocity tracking-based task coordinate frame approach for contouring control of biaxial motion systems
Saarakkala Identification and speed control design of resonating mechanical systems in electric drives
CN111752213B (en) Electromechanical coupling analysis method applied to high-speed high-precision direct-drive feeding system
M Raafat et al. Survey on robust control of precision positioning systems
Väisänen Simulation of Electromechanically Actuated Boom

Legal Events

Date Code Title Description
PB01 Publication
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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200619