WO2015092917A1 - リニアモータの推力定数導出方法および移動制御方法、ならびに、リニアモータの推力定数導出装置および移動制御装置 - Google Patents

リニアモータの推力定数導出方法および移動制御方法、ならびに、リニアモータの推力定数導出装置および移動制御装置 Download PDF

Info

Publication number
WO2015092917A1
WO2015092917A1 PCT/JP2013/084269 JP2013084269W WO2015092917A1 WO 2015092917 A1 WO2015092917 A1 WO 2015092917A1 JP 2013084269 W JP2013084269 W JP 2013084269W WO 2015092917 A1 WO2015092917 A1 WO 2015092917A1
Authority
WO
WIPO (PCT)
Prior art keywords
thrust
thrust constant
constant
deriving
moving body
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.)
Ceased
Application number
PCT/JP2013/084269
Other languages
English (en)
French (fr)
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.)
Fuji Corp
Original Assignee
Fuji Machine Manufacturing Co Ltd
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 Fuji Machine Manufacturing Co Ltd filed Critical Fuji Machine Manufacturing Co Ltd
Priority to PCT/JP2013/084269 priority Critical patent/WO2015092917A1/ja
Priority to CN201380081715.8A priority patent/CN105830321B/zh
Priority to JP2015553301A priority patent/JP6313335B2/ja
Priority to US15/102,543 priority patent/US10097127B2/en
Priority to EP13899693.9A priority patent/EP3086468B1/en
Publication of WO2015092917A1 publication Critical patent/WO2015092917A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/06Linear motors
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/06Linear motors
    • H02P25/064Linear motors of the synchronous type
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02KDYNAMO-ELECTRIC MACHINES
    • H02K41/00Propulsion systems in which a rigid body is moved along a path due to dynamo-electric interaction between the body and a magnetic field travelling along the path
    • H02K41/02Linear motors; Sectional motors
    • H02K41/03Synchronous motors; Motors moving step by step; Reluctance motors
    • H02K41/031Synchronous motors; Motors moving step by step; Reluctance motors of the permanent magnet type
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/06Linear motors
    • H02P25/062Linear motors of the induction type
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/006Controlling linear motors
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P2205/00Indexing scheme relating to controlling arrangements characterised by the control loops
    • H02P2205/05Torque loop, i.e. comparison of the motor torque with a torque reference

