CN102141794A - Continuous locus intersegment connecting speed processing method for numerical control system - Google Patents

Continuous locus intersegment connecting speed processing method for numerical control system Download PDF

Info

Publication number
CN102141794A
CN102141794A CN 201110045938 CN201110045938A CN102141794A CN 102141794 A CN102141794 A CN 102141794A CN 201110045938 CN201110045938 CN 201110045938 CN 201110045938 A CN201110045938 A CN 201110045938A CN 102141794 A CN102141794 A CN 102141794A
Authority
CN
China
Prior art keywords
adjacent
speed
path length
track
sections tracks
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN 201110045938
Other languages
Chinese (zh)
Other versions
CN102141794B (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.)
Guangxi University
Original Assignee
Guangxi University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangxi University filed Critical Guangxi University
Priority to CN2011100459381A priority Critical patent/CN102141794B/en
Publication of CN102141794A publication Critical patent/CN102141794A/en
Application granted granted Critical
Publication of CN102141794B publication Critical patent/CN102141794B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Numerical Control (AREA)

Abstract

The invention discloses a continuous locus intersegment connecting speed processing method for a numerical control system, which comprises the following steps of: determining an adjacent two-segment locus connecting model, and acquiring the movement parameter of the adjacent two-segment locus to be processed; calculating the path length of the adjacent two-segment locus to be processed; calculating the speed of a connecting point of the adjacent two-segment locus under the conditions of acceleration, arch height error, instruction feeding speed and corner constraint; calculating the speed of the connecting point of the adjacent two-segment locus satisfying the condition S curve acceleration and deceleration control-based path length constraint of the front segment locus in the adjacent two-segment locus; and calculating the speed of the connecting point of the adjacent two-segment locus satisfying the condition S curve acceleration and deceleration control-based path length constraint of the rear segment locus in the adjacent two-segment locus, wherein the calculated speed is used as a final value of the adjacent two-segment locus connecting speed. The invention provides the continuous locus intersegment connecting speed processing method for numerical control equipment, so that the processing locus of the numerical control equipment is continuous and the speed is smooth.

Description

