CN112327954B - High-precision positioning method for linear motor controlled by asymmetric S-shaped speed curve - Google Patents

High-precision positioning method for linear motor controlled by asymmetric S-shaped speed curve Download PDF

Info

Publication number
CN112327954B
CN112327954B CN202011273554.0A CN202011273554A CN112327954B CN 112327954 B CN112327954 B CN 112327954B CN 202011273554 A CN202011273554 A CN 202011273554A CN 112327954 B CN112327954 B CN 112327954B
Authority
CN
China
Prior art keywords
time
acceleration
point
speed curve
speed
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.)
Active
Application number
CN202011273554.0A
Other languages
Chinese (zh)
Other versions
CN112327954A (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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN202011273554.0A priority Critical patent/CN112327954B/en
Publication of CN112327954A publication Critical patent/CN112327954A/en
Application granted granted Critical
Publication of CN112327954B publication Critical patent/CN112327954B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D13/00Control of linear speed; Control of angular speed; Control of acceleration or deceleration, e.g. of a prime mover
    • G05D13/62Control of linear speed; Control of angular speed; Control of acceleration or deceleration, e.g. of a prime mover characterised by the use of electric means, e.g. use of a tachometric dynamo, use of a transducer converting an electric value into a displacement

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

The invention discloses a high-precision positioning method for a linear motor controlled by an asymmetric S-shaped speed curve, belongs to the field of high-precision high-speed motion control, and relates to the high-precision positioning method for the linear motor controlled by the asymmetric S-shaped speed curve. According to the method, according to the performance constraint of a linear motor, a scale factor is defined and the shape characteristic of a speed curve of an acceleration section and a deceleration section is set, so that a proper motion control parameter is determined; calculating the matching relation of the S-shaped speed curve time period, and updating the S-shaped speed curve time period based on the point-to-point running distance constraint; according to the relation between the acceleration J and the time t of the S-shaped speed curve, solving by adopting a successive integral method to obtain expressions of the acceleration and the time a-t, the speed and the time v-t and the displacement and the time S-t; and dispersing the expression according to the PLC scanning period to obtain a motion instruction, transmitting the motion instruction to the driver in real time, and driving the motor to complete the action according to the set instruction, so that the stability of the operation process is ensured.

Description

High-precision positioning method for linear motor controlled by asymmetric S-shaped speed curve
Technical Field
The invention belongs to the field of high-precision high-speed motion control, and relates to a high-precision positioning method for a linear motor controlled by an asymmetric S-shaped speed curve.
Background
With the rapid development of the electronic packaging industry, high-speed packaging equipment is more and more, and advanced packaging equipment represented by a high-speed chip mounter not only puts higher requirements on precision, but also has high speed as an important index. The high-speed operation of the high-speed chip mounter from a chip taking point to a chip mounting point requires that the linear motor is frequently started and stopped at a large acceleration in a very short time, so that the vibration of equipment is caused, the chip mounting precision is reduced, and the positioning precision of point-to-point motion under the influence of high-speed operation is mainly reflected. Therefore, in the operation process, a reasonable motion curve planning algorithm and a high-precision motion control technology play a crucial role in coordinating the contradiction between the high speed and the high precision of the paster. However, most of the current drivers are packaged with motion control modules, which have low openness and are difficult to realize complex motion control forms. In order to effectively improve the impact and vibration generated by frequent starting and stopping of the linear motor in the point-to-point high-speed running process, the research of a reasonable motion control method becomes the key for realizing a high-speed high-precision motion control system.
The patent of Zhang Cheng et al, "dynamic tracking method, system and device of moving object based on S-curve acceleration and deceleration", is published under the number CN 107671859A. The patent provides a moving target dynamic tracking method based on S curve acceleration and deceleration, which is mainly applied to mechanical arm shutdown control and has lower requirements on the positioning precision of the tail end; wang King Kong et al patent "Flexible acceleration and deceleration control method and System of machine", patent publication No. CN 109656200A. The patent proposes that the S-shaped acceleration and deceleration curve is adopted to avoid the impact generated during acceleration and deceleration, the acceleration and deceleration processes are symmetrical in theoretical research, and the derivation is not carried out on the asymmetrical acceleration and deceleration process, so that the high-precision positioning of the tail end cannot be met.
Disclosure of Invention
Aiming at the defects of the prior art, the invention discloses a high-precision positioning method for a linear motor controlled by an asymmetric S-shaped speed curve. The method can freely define an S-shaped speed curve, and in the process of setting the motion parameters, the motor performance is taken as constraint, the curve shape characteristics of the acceleration and deceleration section are set by a scale factor, and the appropriate motion control parameters are determined. And dividing the S-shaped speed curve into seven time period controls, calculating a time period matching relation by combining the scale factor, and updating the time period according to the relation between the point-to-point operation displacement and the critical displacement. The continuous expressions of a-t, v-t and x-t are solved by an integral method, the continuous expressions are discretized according to a PLC scanning period, the obtained motion instructions are transmitted to a driver in real time, and the linear motor is driven to complete corresponding actions according to the instructions, so that the running stability is ensured, and the positioning accuracy of the linear motor is improved. The method realizes complex motion control, realizes the stability of the operation process by reasonably planning motion parameters, improves the point-to-point motion precision, and provides theoretical and technical support for high-speed and high-precision motion control.
The technical scheme of the invention is a high-precision positioning method of a linear motor controlled by an asymmetric S-shaped speed curve, which defines a scale factor and sets the shape characteristics of the speed curve of an acceleration section and a deceleration section according to the performance constraint of the linear motor, thereby determining proper motion control parameters; calculating the matching relation of the S-shaped speed curve time period, and updating the S-shaped speed curve time period based on the point-to-point running distance constraint; according to the relation between the acceleration J and the time t of the S-shaped speed curve, solving by adopting a successive integral method to obtain expressions of the acceleration and the time a-t, the speed and the time v-t and the displacement and the time S-t; and dispersing the expression according to the PLC scanning period to obtain a motion instruction, transmitting the motion instruction to the driver in real time, and driving the motor to complete the action according to the set instruction, so that the stability of the operation process is ensured. The method comprises the following specific steps:
step 1: determining high-speed motion control parameters under the constraint of linear motor performance;
for point-to-point motion of a high-speed chip mounter, the initial speed is 0, and in order to reduce impact caused by frequent starting and stopping of a linear motor, a traditional T-shaped speed curve is replaced by a common S-shaped speed curve. During the point-to-point acceleration and deceleration operation process, the acceleration section I, the constant speed section II and the deceleration section III can be divided. Recording the running time of the acceleration section I as taThe running time of the uniform speed section II is tconWith deceleration III for a running time td. In order to facilitate the analysis of the functional expression of the S-shaped speed curve, the acceleration section I is operated for a time taSubdivided into acceleration periods t of rising acceleration1Acceleration period t during which acceleration is constant2And acceleration period t of acceleration drop3Deceleration time tdSubdivided into deceleration periods t of rising acceleration5Deceleration time period t with constant acceleration6And a deceleration period t during which the acceleration decreases7Then there is ta=t1+t2+t3And td=t5+t6+t7. In order to represent the change from the S-shaped speed curve to the T-shaped speed curve, two scaling factors 0 & ltalpha & gt & lt1 & gt and 0 & ltbeta & gt & lt1 & gt are set, the two scaling factors respectively represent the percentage of S-shaped functions of the acceleration section and the deceleration section, and the smaller the value of the scaling factor is, the closer the scaling factor is to the ideal S-shaped curve.
Thus, the seven-segment local time period t of the S-shaped speed curve can be set by the scale factoriI is the relationship between 1, …,7 and the acceleration and deceleration segment time:
Figure GDA0003318912000000031
when α ═ β ═ 0, a pure S-type curve is shown, and when α ═ β ═ 1, a T-type speed curve is shown.
In practical application, parameters of the linear motor in the operation process from the point A to the point B are set according to actual working conditions. Setting the maximum speed V and the acceleration time taAnd deceleration time tdThe maximum jerk J required by the motor in actual operation can be calculated by equations (2) to (3)a,maxAnd maximum deceleration Jd,maxMaximum acceleration aa,maxMaximum deceleration ad,maxComprises the following steps:
Figure GDA0003318912000000032
Figure GDA0003318912000000041
for a given servo linear motor, the maximum acceleration a which can be achieved by the motor can be known according to motor parametersmaxWith maximum jerk Jmax. In order to meet the motor performance, the parameters (maximum speed V, acceleration period time t) set in actual operationaAnd a deceleration period time td) The jerk and acceleration calculated by the equations (2) to (3) satisfy the following conditions:
Figure GDA0003318912000000042
if the conditions are not met, the set initial parameters are invalid, and at the moment, the initial parameters need to be reset to enable the initial parameters not to exceed the maximum performance requirement of the motor.
Step 2: s-shaped speed curve time distribution and updating based on distance constraint
From step 1, the time distribution relationship between the acceleration stage I and the deceleration stage III of the S-shaped velocity curve can be obtained, so as to express the S-shaped velocity curveLine, to determine the time interval t of the uniform velocity section4. However, due to the limitation of the distance between two points, the maximum operating speed cannot be guaranteed to reach the set value V in the actual point-to-point acceleration and deceleration operation process. Based on this, it is necessary to determine the time allocation relationship according to the relationship between the set parameter and the given distance. The calculation method is as follows:
first assume t4=tconThe critical displacement s required to achieve the set maximum speed V during operation from point a to point B at this time can be determined according to equation (5) at 00
Figure GDA0003318912000000043
When the running distance s from the point A to the point B satisfies s & gt s0When the temperature of the water is higher than the set temperature,
Figure GDA0003318912000000051
the running distance s > s from the point A to the point B can be obtained at the moment0Time distribution relation of S-shaped speed curve in process [ T ]]=[t1t2t3t4t5t6t7]。
When the running distance s from the point A to the point B meets the condition that s is less than s0When the speed curve is in S-shaped speed curve displacement and time relation, S is equal to S (t)3) Therefore, the actual operation displacement s and the critical displacement s can be determined0With respect to time tiThe proportionality coefficient ε:
Figure GDA0003318912000000052
the seven-segment time t of the S-shaped speed curve can be calculated by the above formulaiThe updating is as follows:
Figure GDA0003318912000000053
thus, the running distance s from the point A to the point B is less than s0Time distribution relation of S-shaped speed curve in process [ T ]]=[t′1t′2t′3t′4t′5t′6t′7]。
And step 3: establishing an S-shaped speed curve expression;
and analyzing in the step 2 to obtain seven-segment time distribution relation of the S-shaped speed curve under different displacement constraints. Then, an expression between the displacement and the time of the S-shaped velocity curve needs to be determined, and considering that the S-shaped velocity curve has a segmented first-order acceleration, the expression between the acceleration and the time can be obtained as follows:
Figure GDA0003318912000000061
wherein, TiWhere i is 1, …,7 denotes the time coordinate, T1=t1,T2=T1+t2,T3=T2+t3,T4=T3+t4,T5=T4+t5,T6=T5+t6,T7=T6+t7
The relationship between the acceleration and the time is integrated to obtain the relationship between the speed and the time, and the relationship between the displacement and the time is obtained by integrating again, so that the values of the speed v and the displacement S corresponding to the S-shaped function at different moments t can be determined:
Figure GDA0003318912000000062
Figure GDA0003318912000000071
equations (9) - (11) give S-shaped speed curves a-t, v-t and S-t expressions, and accordingly real-time point position control of the linear motor in the operation process from the point A to the point B can be achieved.
And 4, step 4: the method is realized by the rapid point-to-point motion control of a linear motor constrained by a PLC scanning period;
from step 3, the continuous expressions for s-t, v-t and a-t can be obtained, for any given time tiThe position s corresponding to the linear motor at the moment can be calculated according to the continuous equationiV speed of movementiAnd acceleration aiInformation, for convenience of presentation, is denoted as [ t ]i,si,vi,ai]. In practical application, the PLC has the minimum scanning period tsLimitation, therefore, for the moving process from point A to point B, the continuous operation time T needs to be according to the scanning period TsThe integer multiples of which are discretized. When the PLC program is used for calculation, because the time parameter in the PLC can be expressed as the number of scanning periods and the data type is an integer, when the time variable is repeatedly iterated according to t ═ t +1, the scanning period t is expressedsAnd (3) superposition. In order to ensure the consistency of the calculation results, converting the time periods in the formulas (6) and (8) into the number of scanning cycles, and rounding up in a rounding mode;
Figure GDA0003318912000000072
in the formula, n _ tiRepresents a time period tiThe corresponding number of scan cycles, Round, represents the rounding function according to the rounding rule.
Then according to the scanning period tsAnd (3) updating the time distribution relation obtained in the formulas (6) and (8):
Figure GDA0003318912000000081
thereby obtaining an updated time period
Figure GDA0003318912000000082
Successive a-t equations, v-t equations and s-t equations can be calculated according to equations (9) - (11).
On the basis, the operation period T is dispersed into T ═ T according to the minimum scanning periods 2ts 3ts … T]TIn this case, the corresponding information [ t, s, v, a ] of the linear motor in each scanning cycle can be obtained from equations (9) to (11)]. The obtained displacement velocity acceleration [ s, v, a ] in each scanning period]The motion instruction is sent to a motor driver in real time through an EtherCAT network through a PLC control program, and the driver drives the linear motor to complete the rapid point-to-point motion control of the linear motor based on the S-shaped speed curve.
The invention has the obvious effects and benefits that an asymmetric S-shaped speed curve capable of customizing time parameters is provided and applied to point-to-point motion control of the linear motor, so that the problem of overlarge motor motion impact caused by unreasonable speed planning can be effectively reduced, and the method has important significance for improving the positioning precision of the motor. In the process of setting the motion parameters, the performance of a motor is taken as constraint, the curve shape characteristics of the acceleration curve of the acceleration and deceleration section are set by a scale factor, and appropriate motion control parameters are determined; dividing the S-shaped speed curve into seven time period controls, calculating a time period matching relation by combining the scale factor, and updating the time period according to the relation between the point-to-point operation displacement and the critical displacement; solving continuous expressions of a-t, v-t and x-t by using an integral method; and discretizing the continuous expression according to the PLC scanning period, transmitting the obtained motion instruction to the driver in real time, and driving the linear motor to complete corresponding actions according to the instruction, thereby ensuring the operation stability and improving the positioning precision of the linear motor. The method is suitable for various high-speed high-precision motion control occasions requiring complex motion control, particularly high-speed chip mounters, can improve the impact of a linear motor in the operation process and improve the positioning precision, and thus has important significance for improving the packaging precision and performance of electronic products.
Drawings
FIG. 1 is a general flow chart of a high-precision positioning method for a linear motor controlled by an asymmetric S-shaped speed curve.
FIG. 2 shows that the set travel distance s is 40mm (s < s)0)、s=50mm(s=s0) And s 60mm (s > s)0) Time-programmed displacement-time curve, s050mm is the critical displacement required to reach the maximum velocity V. Wherein, the abscissa represents the operation time in ms, and the ordinate represents the displacement in mm.
FIG. 3 shows that the travel distance s is 40mm (s < s)0)、s=50mm(s=s0) And s 60mm (s > s)0) Time programmed speed time curves. Wherein, the abscissa represents the operation time in ms, and the ordinate represents the speed in mm/s.
FIG. 4 shows that the travel distance s is 40mm (s < s)0)、s=50mm(s=s0) And s 60mm (s > s)0) Time programmed acceleration time curves. Wherein the abscissa represents the operation time in ms, and the ordinate represents the acceleration in mm/s2
Fig. 5 shows a set displacement time curve and an actual displacement time curve when the patch distance s is 50 mm. Wherein the abscissa represents the operation time in units of s, and the ordinate represents the displacement in units of mm.
Fig. 6a) shows a displacement time curve of the driver internal control algorithm at a patch distance s of 50mm, wherein the overall operation time T is 40ms, and the settling time T is 17 ms;
fig. 6b) shows a displacement time curve of the proposed control algorithm for a patch distance s of 50mm, wherein the overall operating time T is 31ms and the settling time T is 17 ms.
Detailed Description
The detailed description of the invention will be made in conjunction with the technical solutions and the accompanying drawings
In the high-speed high-precision packaging process of the semiconductor, reasonable movement speed planning can avoid the situation that equipment generates impact or large vibration due to high-speed movement of a motor rotor so as to influence the chip mounting precision of a high-speed chip mounter, and has important significance for improving the stability of the chip mounting process and improving the chip mounting speed and precision. Based on the method, the high-precision positioning method of the linear motor controlled by the asymmetric S-shaped speed curve is invented, and the flow of the method is shown in the attached drawing 1.
The invention uses TwinCAT2 software in a Beifu industrial personal computer as a soft PLC, according to the method provided by the text, the PLC program transmits position, speed, acceleration and time instructions to an Elmo driver through an EtherCAT network cable, and the driver controls a linear motor to output corresponding instructions, thereby completing the high-speed and high-precision positioning control of the linear motor. The specific implementation process of the invention is explained in detail by taking the example of verifying the effectiveness of the method by MATLAB simulation and realizing point-to-point motion of the linear motor controlled by the S-shaped speed curve in TwinCAT2 by utilizing PLC programming.
Firstly, determining the maximum acceleration a of the linear motor according to the performance parameters of the linear motormaxAnd maximum jerk JmaxSetting the maximum speed V and the acceleration time t in the running process according to the actual running environmentaAnd a deceleration time tdAnd the scale factors alpha and beta, the acceleration time and the deceleration time of which can be further subdivided by formula (1) to obtain a time period distribution relation about the S-shaped speed curve. According to the set operation parameters, the maximum acceleration { J ] required by the linear motor in the acceleration section and the deceleration section under the parameters can be calculated according to the formulas (2) to (3)a,max,Jd,maxAnd maximum acceleration aa,max,ad,max}. If equation (4) is satisfied, the set operating parameters can be implemented, otherwise, the parameters need to be readjusted to satisfy the constraint of equation (4).
Secondly, it is considered that a certain displacement distance s is required in order to reach the set maximum operating speed V. However, in an actual working condition, the actual running distance S is limited by the point-to-point distance, and the maximum running speed cannot be guaranteed to reach the set value, so that the time distribution relationship of each section of the S-shaped speed curve needs to be determined according to the relationship between the set parameter and the actual running distance. Firstly, the critical displacement s required for reaching the maximum speed V can be calculated according to the formula (5)0If the actual running distance satisfies s > s0The time distribution [ T ] of the S-shaped speed curve can be determined according to the formula (6)]=[t1t2t3t4t5t6t7](ii) a If the actual running distance satisfies s < s0In which the S-shaped speed profile can be updated according to equations (7) to (8)Time distribution relation [ T ]]=[t′1t′2t′3t′4t′5t′6t′7]。
And then, according to the time distribution relation obtained in the step 2, by using the expressions of the jerk and the time, and calculating and deducing the expressions of the acceleration and the time, the speed and the time and the displacement and the time through integral calculation, wherein the expressions are shown in the formulas (9) to (11). According to the a-t equation, the v-t equation and the S-t equation, the operation process of the linear motor from the point A to the point B can be controlled in real time by the S-shaped speed curve.
And finally, rounding the S-shaped speed curve time distribution relation based on the scanning period in a rounding mode by using equations (12) - (13), calculating to obtain corresponding acceleration, speed and displacement in each scanning period by using equations (9) - (11), and inputting the acceleration, speed and displacement to an Elmo driver through an EtherCAT network cable, so that the linear motor is driven to reach an instruction position according to the instruction speed, and the point-to-point high-precision positioning control of the linear motor controlled by the asymmetric S-shaped speed curve is realized.
Setting the maximum speed V in the operation process to 2500mm/s and the acceleration time t according to the actual operation conditiona20ms, deceleration time tdFor example, when the scaling factor α is set to 0.2 and β is set to 0.1, t can be calculated1=t3=8ms,t2=4ms,t5=t7=9ms,t6At this time, the critical distance s required to reach the set maximum speed V of 2500mm/s can be calculated according to the formula (5) at 2ms050 mm. In the actual chip mounter chip mounting process, the displacements from the chip taking points to the chip mounting points are different from 40-80mm, the chip displacements are 40mm, 50mm and 60mm are taken as examples, the time distribution relation calculated above is recalculated and updated according to the formulas (6) - (8) in the step 2, and the acceleration-time curve, the speed-time curve and the position-time curve in the motor operation process are obtained through simulation in the MATLAB by the formulas (9) - (11) and are shown in the figures 2, 3 and 4.
Then, a PLC control program is established in TwinCAT2 by the method, the scanning period of the system is set to be 0.5ms, S-shaped speed curve time distribution is updated by a rounding method according to formulas (12) - (13), every 0.5ms, the PLC program calculates a group of motor motion information [ t, S, v, a ] according to formulas (9) - (11) and transmits the motor motion information [ t, S, v, a ] to a driver, the driver sends out an instruction to control the linear motor to realize accurate positioning of a chip mounting process according to corresponding motion, an actual operation curve fed back by a position encoder is shown in figure 5, and the actual operation can well track a set value.
Meanwhile, in order to verify the effectiveness of the method, taking the patch distance of 50mm as an example, the algorithm proposed herein is compared with the driver internal control algorithm, and the experimental result is shown in fig. 6, fig. 6a) shows that the driver internal control algorithm shifts a time curve when the patch distance s is 50mm, wherein the set operation time T is 40ms, and the stable adjustment time T is 17 ms. Fig. 6b) shows a displacement time curve of the proposed control algorithm for a patch distance s of 50mm, wherein the set operating time T is 31ms and the settling time T is 17 ms. The experimental result shows that under the same positioning precision, the method can improve the operation set period on the premise of ensuring the positioning period, thereby shortening the operation time of the whole patch compared with the internal algorithm of the driver.
Experimental results show that the high-precision positioning method for the linear motor controlled by the asymmetric S-shaped speed curve can meet the high-precision positioning requirement of the linear motor in high-speed operation, shortens positioning time on the premise of ensuring positioning precision, and provides a guiding function for actually improving high-speed chip mounting efficiency in IC packaging.