Definitions

  • the present invention relates to a linear motor in which a moving body moves along a track member by a thrust generated between a coil and a magnet through which a current is passed. More specifically, the present invention relates to a method and apparatus for deriving a thrust constant that varies depending on the position of the moving body, and a method and apparatus for controlling the movement of the moving body using the derived thrust constant.
  • the thrust constant is an index representing a generation ratio of thrust with respect to a current supplied to the coil, and its unit is represented by N / A (Newton / Ampere).
  • N / A Newton / Ampere
  • the multi-phase linear motor driving device has a reference table that stores information on the amount of thrust change at each position of the mover (moving body), and refers to the table that sequentially references the reference table during driving. It is characterized by having a means.
  • claim 2 of Patent Document 1 discloses a mode in which a linear motor is driven at a substantially constant speed and a reference table is created based on a time series of state quantities obtained at that time. As a result, the linear motor can be driven while sequentially referring to the information on the thrust change at each position in the drive direction, and the drive control is performed with stable and sufficiently reduced control deviation.
  • Patent Documents 2 to 4 have commonality in that the thrust ripple is expressed by a Fourier series.
  • a thrust ripple measuring device for a linear motor disclosed in Patent Literature 2 includes a phase calculator that calculates a phase of a Fourier fundamental wave based on the position of the linear motor, and a parameter determiner that determines a ripple parameter based on the phase and the thrust command. It consists of Thereby, it is said that the thrust ripple can be measured with high accuracy by controlling the linear motor to move at a constant speed and calculating the ripple parameter based on the thrust command and the position of the linear motor. Further, the thrust ripple compensation devices of Patent Literature 3 and Patent Literature 4 disclose devices that perform movement control using the thrust ripple measured in Patent Literature 2.
  • the amount of change in the thrust constant is derived based on the actual measurement result when the moving body is moved at a constant speed, and can cope with individual differences and position dependency of the linear motor. This is preferable.
  • derivation based only on actual measurement at a certain constant speed does not always provide a thrust constant with sufficient accuracy.
  • the moving body is moved at a constant acceleration, and the relationship between the thrust generated at that time and the current flowing in the coil is obtained.
  • the frictional resistance between the moving body and the track member becomes an error factor in deriving the thrust constant depending on the position under any actual measurement conditions of constant speed and constant acceleration.
  • the magnitude of the frictional resistance often changes nonlinearly depending on the speed of the moving body. Therefore, it is preferable to quantitatively grasp the influence of the frictional resistance and reflect it in the derivation of the thrust constant.
  • the present invention has been made in view of the above-described problems of the background art, and is a linear motor thrust constant derivation method and thrust constant derivation in which measurement is performed under a plurality of actual measurement conditions and the thrust constant derivation accuracy is improved as compared with the prior art.
  • Providing a device is a problem to be solved. Furthermore, it is an object to be solved to provide a linear motor movement control method and a movement control device that use a derived thrust constant to improve the accuracy of movement control of a moving body as compared with the conventional technique.
  • the invention of the linear motor thrust constant deriving method according to claim 1 for solving the above-mentioned problems is provided with a track member extending in the moving direction and having one of a magnet and a coil, and movably mounted on the track member.
  • a linear motor that includes a moving body having the other of the magnet and the coil, and generates a thrust in a moving direction between the magnet and a current flowing through the coil, and represents a generation ratio of the thrust with respect to the current
  • a method for deriving a thrust constant wherein an average thrust constant, which is an average thrust constant in a long movement section on the track member, is derived based on actual measurement, and a plurality of local thrusts on the track member are derived.
  • a local thrust constant, which is a constant is derived based on actual measurements, and a thrust constant at an arbitrary position on the track member is derived based on the average thrust constant and the local thrust constant at each location.
  • an average thrust constant can be derived in a long moving section on the track member, and local thrust constants can be derived at a plurality of locations on the track member.
  • the frictional resistance between a mobile body and a track member can be grasped
  • a linear motor movement control method is the thrust at an arbitrary position on the track member derived by the thrust constant deriving method of the linear motor device according to any one of the first to sixth aspects.
  • a thrust constant storage step for storing a thrust constant table in which constants and position coordinates are associated; and when the moving body is controlled to move, a current position coordinate on the track member of the moving body is detected, and the thrust constant
  • a thrust control step for controlling a current supplied to the coil using a thrust constant corresponding to the current position coordinates of the table.
  • the position dependence of the thrust constant with higher accuracy than the prior art is stored in the thrust constant table, and movement control can be performed using the thrust constant corresponding to the current position coordinates of the moving body. Therefore, the accuracy of movement control of the moving body is improved as compared with the conventional case.
  • An invention of a thrust constant deriving device for a linear motor according to claim 8 includes a track member extending in a moving direction and having one of a magnet and a coil, and the magnet and the coil mounted on the track member so as to be movable.
  • a thrust constant representing a generation ratio of the thrust with respect to the current is derived.
  • a device for deriving an average thrust constant which is an average thrust constant within a long movement section on the track member, based on actual measurement, and local thrust constants at a plurality of locations on the track member
  • the invention of the method for deriving the thrust constant of the linear motor according to claim 1 can be implemented as the invention of the device, and the effect of the invention of the device is the same as that of claim 1.
  • the linear motor movement control device is the linear motor thrust constant deriving method according to any one of claims 1 to 6, or the linear motor thrust constant deriving device according to claim 8.
  • a thrust constant storage means for storing a thrust constant table in which a thrust constant at an arbitrary position on the track member and position coordinates associated with each other are stored, and when controlling the movement of the moving body,
  • Thrust control means for detecting a current position coordinate on the track member and controlling a current supplied to the coil by using a thrust constant corresponding to the current position coordinate of the thrust constant table.
  • the invention of the linear motor movement control method according to claim 7 can also be implemented as an invention of the device, and the invention of the device has the same effect as that of claim 7.
  • FIG. 15 is a diagram of an example of an actual measurement waveform of a position coordinate of a moving object corresponding to the alternating current of FIG.
  • FIG. 1 is a plan view illustrating a basic configuration example of the linear motor 1.
  • the linear motor 1 includes a track member 2 and a moving body 4.
  • the track member 2 has two long plates 21 and 22 arranged in parallel.
  • the long plates 21 and 22 extend in the left-right direction of FIG.
  • a plurality of permanent magnets 3 are arranged in the moving direction at a substantially constant arrangement pitch P on the facing surfaces of the long plates 21 and 22.
  • FIG. 1 ten permanent magnets 3 are shown, and a larger number of permanent magnets 3 are arranged in a row before and after the omitted moving direction.
  • the moving body 4 is formed in a substantially rectangular parallelepiped shape that is long in the moving direction.
  • An unillustrated mounting portion of the moving body 4 is engaged with an unillustrated mounting portion of the track member 2.
  • the moving body 4 is mounted so as to be movable between the two long plates 21 and 22 of the track member 2.
  • the movable body 4 and the permanent magnets 3 on both sides of the track member 2 are separated by an approximately equal gap length G.
  • the moving body 4 has three coils 5 that are spaced apart from each other by an arrangement pitch P in the moving direction. Each coil 5 moves relatively while facing the permanent magnets 3 on both sides.
  • the moving body 4 further includes a position detection unit, a movement control unit, and a power supply unit (not shown).
  • the position detection unit detects a current position coordinate Xnow that represents the current position of the moving body 4 on the track member 1.
  • the movement control unit variably controls the magnitude and direction of the current I flowing from the power supply unit to the coil 5 based on the detected current position coordinates Xnow and the movement command. Thereby, a thrust F in the moving direction is generated between the current magnetic field formed by the coil 5 of the moving body 4 and the permanent magnetic field formed by the permanent magnet 3 of the track member 2.
  • the moving body 4 may have four or more coils 5.
  • the coil 5 may be arranged in the track member 2, and the moving body 4 may have the permanent magnet 3.
  • the magnitude of the thrust F varies depending on the strength of the permanent magnet 3 and the magnitude of the current I flowing through the coil 5.
  • the thrust constant N (unit: N / A (Newton / Ampere)) represents the generation ratio of the thrust F to the current I.
  • the thrust F generated between the track member 2 and the moving body 4 when the current of 1 A is supplied to the coil 5 becomes the thrust constant N.
  • the thrust constant N varies depending on the position on the track member 2 due to the change factors of the permanent magnet 3 and the track member 2 illustrated in FIGS. Further, since there is actually a frictional resistance R, it is not possible to use all of the thrust F for accelerating the moving body 4.
  • FIG. 1 also serves to schematically explain individual differences in the strength of the permanent magnet 3 that is a factor for changing the thrust constant N.
  • the ten permanent magnets 3 illustrated in FIG. 1 have individual differences in the strength of the magnetic force indicated by parenthesized percentage values. In the example of FIG. 1, when the average value of the strength of the magnetic force is ⁇ 0%, the individual difference varies in the range of ⁇ 3% to + 3%.
  • the thrust constant N increases near the relatively strong permanent magnet 3 and decreases near the relatively weak permanent magnet 3.
  • FIG. 2 is a diagram for schematically explaining the individual difference in the thickness of the permanent magnet 3 which is a change factor of the thrust constant N.
  • the thickness T ⁇ b> 2 of the second permanent magnet 32 is larger than the thickness T ⁇ b> 1 of the first permanent magnet 31.
  • the gap length G2 between the second permanent magnet 32 and the moving body 4 is smaller than the gap length G1 between the first permanent magnet 31 and the moving body 4. Therefore, the thrust constant N decreases in the vicinity of the first permanent magnet 31 having the large gap length G1, and increases in the vicinity of the second permanent magnet 32 having the small gap length G2.
  • FIG. 3 is a diagram schematically illustrating an individual difference in the adhesive thickness of the permanent magnet 3 that is a change factor of the thrust constant N.
  • the thickness of the adhesive 36 that fixes the fourth permanent magnet 35 to the long plate 21 as compared to the thickness T3 of the adhesive 34 that fixes the third permanent magnet 33 to the long plate 21. T4 is larger.
  • the thrust constant N decreases in the vicinity of the third permanent magnet 33 where the thickness T3 of the adhesive 34 is small, and increases in the vicinity of the fourth permanent magnet 35 where the thickness T4 of the adhesive 36 is large.
  • FIG. 4 is a diagram for schematically explaining the parallelism of the long rails 21 and 22 which is a change factor of the thrust constant N.
  • the distance between the two long rails 21 and 22 changes at one end and the other end. Accordingly, the gap lengths G5 and G6 between the permanent magnet 3 and the moving body 4 change depending on the position, and the thrust constant N changes corresponding to the gap lengths G5 and G6.
  • FIG. 5 is a diagram schematically illustrating the straightness of the long rails 21 and 22 that is a factor for changing the thrust constant N.
  • the distance between the two long rails 21 and 22 varies depending on the position of the two long rails 21 and 22. Accordingly, the gap lengths G7 and G8 between the permanent magnet 3 and the moving body 4 change depending on the position, and the thrust constant N changes corresponding to the gap lengths G7 and G8.
  • the thrust constant N is actually related to individual differences such as the strength and size of the plurality of permanent magnets 3 facing the three coils 5 of the moving body 4 and the gap length G.
  • Newton's equation of motion (second law) is expressed as “the thrust acting on the object is equal to the product of the mass of the object and the generated acceleration”.
  • the entire thrust F is not used for accelerating the moving body 4, but a part thereof is used to compensate against the frictional resistance R between the partial track member 2 and the moving body 4. .
  • the friction resistance R often changes nonlinearly depending on the speed of the moving body 4 and the like, making it difficult to derive an accurate thrust constant N.
  • FIG. 6 is a flowchart illustrating a method for deriving a thrust constant of the linear motor according to the first embodiment.
  • an average value deriving step S1 a local value deriving step S2, a local value correcting step S3, and an interpolation deriving step S4 are sequentially performed.
  • FIG. 7 is a diagram schematically illustrating a method for deriving a thrust constant of the linear motor according to the first embodiment.
  • the horizontal axis represents the position coordinate X on the track member 2
  • the vertical axis represents the thrust constant N.
  • the moving body 4 is moved over a long moving section on the track member 2, and the average thrust constant Nav1 is derived based on actual measurement.
  • an average thrust constant Nav1 representing an average thrust constant in the movement section from which the influence of the frictional resistance R is removed can be derived.
  • FIG. 7 plots the six local thrust constants N1 to N6 of the position coordinates X1 to X6 derived by executing the detailed procedure described in detail later.
  • the local thrust constants N1 to N6 are affected by the local change factors described with reference to FIGS. 1 to 5 and change depending on the position coordinates X1 to X6.
  • an average value Nav2 of the six local thrust constants N1 to N6 is calculated, and a difference ⁇ N from the average thrust constant Nav1 is calculated.
  • the difference ⁇ N illustrated in FIG. 7 can be considered as an error amount of the local thrust constants N1 to N6 that are generally caused by the influence of the frictional resistance R. Accordingly, the six local thrust constants N1 to N6 are corrected by the difference ⁇ N to obtain corrected thrust constants N1c to N6c at the respective locations.
  • N1c N1- ⁇ N
  • N2c N2- ⁇ N
  • ..., N6c N6- ⁇ N.
  • the thrust constants N (X) of arbitrary position coordinates X on the track member 2 are derived by interpolating the six corrected thrust constants N1c to N6c.
  • various known interpolation formulas such as a high-order polynomial can be used.
  • the thrust constant N (X) is represented by a smooth interpolation curve that passes through six corrected thrust constants N1c to N6c in the example of FIG.
  • FIG. 8 is a flowchart for explaining the detailed procedure of the average value deriving step S1.
  • a constant speed actual measurement procedure J1 a constant acceleration actual measurement procedure J2, and an average value calculation procedure J3 are sequentially performed.
  • FIG. 9 is a diagram illustrating a moving section in which the moving body 4 is moved when actual measurement is performed in the average value deriving step S1.
  • the moving section is long.
  • the moving section may be a full length section L in which the moving body 4 can move on the track member 2 or a section obtained by subtracting the minimum start acceleration section Ls and the stop deceleration section Le from the full length section L. preferable.
  • FIG. 10 is an example of measured waveforms of the friction compensation current If and the speed V in the constant speed actual measurement procedure J1 in the average value deriving step S1.
  • the horizontal axis represents time t
  • the vertical axis represents the current I flowing through the coil 5 and the speed V of the moving body 4.
  • the direction of the current I flowing through the coil 5 is controlled to be switched. Therefore, the current I that generates the thrust F in the forward direction of the moving body 4 is expressed as a positive value, and the current I that generates the thrust F in the backward direction of the moving body 4 is expressed as a negative value.
  • the current I in the time period from the time t2 to the time t3 maintaining the speed V is the friction compensation current If. If the friction resistance R does not exist, the friction compensation current If becomes zero. However, since the friction resistance R actually exists, a certain friction compensation current If is required. As illustrated in FIG. 9, the friction compensation current If becomes substantially constant within the movement section. If there is still some variation, the average current value If (avr) illustrated by the one-dot chain line in FIG. 9 is set as the friction compensation current If.
  • the constant speed V is set to a plurality (n types) of speeds Vf1, Vf2,..., Vfn, and the friction compensation currents If1, If2,.
  • the friction compensation currents If1, If2,..., Ifn vary depending on the speeds Vf1, Vf2,.
  • the obtained actual measurement results are stored in a memory or the like as a friction compensation table in which the speeds Vf1, Vf2,..., Vfn and the friction compensation currents If1, If2,. Further, in the friction compensation table, a continuous relationship can be obtained in advance by interpolating a discrete relationship between the speeds Vf1, Vf2,..., Vf and the friction compensation currents If1, If2,. preferable.
  • the moving body 4 is moved at a constant acceleration A over the moving section, and the required acceleration current Ia and the speed V of the moving body 4 required at that time are actually measured.
  • This actual measurement can be performed continuously within the movement section, or can be performed at multiple points even if it is discrete.
  • the required acceleration current Ia and the speed V change according to the position coordinate X in the movement section. Therefore, the actual measurement result can be represented by a graph in which the horizontal axis is time t or the position coordinate X, and the vertical axis is acceleration required current Ia and speed V.
  • a part of the actually measured acceleration required current Ia is used for generating a thrust F for accelerating the moving body 4, and the remaining part is used for compensating the frictional resistance R against it.
  • FIG. 11 shows an example of a measured waveform of the acceleration required current Ia in the constant acceleration actual measurement procedure J2 in the average value deriving step S1, and a waveform example of the acceleration net current Im calculated in the average value calculation procedure J3.
  • the horizontal axis represents time t
  • the vertical axis represents acceleration required current Ia and acceleration net current Im.
  • the constant acceleration A of the moving body 4 is maintained in the first half of the moving section.
  • the acceleration required current Ia during this period is a positive value and gradually increases with some variation.
  • the constant deceleration -A of the moving body 4 is maintained in the latter half of the moving section.
  • the acceleration required current Ia during this period is a negative value and gradually increases in absolute value with some variation.
  • the subsequent arithmetic processing is performed continuously or at discrete points in the movement section.
  • the friction compensation current If corresponding to the velocity V of the moving body 4 measured in the constant acceleration actual measurement procedure J2 is acquired from the friction compensation table.
  • the acceleration net current Im is calculated by subtracting the friction compensation current If acquired from the acceleration required current Ia. Therefore, the acceleration net current Im is a net current used for generating the thrust F.
  • the acceleration net current Im has a substantially constant value in each of the first half and the second half of the movement section.
  • the average thrust constant Nav1 is calculated based on the thrust F obtained from the acceleration A and the deceleration -A and the calculated acceleration net current Im using Newton's equation of motion (second law).
  • the average thrust constant Nav1 is substantially constant within the movement section.
  • the average acceleration net current in the movement section is adopted to calculate the average thrust constant Nav1.
  • FIG. 12 is a flowchart for explaining the detailed procedure of the local value deriving step S2.
  • the vibration measurement procedure J12, the transfer function derivation procedure J13, the model fitting procedure J14, and the local value calculation procedure J15 are performed at a plurality of measurement locations on the track member 2, respectively.
  • a plurality of actually measured points are set in the procedure J11 of FIG.
  • FIG. 13 is a diagram for explaining an example of setting a plurality of actual measurement locations on the track member 2 for deriving the local thrust constants N1 to N6.
  • three coils are arranged in a range of a predetermined length LM in the moving direction of the moving body 4.
  • the actual measurement points for deriving the local thrust constants N1 to N6 are set on the track member 2 at intervals LD smaller than the predetermined length LM.
  • FIG. 13 shows that six actually measured positions corresponding to the position coordinates X1 to X6 in FIG. 7 are set at the interval LD.
  • the intervals LD do not necessarily have to be equal. Further, an interval narrower than the interval LD shown in FIG. 13 may be set.
  • the interval LD is set wider than the predetermined length LM, the permanent magnet 3 that does not face the coil 5 is generated in the local value deriving step S2 and becomes a blind spot. In this case, even if the change factor of the permanent magnet 3 that does not face the coil 5 is extremely small, it is not reflected in the local thrust constants N1 to N6.
  • the moving body 4 is moved to the measurement location of the set position coordinate X1.
  • the AC current Iac and the position coordinate X of the moving body 4 are measured while energizing the coil 5 with an AC current Iac having a variable frequency f and vibrating the moving body 4 back and forth.
  • a high-speed sine wave sweep method is used in which the alternating current Iac energizing the coil 5 is a sine wave current and the frequency f is gradually changed over time. Thereby, the measured waveform examples shown in FIGS. 14 and 15 are obtained.
  • FIG. 14 is a diagram of an example of a measured waveform of the alternating current Iac when the high-speed sine wave sweep method is used in the excitation measurement procedure J12 in the local value derivation step S2.
  • FIG. 15 is a diagram of an example of an actually measured waveform of the position coordinate X of the moving body 4 corresponding to the alternating current Iac of FIG. 14 and 15 are measured with a common sweep time.
  • the frequency f gradually increases with time, and has a substantially constant amplitude except for the initial low frequency region.
  • the vibration waveform of the position coordinate X of the moving body is actually measured.
  • the vibration waveform indicates that the moving body 4 is vibrated minutely by being vibrated before and after the moving direction.
  • a transfer function G (f) in the frequency domain is derived based on the temporal change of the alternating current Iac and the temporal change of the position coordinate X of the moving body 4.
  • a time-dependent change in the position coordinate X of the moving body 4 shown in FIG. 15 is Fourier transformed to derive a frequency domain transfer function G (f).
  • a gain function and a phase function are obtained as a transfer function G (f) by Fourier transform, but only the gain function is used.
  • the gain function represents a characteristic in which the amplitude of the minute vibration of the moving body 4 depends on the frequency f.
  • FIG. 16 is a diagram illustrating the frequency domain transfer function G (f) derived in the transfer function deriving procedure J13 in the local value deriving step S2.
  • the horizontal axis indicates the frequency f
  • the vertical axis indicates the gain g (decibel value).
  • the gain g decreases almost smoothly as the frequency f increases, and the decreasing tendency slows down in a relatively high frequency band from a certain frequency f1 to the frequency f2.
  • the moving body 4 is modeled by fitting a physical model in which the mass M of the moving body 4 and the friction condition during movement are variably modeled to the transfer function G (f) in the frequency domain.
  • Estimated mass M is estimated.
  • a model using a Laplace transform method represented by the following equation is used as a physical model.
  • Is a constant K representing The nonlinearity of the frictional resistance R is modeled by the constant D and the constant K.
  • FIG. 17 is a diagram illustrating an example of the transfer function Gm (s) adjusted in the model matching procedure J14 in the local value derivation step S2.
  • the modeled transfer function Gm (s) is indicated by a broken line and is adjusted to the transfer function G (f).
  • the correction magnification C1 is obtained by dividing the measured mass M0 or the design mass M0 of the moving body 4 by the modeled mass M.
  • the modeled mass M is smaller than the actually measured mass M0, which means that the amount of movement larger than the design value has been actually measured.
  • the thrust F and the local thrust constant N1 are larger than the design values.
  • the fact that the modeled mass M is larger than the actually measured mass M0 means that a movement amount smaller than the design value has been actually measured. That is, it means that the thrust F and the local thrust constant N1 are smaller than the design values. Therefore, the idea of the correction magnification C1 is appropriate.
  • the local thrust constant N1 of the position coordinate X1 can be calculated by multiplying the correction magnification C1 by the thrust constant design value Ndn.
  • next procedure J16 it is determined whether or not it is the last actual measurement location. If it is not the last, it moves to the next actual measurement location and returns to the excitation actual measurement procedure J12.
  • the vibration measurement procedure J12 to the procedure J16 are repeated at six positions of the position coordinates X1 to X6. Thereby, the correction magnifications C1 to C6 and the local thrust constants N1 to N6 can be calculated for the respective position coordinates X1 to X6.
  • step J16 if it is the last measured location, the local value detection step S2 is terminated.
  • a method other than the high-speed sine wave sweep method can also be used in the vibration measurement procedure J12.
  • the amplitude of the minute vibration of the moving body 4 at each frequency f may be measured by changing the frequency f of the alternating current Iac stepwise by the step sweep method.
  • the transfer function derivation procedure J13 the relationship between the frequency f and the amplitude of the minute vibration is obtained instead of the Fourier transform.
  • the local value deriving step S2 can be modified as appropriate.
  • the thrust constant N (X) of the arbitrary position coordinate X on the track member 2 derived in the first embodiment can be used in the linear motor movement control method of the second embodiment.
  • the method for deriving the thrust constant of the linear motor of the first embodiment can be applied to a product inspection method when the linear motor 1 is manufactured. Specifically, for example, the pass / fail of the performance of the linear motor 1 can be determined by providing a limit determination value in the change range of the thrust constant N (X) derived by actual measurement.
  • the linear motor thrust constant deriving method includes a track member 2 having a permanent magnet 3 extending in a moving direction, and a moving body 4 having a coil 5 mounted on the track member 2 so as to be movable.
  • the thrust constant N representing the generation ratio of the thrust F with respect to the current I is derived.
  • an average thrust constant Nav1 which is an average thrust constant within a long movement section on the track member 2, is derived based on actual measurement, and the local positions at a plurality of locations (position coordinates X1 to X6) on the track member 2 are derived.
  • Local thrust constants N1 to N6 which are thrust constants, are derived based on actual measurements, respectively, and an arbitrary position (position coordinate X) on the track member 2 is calculated based on the average thrust constant Nav1 and the local thrust constants N1 to N6 at each location.
  • Thrust constant N ( ) To derive.
  • the average thrust constant Nav1 can be derived in a long moving section on the track member 2, and local thrust constants N1 to N6 are derived at a plurality of locations on the track member 2, respectively. it can. And by comparing and collating both, the frictional resistance R between the mobile body 4 and the track member 2 can be grasped
  • a local value correction step S3 for calculating a difference ⁇ N between Nav2 and the average thrust constant Nav1, and correcting the local thrust constants N1 to N6 at each location by a difference ⁇ N to obtain corrected thrust constants N1c to N6c at each location
  • the local value deriving step S2 applies an AC current Iac having a variable frequency to the coil 5 and vibrates the moving body 4 back and forth in the moving direction. Based on the vibration measurement procedure J12 for actually measuring the current Iac and the position information (position coordinates X) of the moving body 4, the temporal change of the alternating current Iac and the time change of the position information (position coordinates X) of the moving body 4.
  • a correction value C1 to C6 is obtained by dividing by the modeled mass M, and a local value calculation procedure J15 for calculating the local thrust constants N1 to N6 by multiplying the correction magnifications C1 to C6 by the thrust constant design value Ndn; Are performed at a plurality of locations (position coordinates X1 to X6) on the track member 2, respectively.
  • each location (position The local thrust constants N1 to N6 of the coordinates X1 to X6) can be derived with high accuracy.
  • the linear motor thrust constant deriving method of the first embodiment is a high-speed sine wave in which the alternating current Iac energized to the coil 5 is a sine wave current and the frequency f is gradually changed in the vibration measurement procedure J12.
  • the transfer function derivation procedure J13 the temporal change of the position information (position coordinates X) of the moving body 4 is Fourier transformed to derive the frequency domain transfer function G (f).
  • the local thrust constants N1 to N6 can be derived efficiently in a short time by using the high-speed sine wave sweep method.
  • the average value deriving step S1 is a friction compensation current If1 required for compensating the frictional resistance R by moving the moving body 4 at a constant speed over the moving section.
  • ⁇ Ifn is actually measured at a plurality of speeds V1 to Vn, a constant speed measurement procedure J1 for storing a friction compensation table in which the speeds V1 to Vn and the friction compensation currents If1 to Ifn are associated with each other, and a moving section.
  • the friction compensation current If corresponding to the velocity V of the moving body 4 actually measured in step S1 is acquired from the friction compensation table, and the acceleration net current Im obtained by subtracting the friction compensation current If acquired from the acceleration required current Ia.
  • the acceleration net current Im obtained by subtracting the friction compensation current If acquired from the acceleration required current Ia.
  • the average thrust constant Nav1 from which the influence is removed can be derived with extremely high accuracy. Accordingly, the accuracy of deriving the thrust constant N (X) is further enhanced.
  • a plurality of coils 5 are arranged in a range of a predetermined length LM in the moving direction of the moving body 4, and a trajectory for deriving the local thrust constants N1 to N6.
  • a plurality of positions (position coordinates X1 to X6) on the member 2 are set at intervals LD smaller than the predetermined length LM.
  • FIG. 18 is a flowchart for explaining a linear motor movement control method according to the second embodiment. Arithmetic processing and control based on this flowchart are mainly handled by the movement control unit mounted on the moving body 4.
  • the linear motor thrust constant deriving method of the first embodiment is performed to derive a thrust constant N (X) of an arbitrary position coordinate X on the track member 2.
  • a thrust constant table in which the thrust constant N (X) is associated with the position coordinate X is stored in a memory or the like.
  • the thrust constant deriving step S11 and the thrust constant storing step S12 are performed only in a limited case such as immediately after manufacturing the linear motor 1, during repair, or during periodic inspection.
  • the target position setting step S13 and subsequent steps are always repeated when the linear motor 1 is used.
  • target position coordinates Xtgt for moving the moving body 4 are set.
  • the target position coordinate Xtgt may be commanded from the outside, or the movement control unit may determine it by calculation processing.
  • the movement control unit acquires information on the current position coordinates Xnow of the moving body 4 from the position detection unit.
  • the condition of step S15 is not normally satisfied, and the routine proceeds to thrust constant acquisition step S16.
  • a thrust constant N (Xnow) corresponding to the current position coordinate Xnow of the thrust constant table is acquired.
  • the magnitude of the current Inow flowing through the coil 5 is calculated. This calculation can be performed by appropriately applying a known technique using target position coordinates Xtgt, current position coordinates Xnow, thrust constant N (Xnow), and the like as calculation parameters.
  • the calculated current Inow is controlled to actually flow. Thereafter, the process returns to the current position coordinate detection step S14.
  • the current position coordinate detection step S14 to the current control step S18 correspond to the thrust control step of the present invention.
  • the linear motor movement control method of the second embodiment includes a thrust constant N (X) of an arbitrary position coordinate X on the track member 2 derived by the thrust constant deriving method of the linear motor device described in the first embodiment.
  • a thrust constant storage step for storing a thrust constant table associated with the position coordinates X, and when the moving body 4 is controlled to move, a current position coordinate Xnow on the track member 2 of the moving body 4 is detected, and the thrust constant table
  • a thrust control step for controlling a current Inow to be passed through the coil 5 using a thrust constant N (Xnow) corresponding to the current position coordinate Xnow.
  • a thrust constant N (X) with higher accuracy than the prior art is stored in the thrust constant table, and movement control is performed using the thrust constant N (Xnow) corresponding to the current position coordinate Xnow of the moving body 4. Can be done. Therefore, the accuracy of the movement control of the moving body 4 is improved as compared with the prior art.
  • the linear motor thrust constant derivation method of the first embodiment and the linear motor movement control method of the second embodiment can both be implemented as devices.
  • the effect when implemented as an apparatus is the same as that of the first and second embodiments.
  • Various other applications and modifications are possible for the present invention.

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Electromagnetism (AREA)
  • Control Of Linear Motors (AREA)

