KR100664681B1 - The direct interpolation system and method using sampled data of the robot path described with parametric curve - Google Patents
The direct interpolation system and method using sampled data of the robot path described with parametric curve Download PDFInfo
- Publication number
- KR100664681B1 KR100664681B1 KR1020050079861A KR20050079861A KR100664681B1 KR 100664681 B1 KR100664681 B1 KR 100664681B1 KR 1020050079861 A KR1020050079861 A KR 1020050079861A KR 20050079861 A KR20050079861 A KR 20050079861A KR 100664681 B1 KR100664681 B1 KR 100664681B1
- Authority
- KR
- South Korea
- Prior art keywords
- length
- parameter
- curve
- interpolation
- parametric
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/41—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by interpolation, e.g. the computation of intermediate points between programmed end points to define the path to be followed and the rate of travel along that path
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1656—Programme controls characterised by programming, planning systems for manipulators
- B25J9/1664—Programme controls characterised by programming, planning systems for manipulators characterised by motion, path, trajectory planning
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/34—Director, elements to supervisory
- G05B2219/34119—Function generator, filter after interpolator to control position error
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/34—Director, elements to supervisory
- G05B2219/34169—Coarse interpolator, path calculator delivers position, speed, acceleration blocks
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/49—Nc machine tool, till multiple
- G05B2219/49199—For non linear interpolation movement
Abstract
Description
도 1은 본 발명에 의한 로봇의 경로 계획 시스템 구성도.1 is a schematic diagram of a path planning system of a robot according to the present invention;
도 2a 및 도 2b는 본 발명에 의한 커브의 길이를 구하는 알고리즘 순서도2A and 2B are algorithm flow charts for calculating the length of a curve according to the present invention.
<도면의 주요부분에 대한 부호의 설명><Description of the symbols for the main parts of the drawings>
10 : 오퍼레이터 입력부 20 : 경로계획부10: operator input unit 20: route planning unit
21 : 매개변수 곡선 생성부 22 : 길이 계산부21: parameter curve generation unit 22: length calculation unit
23 : 길이-파라메타 테이블 24 : 보간부23 length-parameter table 24 interpolation unit
25 : 좌표 검출부 26 : 이동시간 계산부 25: coordinate detection unit 26: travel time calculation unit
27 : 속력 프로파일 생성부 28 : 주행거리 계산부27: speed profile generation unit 28: driving distance calculation unit
29 : 조인트 각도 레퍼런스 생성부 30 : 제어부29: joint angle reference generation unit 30: control unit
본 발명은 매개변수곡선으로 기술된 로봇 경로의 샘플 테이터를 이용한 직접 보간 시스템 및 방법에 관한 것으로서, 특히 로봇의 경로를 계획하는 단계에서 매개변수 곡선(Parametric Curve)의 길이(Length)를 구하는 방법과 주행거리에 해당하는 파라미터를 구하는 방법에 관한 것이다.The present invention relates to a direct interpolation system and method using sample data of a robot path described by a parametric curve, and more particularly, to a method for obtaining a length of a parametric curve in a step of planning a robot path. The present invention relates to a method for obtaining a parameter corresponding to a mileage.
현대의 캐드 시스템은 2D 혹은 3D 매개변수 곡선을 이용한 디자인 툴을 제공한다. CNC 머신이나 산업용 로봇(industrial robot)은 디자인된 경로 데이터(path data)로부터 머신이 움직일 궤적(trajectory)를 생성해야 한다. 특히 고속(High-speed), 고 정밀 기계가공(high precision machining)을 위해서는 생성된 궤적의 속도 변동(speed fluctuation) 및 윤곽 에러(contour error)가 최소화 되어야한다.Modern CAD systems provide design tools using 2D or 3D parametric curves. CNC machines or industrial robots need to generate trajectory for the machine to move from the designed path data. In particular, for high-speed and high precision machining, the speed fluctuation and contour error of the generated trajectories should be minimized.
이를 위해서 다양한 보간(Interpolation) 기법이 개발되었는데, 크게 두 가지 카테고리로 나눌 수 있다. 하나는 경로(path)를 잘 표현할 수 있는 샘플데이터(sampled data)를 잇는 라인(line) 혹은 스플라인(spline)곡선 등의 근사 곡선(approximation curve)을 이용하여 보간 좌표를 얻는 방법이고, 나머지 하나는 실시간으로 파라메타를 업데이트(update)하면서 매개변수 곡선으로부터 직접보간 좌표를 얻는 방법이다. (이후, 전자의 방법을 경로근사화 방법(Path Approximation Method), 후자의 방법을 실시간 파라메타 업데이트 방법(Realtime Parameter Update Method)라고 칭함.)To this end, various interpolation techniques have been developed, which can be divided into two categories. One method is to obtain interpolation coordinates using an approximation curve such as a line or spline curve connecting sampled data that can express a path well. This is a method of obtaining interpolation coordinates directly from a parameter curve while updating parameters in real time. (Hereinafter, the former method is called Path Approximation Method and the latter method is called Realtime Parameter Update Method.)
초창기 CNC 머신(machine)은 매개변수로 표현된 경로를 직접 처리하기보다는 CAD데이터로부터 경로를 잘 보간할 수 있는 샘플 포인트를 뽑는 전처리 과정을 거친 후 그 샘플 포인트를 보간하는 방법을 사용하였다. 샘플 포인트를 라인 혹은 스플라인으로 보간함으로써 쉽게 보간의 구현이 가능했기 때문이다. 이런 까닭으로 현재도 많은 CNC 머신에서 샘플 데이터를 이용한 경로 보간 기법을 많이 쓴다. 선형적 보간(linear interpolation)은 경로를 아주 작은 선형 블록(linear block)으로 나누어 처리한다. 이 방법은 구현이 쉽고 간단하지만 다른 방법들에 비해 윤곽 에러(contour error) 및 이송-률 변동(feed-rate fluction)이 심하다. 스플라인 보간(spline interpolator)은 각 블록을 부드러운 스플라인 곡선으로 보간한다. 스플라인 보간에 대해서는 최근까지 많은 연구가 이뤄지고 있으며 선형 보간에 비하여 매우 훌륭한 성능을 보인다.Early CNC machines used a method of interpolating the sample points after preprocessing the sample points to interpolate the path from CAD data, rather than directly processing the parameterized paths. This is because interpolation can be easily implemented by interpolating sample points to lines or splines. For this reason, many CNC machines still use path interpolation using sample data. Linear interpolation divides the path into very small linear blocks. This method is easy and simple to implement, but has a greater contour error and feed-rate fluction than other methods. Spline interpolator interpolates each block into a smooth spline curve. Spline interpolation has been studied until recently and shows much better performance than linear interpolation.
그러나 스플라인곡선이 매개변수 곡선을 완벽하게 나타내지 못함에 의해서 생기는 윤곽 에러(contour error)는 스플라인 보간의 피할 수 없는 단점이다. 물론 나트 포인트(knot point)를 증가시키면 원하는 수준의 보간 곡선을 얻을 수 있지만 메모리와 계산량이 늘어나므로 바람직하지 못하다.However, the contour error caused by the spline curve not completely representing the parametric curve is an inevitable disadvantage of spline interpolation. Of course, increasing the knot point will yield the desired level of interpolation curve, but this is not desirable because of the increased memory and computation.
10여년 전부터는 실시간 파라메트릭 보간(realtime parametric interpolator)에 대한 연구가 많이 이뤄지고 있다. 실시간 파라메타 업데이트 방식(Real-time Parameter Update Method)은 샘플 포인트를 사용하지 않으므로 부가적인 메모리의 사용이 없을 뿐 아니라, 매개변수 곡선(parametric curve)으로부터 직접 보간 좌표를 얻으므로 결과적으로도 경로 보간 방법(path approximation method)에 비하여 매우 좋은 성능을 보인다. 그런데 B-스플라인(B-spline)등의 일반적인 매개변수 곡선(parametric curve)의 경우 주행거리에 해당하는 파라메터(parameter)를 얻어내는 식이 크로스 폼 솔루션(closed form solution)으로 구해지지 않는다. 실시간 업데이트(real-time update)에서는 매 샘플링 시간마다 다음 시 간까지의 주행거리의 증분이 원하는 이송-비(feed-rate)가 되도록 실시간으로 파라메터를 업데이트한다.More than a decade ago, research on real-time parametric interpolators has been conducted. Since the Real-time Parameter Update Method does not use sample points, there is no additional memory used, and the interpolation coordinates are obtained directly from the parametric curves. Compared to the path approximation method), the performance is very good. However, in the case of a general parametric curve such as B-spline, a formula for obtaining a parameter corresponding to a mileage cannot be obtained as a closed form solution. In a real-time update, every sampling time is updated in real time so that the increment of the mileage to the next time is the desired feed-rate.
즉, 경로를 완전히 주행하기 전까지는 경로의 길이를 알 수 없다. 따라서 실시간으로 남은 주행 거리를 추정하고 이송-비(feed-rate)를 조정하는 알고리즘이 부가적으로 필요하다. 게다가 길이 정보의 부재는 정확한 작업(task) 실행시간의 예측이 불가능하게 하므로 관련된 작업 프로그래밍(task programming)의 제약을 가져올 수 있다.That is, the length of the route is unknown until the route is completely driven. Therefore, an algorithm for estimating the remaining travel distance and adjusting the feed-rate in real time is additionally needed. In addition, the absence of length information makes it impossible to predict the exact task execution time, which can lead to constraints of the associated task programming.
본 발명은 경로 보간 방식(Path Approximation Method)의 결점인 윤곽 에러(contour error) 문제와 실시간 파라메타 업데이트 방식(Real-time Parameter Update Method)의 실시간 남은 거리 추정 문제를 해결할 수 있도록 두 방식의 단점을 해결하면서 동시에 장점을 포용할 수 있는 새로운 보간 방법(interpolation method)(이하, 파라메타 보간 방식(Parameter Approximation Method)이라 칭함) 및 그 시스템에 관한 것이다.The present invention solves the shortcomings of the two methods to solve the contour error problem, which is a drawback of the path interpolation method and the real-time remaining distance estimation problem of the real-time parameter update method. The present invention relates to a new interpolation method (hereinafter, referred to as a parameter interpolation method) and a system that can simultaneously take advantage of advantages.
파라메타 보간 방식(Parameter Approximation Method)은 경로 보간 방식(Path Approximation Method)에서처럼 샘플 포인트를 사용하여 각 샘플 포인트에서의 주행거리와 파라메터를 계산하지만, 보간좌표를 얻는 과정에서는 실시간 파라메타 업데이터 방식(Real-time Parameter Update Method)에서처럼 매개변수 곡선(parametric curve)로부터 직접 얻는다. 이때 샘플 포인트 사이의 주행거리에 해당 하는 파라미터는 샘플 포인트에서의 경계선 조건(boundary condition)을 이용하여 5차 방정식의 다항식으로 근사화하여 구한다.The Parameter Approximation Method calculates the mileage and parameters at each sample point using the sample points as in the Path Approximation Method.However, in the process of obtaining the interpolation coordinates, the real-time parameter update method is used. Directly from the parametric curve as in the Parameter Update Method. At this time, the parameter corresponding to the distance traveled between the sample points is obtained by approximating the polynomial of the fifth-order equation using the boundary condition at the sample point.
즉, 경로 보간 방식에서는 샘플 포인트 사이의 경로를 보간하여 보간된 커브(approximated curve)로부터 보간좌표를 구했다면, 파라메타 보간 방식은 샘플 포인트 사이의 주행거리와 파라메터 사이의 관계를 근사화하여 근사화된 파라메타를 이용하여 매개변수 곡선으로부터 직접 보간좌표를 얻는다. 이를 통해 경로 보간 방식에 의한 윤곽 에러와 실시간 파라메터 업데이트 방식에 의한 실시간 남은 주행거리 문제를 해결할 수 있다.That is, in the path interpolation method, if the interpolation coordinates are obtained from the interpolated curve by interpolating the path between sample points, the parameter interpolation method approximates the relationship between the mileage and the parameter between the sample points and approximates the parameter. To obtain the interpolation coordinates directly from the parametric curves. Through this, it is possible to solve the contour error caused by the path interpolation method and the remaining distance traveled by the real time parameter update method.
본 발명은 CNC머신이나 산업용 로보트(industrial robot)의 이동 경로를 계획하는 방법에 관한 것으로 경로의 곡률(curvature)을 고려하고 매개변수 곡선(parametric curve)의 길이를 구하고, 길이-파라메터 테이블(length to parameter table) 참조와 5차 다항식 보간(quintic polynomial interpolation)을 이용해서 주행거리에 대응하는 파라메타를 찾는 것을 특징으로 한다.The present invention relates to a method for planning a moving path of a CNC machine or an industrial robot, taking into account the curvature of the path, obtaining a length of a parametric curve, and measuring a length-parameter table. parameter table) and using a quintic polynomial interpolation to find the parameter corresponding to the mileage.
도 1은 본 발명을 설명하기 위한 로봇의 경로계획 시스템 구성도이다.1 is a configuration diagram of a route planning system of a robot for explaining the present invention.
도 1에 도시된 바와 같이, 로봇의 이동 경로계획 단계는 입력부(10)으로부터 매개변수 곡선의 종류와 경유점 등을 입력 받아 해당하는 파라메타 경로(parametric path)를 생성하는 파라메타 커브 생성부(21)와, 해당 커브의 전체 길이를 계산하는 길이 계산부(22)와, 샘플 포인트에서의 파라메타와 그 포인트까의 주행거리를 맵핑(mapping)하는 길이-파라메타 테이블(Length to Parameter table)(23)과, 이동거리 계산부(28)로부터 입력받은 이동거리에 해당하는 파라메타를 길이-파라메타 테이블(23)의 테이블과 5차 다항식 보간법(quintic polynomial interpolation)을 이용하여 구하는 보간(Interporation)부(24)와, 보간부(24)에서 입력받은 파라메타의 직교(Cartesian) 좌표의 위치, 속도, 가속도를 파라메타 커브 생성부(21)에서 생성한 파라메타 커브로부터 직접 구하는 보간좌표 계산부(25)와, 로봇 제어에 필요한 레퍼런스(reference)를 획득할 주기를 설정하는 이동시간계산부(26), 상기 길이계산부(22)에서 계산된 길이를 이용하여 속력 프로파일을 생성하는 속력 프로파일(speed Profile)생성부(27)과, 길이계산부(26)에서 생성한 샘플링 시간(sampling time)마다 속력 프로파일 생성부(27)의 속력 프로파일을 참조하여 주행거리를 계산하는 이동거리계산부(28)과, 상기 보간좌표 생성부(25)에서 생성한 좌표를 모터제어에 필요한 조인트 각도 레퍼런스(joint angle reference)로 변환하는 레퍼런스(reference) 생성부(29)로 구성된다.As shown in FIG. 1, in the movement path planning step of the robot, the
이와 같이 구성되는 시스템은 기본적인 로봇 제어 시스템으로서, 오퍼레이터 입력부(10)는 운전자가 입력한 모션명령을 파라메터 커브의 종류, 경유점, 포인트개수 등으로 정형화 하는 것이고, 경로계획부(20)는 입력정보에 의해 로봇이 이동해야할 경로와 속력 프로파일을 생성하여 경로계획을 하게 된다. 그리고 제어부(30)는 경로계획부(20)에서 입력 받은 조인트 각도 레퍼런스(joint angle reference)를 잘 추종할 수 있게끔 제어하게 된다.The system configured as described above is a basic robot control system, in which the
상기 경로계획부(20)는 제어부(30)를 포함하는 콘트롤러내에 프로그램 모듈 로서 설치되는 것이며, 본 발명에서는 경로계획부(20)의 경로계획 중에 길이 계산부(22)에서 윤곽 에러(contour error)를 줄이는 커브의 길이를 생성하는 방법과 보간부(24)에서 주행거리에 해당하는 파라메타를 생성하는 방법에 주 요지가 있다.The
먼저, 파라메타 커브의 길이를 계산하는 방법에 대해 살펴보면 다음과 같다. First, the method of calculating the length of the parameter curve is as follows.
커브 길이는 기존의 경로 근사화 방법(Path approximation method)과 마찬가지로 샘플 데이터를 근사화(approximation)하여 계산한다. 하지만 커브를 등간격 샘플링하고, 샘플포인트 사이의 현 길이(chord length)를 근사화하여 커브의 길이를 구하는 기존의 방법 대신 곡률 차 의존 샘플링(curvature-difference dependent sampling)을 하고, 각 샘플 포인트 사이의 현의 길이(chordal-length)와 샘플 포인트에서의 곡률을 이용하여 그 사이의 호의 길이(arc-length)를 근사화하여 커브의 길이를 구하는 새로운 방법을 사용하였다.The curve length is calculated by approximating the sample data as in the path approximation method. However, instead of the conventional method of equally sampling the curve, approximating the chord length between sample points to find the length of the curve, perform curvature-difference dependent sampling, and chord between each sample point. Using the chordal-length and the curvature at the sample point, we used a new method to approximate the arc-length between them to find the length of the curve.
두 샘플 포인트 와 사이의 실제 길이를 , 현의 길이(chordal-length)를 , 추정된 길이를 라고 하자.Two sample points Wow The actual length between , The chordal-length , The estimated length Let's say
원호 같이 곡률이 일정한 곡선의 경우 현의 길이와 실제 길이는 곡률 를 이용하여 다음과 같이 표현가능하다.For curves with constant curvature, such as arcs, the length and actual length of the strings Can be expressed as follows.
, ,
만약 두 샘플 포인트 사이의 곡률이 일정하지 않다면 두 점 사이의 곡률의 대표값을 사용하여 추정할 수 밖에 없다. 두 점 사이에서의 곡률 변화가 단일변화(monotonic)하다고 가정하고, 곡률 대표값을 라 두자. 그러면 추정된 호의 길이(arc-length)는 다음과 같이 주어진다.If the curvature between two sample points is not constant, it can only be estimated using the representative value of the curvature between the two points. Assuming that the change in curvature between two points is monotonic, Let's do it. Then the estimated arc length is given by
여기서 추정오차와 윤곽 길이(chordal length) 및 곡률(curvature)와의 관계를 살펴볼 필요가 있다. 계산의 편의를 위해 실제길이 와 와의 관계가 와 같다고 가정하자.Here, it is necessary to examine the relationship between the estimation error, the chordal length, and the curvature. Actual length for convenience of calculation Wow Relationship with Assume that
( 단, )( only, )
그러면 길이추정 오차 바운드를 다음과 같이 얻을 수 있다.Then we can get the length estimation error bound as
= =
여기서 은 관계 길이 추정 에러 경계(relative length estimation error bound)로 부르며, 다음과 같이 풀어 쓸 수 있다.here Is called relative length estimation error bound, and can be solved as follows.
이제 샘플링 포인트의 결정에 대해서 생각해 보자. 이송-률 변동(feed-rate fluctuation)을 고려하여 샘플링하려고 한다면 두 샘플링 포인트 사이의 현의 길이(chordal length) 뿐만 아니라 현(chord)의 양 끝에서의 곡률을 고려하여 결정해야함을 알 수 있다. 그런데 현재 샘플링 포인트(current sampling point)에서 조건식을 만족하는 최적의 다음 샘플링 포인트를 찾는 것은 쉽지 않다. 따라서 보다 효율적인 구현을 위해서 1차로 간격으로 샘플링을 하고, 필요한 구간에 대해서는 재귀적으로 바이섹션(bisection)하는 방식을 제안한다. 이때 은 로보트(robot)가 트래킹(tracking)해야 할 최대 곡률(maximum curvature)을 고려하여 결정되어야 한다.Now consider the determination of the sampling point. If you want to sample with feed-rate fluctuation in mind, the chordal length between the two sampling points In addition, it can be seen that the decision should be made considering the curvature at both ends of the chord. However, finding the next optimal sampling point that satisfies the conditional expression at the current sampling point is not easy. Therefore, for a more efficient implementation, We propose a method of sampling at intervals and recursively bisectioning the necessary sections. At this time Must be determined in consideration of the maximum curvature that the robot should track.
도 2a 및 도 2b는 본 발명을 이상에서 설명한 방법을 파라메트릭 커브(parametric curve)의 길이를 구하는 알고리즘의 순서도이다. 2A and 2B are flowcharts of an algorithm for obtaining the length of a parametric curve in the method described above.
이에 도시된 바와 같이, 샘플링 포인트마다 그 포인트에서의 파라메터와 시작점에 그 포인트까지의 커브의 길이를 저장할 테이블을 초기화하는 단계(S20)와; 을 이용하여 다음 샘플 포인트 간격 를 결정하는 단계(S30)와;As shown therein, initializing a table for storing the parameters of the point and the length of the curve to the point at the starting point for each sampling point (S20); Next sample point interval using Determining (S30);
두 샘플링 포인트에서의 곡률 와 두 샘플링 포인트 사이의 대표 곡률 를 계산하는 단계(S41)와; 두 샘플링 포인트 사이 현의 길이(chordal lengt)h 를 계산하는 단계(S42)와; 길이 추정 오차 바운드 을 계산하는 단계(S43)와; 길이 추정 오차 바운드 이 주어진 허용 바운드(bound)보다 크거나 파라메타의 단일변화에 대한 가정이 깨지는 경우를 체크하는 단계(S44)와; (S44)에서의 조건이 거짓(false)이면 호의 길이(Arc-Length)를 계산하고 테이블에 오름차순으로 저장하는 단계(S45)와; (S44)에서의 조건이 참(true)이면 샘플링 포인트 간격을 바이섹션(bisection)하는 단계(S46)와; 샘플링 포인트 u를 업데이트() 하는 단계(S50)와; 종단 조건 를 체크하여 거짓(false)이면 (S30)부터 다시 수행하고 참(ture)이면 종료(S70)하는 단계(S60)을 수행하여 파라메트릭 곡률(Parametric Curve)의 전체길이와 길이-파라메타 테이블(length to parameter table)을 업데이트 한다.Curvature at two sampling points Representative curvature between and two sampling points Calculating (S41); Length of chord between two sampling points (chordal lengt) h Calculating (S42); Length estimation error bound Calculating (S43); Length estimation error bound Given allowable bound Checking a case in which an assumption about a single change of a parameter is larger or larger is broken (S44); If the condition at S44 is false, calculating a length of the arc Arc-Length and storing it in the table in ascending order (S45); Bisecting the sampling point interval if the condition in S44 is true; Update sampling point u ( Step (S50); Termination Condition If it is false, perform the operation again from (S30) if it is false, and finish (S70) if it is true (S70) to perform the full length and length-parameter table of the parametric curvature (Parametric Curve). update the parameter table).
다음으로, 주행거리에 해당하는 파라메타를 생성하는 방법에 대해 살펴보면 다음과 같다. Next, a method of generating a parameter corresponding to a driving distance is as follows.
보간좌표를 얻는 보간단계에서는 주행거리 에 해당하는 파라미터 를 찾아내야 한다. 일단 를 계산해 내면 그 시간에서의 위치, 속도 및 가속도를 구할 수 있다. 우리는 먼저 가 어느 샘플 포인트 사이에 위치하는지 알아야한다. 테이블의 데이터는 모두 오름차순으로 저장되어 있으므로, 가 속한 구간을 찾는 것은 다음 부등식을 만족하는 를 찾음으로서 해결할 수 있다.In the interpolation step to obtain interpolation coordinates Corresponding parameter Should be found. First We can calculate the position, velocity and acceleration at that time. We first You need to know which sample point is located between. Since all data in the table is stored in ascending order, Finding the interval to which it belongs satisfies the following inequality This can be solved by finding
테이블이 오름차순으로 정렬되어 있으므로 일반적으로 바이너리 서치(binary search)를 통하여 찾을 수 있다. 일단 가 속한 구간을 찾으면 그 구간의 양 끝점에서 파라미터의 길이에 대한 1차 미분과 2차 미분 정보를 이용하여 파라미터와 길의와의 관계를 5차 방정식의 다항식으로 근사화하여 에 해당하는 를 구한다. 파라메타 는 의 5차 함수를 이용하여 근사화(approximation)할 수 있다.Since the tables are sorted in ascending order, they can usually be found by binary search. First After finding the section to which the term belongs, the relation between the parameter and the path is approximated by the polynomial of the fifth order equation using the first and second derivatives of the parameter lengths at both ends of the interval. Equivalent to Obtain Parameter Is Approximation can be made using the 5th order function of.
, ,
위 식의 각 계수는 다음의 여섯가지 경계조건(boundary condition)으로부터 구해진다.Each coefficient of the above equation is obtained from the following six boundary conditions.
이상에서 상세히 설명한 바와 같이 본 발명에 의하면 일정한 간격의 샘플링 포인트를 길이 추정 오차 바운드 가 허용 기준 바운드 를 넘으면 그 구간을 둘로 나누는(bisection) 방법으로 재조정하여 샘플링 포인트 사이의 호 길이를 구하기 때문에 파라메트릭 커브를 등간격으로 샘플링하여 커브의 길이를 구하는 기존의 방법보다 윤곽 에러(contour error)를 줄일 수 있다. As described in detail above, according to the present invention, a sampling interval of a constant interval is bound to a length estimation error. Bound Tolerance If it exceeds, it re-adjusts the section by bisection to find the arc length between sampling points, so it is possible to reduce the contour error compared to the conventional method of sampling the parametric curves at equal intervals. have.
또한, 에 해당하는 를 5차 다항식으로 근사화하여 구하기 때문에 선형으로 보간하는 방법보다 근사화 오차를 줄일 수 있다.Also, Equivalent to Since it is approximated by the fifth-order polynomial, the approximation error can be reduced more than the linear interpolation method.
본 발명을 CNC 머신이나 산업용 로보트의 이동 경로를 계획하는 단계에 적용하면 더욱 정교한 경로의 트래킹 성능 높일 수 있는 효과가 있다.Applying the present invention to the step of planning the movement path of the CNC machine or industrial robot has the effect of increasing the tracking performance of the more sophisticated path.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020050079861A KR100664681B1 (en) | 2005-08-30 | 2005-08-30 | The direct interpolation system and method using sampled data of the robot path described with parametric curve |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020050079861A KR100664681B1 (en) | 2005-08-30 | 2005-08-30 | The direct interpolation system and method using sampled data of the robot path described with parametric curve |
Publications (1)
Publication Number | Publication Date |
---|---|
KR100664681B1 true KR100664681B1 (en) | 2007-01-04 |
Family
ID=37866897
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020050079861A KR100664681B1 (en) | 2005-08-30 | 2005-08-30 | The direct interpolation system and method using sampled data of the robot path described with parametric curve |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR100664681B1 (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109227547A (en) * | 2018-11-01 | 2019-01-18 | 珠海格力电器股份有限公司 | A kind of progress control method of robot, device, storage medium and robot |
CN109968357A (en) * | 2019-03-27 | 2019-07-05 | 杭州电子科技大学 | The quick planing method of transient motion between a kind of straight path |
KR20200003480A (en) * | 2018-07-02 | 2020-01-10 | 한국전자통신연구원 | Apparatus and method of producing data for automatic perforation |
CN113835423A (en) * | 2020-06-08 | 2021-12-24 | 苏州艾吉威机器人有限公司 | B-spline curve-based path planning method, motion control method and system |
CN114265367A (en) * | 2021-12-31 | 2022-04-01 | 杭州电子科技大学 | High-precision parameter path interpolation method and system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR19980054690A (en) * | 1996-12-27 | 1998-09-25 | 추호석 | Circular interpolation method and device for articulated robot |
KR19990059974A (en) * | 1997-12-31 | 1999-07-26 | 윤종용 | Path planning device and path planning method for robot control |
JPH11219350A (en) | 1998-02-03 | 1999-08-10 | Denso Corp | Interpolation value arithmetic unit, electronic controller and interpolation data table adaptation device |
KR20010035959A (en) * | 1999-10-05 | 2001-05-07 | 윤덕용 | 3-dimensional linear and circular interpolation method using the intersection criteria of two surfaces |
-
2005
- 2005-08-30 KR KR1020050079861A patent/KR100664681B1/en not_active IP Right Cessation
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR19980054690A (en) * | 1996-12-27 | 1998-09-25 | 추호석 | Circular interpolation method and device for articulated robot |
KR19990059974A (en) * | 1997-12-31 | 1999-07-26 | 윤종용 | Path planning device and path planning method for robot control |
JPH11219350A (en) | 1998-02-03 | 1999-08-10 | Denso Corp | Interpolation value arithmetic unit, electronic controller and interpolation data table adaptation device |
KR20010035959A (en) * | 1999-10-05 | 2001-05-07 | 윤덕용 | 3-dimensional linear and circular interpolation method using the intersection criteria of two surfaces |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20200003480A (en) * | 2018-07-02 | 2020-01-10 | 한국전자통신연구원 | Apparatus and method of producing data for automatic perforation |
KR102534879B1 (en) * | 2018-07-02 | 2023-05-22 | 한국전자통신연구원 | Apparatus and method of producing data for automatic perforation |
CN109227547A (en) * | 2018-11-01 | 2019-01-18 | 珠海格力电器股份有限公司 | A kind of progress control method of robot, device, storage medium and robot |
CN109227547B (en) * | 2018-11-01 | 2020-08-04 | 珠海格力电器股份有限公司 | Robot operation control method and device, storage medium and robot |
CN109968357A (en) * | 2019-03-27 | 2019-07-05 | 杭州电子科技大学 | The quick planing method of transient motion between a kind of straight path |
CN109968357B (en) * | 2019-03-27 | 2022-04-05 | 杭州电子科技大学 | Rapid planning method for transition motion between linear tracks |
CN113835423A (en) * | 2020-06-08 | 2021-12-24 | 苏州艾吉威机器人有限公司 | B-spline curve-based path planning method, motion control method and system |
CN113835423B (en) * | 2020-06-08 | 2024-03-26 | 苏州艾吉威机器人有限公司 | Path planning method, motion control method and system based on B spline curve |
CN114265367A (en) * | 2021-12-31 | 2022-04-01 | 杭州电子科技大学 | High-precision parameter path interpolation method and system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101493687B (en) | Real time forward looking whole-process acceleration and deceleration controlled NURBS curve self-adapting subsection interpolation method | |
US8774965B2 (en) | Method and device for controlling a manipulator | |
Pateloup et al. | Bspline approximation of circle arc and straight line for pocket machining | |
Nam et al. | A study on a generalized parametric interpolator with real-time jerk-limited acceleration | |
Lin | Real-time surface interpolator for 3-D parametric surface machining on 3-axis machine tools | |
KR100743143B1 (en) | apparatus for transition trajectory generation between non-smoothly connected paths of Industrial Robots | |
CN110865610B (en) | Tool path interpolation and speed planning method based on machine tool vibration suppression | |
KR100664681B1 (en) | The direct interpolation system and method using sampled data of the robot path described with parametric curve | |
US4683543A (en) | Time-based interpolation control of a robot | |
CN109765887A (en) | A kind of automatic Pilot control method | |
US10416613B2 (en) | Machining time estimating apparatus | |
CN104331025A (en) | Velocity look-ahead preprocessing method oriented to efficient processing of tiny line segments | |
CN110147077B (en) | Cosine interpolation method of industrial robot under space linear operation | |
Ye et al. | Interpolation of continuous micro line segment trajectories based on look-ahead algorithm in high-speed machining | |
JP2009053926A (en) | Route-planning device and route-planning method | |
JPH0969006A (en) | Curve interpolating method for speed control at connecting operation of robot | |
CN112549027B (en) | Cartesian space and joint space curve smooth transition method and device | |
CN109313429A (en) | S type speed planning method, apparatus, system, robot and numerically-controlled machine tool | |
CN106647636A (en) | Numerical controller performing speed control with curvature and curvature change amount | |
JP2790643B2 (en) | Numerical control unit | |
CN108513636B (en) | Robot motion control method and related device | |
CN115202293B (en) | Two-section type speed planning method for industrial robot | |
CN114002996B (en) | C3 continuous five-axis path switching fairing method for hybrid robot | |
JP2001051708A (en) | Numerical control method | |
CN113156962B (en) | Motion control method, motion control device, robot and storage medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A201 | Request for examination | ||
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant | ||
FPAY | Annual fee payment |
Payment date: 20121224 Year of fee payment: 7 |
|
FPAY | Annual fee payment |
Payment date: 20131125 Year of fee payment: 8 |
|
FPAY | Annual fee payment |
Payment date: 20141119 Year of fee payment: 9 |
|
FPAY | Annual fee payment |
Payment date: 20151119 Year of fee payment: 10 |
|
FPAY | Annual fee payment |
Payment date: 20161114 Year of fee payment: 11 |
|
FPAY | Annual fee payment |
Payment date: 20171129 Year of fee payment: 12 |
|
LAPS | Lapse due to unpaid annual fee |