Claims (1)

1. A high-precision positioning method for a linear motor controlled by an asymmetric S-shaped speed curve is characterized in that according to the performance constraint of the linear motor, a scale factor is defined and the shape characteristic of the speed curve of an acceleration section and a deceleration section is set, so that a proper motion control parameter is determined; calculating the matching relation of the S-shaped speed curve time period, and updating the S-shaped speed curve time period based on the point-to-point running distance constraint; according to the relation between the acceleration J and the time t of the S-shaped speed curve, solving by adopting a successive integral method to obtain expressions of the acceleration and the time a-t, the speed and the time v-t and the displacement and the time S-t; the motion instruction obtained by dispersing the expression is transmitted to the driver in real time according to the PLC scanning period, and the driving motor completes the action according to the set instruction, so that the stability of the operation process is ensured;
the method comprises the following specific steps:
step 1: high-speed motion control parameter determination under linear motor performance constraint
For point-to-point motion of a high-speed chip mounter, the initial speed is 0, and in order to reduce impact caused by frequent starting and stopping of a linear motor, a traditional T-shaped speed curve is replaced by a common S-shaped speed curve; dividing the acceleration section I, the constant speed section II and the deceleration section III in the point-to-point acceleration and deceleration running process; recording the running time of the acceleration section I as taThe running time of the uniform speed section II is tconWith deceleration III for a running time td(ii) a In order to facilitate the analysis of the functional expression of the S-shaped speed curve, the acceleration section I is operated for a time taSubdivided into acceleration periods t of rising acceleration1Acceleration period t during which acceleration is constant2And acceleration period t of acceleration drop3Deceleration time tdSubdivided into deceleration periods t of rising acceleration5Deceleration time period t with constant acceleration6And a deceleration period t during which the acceleration decreases7Then there is ta=t1+t2+t3And td=t5+t6+t7(ii) a In order to represent the change from the S-shaped speed curve to the T-shaped speed curve, two scaling factors of which the alpha is more than or equal to 0 and less than or equal to 1 and the beta is more than or equal to 0 and less than or equal to 1 are set to respectively represent the percentage of the S-shaped function of the acceleration section and the deceleration section, and the smaller the value of the scaling factor is, the closer the scaling factor is to the ideal S-shaped curve;
setting seven-segment local time segment t of S-shaped speed curve by scale factoriAnd acceleration and deceleration segment time, wherein i 1.., 7:
Figure FDA0003180049990000021
when α ═ β ═ 0, a pure S-type curve is shown, and when α ═ β ═ 1, a T-type speed curve is shown;
setting the linear motor slave according to the actual working conditionParameters during the operation from point a to point B: setting the maximum speed V and the acceleration time taAnd deceleration time tdCalculating the maximum acceleration J required by the motor in actual operation through formulas (2) to (3)a,maxAnd maximum deceleration Jd,maxMaximum acceleration aa,maxMaximum deceleration ad,maxComprises the following steps:
Figure FDA0003180049990000022
Figure FDA0003180049990000023
for a given servo linear motor, the maximum acceleration a which can be achieved by the motor can be known according to motor parametersmaxWith maximum jerk Jmax(ii) a In order to meet the motor performance, the parameters set in actual operation are as follows: maximum speed V, acceleration period time taAnd a deceleration period time tdThe jerk and acceleration calculated by the equations (2) to (3) satisfy the following conditions:
Figure FDA0003180049990000024
if the conditions are not met, the set initial parameters are invalid, and at the moment, the initial parameters need to be reset to ensure that the initial parameters do not exceed the maximum performance requirement of the motor;
step 2: s-shaped speed curve time distribution and updating based on distance constraint;
from step 1, a time distribution relation between an acceleration section I and a deceleration section III of the S-shaped speed curve is obtained, and in order to express the S-shaped speed curve, a time interval t of a constant speed section needs to be determined4(ii) a However, limited by the distance between two points, the maximum speed of operation cannot be guaranteed to reach the set value V in the actual point-to-point acceleration and deceleration operation process; based on this, it is necessary to determine the time point according to the relationship between the set parameter and the given distanceMatching relation;
first assume t4=tconThe critical displacement s required to achieve the set maximum speed V during the operation from point a to point B at this time is determined according to equation (5) at 00
Figure FDA0003180049990000031
When the running distance s from the point A to the point B satisfies s & gt s0When the temperature of the water is higher than the set temperature,
Figure FDA0003180049990000032
obtaining the running distance s > s from the point A to the point B0In the process, the time distribution relation of S-shaped speed curve [ T]=[t1t2t3t4t5t6t7];
When the running distance s from the point A to the point B meets the condition that s is less than s0When the speed curve is in S-shaped speed curve displacement and time relation, S is equal to S (t)3) Therefore, based on the actual operating displacement s and the critical displacement s0With respect to time tiThe proportionality coefficient ε:
Figure FDA0003180049990000033
calculating seven time t of S-shaped speed curve by the above formulaiThe updating is as follows:
Figure FDA0003180049990000041
the running distance s from the point A to the point B is less than s0In the process, the time distribution relation of S-shaped speed curve [ T]=[t′1t′2t′3t′4t′5t′6t′7];
And step 3: establishing an S-shaped speed curve expression;
analyzing and obtaining seven-segment time distribution relations of S-shaped speed curves under different displacement constraints by the step 2; then, an expression between the displacement and the time of the S-shaped velocity curve needs to be determined, and the expression between the acceleration and the time is obtained by considering that the S-shaped velocity curve has a segmented first-order acceleration:
Figure FDA0003180049990000042
wherein, Ti1, 7 denotes a time coordinate, T1=t1,T2=T1+t2,T3=T2+t3,T4=T3+t4,T5=T4+t5,T6=T5+t6,T7=T6+t7
Integrating the relation between the acceleration and the time to obtain the relation between the speed and the time; the relationship between the displacement and the time can be obtained by integrating again, so that the values of the speed v and the displacement S corresponding to the sigmoid function at different time t can be determined:
Figure FDA0003180049990000051
Figure FDA0003180049990000052
formulas (9) - (11) give S-shaped speed curves a-t, v-t and S-t expressions, so that the real-time point location control of the linear motor in the operation process from the point A to the point B is realized;
and 4, step 4: the method is realized by the rapid point-to-point motion control of a linear motor constrained by a PLC scanning period;
the successive expressions of s-t, v-t and a-t are obtained from step 3, for any given timeMoment tiThe position s corresponding to the linear motor at the moment can be calculated according to the continuous equationiV speed of movementiAnd acceleration aiInformation, for convenience of presentation, is denoted as [ t ]i,si,vi,ai](ii) a In practical application, the PLC has the minimum scanning period tsLimitation, therefore, for the moving process from point A to point B, the continuous operation time T needs to be according to the scanning period TsThe integral multiples of the integer are dispersed; when the PLC program is used for calculation, because the time parameter in the PLC is represented by the number of scanning periods and the data type is an integer, when the time variable is repeatedly iterated according to t ═ t +1, the scanning period t is representedsSuperposition of (2); in order to ensure the consistency of the calculation results, converting the time periods in the formulas (6) and (8) into the number of scanning cycles, and rounding up in a rounding mode;
Figure FDA0003180049990000061
in the formula, n _ tiRepresents a time period tiCorresponding scanning period number, Round represents an integer function according to a rounding rule;
then according to the scanning period tsAnd (3) updating the time distribution relation obtained in the formulas (6) and (8):
Figure FDA0003180049990000062
thereby obtaining an updated time period
Figure FDA0003180049990000063
Calculating to obtain a continuous a-t equation, a v-t equation and an s-t equation according to the formulas (9) to (11);
on the basis, the operation period T is dispersed into T ═ T according to the minimum scanning periods 2ts 3ts…T]TAt this time, the linear motor is obtained by equations (9) to (11) at each timeCorresponding information [ t, s, v, a ] for each scanning cycle](ii) a The obtained displacement velocity acceleration [ s, v, a ] in each scanning period]The motion instruction is sent to a motor driver in real time through an EtherCAT network through a PLC control program, and the driver drives the linear motor to complete the rapid point-to-point motion control of the linear motor based on the S-shaped speed curve.
CN202011273554.0A 2020-11-14 2020-11-14 High-precision positioning method for linear motor controlled by asymmetric S-shaped speed curve Active CN112327954B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011273554.0A CN112327954B (en) 2020-11-14 2020-11-14 High-precision positioning method for linear motor controlled by asymmetric S-shaped speed curve

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011273554.0A CN112327954B (en) 2020-11-14 2020-11-14 High-precision positioning method for linear motor controlled by asymmetric S-shaped speed curve