Abstract

 本発明は、移動方向に延在して磁石を有する軌道部材と、軌道部材に移動可能に装架されてコイルを有する移動体とを備え、コイルに電流を通電して磁石との間に移動方向の推力を発生するリニアモータにおいて、電流に対する推力の発生比率を表す推力定数を導出する方法であって、軌道部材上の長い移動区間内の平均的な推力定数である平均推力定数を実測に基づいて導出し、軌道部材上の複数箇所の局所的な推力定数である局所推力定数をそれぞれ実測に基づいて導出し、平均推力定数および各箇所の局所推力定数に基づいて軌道部材上の任意の位置の推力定数を導出する。これによれば、複数の実測条件での実測に基づき推力定数の導出精度を従来よりも高めることができる。

Description

リニアモータの推力定数導出方法および移動制御方法、ならびに、リニアモータの推力定数導出装置および移動制御装置
 本発明は、電流を通電したコイルと磁石との間に発生する推力により移動体が軌道部材に沿って移動するリニアモータに関する。より詳細には、移動体の位置に依存して変化する推力定数を導出する方法および装置、ならびに、導出した推力定数を用いて移動体の移動を制御する方法および装置に関する。
 多数の電子部品が装着された基板を生産する設備として、半田印刷機、部品実装機、リフロー炉、基板検査機などがある。これらの設備を基板搬送装置で連結して基板生産ラインを構築することが一般的になっている。部品実装機や基板検査機では、実装ヘッドや検査ヘッドの駆動装置として、従来から送りねじ駆動装置が用いられてきた。近年では、ヘッドの高速移動および高精度位置制御への要求が高まり、駆動装置としてリニアモータが用いられるようになってきている。リニアモータ装置は、部品実装機や基板検査機の用途に限定されず、直進可動部を有する各種の産業機械に広く用いられている。
 この種のリニアモータでは、モータ特性の一項目である推力定数の位置依存性が制御特性に影響を及ぼすことが判っている。推力定数とは、コイルに通電された電流に対する推力の発生比率を表す指標であり、単位はN/A(ニュートン/アンペア)で表される。例えば、或るリニアモータ制御システムでは、軌道部材上の移動体の位置に依存して推力定数が10%程度低下すると、移動制御時に移動体が目標位置に到達するまでの所要時間が10ms程度長引くことが確認されている。
 推力定数が位置に依存して変化する変化要因として、個々の磁石の強さおよび大きさの個体差や、列設される磁石の配設位置のばらつきなどが考えられる。これらの変化要因は、製造工程での部品管理や組付作業管理を厳密化することでその影響を低減できる。しかしながら、管理の厳密化は、コストアップのデメリットに直結している。このような状況下で、推力定数の変化量を実測に基づいて導出し、移動制御時に制御パラメータを可変に調整して制御性能を高める技術が一般的に用いられている。推力定数を実測に基づいて導出し、移動制御に用いる技術例が特許文献1~4に開示されている。
 特許文献1の請求項1の多相リニアモータ駆動装置は、可動子(移動体)の各位置における推力変化分の情報を記憶した参照テーブルを有し、駆動時に参照テーブルを逐次参照するテーブル参照手段を具備することを特徴としている。さらに、特許文献1の請求項2には、リニアモータをほぼ一定速度で駆動させ、そのときに得られる状態量の時系列に基づいて参照テーブルを作成する態様が開示されている。これにより、駆動方向の各位置における推力変化分の情報を逐次参照しながらリニアモータを駆動でき、安定でかつ制御偏差が十分低減された駆動制御が行われる、とされている。
 また、特許文献2~4の技術は、推力リップルをフーリエ級数で表現する点で共通性を有している。特許文献2のリニアモータの推力リップル測定装置は、リニアモータの位置に基づいてフーリエ基本波の位相を算出する位相計算器と、前記位相および推力指令に基づいてリップルパラメータを決定するパラメータ決定器とで構成されている。これにより、リニアモータを一定速度で移動させるように制御し、推力指令およびリニアモータの位置に基づいてリップルパラメータを計算することで、高精度に推力リップルを測定できる、とされている。さらに、特許文献3および特許文献4の推力リップル補償装置は、特許文献2で測定された推力リップルを用いて移動制御を行う装置を開示している。