A kind of intersegmental linking velocity process of continuous path method that is used for digital control system
Technical field
The present invention relates to a kind of intersegmental linking velocity process of continuous path method that is used for digital control system in the method for planning track in the digital control system field, particularly numerically-controlled machine and the robot motion's control technology.
Background technology
In numerically-controlled machine processing and robot motion's control, trajectory planning occupies an important position.In machine tooling, trajectory planning mainly is that the speed of lathe is planned.For the numerically-controlled machine of continuous cutting, the control Significance of speed of feed is great.It not only directly influences the roughness and the precision of processing parts, and closely related with life-span, the production efficiency of cutter and lathe.In robot motion control, the robot operational process requires action continuously smooth and avoiding obstacles smoothly.
In practice, the processing of numerically-controlled machine and the motion of robot generally are made up of the multistage movement locus, if every section track to be processed is all adopted from the static target velocity that accelerates to, decelerate to zero method for planning track again, to reduce average speed of feed, influence working (machining) efficiency; And because the frequent start-stop of motor and acceleration and deceleration running, cause not only that surface of the work is coarse, precision is low, and the expanded motor load, electrical machinery life reduced.
At present, for addressing the above problem main employing continuous path planing method.In the continuous path planing method, acceleration and deceleration control method commonly used mainly is S curve acceleration and deceleration control methods.But the algorithm of the intersegmental linking speed of existing continuous path, not considering will be based on the path of adjacent two sections tracks of S curve acceleration and deceleration control as constraint condition, the actual feed that causes being connected the some place at adjacent two sections tracks does not reach precalculated speed, linking speed between adjacent two sections orbit segments produces sudden change, and speed of feed can not link up smoothly between adjacent two sections orbit segments; When part processing was finished, tool speed can not in time drop to zero, produces and cuts.
Summary of the invention
Purpose of the present invention aims to provide a kind of intersegmental linking velocity process of continuous path method that is used for digital control system, make that being connected an actual feed at place at adjacent two sections tracks can arrive precalculated speed, can link up smoothly between promptly adjacent two sections orbit segments.Digital control system does not need to thresh last spot speed accessibility judgement and adjustment when doing S curve acceleration and deceleration control, thereby improves digital control system acceleration and deceleration control efficient.
For realizing above-mentioned target, the technical solution adopted in the present invention is as follows:
Digital control system in the numerical control device is deciphered processing to the numerical control program in machine code, and store decode results into the decoding buffer zone, utilize the program segment information of deciphering in the buffer zone to the intersegmental linking velocity process of continuous path to be processed, make that the numerical control device machining locus is continuous, rate smoothing, the intersegmental linking velocity process of described continuous path method may further comprise the steps:
(1) determines that adjacent two sections tracks are connected model, and obtain the kinematic parameter of adjacent two sections tracks to be processed;
(2) calculate the path of adjacent two sections tracks to be processed;
(3) speed of the linking point of the adjacent two section tracks of calculating under acceleration, bow high level error, instruction speed of feed and turning constraint condition;
(4) calculate to satisfy speed based on the linking point of adjacent two sections tracks of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control;
(5) calculate to satisfy speed based on the linking point of adjacent two sections tracks of the path length constraint condition of the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control, and as the end value of linking speed between adjacent two sections orbit segments.
The path of described adjacent two sections tracks to be processed comprises straight line path length and circular arc path length; Described circular arc path length calculation step comprises: (2.1) calculate the positive dirction angle in the horizontal axis of circular arc initial point place principal normal unit vector and workpiece coordinate system, (2.2) calculating is at the positive dirction angle of the horizontal axis of circular arc end point place's principal normal unit vector and workpiece coordinate system, and (2.3) calculate circular arc path length;
(2.1) calculate in circular arc initial point place principal normal unit vector
Figure BSA00000439941800021
Positive dirction angle α with the horizontal axis of workpiece coordinate system I_s:
sin α i _ s = y i _ c - y i R i - - - ( 1 )
cos α i _ s = x i _ c - x i R i - - - ( 2 )
Wherein: i=1,2,3...N, N represent orbit segment number to be processed; (x i, y i) be i section Track Initiation point coordinate; (x I_c, y I_c) be the central coordinate of circle of circular arc when i section track is circular arc; R iBe the radius of circular arc when i section track is circular arc, get α by formula (1) and (2) I_sComputing formula be:
Figure BSA00000439941800024
(2.2) calculate principal normal unit vector at some place, circular arc end
Figure BSA00000439941800025
Positive dirction angle α with the horizontal axis of workpiece coordinate system I_e:
sin α i _ e = y i _ c - y i + 1 R i - - - ( 4 )
cos α i _ e = x i _ c - x i + 1 R i - - - ( 5 )
Wherein: (x I+1, y I+1) be i section track end point coordinate;
Get α by formula (4) and (5) I_eComputing formula be:
Figure BSA00000439941800028
(2.3) calculate circular arc path length S i:
Get S by formula (3) and (6) iComputing formula be:
Described path length constraint condition based on the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control is inequality group (8), (9) and (10):
| V i , T - V i _ s | &GreaterEqual; a max 2 J max S i > 2 V i _ s a max J max + a max 3 J max 2 S i < J max ( V i , T 2 - V i _ s 2 ) + a max 2 ( V i , T + V i _ s ) 2 J max a max - - - ( 8 )
| V i , T - V i _ s | < a max 2 J max S i < ( V i , T + V i _ s ) | V i , T - V i _ s | J max - - - ( 9 )
| V i , T - V i _ s | &GreaterEqual; a max 2 J max S i &le; 2 V i _ s a max J max + a max 3 J max 2 - - - ( 10 )
Wherein: V I_sBe the initial point speed of i section track; a MaxPeak acceleration for the digital control system permission; J MaxMaximum acceleration for the digital control system permission; S iIt is the path of i section track; V I, TSpeed for the linking point of adjacent two sections tracks of under acceleration, bow high level error, instruction speed of feed and turning constraint condition, calculating;
When inequality group (8), (9) and (10) when all being false, the speed V of the linking point of the adjacent two sections tracks that calculate is described under acceleration, bow high level error, instruction speed of feed and turning constraint condition I, TCan satisfy path length constraint condition, satisfy speed V based on the linking point of adjacent two sections tracks of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control based on the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control I, T_fEqual V I, T
When inequality group (8), (9) and (10) by and only form immediately by one, the speed V of the linking point of the adjacent two sections tracks that calculate under acceleration, bow high level error, instruction speed of feed and turning constraint condition is described I, TCan not satisfy path length constraint condition, satisfy speed V based on the linking point of adjacent two sections tracks of the path length constraint of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control based on the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control I, T_fComputing method are: when inequality group (8) is set up, and V I, T_fCalculate according to formula (11); When inequality group (9) or (10) establishment, V I, T_fCalculate according to formula (12) and (13);
V i , T _ f = a max 4 - 4 J max ( a max 2 V i _ s - J max V i _ s 2 - 2 J max a max S i ) - a max 2 2 J max - - - ( 11 )
V i , T _ f = V i _ s + J max T jerk 2 - - - ( 12 )
J max T jerk 3 + 2 V i , T _ f - S i = 0 - - - ( 13 )
Wherein: T JerkFor adding the time of accelerating sections;
Described path length constraint condition based on the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control is inequality group (14), (15) and (16):
V i , T _ f &GreaterEqual; a max 2 J max S i + 1 > 2 V i , T _ f a max J max + a max 3 J max 2 S i + 1 < a max 2 V i , T _ f - J max V i , T _ f 2 2 J max a max - - - ( 14 )
V i , T _ f < a max 2 J max S i + 1 < V i , T _ f V i , T _ f J max - - - ( 15 )
V i , T _ f &GreaterEqual; a max 2 J max S i + 1 &le; 2 V i , T _ f a max J max + a max 3 J max 2 - - - ( 16 )
Wherein: S I+1It is the path of i+1 section track; V I, T_fFor satisfying speed based on the linking point of adjacent two sections tracks of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control;
When (14), (15) and (16) when all being false, the speed V that satisfies based on the linking point of adjacent two sections tracks of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control is described I, T_fCan satisfy path length constraint condition, satisfy speed V based on the linking point of adjacent two sections tracks of the path length constraint condition of the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control based on the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control I, T_lEqual V I, T_f
When inequality group (14), (15) and (16) by and only form immediately by one, the speed V that satisfies based on the linking point of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control is described I, T_fCan not satisfy path length constraint condition, satisfy speed V based on the linking point of adjacent two sections tracks of the path length constraint condition of the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control based on the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control I, T_lComputing method are: when inequality group (14) is set up, and V I, T_lCalculate according to formula (17); When inequality group (15) or (16) establishment, V I, T_lCalculate according to formula formula (18);
V i , T _ l = a max 2 - a max 4 - 8 J max 2 S i + 1 2 J max - - - ( 17 )
V i , T _ l = 3 J max S i + 1 2 - - - ( 18 )
The present invention compared with prior art has the following advantages and effect:
(1) will be based on the path of adjacent two sections tracks of S curve acceleration and deceleration control as constraint condition
Comprise based on the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control with based on the path length constraint condition of the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control based on the path length constraint condition of adjacent two sections tracks of S curve acceleration and deceleration control, make that being connected an actual feed at place at adjacent two sections tracks can arrive precalculated speed, linking speed between adjacent two sections orbit segments can not produce sudden change, be to link up smoothly between adjacent two sections orbit segments, make the numerical control device machining locus continuous, rate smoothing, and when part processing is finished, tool speed can in time drop to zero, prevents from that cutter from crossing to cut.
(2) improve digital control system acceleration and deceleration control efficient
This disposal route can guarantee that S curve acceleration and deceleration control input parameter is reasonable, and digital control system does not need to thresh last spot speed accessibility judgement and adjustment when doing S curve acceleration and deceleration control, thereby improves digital control system acceleration and deceleration control efficient.
(3) have wide range of applications
This disposal route also can be applied to nurbs curve subsection interpolation algorithm, improves the pre-service efficient of the prediction pretreatment module in the nurbs curve subsection interpolation algorithm.
Description of drawings
Fig. 1 is the intersegmental linking velocity process of a continuous path of the present invention method overview flow chart
Fig. 2 (a) is that straight path and straight path are connected the model synoptic diagram
Fig. 2 (b) is that straight path and arc track are connected the model synoptic diagram
Fig. 2 (c) is that arc track and straight path are connected the model synoptic diagram
Fig. 2 (d) is that arc track and arc track are connected the model synoptic diagram
Fig. 3 is the path calculation flow chart
Fig. 4 is the speed calculation process flow diagram of the linking point of the adjacent two sections tracks under acceleration, bow high level error, instruction speed of feed and turning constraint condition
Embodiment
The invention will be further described below in conjunction with accompanying drawing.
Shown in the intersegmental linking velocity process of accompanying drawing 1 continuous path of the present invention method overview flow chart, digital control system in the numerical control device is deciphered processing to the numerical control program in machine code, and store decode results into the decoding buffer zone, utilize the program segment information of deciphering in the buffer zone to the intersegmental linking velocity process of continuous path to be processed, the intersegmental linking velocity process of described continuous path may further comprise the steps:
(1) determines that adjacent two sections tracks are connected model, and obtain the kinematic parameter of adjacent two sections tracks to be processed;
Adjacent two sections tracks of the present invention are connected model and comprise following four kinds: straight path shown in the accompanying drawing 2 (a) and straight path are connected model; Straight path shown in the accompanying drawing 2 (b) and arc track are connected model; Arc track shown in the accompanying drawing 2 (c) and straight path are connected model; Arc track shown in the accompanying drawing 2 (d) and arc track are connected model.After determining that pending adjacent two sections tracks are connected model, obtain the kinematic parameter of adjacent two sections tracks to be processed, kinematic parameter of the present invention comprises the path of instruction speed of feed, track.
Repeating step (1), all program segments dispose in the decoding buffer zone.
(2) calculate the path of adjacent two sections tracks to be processed, the calculation procedure flow process as shown in Figure 3;
When i section track is a straight-line segment,
S i = ( x i + 1 - x i ) 2 + ( y i + 1 - y i ) 2 - - - ( 1 )
Wherein: i=1,2,3...N, N represent orbit segment number to be processed; (x i, y i) be i section Track Initiation point coordinate; (x I+1, y I+1) be i section track end point coordinate;
When i section track is an arc section, circular arc path length calculation step comprises: (2.1) calculate the positive dirction angle in the horizontal axis of circular arc initial point place principal normal unit vector and workpiece coordinate system, (2.2) calculating is at the positive dirction angle of the horizontal axis of circular arc end point place's principal normal unit vector and workpiece coordinate system, and (2.3) calculate circular arc path length;
(2.1) calculate in circular arc initial point place principal normal unit vector
Figure BSA00000439941800062
Positive dirction angle α with the horizontal axis of workpiece coordinate system I_s:
sin &alpha; i _ s = y i _ c - y i R i - - - ( 2 )
cos &alpha; i _ s = x i _ c - x i R i - - - ( 3 )
Wherein: i=1,2,3...N, N represent orbit segment number to be processed; (x i, y i) be i section Track Initiation point coordinate; (x I_c, y I_c) be the central coordinate of circle of circular arc when i section track is circular arc; R iBe the radius of circular arc when i section track is circular arc, get α by formula (2) and (3) I_sComputing formula be:
(2.2) calculate principal normal unit vector at some place, circular arc end
Figure BSA00000439941800066
Positive dirction angle α with the horizontal axis of workpiece coordinate system I_e:
sin &alpha; i _ e = y i _ c - y i + 1 R i - - - ( 5 )
cos &alpha; i _ e = x i _ c - x i + 1 R i - - - ( 6 )
Wherein: (x I+1, y I+1) be i section track end point coordinate;
Get α by formula (5) and (6) I_eComputing formula be:
Figure BSA00000439941800072
(2.3) calculate circular arc path length S i:
Get S by formula (4) and (7) iComputing formula be:
Figure BSA00000439941800073
Repeating step (2), all program segments dispose in the decoding buffer zone.
(3) the speed V of the linking point of the adjacent two section tracks of calculating under acceleration, bow high level error, instruction speed of feed and turning constraint condition I, T, calculation process as shown in Figure 4;
(3.1) calculating is in the unit vector of Track Initiation point place tangent line Sine value sin θ with the positive dirction angle of the horizontal axis of workpiece coordinate system I_sWith cosine value cos θ I_s
When i section track is a straight-line segment, shown in accompanying drawing 2 (a) and accompanying drawing 2 (b),
cos &theta; i _ s = y i S i - - - ( 9 )
sin &theta; i _ s = x i S i - - - ( 10 )
When i section track is an arc section, shown in accompanying drawing 2 (c) and accompanying drawing 2 (d),
Figure BSA00000439941800077
Figure BSA00000439941800078
Wherein: θ I_sBe unit vector at Track Initiation point place tangent line
Figure BSA00000439941800079
Positive dirction angle with the horizontal axis of workpiece coordinate system; (x I_c, y I_c) be the central coordinate of circle of circular arc; R iArc radius for circular arc;
(3.2) calculating is in the unit vector of some place, track end tangent line
Figure BSA00000439941800081
Sine value sin θ with the positive dirction angle of the horizontal axis of workpiece coordinate system I_eWith cosine value cos θ I_e
When i section track is a straight-line segment, shown in accompanying drawing 2 (a) and accompanying drawing 2 (b),
cos &theta; i _ e = y i + 1 S i - - - ( 13 )
sin &theta; i _ e = x i + 1 S i - - - ( 14 )
When i section track is an arc section, shown in accompanying drawing 2 (c) and accompanying drawing 2 (d),
sin &theta; i _ e = y i _ c - y i + 1 R i - - - ( 15 )
Figure BSA00000439941800085
Wherein: θ I_eBe unit vector at some place, track end tangent line
Figure BSA00000439941800086
Positive dirction angle with the horizontal axis of workpiece coordinate system
(3.3) the linking contact speed of calculating under acceleration and the adjacent two sections track angles constraint is V I, A
V i , A = min { a max _ x T | cos &theta; i + 1 _ s - cos &theta; i _ e | , a max _ y T | sin &theta; i + 1 _ s - sin &theta; i _ e | } - - - ( 17 )
Wherein: a Max_xBe digital control system x axle acceleration; a Max_yBe digital control system y axle acceleration; T is the sampling period; θ I+1_sBe unit vector at Track Initiation point place tangent line
Figure BSA00000439941800088
Positive dirction angle with the horizontal axis of workpiece coordinate system;
(3.4) the linking contact speed of calculating under the constraint of bow high level error is V I, G
V i , G = min { 8 R i E T , 8 R i + 1 E T } - - - ( 18 )
Wherein: E is the bow high level error that system allows
(3.5) get V I, A, V I, G, i section track instruction speed of feed V I, FWith i+1 section track instruction speed of feed V I+1, FFour minimum value is as the speed V of the linking point of the adjacent two sections tracks under acceleration, bow high level error, instruction speed of feed and turning constraint condition I, T, promptly have:
V i,T=min{V i,A,V i,G,V i,F,V i+1,F}(19)
(4) calculate to satisfy speed based on the linking point of adjacent two sections tracks of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control;
Described path length constraint condition based on the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control is inequality group (20), (21) and (22):
| V i , T - V i _ s | &GreaterEqual; a max 2 J max S i > 2 V i _ s a max J max + a max 3 J max 2 S i < J max ( V i , T 2 - V i _ s 2 ) + a max 2 ( V i , T + V i _ s ) 2 J max a max - - - ( 20 )
| V i , T - V i _ s | < a max 2 J max S i < ( V i , T + V i _ s ) | V i , T - V i _ s | J max - - - ( 21 )
| V i , T - V i _ s | &GreaterEqual; a max 2 J max S i &le; 2 V i _ s a max J max + a max 3 J max 2 - - - ( 22 )
Wherein: V I_sBe the initial point speed of i section track; a MaxPeak acceleration for the digital control system permission; J MaxMaximum acceleration for the digital control system permission; S iIt is the path of i section track; V I, TSpeed for the linking point of adjacent two sections tracks of under acceleration, bow high level error, instruction speed of feed and turning constraint condition, calculating;
When inequality group (20), (21) and (22) when all being false, the speed V of the linking point of the adjacent two sections tracks that calculate is described under acceleration, bow high level error, instruction speed of feed and turning constraint condition I, TCan satisfy path length constraint condition, satisfy speed V based on the linking point of adjacent two sections tracks of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control based on the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control I, T_fEqual V I, T
When inequality group (20), (21) and (22) by and only form immediately by one, the speed V of the linking point of the adjacent two sections tracks that calculate under acceleration, bow high level error, instruction speed of feed and turning constraint condition is described I, TCan not satisfy path length constraint condition, satisfy speed V based on the linking point of adjacent two sections tracks of the path length constraint of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control based on the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control I, T_fComputing method are: when inequality group (20) is set up, and V I, T_fCalculate according to formula (23); When inequality group (21) or (22) establishment, V I, T_fCalculate according to formula (24) and (25);
V i , T _ f = a max 4 - 4 J max ( a max 2 V i _ s - J max V i _ s 2 - 2 J max a max S i ) - a max 2 2 J max - - - ( 23 )
V i , T _ f = V i _ s + J max T jerk 2 - - - ( 24 )
J max T jerk 3 + 2 V i , T _ f - S i = 0 - - - ( 25 )
Wherein: T JerkFor adding the time of accelerating sections;
(5) calculate to satisfy speed based on the linking point of adjacent two sections tracks of the path length constraint condition of the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control, and as the end value of linking speed between adjacent two sections orbit segments.
Described path length constraint condition based on the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control is inequality group (26), (27) and (28):
V i , T _ f &GreaterEqual; a max 2 J max S i + 1 > 2 V i , T _ f a max J max + a max 3 J max 2 S i + 1 < a max 2 V i , T _ f - J max V i , T _ f 2 2 J max a max - - - ( 26 )
V i , T _ f < a max 2 J max S i + 1 < V i , T _ f V i , T _ f J max - - - ( 27 )
V i , T _ f &GreaterEqual; a max 2 J max S i + 1 &le; 2 V i , T _ f a max J max + a max 3 J max 2 - - - ( 28 )
Wherein: S I+1It is the path of i+1 section track; V I, T_fFor satisfying speed based on the linking point of adjacent two sections tracks of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control;
When (26), (27) and (28) when all being false, the speed V that satisfies based on the linking point of adjacent two sections tracks of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control is described I, T_fCan satisfy path length constraint condition, satisfy speed V based on the linking point of adjacent two sections tracks of the path length constraint condition of the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control based on the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control I, T_lEqual V I, T_f
When inequality group (26), (27) and (28) by and only form immediately by one, the speed V that satisfies based on the linking point of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control is described I, T_fCan not satisfy path length constraint condition, satisfy speed V based on the linking point of adjacent two sections tracks of the path length constraint condition of the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control based on the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control I, T_lComputing method are: when inequality group (26) is set up, and V I, T_lCalculate according to formula (29); When inequality group (27) or (28) establishment, V I, T_lCalculate according to formula formula (30);
V i , T _ l = a max 2 - a max 4 - 8 J max 2 S i + 1 2 J max - - - ( 29 )
V i , T _ l = 3 J max S i + 1 2 - - - ( 30 )
Repeating step (3) is to (5), and all program segments dispose in the decoding buffer zone.
Explanation is that a kind of intersegmental linking velocity process of continuous path method that is used for digital control system of the present invention is not limited to the foregoing description at last, can also make various modifications, conversion and distortion.Therefore, instructions and accompanying drawing are regarded in an illustrative, rather than a restrictive.Every foundation technical scheme of the present invention is made amendment, modification or equivalent variations, and does not break away from the thought and the scope of technical solution of the present invention, and it all should be encompassed in the middle of the claim scope of the present invention.