Publications (2)

Publication Number Publication Date
CN112327954A CN112327954A (en) 2021-02-05
CN112327954B true CN112327954B (en) 2021-11-26

Family

ID=74318626

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011273554.0A Active CN112327954B (en) 2020-11-14 2020-11-14 High-precision positioning method for linear motor controlled by asymmetric S-shaped speed curve

Country Status (1)

Country Link
CN (1) CN112327954B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113253772B (en) * 2021-07-02 2021-10-12 武汉市正弦电气技术有限公司 Servo positioning system and method based on error estimation and compensation
CN114326842B (en) * 2021-11-25 2024-03-01 广州极飞科技股份有限公司 Control method and device of driving device and control method of unmanned equipment
CN116483026B (en) * 2023-06-20 2023-09-22 英孚康(浙江)工业技术有限公司 Multi-mode multi-type asymmetric S-shaped flexible speed curve bidirectional rapid planning method
CN116954280A (en) * 2023-09-21 2023-10-27 北京一控软件技术有限公司 Speed control system of metal rolling post-treatment process line

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105964705A (en) * 2016-07-08 2016-09-28 中国重型机械研究院股份公司 Feed-in synchronization mechanism for large cold rolling pipe mills and synchronization method thereof
CN208678318U (en) * 2018-07-31 2019-04-02 大连理工大学 Intelligent badminton service robot based on crank guide rod mechanism

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB863125A (en) * 1956-11-15 1961-03-15 Inductosyn Ltd Numerically controlled curve interpolating machine
US5453672A (en) * 1991-12-31 1995-09-26 Avitan; Isaac Regulation system for decoupled efficiency optimized operation of DC traction motors
CN101009455B (en) * 2006-12-12 2011-05-11 大连交通大学 Long range big-push force permanent-magnet brushless linear DC motor
US7479751B2 (en) * 2007-01-29 2009-01-20 Rockwell Automation Technologies, Inc. Elimination of unintended velocity reversals in s-curve velocity profiles
CN101477354B (en) * 2009-01-19 2010-08-04 东元总合科技(杭州)有限公司 Position S type instruction generation method
WO2011110289A1 (en) * 2010-03-12 2011-09-15 Sew-Eurodrive Gmbh & Co. Kg Method for determining an optimized desired curve profile in a multi-shaft machine for the on-the-fly processing of material to be conveyed
CN102141794B (en) * 2011-02-25 2012-11-21 广西大学 Continuous locus intersegment connecting speed processing method for numerical control system
CN102664585B (en) * 2012-04-13 2015-01-14 深圳众为兴技术股份有限公司 Precise control device for acceleration/deceleration of motor
CN103558808B (en) * 2013-09-28 2015-12-09 大连理工大学 The kinematics control method that complex-curved five-shaft numerical control processing cutter is vowed
CN104898577B (en) * 2015-04-08 2020-02-11 北京配天技术有限公司 S-shaped curve speed planning method and device, numerical control method and numerical control machine tool
CN106444635B (en) * 2015-08-04 2019-09-13 深圳市雷赛智能控制股份有限公司 A kind of asymmetric S curve Acceleration-deceleration Control Method and device of kinetic control system
CN107850887B (en) * 2015-12-28 2020-06-02 深圳配天智能技术研究院有限公司 S-shaped curve planning method and device and numerical control machine tool
CN106168790B (en) * 2016-02-29 2020-08-18 华南理工大学 S-shaped acceleration and deceleration control method for changing target speed and position on line
CN105892402B (en) * 2016-06-17 2018-06-15 浙江理工大学 Mechanical arm point-to-point motion control method
CN106800040B (en) * 2017-02-24 2022-10-21 南京航空航天大学 Automobile electric control composite steering system and multi-objective optimization method thereof
CN108287527B (en) * 2017-07-18 2019-02-01 山东大学 A kind of improvement S curve Acceleration-deceleration Control Method based on trigonometric function
CN107825424B (en) * 2017-10-31 2020-05-26 河北工业大学 Asymmetric S-shaped trajectory planning method for reducing residual vibration of high-speed manipulator
CN108958173B (en) * 2018-07-27 2021-07-09 山东大学 S curve acceleration and deceleration planning method based on trapezoidal solution and under any displacement speed
CN111381565B (en) * 2018-12-27 2023-02-21 沈阳中科数控技术股份有限公司 Acceleration recursive calculation limiting method for cubic polynomial S curve acceleration and deceleration
CN110286653A (en) * 2019-06-14 2019-09-27 杭州爱科科技股份有限公司 Speed calculation method for arbitrary curve movement S feed speed control
CN110515302B (en) * 2019-08-13 2022-04-15 西安理工大学 Reverse step self-adaptive neural network control method for hydraulic position servo system
CN111158318B (en) * 2020-01-16 2022-10-18 江南大学 Flexible acceleration and deceleration planning method for asymmetric quartic curve
CN111245332B (en) * 2020-01-18 2022-05-24 南京岸鸣智能科技有限公司 Motor speed control method and control system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105964705A (en) * 2016-07-08 2016-09-28 中国重型机械研究院股份公司 Feed-in synchronization mechanism for large cold rolling pipe mills and synchronization method thereof
CN208678318U (en) * 2018-07-31 2019-04-02 大连理工大学 Intelligent badminton service robot based on crank guide rod mechanism