特開平9-65511号公報 特開2009-159741号公報 特開2009-159751号公報 特開2010-130854号公報
 ところで、特許文献1および特許文献2では、移動体を一定速度で移動させたときの実測結果に基づいて推力定数の変化分を導出しており、リニアモータの個体差および位置依存性に対応できる点で好ましい。しかしながら、或る一定速度の実測のみによる導出では、必ずしも十分な精度の推力定数が得られるとは限らない。また、一般的な推力定数の導出方法では、移動体を一定加速度で移動させ、そのときに発生している推力と、コイルに流れている電流との関係を求める。ここで、一定速度および一定加速度のいずれの実測条件においても、移動体と軌道部材との間の摩擦抵抗は、位置に依存する推力定数を導出する際の誤差要因になる。摩擦抵抗の大きさは、移動体の速度に依存して非線形に変化する場合が多い。したがって、摩擦抵抗の影響を定量的に把握して、推力定数の導出に反映することが好ましい。
 本発明は、上記背景技術の問題点に鑑みてなされたものであり、複数の実測条件での実測を行い推力定数の導出精度を従来よりも高めたリニアモータの推力定数導出方法および推力定数導出装置を提供することを解決すべき課題とする。さらには、導出した推力定数を用いることにより移動体の移動制御の精度を従来よりも高めたリニアモータの移動制御方法および移動制御装置を提供することを解決すべき課題とする。
 上記課題を解決する請求項1に係るリニアモータの推力定数導出方法の発明は、移動方向に延在して磁石およびコイルの一方を有する軌道部材と、前記軌道部材に移動可能に装架されて前記磁石および前記コイルの他方を有する移動体とを備え、前記コイルに電流を通電して前記磁石との間に移動方向の推力を発生するリニアモータにおいて、前記電流に対する前記推力の発生比率を表す推力定数を導出する方法であって、前記軌道部材上の長い移動区間内の平均的な推力定数である平均推力定数を実測に基づいて導出し、前記軌道部材上の複数箇所の局所的な推力定数である局所推力定数をそれぞれ実測に基づいて導出し、前記平均推力定数および各箇所の前記局所推力定数に基づいて前記軌道部材上の任意の位置の推力定数を導出する。
 これによれば、複数の実測条件で実測を行うことにより、軌道部材上の長い移動区間で平均推力定数を導出でき、軌道部材上の複数箇所でそれぞれ局所推力定数を導出できる。そして、両者を比較照合することにより移動体と軌道部材との間の摩擦抵抗を定量的に把握して、その影響を低減することができる。したがって、複数の実測条件での実測に基づく軌道部材上の任意の位置の推力定数は、単一の実測条件に基づく従来技術と比較して導出精度が高められる。
 請求項7に係るリニアモータの移動制御方法の発明は、請求項1~6のいずれか一項に記載のリニアモータ装置の推力定数導出方法によって導出された前記軌道部材上の任意の位置の推力定数と位置座標とを対応付けた推力定数テーブルを記憶する推力定数記憶ステップと、前記移動体を移動制御するときに、前記移動体の前記軌道部材上の現在位置座標を検出し、前記推力定数テーブルの前記現在位置座標に対応する推力定数を用いて前記コイルに通電する電流を制御する推力制御ステップと、を有する。
 これによれば、従来技術よりも精度の高い推力定数の位置依存性が推力定数テーブルに記憶されており、移動体の現在位置座標に対応する推力定数を用いて移動制御を行える。したがって、移動体の移動制御の精度が従来よりも高められる。
 請求項8に係るリニアモータの推力定数導出装置の発明は、移動方向に延在して磁石およびコイルの一方を有する軌道部材と、前記軌道部材に移動可能に装架されて前記磁石および前記コイルの他方を有する移動体とを備え、前記コイルに電流を通電して前記磁石との間に移動方向の推力を発生するリニアモータにおいて、前記電流に対する前記推力の発生比率を表す推力定数を導出する装置であって、前記軌道部材上の長い移動区間内の平均的な推力定数である平均推力定数を実測に基づいて導出する手段と、前記軌道部材上の複数箇所の局所的な推力定数である局所推力定数をそれぞれ実測に基づいて導出する手段と、前記平均推力定数および各箇所の前記局所推力定数に基づいて前記軌道部材上の任意の位置の推力定数を導出する手段と、を有する。
 これによれば、請求項1に係るリニアモータの推力定数導出方法の発明は装置の発明として実施することもでき、装置の発明でも請求項1と同様の効果が生じる。
 請求項9に係るリニアモータの移動制御装置の発明は、請求項1~6のいずれか一項に記載のリニアモータの推力定数導出方法、あるいは請求項8に記載のリニアモータの推力定数導出装置によって導出された前記軌道部材上の任意の位置の推力定数と位置座標とを対応付けた推力定数テーブルを記憶する推力定数記憶手段と、前記移動体を移動制御するときに、前記移動体の前記軌道部材上の現在位置座標を検出し、前記推力定数テーブルの前記現在位置座標に対応する推力定数を用いて前記コイルに通電する電流を制御する推力制御手段と、を有する。
 これによれば、請求項7に係るリニアモータの移動制御方法の発明は装置の発明として実施することもでき、装置の発明でも請求項7と同様の効果が生じる。