Claims (4)

1. the intersegmental linking velocity process of continuous path method that is used for digital control system, digital control system in the numerical control device is deciphered processing to the numerical control program in machine code, and store decode results into the decoding buffer zone, it is characterized in that, utilize the program segment information of deciphering in the buffer zone to the intersegmental linking velocity process of continuous path to be processed, make that the numerical control device machining locus is continuous, rate smoothing, the intersegmental linking velocity process of described continuous path method may further comprise the steps:
(1) determines that adjacent two sections tracks are connected model, and obtain the kinematic parameter of adjacent two sections tracks to be processed;
(2) calculate the path of adjacent two sections tracks to be processed;
(3) speed of the linking point of the adjacent two section tracks of calculating under acceleration, bow high level error, instruction speed of feed and turning constraint condition;
(4) calculate to satisfy speed based on the linking point of adjacent two sections tracks of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control;
(5) calculate to satisfy speed based on the linking point of adjacent two sections tracks of the path length constraint condition of the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control, and as the end value of linking speed between adjacent two sections orbit segments.
2. a kind of intersegmental linking velocity process of continuous path method that is used for digital control system according to claim 1 is characterized in that, the path of described adjacent two sections tracks to be processed comprises straight line path length and circular arc path length; Described circular arc path length calculation step comprises: (2.1) calculate the positive dirction angle in the horizontal axis of circular arc initial point place principal normal unit vector and workpiece coordinate system, (2.2) calculating is at the positive dirction angle of the horizontal axis of circular arc end point place's principal normal unit vector and workpiece coordinate system, and (2.3) calculate circular arc path length;
(2.1) calculate in circular arc initial point place principal normal unit vector
Figure FSA00000439941700011
Positive dirction angle α with the horizontal axis of workpiece coordinate system I_s:
Figure FSA00000439941700012
Figure FSA00000439941700013
Wherein: i=1,2,3...N, N represent orbit segment number to be processed; (x i, y i) be i section Track Initiation point coordinate; (x I_c, y I_c) be the central coordinate of circle of circular arc when i section track is circular arc; R iBe the radius of circular arc when i section track is circular arc, get α by formula (1) and (2) I_sComputing formula be:
(2.2) calculate principal normal unit vector at some place, circular arc end
Figure FSA00000439941700015
Positive dirction angle α with the horizontal axis of workpiece coordinate system I_e:
Figure FSA00000439941700022
Wherein: (x I+1, y I+1) be i section track end point coordinate;
Get α by formula (4) and (5) I_eComputing formula be:
Figure FSA00000439941700023
(2.3) calculate circular arc path length S i:
Get S by formula (3) and (6) iComputing formula be:
Figure FSA00000439941700024
3. a kind of intersegmental linking velocity process of continuous path method that is used for digital control system according to claim 1, it is characterized in that described path length constraint condition based on the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control is inequality group (8), (9) and (10):
Figure FSA00000439941700025
Figure FSA00000439941700026
Figure FSA00000439941700027
Wherein: V I_sBe the initial point speed of i section track; a MaxPeak acceleration for the digital control system permission; J MaxMaximum acceleration for the digital control system permission; S ' iIt is the path of i section track; V I, TSpeed for the linking point of adjacent two sections tracks of under acceleration, bow high level error, instruction speed of feed and turning constraint condition, calculating;
When inequality group (8), (9) and (10) when all being false, the speed V of the linking point of the adjacent two sections tracks that calculate is described under acceleration, bow high level error, instruction speed of feed and turning constraint condition I, TCan satisfy path length constraint condition, satisfy speed V based on the linking point of adjacent two sections tracks of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control based on the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control I, T_fEqual V I, T
When inequality group (8), (9) and (10) by and only form immediately by one, the speed V of the linking point of the adjacent two sections tracks that calculate under acceleration, bow high level error, instruction speed of feed and turning constraint condition is described I, TCan not satisfy path length constraint condition, satisfy speed V based on the linking point of adjacent two sections tracks of the path length constraint of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control based on the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control I, T_fComputing method are: when inequality group (8) is set up, and V I, T_fCalculate according to formula (11); When inequality group (9) or (10) establishment, V I, T_fCalculate according to formula (12) and (13);
Figure FSA00000439941700031
Figure FSA00000439941700033
Wherein: T JerkFor adding the time of accelerating sections.
4. a kind of intersegmental linking velocity process of continuous path method that is used for digital control system according to claim 1, it is characterized in that described path length constraint condition based on the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control is inequality group (14), (15) and (16):
Figure FSA00000439941700034
Figure FSA00000439941700035
Figure FSA00000439941700041
Wherein: S I+1It is the path of i+1 section track; V I, T_fFor satisfying speed based on the linking point of adjacent two sections tracks of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control;
When (14), (15) and (16) when all being false, the speed V that satisfies based on the linking point of adjacent two sections tracks of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control is described I, T_fCan satisfy path length constraint condition, satisfy speed V based on the linking point of adjacent two sections tracks of the path length constraint condition of the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control based on the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control I, T_fEqual V I, T_f
When inequality group (14), (15) and (16) by and only form immediately by one, the speed V that satisfies based on the linking point of the path length constraint condition of the leading portion track in adjacent two sections tracks of S curve acceleration and deceleration control is described I, T_fCan not satisfy path length constraint condition, satisfy speed V based on the linking point of adjacent two sections tracks of the path length constraint condition of the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control based on the back segment track in adjacent two sections tracks of S curve acceleration and deceleration control I, T_fComputing method are: when inequality group (14) is set up, and V I, T_fCalculate according to formula (17); When inequality group (15) or (16) establishment, V I, T_fCalculate according to formula formula (18);
Figure FSA00000439941700042
Figure FSA00000439941700043
CN2011100459381A 2011-02-25 2011-02-25 Continuous locus intersegment connecting speed processing method for numerical control system Expired - Fee Related CN102141794B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011100459381A CN102141794B (en) 2011-02-25 2011-02-25 Continuous locus intersegment connecting speed processing method for numerical control system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011100459381A CN102141794B (en) 2011-02-25 2011-02-25 Continuous locus intersegment connecting speed processing method for numerical control system