Also Published As

Publication number Publication date
CN112327954A (en) 2021-02-05

Similar Documents

Publication Publication Date Title
CN112327954B (en) High-precision positioning method for linear motor controlled by asymmetric S-shaped speed curve
CN103064344B (en) Non uniform rational B spline (NURBS) interpolation based speed smooth control method
CN101493687B (en) Real time forward looking whole-process acceleration and deceleration controlled NURBS curve self-adapting subsection interpolation method
Fang et al. An approach for jerk-continuous trajectory generation of robotic manipulators with kinematical constraints
CN103676787B (en) A kind of center of circle model space circular interpolation method for kinetic control system
CN104281099A (en) NURBS direct interpolation method and device with processing characteristics considered
CN110083160B (en) Robot trajectory planning method based on deep learning
CN111158318B (en) Flexible acceleration and deceleration planning method for asymmetric quartic curve
CN113253677B (en) Robot motion control method combining speed optimization and feedforward compensation
CN102354146B (en) Motion control system and position control method thereof
CN112433507B (en) LSO-LSSVM (least squares support vector machine) -based five-axis numerical control machine tool thermal error comprehensive modeling method
CN111459016A (en) Method for improving tracking control precision of cutting machine trajectory profile
CN203396032U (en) Room temperature control device based on fuzzy self-adaption PID (proportion integration differentiation)
CN114077226A (en) S-shaped curve speed planning method, control terminal and computer readable storage medium
CN103809521A (en) Spline curve interpolation method based on secant method
CN108099908A (en) A kind of vehicle adaptive cruise optimal control computational methods
Xu et al. An interpolation method based on adaptive smooth feedrate scheduling and parameter increment compensation for NURBS curve
CN103345195A (en) Independent control method and system for steel cord of double twisting strander
CN105329788A (en) Uniform velocity amplitude-varying control method for jib crane
CN105824270A (en) Multi-shaft motion controller and control method of closed-loop PWM pulse
CN111722591B (en) High-precision linkage interpolation method for trademark die-cutting machine
CN112549027A (en) Cartesian space and joint space curve smooth transition method and device
CN110879568B (en) Motion control method of multi-axis linkage economic numerical control system
Liu et al. Adaptive positioning control of servomotor based on model identification
CN112462608A (en) Discrete sliding mode track and speed tracking control method for high-speed train

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