リニアモータの基本構成例を説明した平面図であり、推力定数の変化要因である永久磁石の強さの個体差を模式的に説明する図を兼ねている。 推力定数の変化要因である永久磁石の厚みの個体差を模式的に説明する図である。 推力定数の変化要因である永久磁石の接着剤厚みの個体差を模式的に説明する図である。 推力定数の変化要因である長尺レールの平行度を模式的に説明する図である。 推力定数の変化要因である長尺レールの直線度を模式的に説明する図である。 第1実施形態のリニアモータの推力定数導出方法を説明するフロー図である。 第1実施形態のリニアモータの推力定数導出方法を図式的に説明する図である。 平均値導出ステップの詳細手順を説明するフロー図である。 平均値導出ステップで実測を行うときに、移動体を移動させる移動区間を説明する図である。 平均値導出ステップ中の定速実測手順における摩擦補償電流および速度の実測波形例である。 平均値導出ステップ中の一定加速度実測手順における加速所要電流の実測波形例、ならびに平均値算出手順で算出した加速正味電流の波形例である。 局所値導出ステップの詳細手順を説明するフロー図である。 局所推力定数を導出する軌道部材上の複数の実測箇所の設定例を説明する図である。 局所値導出ステップ中の加振実測手順で高速正弦波掃引法を用いたときの交流電流の実測波形例の図である。 加振実測手順において、図14の交流電流に対応した移動体の位置座標の実測波形例の図である。 局所値導出ステップ中の伝達関数導出手順で導出された周波数領域の伝達関数を例示した図である。 局所値導出ステップ中のモデル合わせ込み手順で、合わせ込まれた伝達関数の例を示した図である。 第2実施形態のリニアモータの移動制御方法を説明するフロー図である。
 本発明の第1実施形態のリニアモータの推力定数導出方法について、図1~図17を参考にして説明する。まず、第1実施形態に用いるリニアモータ1の基本構成例について説明する。図1は、リニアモータ1の基本構成例を説明した平面図である。リニアモータ1は、軌道部材2および移動体4からなる。
 軌道部材2は、平行配置された2本の長尺プレート21、22を有している。長尺プレート21、22は、移動体4の移動方向となる図2の左右方向に延在している。各長尺プレート21、22の向かい合う面には、複数の永久磁石3が概ね一定の配設ピッチPで移動方向に列設されている。図1には10個の永久磁石3が示されており、省略された移動方向前後には、さらに多数の永久磁石3が列設されている。
 移動体4は、移動方向に長い略直方体形状に形成されている。移動体4の図略の被装架部は、軌道部材2の図略の装架部に係合されている。これにより、移動体4は、軌道部材2の2枚の長尺プレート21、22の間に移動可能に装架される。移動体4と軌道部材2の両側の永久磁石3との間は、概ね等しいギャップ長Gだけ離れている。移動体4は、移動方向に概ね配設ピッチPだけ離隔して配置された3個のコイル5を有している。各コイル5は、両側の永久磁石3に対向しながら相対移動する。
 移動体4は、さらに、図略の位置検出部、移動制御部、および電源部を搭載している。位置検出部は、軌道部材1上における移動体4の現在位置を表す現在位置座標Xnowを検出する。移動制御部は、検出された現在位置座標Xnowおよび移動指令に基づいて、電源部からコイル5に流れる電流Iの大きさおよび方向を可変に制御する。これにより、移動体4のコイル5が形成する電流磁界と、軌道部材2の永久磁石3が形成する永久磁界との間に、移動方向の推力Fが発生する。
 軌道部材2および移動体4の構成は、上述に限定されず、公知の各種技術を応用できる。例えば、移動体4が4個以上のコイル5を有していてもよい。また例えば、軌道部材2にコイル5が列設され、移動体4が永久磁石3を有していてもよい。いずれの構成においても、推力Fの大きさは、永久磁石3の強さおよびコイル5に流れる電流Iの大きさに依存して変化する。
 ここで、推力定数N(単位:N/A(ニュートン/アンペア))は、電流Iに対する推力Fの発生比率を表す。換言すれば、コイル5に1Aの電流を通電したときに軌道部材2と移動体4との間に発生する推力Fが推力定数Nとなる。各永久磁石3が等しい永久磁界を生成し、かつ軌道部材2と移動体4との間に摩擦抵抗Rが無い理想条件において、推力Fは、現在位置座標Xnowに依存しない一定値になる。しかしながら、実際には図1~図5に例示される永久磁石3および軌道部材2の変化要因により、推力定数Nは軌道部材2上の位置に依存して変化する。また、実際には摩擦抵抗Rが有るので、推力Fの全てを移動体4の加速に使用することはできない。
 図1は、推力定数Nの変化要因である永久磁石3の強さの個体差を模式的に説明する図を兼ねている。図1に例示された10個の永久磁石3は、かっこ付きのパーセント値で示された磁力の強さの個体差を有している。図1の例で、磁力の強さの平均値を±0%としたときに、個体差はー3%~+3%の範囲にばらついている。推力定数Nは、相対的に磁力の強い永久磁石3の近傍で増加し、相対的に磁力の弱い永久磁石3の近傍で減少する。
 図2は、推力定数Nの変化要因である永久磁石3の厚みの個体差を模式的に説明する図である。図2の例では、第1の永久磁石31の厚みT1に比較して、第2の永久磁石32の厚みT2の方が大きい。このため、第1の永久磁石31と移動体4とのギャップ長G1よりも、第2の永久磁石32と移動体4とのギャップ長G2のほうが小さくなる。したがって、推力定数Nは、ギャップ長G1が大きな第1の永久磁石31の近傍で減少し、ギャップ長G2が小さな第2の永久磁石32の近傍で増加する。
 図3は、推力定数Nの変化要因である永久磁石3の接着剤厚みの個体差を模式的に説明する図である。図3の例では、第3の永久磁石33を長尺プレート21に固定する接着剤34の厚みT3に比較して、第4の永久磁石35を長尺プレート21に固定する接着剤36の厚みT4の方が大きい。このため、2個の永久磁石33、35の厚みが同じであっても、第3の永久磁石33側のギャップ長G3よりも、第4の永久磁石35側のギャップ長G4のほうが小さくなる。したがって、推力定数Nは、接着剤34の厚みT3が小さな第3の永久磁石33の近傍で減少し、接着剤36の厚みT4が大きな第4の永久磁石35の近傍で増加する。
 図4は、推力定数Nの変化要因である長尺レール21、22の平行度を模式的に説明する図である。通常、2本の長尺レール21、22の組み付け精度に誤差が有るため、平行度が不足するおそれが皆無ではない。この場合、図4に示されるように、2本の長尺レール21、22の一端および他端で相互間距離が変化する。これにより、位置に依存して永久磁石3と移動体4とのギャップ長G5、G6が変化し、ギャップ長G5、G6に対応して推力定数Nが変化する。
 図5は、推力定数Nの変化要因である長尺レール21、22の直線度を模式的に説明する図である。通常、2本の長尺レール21、22の加工精度に誤差が有るため、少なくとも1本の長尺レール21、22で直線度が不足するおそれが皆無ではない。この場合、図5に示されるように、2本の長尺レール21、22の位置によって相互間距離が変化する。これにより、位置に依存して永久磁石3と移動体4とのギャップ長G7、G8が変化し、ギャップ長G7、G8に対応して推力定数Nが変化する。
 図1~図5において、実際には、移動体4の3個のコイル5が対向する複数の永久磁石3の強さや大きさなどの個体差、およびギャップ長Gが関連しあって推力定数Nを変化させる。実測に基づいて推力定数Nを導出するために、「物体に作用している推力は、物体の質量と発生している加速度との積に等しい。」というニュートンの運動方程式(第2法則)を利用する。ここで、推力Fの全部は移動体4の加速に使用されずに、一部は部軌道部材2と移動体4との間の摩擦抵抗Rに対抗してこれを補償するために使用される。摩擦抵抗Rは、その大きさが移動体4の速度などに依存して非線形に変化する場合が多く、正確な推力定数Nの導出を難しくしている。
 そこで、第1実施形態のリニアモータの推力定数導出方法では、2種類の実測条件で実測を行い、実測結果を組み合わせて軌道部材2上の任意の位置の推力定数Nを精度よく導出する。図6は、第1実施形態のリニアモータの推力定数導出方法を説明するフロー図である。第1実施形態のリニアモータの推力定数導出方法は、平均値導出ステップS1、局所値導出ステップS2、局所値補正ステップS3、および補間導出ステップS4を順番に実施する。第1実施形態のリニアモータの推力定数導出方法は、移動体4の位置検出部および移動制御部に外部のディジタル計測装置やデータ処理装置などを適宜連携させて実施することができる。図7は、第1実施形態のリニアモータの推力定数導出方法を図式的に説明する図である。図7で、横軸は軌道部材2上の位置座標X、縦軸は推力定数Nを表している。
 まず、図6および図7を参考にして、推力定数導出方法の全体手順について説明する。図6の平均値導出ステップS1で、軌道部材2上の長い移動区間に渡って移動体4を移動させ、平均推力定数Nav1を実測に基づいて導出する。後に詳述する詳細手順を実施することより、摩擦抵抗Rの影響を除去した移動区間内の平均的な推力定数を表す平均推力定数Nav1を導出できる。
 次に、局所値導出ステップS2で、軌道部材2上の複数箇所で移動体4を移動方向に微小移動させて、各箇所の局所推力定数をそれぞれ導出する。図7には、後に詳述する詳細手順を実施して導出した位置座標X1~X6の6箇所の局所推力定数N1~N6がプロットされている。局所推力定数N1~N6は、図1~図5で説明した局所的な変化要因の影響を受け、位置座標X1~X6に依存して変化している。
 次に、局所値補正ステップS3で、6箇所の局所推力定数N1~N6の平均値Nav2を算出し、平均推力定数Nav1との差分ΔNを算出する。図7に例示される差分ΔNは、概ね摩擦抵抗Rの影響に起因する局所推力定数N1~N6の誤差分と考えることができる。したがって、6箇所の局所推力定数N1~N6を差分ΔNだけ補正して各箇所の補正推力定数N1c~N6cとする。図7の例では、N1c=N1-ΔN、N2c=N2-ΔN、……、N6c=N6-ΔNとする。
 次に、補間導出ステップS4で、6箇所の補正推力定数N1c~N6cを補間することにより、軌道部材2上の任意の位置座標Xの推力定数N(X)を導出する。補間の方法としては、公知の各種補間式、例えば高次多項式などを用いることができる。これにより、摩擦抵抗Rの影響を除外しつつ局所的な変化要因の影響を反映した精度の高い推力定数N(X)を導出できる。推力定数N(X)は、図7の例では6箇所の補正推力定数N1c~N6cを通る滑らかな補間曲線で表されている。
 次に、平均値導出ステップS1の詳細手順について説明する。図8は、平均値導出ステップS1の詳細手順を説明するフロー図である。平均値導出ステップS1では、定速実測手順J1、一定加速度実測手順J2、および平均値算出手順J3を順番に実施する。図9は、平均値導出ステップS1で実測を行うときに、移動体4を移動させる移動区間を説明する図である。平均値導出ステップS1で、移動区間は長いことが好ましい。換言すれば、移動区間は、軌道部材2上を移動体4が移動し得る全長区間L、もしくは、全長区間Lから最小限の始動加速区間Lsおよび停止減速区間Leを差し引いた区間とすることが好ましい。
 定速実測手順J1では、軌道部材2上の移動区間に渡って移動体4を一定の速度Vで移動させて摩擦抵抗Rの補償に要する摩擦補償電流Ifを実測する。図10は、平均値導出ステップS1中の定速実測手順J1における摩擦補償電流Ifおよび速度Vの実測波形例である。図10で、横軸は時間tを表し、縦軸はコイル5に流れる電流Iおよび移動体4の速度Vを表している。移動体4の位置座標Xに応じて、コイル5に流れる電流Iの方向は切り替え制御される。そこで、移動体4の前進方向の推力Fを発生する電流Iを正値で表し、移動体4の後進方向の推力Fを発生する電流Iを負値で表す。
 時刻t1に正値の電流Iが流れ始めると、移動体4は加速しながら始動加速区間Lsを移動する。始動加速区間Lsでは、加速のために大きな正値の電流Iが流れる。時刻t2に移動体4が所定の速度Vに達すると始動加速区間Lsは終了する。そして、以降は速度Vを一定に維持するように電流Iを可変に制御する。速度Vの維持は、時刻t3に移動体4が停止減速区間Leに差しかかるまで継続される。時刻t3以降には負値の大きな電流Iを通電して移動体4を制動し、時刻t4に移動体4は停止する。
 速度Vを維持している時刻t2から時刻t3までの時間帯の電流Iが摩擦補償電流Ifである。仮に摩擦抵抗Rが存在しなければ摩擦補償電流Ifはゼロとなるが、実際には摩擦抵抗Rが存在するので或る摩擦補償電流Ifが必要になる。摩擦補償電流Ifは、図9に例示されるように、移動区間内で概ね一定になる。それでも、多少のばらつきがある場合には、図9に一点鎖線で例示された平均電流値If(avr)を摩擦補償電流Ifとする。
 さらに、一定の速度Vを複数(n種類)の速度Vf1、Vf2、……、Vfnに設定し、それぞれについて摩擦補償電流If1、If2、……、Ifnを実測する。摩擦補償電流If1、If2、……、Ifnは、速度Vf1、Vf2、……、Vfnに依存して変化する。得られた実測結果は、速度Vf1、Vf2、……、Vfnと摩擦補償電流If1、If2、……、Ifnとを対応付けた摩擦補償テーブルにして、メモリなどに記憶する。さらに、摩擦補償テーブルでは、速度Vf1、Vf2、……、Vfと摩擦補償電流If1、If2、……、Ifnとの離散的な関係を補間して、予め連続的な関係を得ておくことが好ましい。
 次の一定加速度実測手順J2では、移動区間に渡って移動体4を一定の加速度Aで移動させ、そのときに必要であった加速所要電流Iaおよび移動体4の速度Vを実測する。この実測は、移動区間内で連続的に行え、あるいは離散的であっても多数点で行える。加速所要電流Iaおよび速度Vは移動区間内の位置座標Xに応じて変化する。したがって、実測結果は、例えば、横軸を時間tまたは位置座標Xとし、縦軸を加速所要電流Iaおよび速度Vとしたグラフで表すことができる。実測された加速所要電流Iaは、一部が移動体4を加速させる推力Fの発生に使用され、残部が摩擦抵抗Rに対抗してこれを補償するために使用される。
 一定加速度実測手順J2では、移動区間の前半で一定の加速度Aとし、移動区間の後半で異符号同大の減速度-Aとする実測条件を採用してもよい。この実測条件で、例えば、図11に例示される実測結果が得られる。図11は、平均値導出ステップS1中の一定加速度実測手順J2における加速所要電流Iaの実測波形例、ならびに平均値算出手順J3で算出した加速正味電流Imの波形例である。図11で、横軸は時間tを表し、縦軸は加速所要電流Iaおよび加速正味電流Imを表している。
 図11の実測波形例において、時刻t11から時刻t12までは、移動区間の前半で移動体4の一定の加速度Aを維持している。この間の加速所要電流Iaは、正値であって多少のばらつきを含みつつ徐々に増加している。また、時刻t12から時刻t13までは、移動区間の後半で移動体4の一定の減速度-Aを維持している。この間の加速所要電流Iaは、負値であって多少のばらつきを含みつつ絶対値が徐々に増加している。
 次の平均値算出手順J3では、移動区間内で連続的に、あるいは離散的な多数点で以降の演算処理を行う。まず、一定加速実測手順J2で実測した移動体4の速度Vに対応する摩擦補償電流Ifを摩擦補償テーブルから取得する。2番目に、加速所要電流Iaから取得した摩擦補償電流Ifを減算して加速正味電流Imを算出する。したがって、加速正味電流Imは、推力Fの発生に使用される正味の電流となる。図11の例で、加速正味電流Imは、移動区間の前半および後半でそれぞれ概ね一定値となっている。3番目に、ニュートンの運動方程式(第2法則)を利用して加速度Aおよび減速度-Aから得られる推力Fと、算出した加速正味電流Imとに基づいて、平均推力定数Nav1を算出する。平均推力定数Nav1は、移動区間内で概ね一定になる。加速正味電流Imのばらつきがある場合には、移動区間内の平均的な加速正味電流を採用して、平均推力定数Nav1を算出する。
 なお、移動体4の往復の移動方向によって特性差が生じる場合には、往動方向と復動方向でそれぞれ実測を行い、往動方向と復動方向とで異なる平均推力定数Nav1を算出してもよい。さらに、定速実測手順J1および一定加速度実測手順J2は、それぞれ複数回の実測を行って平均値を求めることにより、導出精度を高めることもできる。その他にも、平均値導出ステップS1は、適宜変形できる。
 次に、局所値導出ステップS2の詳細手順について説明する。図12は、局所値導出ステップS2の詳細手順を説明するフロー図である。局所値導出ステップS2では、加振実測手順J12、伝達関数導出手順J13、モデル合わせ込み手順J14、および局所値算出手順J15を軌道部材2上の複数の実測箇所でそれぞれ実施する。まず、図12の手順J11で、複数の実測箇所を設定する。
 図13は、局所推力定数N1~N6を導出する軌道部材2上の複数の実測箇所の設定例を説明する図である。図示されるように、移動体4の移動方向の所定長LMの範囲に3個のコイルが列設されている。そして、局所推力定数N1~N6を導出する実測箇所を、所定長LMよりも狭い間隔LDで軌道部材2上に設定する。図13は、図7の位置座標X1~X6に相当する6つの実測箇所が間隔LDで設定されたことを示している。なお、間隔LDは、必ずしも等間隔である必要はない。また、図13に示される間隔LDよりもさらに狭い間隔を設定してもよい。
 このように間隔LDを或る程度狭く設定することで、全ての永久磁石3の変化要因の影響が局所推力定数N1~N6の少なくとも1量に反映される。これにより、推力定数Nの導出精度が確保される。仮に、間隔LDを所定長LMよりも広く設定すると、局所値導出ステップS2でコイル5に対向しない永久磁石3が生じて盲点となる。この場合、コイル5に対向しない永久磁石3の変化要因が著しく小さくとも局所推力定数N1~N6に反映されないので、推力定数Nの導出精度の著しい低下が懸念される。
 次の加振実測手順J12では、まず、設定された位置座標X1の実測箇所に移動体4を移動する。次に、コイル5に周波数f可変の交流電流Iacを通電して移動体4を移動方向前後に加振しながら、交流電流Iacおよび移動体4の位置座標Xを実測する。第1実施形態では、コイル5に通電する交流電流Iacを正弦波電流とし、かつ周波数fを時間的に徐変させる高速正弦波掃引法を用いる。これにより、図14および図15に示される実測波形例が得られる。
 図14は、局所値導出ステップS2中の加振実測手順J12で高速正弦波掃引法を用いたときの交流電流Iacの実測波形例の図である。また、図15は、加振実測手順J12において、図14の交流電流Iacに対応した移動体4の位置座標Xの実測波形例の図である。図14および図15は、共通の掃引時間で実測されている。図14に示される交流電流Iacの波形は、時間の経過とともに周波数fが徐々に増加しており、かつ初期の低周波領域を除いて概ね一定の振幅になっている。これに対応して、図15に示されるように、移動体の位置座標Xの振動波形が実測されている。振動波形は、移動体4が移動方向前後に加振されて微小振動していることを示している。
 次の伝達関数導出手順J13では、交流電流Iacの時間的変化および移動体4の位置座標Xの時間的変化に基づいて、周波数領域の伝達関数G(f)を導出する。第1実施形態では、前提条件として、図14に示された高速正弦波掃引法の交流電流Iacの波形が誤っていないことを確認する。その上で、図15に示された移動体4の位置座標Xの時間的変化をフーリエ変換して周波数領域の伝達関数G(f)を導出する。フーリエ変換により伝達関数G(f)としてゲイン関数および位相関数が得られるが、ゲイン関数のみを使用する。ゲイン関数は、移動体4の微小振動の振幅が周波数fに依存する特性を表している。
 図16は、局所値導出ステップS2中の伝達関数導出手順J13で導出された周波数領域の伝達関数G(f)を例示した図である。図16で、横軸は周波数fを示し、縦軸はゲインg(デシベル値)を示している。ゲインgは、周波数fの増加に伴って概ね滑らかに減少しており、比較的高い或る周波数f1から周波数f2までの周波数帯で減少傾向が鈍化している。
 次のモデル合わせ込み手順J14では、移動体4の質量Mおよび移動時の摩擦条件を可変にモデル化した物理モデルを周波数領域の伝達関数G(f)に合わせ込むことにより、移動体のモデル化された質量Mを推定する。物理モデルとして、次式で表されるラプラス変換法を用いたモデルを使用する。
   モデル化された伝達関数Gm(s)=(X/F)=1/(Ms+Ds+K)