Publications (2)

Publication Number Publication Date
CN102141794A true CN102141794A (en) 2011-08-03
CN102141794B CN102141794B (en) 2012-11-21

Family

ID=44409376

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011100459381A Expired - Fee Related CN102141794B (en) 2011-02-25 2011-02-25 Continuous locus intersegment connecting speed processing method for numerical control system

Country Status (1)

Country Link
CN (1) CN102141794B (en)

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103064344A (en) * 2012-12-25 2013-04-24 广东省自动化研究所 Non uniform rational B spline (NURBS) interpolation based speed smooth control method
CN103135501A (en) * 2013-02-01 2013-06-05 北京配天大富精密机械有限公司 Acceleration and deceleration controlling method based on S-shaped curve and equipment using the same and numerically-controlled machine tool
CN103353761A (en) * 2013-07-12 2013-10-16 北京配天大富精密机械有限公司 Continuous processing method and apparatus of robot and smooth switching method and apparatus
CN103353737A (en) * 2013-07-12 2013-10-16 北京配天大富精密机械有限公司 Continuous processing method and apparatus of robot and smooth switching method and apparatus
CN103365244A (en) * 2013-07-12 2013-10-23 北京配天大富精密机械有限公司 Robot continuous machining method and device and smooth switching method and device
CN104932424A (en) * 2015-07-13 2015-09-23 苏州谷夫道自动化科技有限公司 CNC numerical control equipment tool path corner smoothing processing method
CN105045211A (en) * 2015-07-13 2015-11-11 太原科技大学 Equal-chord-error variable-step tangent interpolation method
CN105500354A (en) * 2016-02-02 2016-04-20 南京埃斯顿机器人工程有限公司 Transitional track planning method applied by industrial robot
CN105892408A (en) * 2015-09-07 2016-08-24 广东省自动化研究所 Bender speed planning method and device
CN106444637A (en) * 2015-08-04 2017-02-22 深圳市雷赛智能控制股份有限公司 Micro line segment dynamic prospect control method and micro line segment dynamic prospect control device
CN108415427A (en) * 2018-02-11 2018-08-17 昆山艾派科技有限公司 The path joining method that robot multistep continuously moves
CN109416529A (en) * 2017-12-19 2019-03-01 深圳配天智能技术研究院有限公司 Digital control system, the device with store function and nurbs curve interpolating method
CN109732219A (en) * 2018-12-05 2019-05-10 大族激光科技产业集团股份有限公司 A kind of rounded corner method and system of laser cutting
CN110450160A (en) * 2019-08-14 2019-11-15 合肥工业大学 A kind of open control method for exempting to program robot
CN111722591A (en) * 2020-06-05 2020-09-29 浙江工业大学 High-precision linkage interpolation method for trademark die-cutting machine
CN112276949A (en) * 2020-10-21 2021-01-29 哈工大机器人(合肥)国际创新研究院 Adjacent joint space-Cartesian space trajectory transition method and device
CN112327954A (en) * 2020-11-14 2021-02-05 大连理工大学 High-precision positioning method for linear motor controlled by asymmetric S-shaped speed curve
CN112650145A (en) * 2020-12-24 2021-04-13 苏州谋迅智能科技有限公司 S-shaped speed curve self-adaptive combination evaluation method
CN112947296A (en) * 2021-03-10 2021-06-11 深圳市泰达智能装备有限公司 Five-axis speed planning method and device, computer equipment and storage medium
CN113272756A (en) * 2020-09-08 2021-08-17 深圳市大疆创新科技有限公司 Holder control method and device, holder and storage medium
CN113341886A (en) * 2021-05-25 2021-09-03 深圳市汇川技术股份有限公司 Smooth feed speed planning method, device and computer readable storage medium
CN113946139A (en) * 2020-07-17 2022-01-18 配天机器人技术有限公司 Speed prediction method of numerical control system, control method of numerical control system and numerical control system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1235126A1 (en) * 2001-02-26 2002-08-28 Hitachi, Ltd. Numerically controlled curved surface machining unit
CN1811632A (en) * 2005-12-19 2006-08-02 深圳信息职业技术学院 A digital control code encoder and method for establishing digital control system software based on the same
CN1967421A (en) * 2005-11-18 2007-05-23 中国科学院沈阳计算技术研究所有限公司 Dynamic forward-looking processing method of small program segment and implementing device
JP2007172325A (en) * 2005-12-22 2007-07-05 Toshiba Mach Co Ltd Method of machining free curve and numerical control device
JP2010146176A (en) * 2008-12-17 2010-07-01 Fanuc Ltd Numerical control device for machine tool having rotating shaft

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1235126A1 (en) * 2001-02-26 2002-08-28 Hitachi, Ltd. Numerically controlled curved surface machining unit
CN1967421A (en) * 2005-11-18 2007-05-23 中国科学院沈阳计算技术研究所有限公司 Dynamic forward-looking processing method of small program segment and implementing device
CN1811632A (en) * 2005-12-19 2006-08-02 深圳信息职业技术学院 A digital control code encoder and method for establishing digital control system software based on the same
JP2007172325A (en) * 2005-12-22 2007-07-05 Toshiba Mach Co Ltd Method of machining free curve and numerical control device
JP2010146176A (en) * 2008-12-17 2010-07-01 Fanuc Ltd Numerical control device for machine tool having rotating shaft

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《中国机械工程》 20100131 潘海鸿等 全程S曲线加减速控制的自适应分段NURBS曲线插补算法 190-195 1-4 第21卷, 第2期 *