ただし、ラプラス演算子s、位置座標X、推力F、モデル化された質量M、摩擦抵抗Rのうち速度Vに比例する成分を表す定数D、および摩擦抵抗Rのうち速度Vに関係しない一定成分を表す定数Kである。定数Dおよび定数Kにより、摩擦抵抗Rの非線形性がモデル化されている。
 ここで、質量M、定数D、および定数Kを可変に調整することで、モデル化された伝達関数Gm(s)を導出された伝達関数G(f)に合わせ込むことができる。図17は、局所値導出ステップS2中のモデル合わせ込み手順J14で、合わせ込まれた伝達関数Gm(s)の例を示した図である。図17で、モデル化された伝達関数Gm(s)は破線で示され、伝達関数G(f)に合わせ込まれている。このとき、ゲインgの大きな周波数f1以下の低周波帯で2つの伝達関数Gm(s)および伝達関数G(f)を良く一致させることが重要である。また、ゲインgの小さな周波数f1~周波数f2の周波数帯では、或る程度の乖離を許容する。これにより、モデル化された質量Mが定まる。
 次の局所値算出手順J15では、移動体4の実測質量M0または設計質量M0をモデル化された質量Mで除算して補正倍率C1を求める。ここで、移動体4の実際の質量は当然ながら不変であるので、モデル化された質量Mが実測質量M0よりも小さめであるということは、設計値よりも大きめの移動量が実測されたことを意味する。つまり、推力Fおよび局所推力定数N1が設計値よりも大きめであることを意味する。逆に、モデル化された質量Mが実測質量M0よりも大きめであるということは、設計値よりも小さめの移動量が実測されたことを意味する。つまり、推力Fおよび局所推力定数N1が設計値よりも小さめであることを意味する。したがって、補正倍率C1の考え方は妥当なものとなる。そして、補正倍率C1を推力定数の設計値Ndnに乗算することで位置座標X1の局所推力定数N1を算出できる。
 次の手順J16では、最後の実測箇所であるか否か判定し、最後でなければ、次の実測箇所に移って加振実測手順J12に戻る。図13の例では、位置座標X1~X6の6箇所で加振実測手順J12から手順J16までを繰り返す。これにより、各位置座標X1~X6についてそれぞれ、補正倍率C1~C6および局所推力定数N1~N6を算出できる。手順J16で、最後の実測箇所であると、局所値検出ステップS2を終了する。
 なお、加振実測手順J12で、高速正弦波掃引法以外の方法を用いることもできる。例えば、ステップ掃引法により交流電流Iacの周波数fを段階的に変化させ、各周波数fにおける移動体4の微小振動の振幅を実測してもよい。この場合、伝達関数導出手順J13で、フーリエ変換に代えて周波数fと微小振動の振幅との関係を求める。その他にも、局所値導出ステップS2は、適宜変形できる。
 第1実施形態で導出された軌道部材2上の任意の位置座標Xの推力定数N(X)は、第2実施形態のリニアモータの移動制御方法に用いることができる。これに限定されず、第1実施形態のリニアモータの推力定数導出方法は、リニアモータ1を製造したときの製品検査方法に応用することができる。具体的に例えば、実測により導出した推力定数N(X)の変化範囲に制限判定値を設けることで、リニアモータ1の性能の合否を判定できる。
 第1実施形態のリニアモータの推力定数導出方法は、移動方向に延在して永久磁石3を有する軌道部材2と、軌道部材2に移動可能に装架されてコイル5を有する移動体4とを備え、コイル5に電流Iを通電して永久磁石3との間に移動方向の推力Fを発生するリニアモータ1において、電流Iに対する推力Fの発生比率を表す推力定数Nを導出する方法であって、軌道部材2上の長い移動区間内の平均的な推力定数である平均推力定数Nav1を実測に基づいて導出し、軌道部材2上の複数箇所(位置座標X1~X6)の局所的な推力定数である局所推力定数N1~N6をそれぞれ実測に基づいて導出し、平均推力定数Nav1および各箇所の局所推力定数N1~N6に基づいて軌道部材2上の任意の位置(位置座標X)の推力定数N(X)を導出する。
 これによれば、複数の実測条件で実測を行うことにより、軌道部材2上の長い移動区間で平均推力定数Nav1を導出でき、軌道部材2上の複数箇所でそれぞれ局所推力定数N1~N6を導出できる。そして、両者を比較照合することにより移動体4と軌道部材2との間の摩擦抵抗Rを定量的に把握して、その影響を低減することができる。したがって、複数の実測で導出された平均推力定数Nav1および各箇所の局所推力定数N1~N6に基づく軌道部材2上の任意の位置(位置座標X)の推力定数N(X)は、単一の実測条件に基づく従来技術と比較して導出精度が高められる。
 さらに、第1実施形態のリニアモータの推力定数導出方法は、軌道部材2上の長い移動区間に渡って移動体4を移動させ、平均推力定数Nav1を導出する平均値導出ステップS1と、軌道部材2上の複数箇所で移動体2を移動方向に微小移動させて、各箇所の局所推力定数N1~N6をそれぞれ導出する局所値導出ステップS2と、各箇所の局所推力定数N1~N6の平均値Nav2と平均推力定数Nav1との差分ΔNを算出し、各箇所の局所推力定数N1~N6を差分ΔNだけ補正して各箇所の補正推力定数N1c~N6cとする局所値補正ステップS3と、各箇所の補正推力定数N1c~N6cを補間することにより、軌道部材2上の任意の位置(位置座標X)の推力定数N(X)を導出する補間導出ステップS4と、を有する。
 これによれば、摩擦抵抗Rの影響を除去した平均推力定数Nav1と摩擦抵抗Rの影響を受ける各箇所の局所推力定数N1~N6の平均値Nav2との差分ΔNを算出するので、摩擦抵抗Rの影響を確実に低減できる。これにより、推力定数N(X)の導出精度は、従来よりも格段に高められる。
 さらに、第1実施形態のリニアモータの推力定数導出方法で、局所値導出ステップS2は、コイル5に周波数可変の交流電流Iacを通電して移動体4を移動方向前後に加振しながら、交流電流Iacおよび移動体4の位置情報(位置座標X)を実測する加振実測手順J12と、交流電流Iacの時間的変化および移動体4の位置情報(位置座標X)の時間的変化に基づいて、周波数領域の伝達関数G(f)を導出する伝達関数導出手順J13と、移動体4の質量Mおよび移動時の摩擦条件(定数Dおよび定数K)を可変にモデル化した物理モデルを周波数領域の伝達関数G(f)に合わせ込むことにより、移動体4のモデル化された質量Mを推定するモデル合わせ込み手順J14と、移動体4の実測質量M0または設計質量M0を移動体4のモデル化された質量Mで除算して補正倍率C1~C6を求め、補正倍率C1~C6を推力定数の設計値Ndnに乗算して局所推力定数N1~N6を算出する局所値算出手順J15と、を軌道部材2上の複数箇所(位置座標X1~X6)でそれぞれ実施する。
 これによれば、永久磁石3の強さおよび大きさの個体差や、列設される永久磁石3の配設位置のばらつきなどの変化要因があっても、軌道部材2上の各箇所(位置座標X1~X6)の局所推力定数N1~N6を精度よく導出できる。
 さらに、第1実施形態のリニアモータの推力定数導出方法は、加振実測手順J12で、コイル5に通電する交流電流Iacを正弦波電流とし、かつ周波数fを時間的に徐変させる高速正弦波掃引法を用い、伝達関数導出手順J13で、移動体4の位置情報(位置座標X)の時間的変化をフーリエ変換して周波数領域の伝達関数G(f)を導出する。
 これによれば、高速正弦波掃引法を用いることにより、短時間で効率よく局所推力定数N1~N6を導出できる。
 また、第1実施形態のリニアモータの推力定数導出方法で、平均値導出ステップS1は、移動区間に渡って移動体4を一定の速度で移動させて摩擦抵抗Rの補償に要する摩擦補償電流If1~Ifnを実測することを複数の速度V1~Vnで行い、速度V1~Vnと摩擦補償電流If1~Ifnとを対応付けた摩擦補償テーブルを記憶する定速実測手順J1と、移動区間に渡って移動体4を一定の加速度Aおよび減速度-Aで移動させ、そのときに必要であった加速所要電流Iaおよび移動体4の速度Vを実測する一定加速実測手順J2と、一定加速実測手順J2で実測した移動体4の速度Vに対応する摩擦補償電流Ifを摩擦補償テーブルから取得し、加速所要電流Iaから取得した摩擦補償電流Ifを減算して求めた加速正味電流Imに基づいて平均推力定数Nav1を算出する平均値算出手順J3と、を含む。
 これによれば、摩擦抵抗Rが速度Vに対して非線形に変化してもその影響を除去した平均推力定数Nav1を極めて高精度に導出できる。したがって、推力定数N(X)の導出精度がさらに一層高められる。
 また、第1実施形態のリニアモータの推力定数導出方法は、移動体4の移動方向の所定長LMの範囲に複数のコイル5が列設されており、局所推力定数N1~N6を導出する軌道部材2上の複数箇所(位置座標X1~X6)を所定長LMよりも狭い間隔LDで設定する。
 これによれば、全ての永久磁石3の変化要因の影響が局所推力定数N1~N6の少なくとも1量に反映されて盲点が生じないので、推力定数N(X)の導出精度がさらに一層高められる。
 次に、第2実施形態のリニアモータの移動制御方法について、図18を参考にして説明する。第2実施形態では、第1実施形態と同じリニアモータ1の移動体4を制御対象として移動制御を行う。図18は、第2実施形態のリニアモータの移動制御方法を説明するフロー図である。このフロー図に基づく演算処理や制御は、主に移動体4に搭載された移動制御部が担当する。
 推力定数導出ステップS11で、第1実施形態のリニアモータの推力定数導出方法を実施して、軌道部材2上の任意の位置座標Xの推力定数N(X)を導出する。次の推力定数記憶ステップS12で、推力定数N(X)と位置座標Xとを対応付けた推力定数テーブルをメモリなどに記憶する。推力定数導出ステップS11および推力定数記憶ステップS12は、リニアモータ1の製造直後や修理時、定期点検時などの限られた場合にだけ実施する。これに対して、目標位置設定ステップS13以降は、リニアモータ1の使用時に常に繰り返して実施する。
 目標位置設定ステップS13で、移動体4を移動させる目標位置座標Xtgtを設定する。目標位置座標Xtgtは、外部から指令されてもよく、移動制御部が演算処理により決定してもよい。次の現在位置座標検出ステップS14で、移動制御部は、位置検出部から移動体4の現在位置座標Xnowの情報を取得する。次のステップS15で、移動体4が目標位置座標Xtgtに到達しているか否か、すなわち、現在位置座標Xnow=目標位置座標Xtgtが成り立つか否かを判定する。目標位置座標Xtgtを設定した直後にはステップS15の条件は通常は成り立たず、推力定数取得ステップS16に進む。
 推力定数取得ステップS16で、推力定数テーブルの現在位置座標Xnowに対応する推力定数N(Xnow)を取得する。次の電流演算ステップS17で、コイル5に流す電流Inowの大きさを演算する。この演算は、目標位置座標Xtgt、現在位置座標Xnow、および推力定数N(Xnow)などを演算パラメータとする公知技術を適宜応用して行うことができる。次の電流制御ステップS18で、演算した電流Inowを実際に流すように制御する。その後、現在位置座標検出ステップS14に戻る。
 移動制御部は、2度目以降の現在位置座標検出ステップS14で、再度現在位置座標Xnowの情報を取得し、ステップS15で、現在位置座標Xnow=目標位置座標Xtgtが成り立つか否かを判定する。ステップS15の条件が成り立たない間、移動制御部は、現在位置座標検出ステップS14から電流制御ステップS18までのループを繰り返して実施する。移動体4が目標位置座標Xtgtに到達すると、ステップS15からステップS19に進み、移動体4を停止させる。次に、移動体4の移動が必要になると、目標位置設定ステップS13に戻り、以降は同様の移動制御を繰り返す。現在位置座標検出ステップS14から電流制御ステップS18までは、本発明の推力制御ステップに相当する。
 第2実施形態のリニアモータの移動制御方法は、第1実施形態に記載のリニアモータ装置の推力定数導出方法によって導出された軌道部材2上の任意の位置座標Xの推力定数N(X)と位置座標Xとを対応付けた推力定数テーブルを記憶する推力定数記憶ステップと、移動体4を移動制御するときに、移動体4の軌道部材2上の現在位置座標Xnowを検出し、推力定数テーブルの現在位置座標Xnowに対応する推力定数N(Xnow)を用いてコイル5に通電する電流Inowを制御する推力制御ステップと、を有する。
 これによれば、従来技術よりも精度の高い推力定数N(X)が推力定数テーブルに記憶されており、移動体4の現在位置座標Xnowに対応する推力定数N(Xnow)を用いて移動制御を行える。したがって、移動体4の移動制御の精度が従来よりも高められる。
 なお、第1実施形態のリニアモータの推力定数導出方法、ならびに、第2実施形態のリニアモータの移動制御方法は、ともに装置として実施することもできる。また、装置として実施した場合の効果は、第1および第2実施形態と同様である。本発明は、その他にも様々な応用や変形が可能である。
  1:リニアモータ  2:軌道部材  21、22:長尺プレート
  3、31、32、33、35:永久磁石  4:移動体  5:コイル
  G、G1~G8:ギャップ長  X、X1~X6:位置座標
  LM:所定長  LD:間隔
  N、N(X):推力定数
  Nav1:平均推力定数  N1~N6:局所推力定数
  Nav2:局所推力定数の平均値  ΔN:差分
  N1c~N6c:補正推力定数
  V:速度
  If:摩擦補償電流  Ia:加速所要電流  Im:加速正味電流
  G(f):周波数領域の伝達関数
  Gm(s):モデル化された伝達関数

Claims (9)

  1.  移動方向に延在して磁石およびコイルの一方を有する軌道部材と、前記軌道部材に移動可能に装架されて前記磁石および前記コイルの他方を有する移動体とを備え、前記コイルに電流を通電して前記磁石との間に移動方向の推力を発生するリニアモータにおいて、前記電流に対する前記推力の発生比率を表す推力定数を導出する方法であって、
     前記軌道部材上の長い移動区間内の平均的な推力定数である平均推力定数を実測に基づいて導出し、前記軌道部材上の複数箇所の局所的な推力定数である局所推力定数をそれぞれ実測に基づいて導出し、前記平均推力定数および各箇所の前記局所推力定数に基づいて前記軌道部材上の任意の位置の推力定数を導出するリニアモータの推力定数導出方法。
  2.  前記軌道部材上の長い移動区間に渡って前記移動体を移動させ、前記平均推力定数を導出する平均値導出ステップと、
     前記軌道部材上の複数箇所で前記移動体を移動方向に微小移動させて、各箇所の前記局所推力定数をそれぞれ導出する局所値導出ステップと、
     各箇所の前記局所推力定数の平均値と前記平均推力定数との差分を算出し、各箇所の前記局所推力定数を前記差分だけ補正して各箇所の補正推力定数とする局所値補正ステップと、
     各箇所の前記補正推力定数を補間することにより、前記軌道部材上の任意の位置の推力定数を導出する補間導出ステップと、
     を有する請求項1にリニアモータの推力定数導出方法。
  3.  前記局所値導出ステップは、
     前記コイルに周波数可変の交流電流を通電して前記移動体を移動方向前後に加振しながら、前記交流電流および前記移動体の位置情報を実測する加振実測手順と、
     前記交流電流の時間的変化および前記移動体の位置情報の時間的変化に基づいて、周波数領域の伝達関数を導出する伝達関数導出手順と、
     移動体の質量および移動時の摩擦条件を可変にモデル化した物理モデルを前記周波数領域の伝達関数に合わせ込むことにより、移動体のモデル化された質量を推定するモデル合わせ込み手順と、
     前記移動体の実測質量または設計質量を前記移動体のモデル化された質量で除算して補正倍率を求め、前記補正倍率を推力定数の設計値に乗算して前記局所推力定数を算出する局所値算出手順と、
     を前記軌道部材上の複数箇所でそれぞれ実施する請求項2に記載のリニアモータの推力定数導出方法。
  4.  前記加振実測手順で、前記コイルに通電する交流電流を正弦波電流とし、かつ前記周波数を時間的に徐変させる高速正弦波掃引法を用い、
     前記伝達関数導出手順で、前記移動体の位置情報の時間的変化をフーリエ変換して前記周波数領域の伝達関数を導出する請求項3に記載のリニアモータの推力定数導出方法。
  5.  前記平均値導出ステップは、
     前記移動区間に渡って前記移動体を一定の速度で移動させて摩擦抵抗の補償に要する摩擦補償電流を実測することを複数の速度で行い、前記速度と前記摩擦補償電流とを対応付けた摩擦補償テーブルを記憶する定速実測手順と、
     前記移動区間に渡って前記移動体を一定の加速度で移動させ、そのときに必要であった加速所要電流および移動体の速度を実測する一定加速実測手順と、
     前記一定加速実測手順で実測した移動体の速度に対応する摩擦補償電流を前記摩擦補償テーブルから取得し、前記加速所要電流から取得した摩擦補償電流を減算して求めた加速正味電流に基づいて前記平均推力定数を算出する平均値算出手順と、
     を含む請求項2~4のいずれか一項記載のリニアモータの推力定数導出方法。
  6.  前記移動体の移動方向の所定長の範囲に複数の磁石またはコイルが列設されており、
     前記局所推力定数を導出する前記軌道部材上の複数箇所を前記所定長よりも狭い間隔で設定する請求項1~5のいずれか一項記載のリニアモータの推力定数導出方法。
  7.  請求項1~6のいずれか一項に記載のリニアモータ装置の推力定数導出方法によって導出された前記軌道部材上の任意の位置の推力定数と位置座標とを対応付けた推力定数テーブルを記憶する推力定数記憶ステップと、
     前記移動体を移動制御するときに、前記移動体の前記軌道部材上の現在位置座標を検出し、前記推力定数テーブルの前記現在位置座標に対応する推力定数を用いて前記コイルに通電する電流を制御する推力制御ステップと、を有するリニアモータの移動制御方法。
  8.  移動方向に延在して磁石およびコイルの一方を有する軌道部材と、前記軌道部材に移動可能に装架されて前記磁石および前記コイルの他方を有する移動体とを備え、前記コイルに電流を通電して前記磁石との間に移動方向の推力を発生するリニアモータにおいて、前記電流に対する前記推力の発生比率を表す推力定数を導出する装置であって、
     前記軌道部材上の長い移動区間内の平均的な推力定数である平均推力定数を実測に基づいて導出する手段と、前記軌道部材上の複数箇所の局所的な推力定数である局所推力定数をそれぞれ実測に基づいて導出する手段と、前記平均推力定数および各箇所の前記局所推力定数に基づいて前記軌道部材上の任意の位置の推力定数を導出する手段と、を有するリニアモータの推力定数導出装置。
  9.  請求項1~6のいずれか一項に記載のリニアモータの推力定数導出方法、あるいは請求項8に記載のリニアモータの推力定数導出装置によって導出された前記軌道部材上の任意の位置の推力定数と位置座標とを対応付けた推力定数テーブルを記憶する推力定数記憶手段と、
     前記移動体を移動制御するときに、前記移動体の前記軌道部材上の現在位置座標を検出し、前記推力定数テーブルの前記現在位置座標に対応する推力定数を用いて前記コイルに通電する電流を制御する推力制御手段と、を有するリニアモータの移動制御装置。