Cited By (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103064344B (en) * 2012-12-25 2014-11-05 广东省自动化研究所 Non uniform rational B spline (NURBS) interpolation based speed smooth control method
CN103064344A (en) * 2012-12-25 2013-04-24 广东省自动化研究所 Non uniform rational B spline (NURBS) interpolation based speed smooth control method
CN103135501A (en) * 2013-02-01 2013-06-05 北京配天大富精密机械有限公司 Acceleration and deceleration controlling method based on S-shaped curve and equipment using the same and numerically-controlled machine tool
CN103135501B (en) * 2013-02-01 2015-10-21 北京配天技术有限公司 Based on the Acceleration-deceleration Control Method of S type curve and device and numerically-controlled machine
CN103353761B (en) * 2013-07-12 2016-01-20 北京配天技术有限公司 Robot continuous process, device and level and smooth forwarding method, device
CN103365244A (en) * 2013-07-12 2013-10-23 北京配天大富精密机械有限公司 Robot continuous machining method and device and smooth switching method and device
CN103353737A (en) * 2013-07-12 2013-10-16 北京配天大富精密机械有限公司 Continuous processing method and apparatus of robot and smooth switching method and apparatus
CN103365244B (en) * 2013-07-12 2015-10-21 北京配天技术有限公司 Robot continuous process, device and level and smooth forwarding method, device
CN103353761A (en) * 2013-07-12 2013-10-16 北京配天大富精密机械有限公司 Continuous processing method and apparatus of robot and smooth switching method and apparatus
CN103353737B (en) * 2013-07-12 2016-08-10 北京配天技术有限公司 Robot continuous process, device and smooth forwarding method, device
CN104932424A (en) * 2015-07-13 2015-09-23 苏州谷夫道自动化科技有限公司 CNC numerical control equipment tool path corner smoothing processing method
CN105045211A (en) * 2015-07-13 2015-11-11 太原科技大学 Equal-chord-error variable-step tangent interpolation method
CN104932424B (en) * 2015-07-13 2019-10-29 苏州谷夫道自动化科技有限公司 CNC numerical control device cutter track corner smoothing processing method
CN105045211B (en) * 2015-07-13 2017-11-07 太原科技大学 One kind etc. bends high error variable step tangent line interpolating method
CN106444637A (en) * 2015-08-04 2017-02-22 深圳市雷赛智能控制股份有限公司 Micro line segment dynamic prospect control method and micro line segment dynamic prospect control device
CN105892408A (en) * 2015-09-07 2016-08-24 广东省自动化研究所 Bender speed planning method and device
CN105892408B (en) * 2015-09-07 2019-08-06 广东省智能制造研究所 A kind of bending machine speed planning method and device
CN105500354A (en) * 2016-02-02 2016-04-20 南京埃斯顿机器人工程有限公司 Transitional track planning method applied by industrial robot
CN109416529A (en) * 2017-12-19 2019-03-01 深圳配天智能技术研究院有限公司 Digital control system, the device with store function and nurbs curve interpolating method
CN109416529B (en) * 2017-12-19 2021-12-07 深圳配天智能技术研究院有限公司 Numerical control system, device with storage function and NURBS curve interpolation method
CN108415427A (en) * 2018-02-11 2018-08-17 昆山艾派科技有限公司 The path joining method that robot multistep continuously moves
CN109732219A (en) * 2018-12-05 2019-05-10 大族激光科技产业集团股份有限公司 A kind of rounded corner method and system of laser cutting
CN109732219B (en) * 2018-12-05 2021-08-27 大族激光科技产业集团股份有限公司 Laser cutting rounding method and system
CN110450160A (en) * 2019-08-14 2019-11-15 合肥工业大学 A kind of open control method for exempting to program robot
CN111722591A (en) * 2020-06-05 2020-09-29 浙江工业大学 High-precision linkage interpolation method for trademark die-cutting machine
CN113946139A (en) * 2020-07-17 2022-01-18 配天机器人技术有限公司 Speed prediction method of numerical control system, control method of numerical control system and numerical control system
CN113272756A (en) * 2020-09-08 2021-08-17 深圳市大疆创新科技有限公司 Holder control method and device, holder and storage medium
WO2022051916A1 (en) * 2020-09-08 2022-03-17 深圳市大疆创新科技有限公司 Gimbal control method and apparatus, gimbal, and storage medium
CN112276949A (en) * 2020-10-21 2021-01-29 哈工大机器人(合肥)国际创新研究院 Adjacent joint space-Cartesian space trajectory transition method and device
CN112327954A (en) * 2020-11-14 2021-02-05 大连理工大学 High-precision positioning method for linear motor controlled by asymmetric S-shaped speed curve
CN112650145A (en) * 2020-12-24 2021-04-13 苏州谋迅智能科技有限公司 S-shaped speed curve self-adaptive combination evaluation method
CN112947296A (en) * 2021-03-10 2021-06-11 深圳市泰达智能装备有限公司 Five-axis speed planning method and device, computer equipment and storage medium
CN112947296B (en) * 2021-03-10 2022-04-26 深圳市泰达智能装备有限公司 Five-axis speed planning method and device, computer equipment and storage medium
CN113341886A (en) * 2021-05-25 2021-09-03 深圳市汇川技术股份有限公司 Smooth feed speed planning method, device and computer readable storage medium

Also Published As

Publication number Publication date
CN102141794B (en) 2012-11-21

Similar Documents

Publication Publication Date Title
CN102141794B (en) Continuous locus intersegment connecting speed processing method for numerical control system
CN102147600B (en) Numerical control interpolation system for real-time generation of curvature-continuous path
CN103135501B (en) Based on the Acceleration-deceleration Control Method of S type curve and device and numerically-controlled machine
CN102419570B (en) Acceleration and deceleration look-ahead control method for high-speed machining of numerical control machine tool
CN100451891C (en) Dynamic forward-looking processing method of small program segment and implementing device
CN101497140B (en) Off-line planning method for cutting feed rate of five-shaft numerical control side milling machining
CN103699056B (en) The little line segment real-time smooth transition interpolation method of high-speed, high precision digital control processing
CN104615084B (en) Machining feed speed optimized tool path curve contour error compensation method
CN101510087B (en) Forward looking self-adapting speed controlling method for high-speed processing tiny line segment
Wang et al. A look-ahead and adaptive speed control algorithm for high-speed CNC equipment
CN103728923A (en) High-efficiency high-stability flexible acceleration and deceleration control method for numerical control system
CN108829045B (en) Optimization method and system for connection speed of continuous micro-linear segment
CN104678899A (en) Curve velocity planning method and device, as well as numerical control machining path data processing method
CN110147077B (en) Cosine interpolation method of industrial robot under space linear operation
US8838265B2 (en) Machine tool, machining method, program and NC data generation device
CN104160617A (en) Motor control device
CN101169647A (en) Mini line segment curve five-axis linkage machining speed smoothing method
CN104281099A (en) NURBS direct interpolation method and device with processing characteristics considered
Zhao et al. Pocketing toolpath optimization for sharp corners
CN106292529A (en) A kind of processing path generation method of lathe
JP2011158982A (en) Controller for machine tool
CN109901518B (en) Method for planning acceleration and deceleration speed of numerical control machine tool under constant force constraint condition
CN102707671A (en) Processing path optimization method applied to machine tool
CN110879569A (en) Prospective control method
CN105629882A (en) Trigonometric function speed planning method used for spline interpolation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20121121

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