PCT/JP2013/084269 2013-12-20 2013-12-20 リニアモータの推力定数導出方法および移動制御方法、ならびに、リニアモータの推力定数導出装置および移動制御装置 Ceased WO2015092917A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
PCT/JP2013/084269 WO2015092917A1 (ja) 2013-12-20 2013-12-20 リニアモータの推力定数導出方法および移動制御方法、ならびに、リニアモータの推力定数導出装置および移動制御装置
CN201380081715.8A CN105830321B (zh) 2013-12-20 2013-12-20 线性马达的推力常数导出方法、移动控制方法、线性马达的推力常数导出装置及移动控制装置
JP2015553301A JP6313335B2 (ja) 2013-12-20 2013-12-20 リニアモータの推力定数導出方法および移動制御方法、ならびに、リニアモータの推力定数導出装置および移動制御装置
US15/102,543 US10097127B2 (en) 2013-12-20 2013-12-20 Thrust constant derivation method and movement control method of linear motor, and thrust constant derivation device and movement control device of linear motor
EP13899693.9A EP3086468B1 (en) 2013-12-20 2013-12-20 Linear motor thrust constant deriving method and movement control method, and linear motor thrust constant deriving device and movement control device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2013/084269 WO2015092917A1 (ja) 2013-12-20 2013-12-20 リニアモータの推力定数導出方法および移動制御方法、ならびに、リニアモータの推力定数導出装置および移動制御装置

Publications (1)

Publication Number Publication Date
WO2015092917A1 true WO2015092917A1 (ja) 2015-06-25

Family

ID=53402309

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2013/084269 Ceased WO2015092917A1 (ja) 2013-12-20 2013-12-20 リニアモータの推力定数導出方法および移動制御方法、ならびに、リニアモータの推力定数導出装置および移動制御装置

Country Status (5)

Country Link
US (1) US10097127B2 (ja)
EP (1) EP3086468B1 (ja)
JP (1) JP6313335B2 (ja)
CN (1) CN105830321B (ja)
WO (1) WO2015092917A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020129854A (ja) * 2019-02-07 2020-08-27 株式会社Fuji キャリア及び搬送システム
WO2026034007A1 (ja) * 2024-08-08 2026-02-12 株式会社日立ハイテク 搬送装置、検体搬送システム、及び被搬送体の搬送方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6020537B2 (ja) * 2014-11-21 2016-11-02 株式会社安川電機 モータ制御装置及びモータ制御方法
JP6212068B2 (ja) * 2015-04-24 2017-10-11 ファナック株式会社 機械の周波数特性をオンラインで取得する機能を有するサーボ制御装置
JP6704705B2 (ja) * 2015-10-22 2020-06-03 キヤノン株式会社 可動磁石型リニアモータ制御システム及びその制御方法
US10243502B1 (en) * 2018-02-22 2019-03-26 Hamilton Sundstrand Corporation System and method for detecting latent faults in a redundant motor application
CN113283152B (zh) * 2021-06-01 2023-05-02 上海艾为电子技术股份有限公司 驱动电信号的调整方法、驱动方法及驱动电路、电子设备

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0965511A (ja) 1995-08-24 1997-03-07 Canon Inc 多相リニアモータ駆動装置
JP2006211873A (ja) * 2005-01-31 2006-08-10 Canon Inc 移動体制御装置及び移動体制御方法
JP2008131762A (ja) * 2006-11-21 2008-06-05 Canon Inc 制御装置及び露光装置
JP2009159741A (ja) 2007-12-27 2009-07-16 Yaskawa Electric Corp リニアモータの推力リップル測定装置およびその測定方法とモータ制御装置
JP2009159751A (ja) 2007-12-27 2009-07-16 Yaskawa Electric Corp リニアモータの推力リップル補償装置とモータ制御装置
JP2010063343A (ja) * 2008-08-06 2010-03-18 Yaskawa Electric Corp モータ制御装置とそのトルクリップル補正方法及びモータ制御システム
JP2010130854A (ja) 2008-11-28 2010-06-10 Yaskawa Electric Corp リニアモータの推力リップル補償装置およびその補償方法
JP2010207038A (ja) * 2009-03-05 2010-09-16 Yokogawa Electric Corp 推力リップル解析方法及び補正方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE69412719T2 (de) 1993-06-24 1999-02-18 Canon K.K., Tokio/Tokyo Steuerung für einen mehrphasigen Motor
EP1262264A4 (en) * 2000-01-12 2004-05-26 Mitsubishi Electric Corp DRAWER, METHOD AND DEVICE FOR CONTROLLING A DRAWER
US6671637B2 (en) * 2002-01-29 2003-12-30 Sumitomo Heavy Industries, Ltd. Thrust ripple measuring apparatus and method in linear motor
DE10236847A1 (de) * 2002-08-08 2004-02-26 Dr. Johannes Heidenhain Gmbh Verfahren zur Bestimmung des Massenträgheitsmomentes eines elektromotorischen Antriebssystems
US7005810B2 (en) * 2003-02-21 2006-02-28 Matsushita Electric Industrial Co., Ltd. Motor driving apparatus
US7672741B2 (en) * 2003-07-24 2010-03-02 Keio University Position/force control device
JP4474151B2 (ja) * 2003-11-28 2010-06-02 キヤノン株式会社 モータおよびそれを用いた露光装置ならびにデバイス製造方法
JP2006026750A (ja) * 2004-07-12 2006-02-02 Koyo Seiko Co Ltd リニアモータの制御方法
JP4739442B2 (ja) * 2009-06-09 2011-08-03 キヤノン株式会社 ステージ装置、及びステージ装置の制御方法、それを用いた露光装置及びデバイスの製造方法
US8344669B1 (en) * 2012-04-02 2013-01-01 Etagen, Inc. Methods and systems for controlling a multiphase electromagnetic machine

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0965511A (ja) 1995-08-24 1997-03-07 Canon Inc 多相リニアモータ駆動装置
JP2006211873A (ja) * 2005-01-31 2006-08-10 Canon Inc 移動体制御装置及び移動体制御方法
JP2008131762A (ja) * 2006-11-21 2008-06-05 Canon Inc 制御装置及び露光装置
JP2009159741A (ja) 2007-12-27 2009-07-16 Yaskawa Electric Corp リニアモータの推力リップル測定装置およびその測定方法とモータ制御装置
JP2009159751A (ja) 2007-12-27 2009-07-16 Yaskawa Electric Corp リニアモータの推力リップル補償装置とモータ制御装置
JP2010063343A (ja) * 2008-08-06 2010-03-18 Yaskawa Electric Corp モータ制御装置とそのトルクリップル補正方法及びモータ制御システム
JP2010130854A (ja) 2008-11-28 2010-06-10 Yaskawa Electric Corp リニアモータの推力リップル補償装置およびその補償方法
JP2010207038A (ja) * 2009-03-05 2010-09-16 Yokogawa Electric Corp 推力リップル解析方法及び補正方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP3086468A4

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020129854A (ja) * 2019-02-07 2020-08-27 株式会社Fuji キャリア及び搬送システム
JP7313836B2 (ja) 2019-02-07 2023-07-25 株式会社Fuji キャリア及び搬送システム
WO2026034007A1 (ja) * 2024-08-08 2026-02-12 株式会社日立ハイテク 搬送装置、検体搬送システム、及び被搬送体の搬送方法

Also Published As

Publication number Publication date
JP6313335B2 (ja) 2018-04-18
JPWO2015092917A1 (ja) 2017-03-16
EP3086468A4 (en) 2017-01-18
EP3086468B1 (en) 2020-04-15
US20160315574A1 (en) 2016-10-27
EP3086468A1 (en) 2016-10-26
CN105830321A (zh) 2016-08-03
US10097127B2 (en) 2018-10-09
CN105830321B (zh) 2019-04-09

Similar Documents

Publication Publication Date Title
JP6313335B2 (ja) リニアモータの推力定数導出方法および移動制御方法、ならびに、リニアモータの推力定数導出装置および移動制御装置
KR101688360B1 (ko) 서보 제어 장치
JP5293977B2 (ja) クレーンの振れ止め制御方法及び振れ止め制御装置
JP2011010533A (ja) モータ制御装置及びモータ制御システム
CN105814506A (zh) 摩擦辨识方法以及摩擦辨识装置
CN105871282A (zh) 一种基于电机转动惯量的控制器pi参数整定方法
Renton et al. Motion control for linear motor feed drives in advanced machine tools
JP2010079845A (ja) 数値制御装置
US10216151B2 (en) Power consumption-amount estimation apparatus
KR102103248B1 (ko) 구동 장치, 위치결정 장치, 리소그래피 장치, 및 물품 제조 방법
JP5441944B2 (ja) モータ制御装置
US9024548B2 (en) Method of controlling linear motor
CN1243627C (zh) 永磁直线同步电机进给装置
CN103684194A (zh) 异步感应直线电机及异步感应直线电机的驱动方法
CN108075703A (zh) 用于降低永磁激励电马达的同步波动的方法
JP2010130854A (ja) リニアモータの推力リップル補償装置およびその補償方法
JP2010038896A (ja) 速度検出装置および位置検出装置および位置決め装置
JP6661676B2 (ja) ロボット制御装置
CN113541561B (zh) 马达的控制方法和马达的控制装置
CN116647155A (zh) 用自动化设备确定运输转子与定子区段之间的气隙的方法
JP2010207038A (ja) 推力リップル解析方法及び補正方法
JP2017005929A (ja) 駆動システム、推力リップルの補償値の決定方法、モータの制御方法、リソグラフィ装置および物品製造方法
JP2021175365A (ja) モータの制御方法及び制御装置
CN107533334B (zh) 定位控制装置
JP5836161B2 (ja) 制御装置および位置決めテーブル装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13899693

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2015553301

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 15102543

Country of ref document: US

REEP Request for entry into the european phase

Ref document number: 2013899693

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2013899693

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE