KR101121474B1 - Dc 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법 - Google Patents

Dc 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법 Download PDF

Info

Publication number
KR101121474B1
KR101121474B1 KR1020100121366A KR20100121366A KR101121474B1 KR 101121474 B1 KR101121474 B1 KR 101121474B1 KR 1020100121366 A KR1020100121366 A KR 1020100121366A KR 20100121366 A KR20100121366 A KR 20100121366A KR 101121474 B1 KR101121474 B1 KR 101121474B1
Authority
KR
South Korea
Prior art keywords
inertia
equation
moment
value
friction coefficient
Prior art date
Application number
KR1020100121366A
Other languages
English (en)
Inventor
이상철
Original Assignee
대성전기공업 주식회사
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 대성전기공업 주식회사 filed Critical 대성전기공업 주식회사
Priority to KR1020100121366A priority Critical patent/KR101121474B1/ko
Application granted granted Critical
Publication of KR101121474B1 publication Critical patent/KR101121474B1/ko

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/14Electronic commutators
    • H02P6/16Circuit arrangements for detecting position
    • H02P6/18Circuit arrangements for detecting position without separate position detecting elements
    • H02P6/182Circuit arrangements for detecting position without separate position detecting elements using back-emf in windings

Landscapes

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

Abstract

본 발명은 DC 서보모터 시스템의 관성모멘트(J)와 마찰계수(B)를 구하는 방법에 관한 것으로서, DC 서보모터 시스템에 다양한 전압(Va)을 입력하고 각 입력에 따른 출력으로서 각속도(wm)와 각가속도(
Figure 112010079216370-pat00018
)를 측정하는 단계와; 상기 DC 서보모터 시스템에 적용되는 전기자 회로의 전압방정식, 유기기전력 관계식, 회전력 관계식 및 운동방정식으로부터 관계식
Figure 112010079216370-pat00019
을 유도하는 단계와(여기서, ξ1 = (Ra/ktφf)?J, ξ2 = (Ra/ktφf)?(B + ktkeφf 2/Ra); 이때, Ra는 전기자 권선의 저항, kt는 토크상수, φf는 계자 자속, ke는 역기전력 상수); 상기 측정된 전압(Va)에 따른 각속도(wm)와 각가속도(
Figure 112010079216370-pat00020
)를 상기 관계식에 대입하여 ξ1과 ξ2를 구하되, 완전최소자승법을 이용하여 정확한 ξ1 값과 ξ2 값을 추정하는 단계; 및 상기 추정된 ξ1 값과 ξ2 값으로부터 관성모멘트(J)와 마찰계수(B)를 구하는 단계를 포함하는 것을 특징으로 한다. 이에 의해, 파라미터의 측정을 위한 별도의 장비를 사용하지 않고도 단순히 입력 데이터와 출력 데이터만을 이용하여 파라미터 값인 관성모멘트와 마찰계수를 구하되 종래의 방법에 비해 좀더 향상된 정확도를 얻을 수 있다.

Description

DC 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법{METHOD FOR DETERMINING MOMENT OF INERTIA AND FRICTION COEFFICIENT OF DC SERVO MOTOR SYSTEM}
본 발명은 DC 서보모터 시스템의 제어를 위한 모델링식의 도출에 필요한 파라미터 값으로서 관성모멘트와 마찰계수를 구하는 방법에 관한 것이다.
일반적으로 DC 서보모터를 위치제어 하거나 속도제어를 하기 위해서는 비례적분 (PI)제어기가 광범위하게 이용되고 있다.
이러한 비례적분 제어기를 설계하기 위해서는 우선 DC 서보모터의 수학적인 모델링 과정과 물리적인 모델링 과정이 선행되어야 한다.
모델링된 최종 시스템의 파라미터 값들로는 전기자 권선의 저항, 전기자 권선의 인덕턴스, 토크상수(즉, 역기전력 상수), 전체시스템의 관성모멘트(J) 그리고 마찰계수(B) 등이 포함된다.
따라서, 정확한 모델링식을 도출해내기 위해서는 각각의 파라미터 값을 개별적으로 측정해야 한다.
일반적으로, 인덕턴스와 임피던스를 측정하기 위해서는 주파수에 따른 전기적인 특성들을 측정할 수 있는 LCR 미터를 이용한다.
하지만, 부하시스템을 구동하는 DC 서보모터의 관성모멘트와 마찰계수를 측정하기란 좀처럼 쉽지 않을 뿐 아니라 측정을 위한 관성모멘트 측정기 등의 장비들도 대체로 고가이므로 비용 문제를 고려하지 않을 수 없다.
한편, 상기 관성모멘트 측정기를 이용하지 않을 경우의 종래 관성모멘트 측정 방법은 주로 최소자승법을 이용하는 것이었다.
하지만, 최소자승법에서는 모든 오차가 측정값에만 한정되어 있다는 가정이 필요하다.
다시 말하면, 아래의 식 (a)와 같은 형태의 방정식에서 오차의 squared norm ||e||2이 최소화 되도록 식을 푸는 것이 목적이 된다.
Ax=b-e ... (a)
이것은 에러값이 b 벡터에만 포함된다고 가정하는 것이다.
하지만, 이런 가정은 시스템 모델링에 있어서의 오차, 표준화의 오차 그리고 계측기기의 오차 등에 대해서도 위 식과 같이 가정하기가 어렵다는 문제가 있다.
즉, 실제 이러한 오차들은 행렬 A에도 포함되어 있기 때문이다.
본 발명은 파라미터의 측정을 위한 별도의 장비를 사용하지 않고도 단순히 입력 데이터와 출력 데이터만을 이용하여 파라미터 값인 관성모멘트와 마찰계수를 구하되 종래의 방법에 비해 좀더 향상된 정확도를 확보할 수 있는 DC 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법을 제공하는 데 있다.
상기 목적을 달성하기 위해 본 발명은 DC 서보모터 시스템의 관성모멘트(J)와 마찰계수(B)를 구하는 방법에 있어서, 상기 DC 서보모터 시스템에 다양한 전압(Va)을 입력하고, 각각의 입력에 따른 상기 DC 서보모터 시스템의 출력으로서 각속도(wm)와 각가속도(
Figure 112010079216370-pat00001
)를 측정하는 단계와; 상기 DC 서보모터 시스템에 적용되는 전기자 회로의 전압방정식, 유기기전력 관계식, 회전력 관계식 및 운동방정식으로부터 아래의 관계식 (Ⅰ)을 유도하는 단계와;
Figure 112010079216370-pat00002
... (Ⅰ)
(여기서, ξ1 = (Ra/ktφf)?J, ξ2 = (Ra/ktφf)?(B + ktkeφf 2/Ra); 이때, Ra는 전기자 권선의 저항, kt는 토크상수, φf는 계자 자속, ke는 역기전력 상수)
상기 측정된 전압(Va)에 따른 각속도(wm)와 각가속도(
Figure 112010079216370-pat00003
)를 상기 관계식 (Ⅰ)에 대입하여 상기 ξ1과 ξ2를 구하되, 완전최소자승법을 이용하여 정확한 ξ1 값과 ξ2 값을 추정하는 단계; 및 상기 추정된 ξ1 값과 ξ2 값으로부터 상기 관성모멘트(J)와 마찰계수(B)를 구하는 단계를 포함하는 것을 특징으로 하는 DC 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법을 제공한다.
상기한 본 발명에 따른 DC 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법에 의하면, DC 서보모터의 파라미터의 측정을 위한 별도의 장비를 사용하지 않고도 단순히 입력 데이터와 출력 데이터만을 이용하여 파라미터 값인 관성모멘트와 마찰계수를 구할 수 있음은 물론 종래 사용하던 방법에 비해 더욱 향상된 정확도를 얻을 수 있다.
영구자석 DC 서보모터 시스템은 전기자 회로, 유기기전력, 회전력 그리고 기계시스템을 나타내는 식으로 표현될 수 있다.
1) 전기자 회로의 전압방정식
전기자 회로의 전압방정식은 아래의 식 (1)과 같이 표현될 수 있다.
Va = Raia + La?(dia/dt) + ea ... (1)
여기서, Va는 입력전압[V], ia는 전기자 권선의 전류[A], La는 전기자 권선의 인덕턴스[H], Ra는 전기자 권선의 저항[Ω], ea는 역기전력[V]을 나타낸다.
2) 회전력(Torque)
계자 자속 φf가 형성되어 있는 곳에 전류 ia가 흐르고 있는 전기자 도체에 발휘되는 DC 모터의 토크식은 아래의 식 (2)와 같다.
T = ktf?ia ... (2)
여기서, kt는 사전에 알고 있는 토크상수이다.
3) 유기기전력
자속 φf가 존재하는 공극 내에 wm[rad/s]의 각속도로 회전하는 전기자 도체에 유기되는 기전력은 아래의 식 (3)과 같다.
ea = kef?wm ... (3)
여기서, ke는 사전에 알고 있는 역기전력 상수이다.
MKS 단위계에서 ke 값은 앞서 설명한 토크상수 kt와 같다.
4) 운동방정식
직류 전동기의 기계 기스템에 대한 운동방정식은 각속도 wm을 이용하여 다음 식 (4)와 같이 표현할 수 있다.
T = J?(dwm/dt) + B?wm ... (4)
여기서, J는 전체 시스템의 관성모멘트[km-m2], B는 마찰력 상수, T는 부하토크[Nm]를 나타낸다.
이상의 1) 내지 4)의 관계식을 통합하는 전체시스템을 유도하기 위해, 전기자 회로의 전압방정식인 상기 식 (1)과 유기기전력을 이용한 상기 식 (3)을 합하면 아래의 식 (5)와 같이 나타낼 수 있다.
Ra?ia + La?(dia/dt) = Va - kef?wm ... (5)
또한, 기계 시스템의 운동방정식인 식 (4)와 토크 방정식인 상기 식 (2)를 합하면 아래의 식 (6)과 같이 나타낼 수 있다.
T = J?(dwm/dt) + B?wm = ktf?ia ... (6)
여기에서, La《 Ra 를 가정하여, 위 두 식 (5)와 (6)을 아래와 같이 전체 시스템으로 모델링할 수 있다.
J?(dwm/dt) + (B + ktkeφf 2/Ra)?wm = (ktφf/Ra)?Va ... (7)
(Ra/ktφf)?J?(dwm/dt) + (Ra/ktφf)?(B + ktkeφf 2/Ra)?wm = Va ... (8)
식 (8)의 간략화를 위해서 아래와 같은 변수 표기법을 사용한다.
Figure 112010079216370-pat00004
... (9)
여기서,
Figure 112010079216370-pat00005
이다.
따라서, 상기 관계식 (9)에 완전최소자승법을 적용하면 ξ1과 ξ2의 값을 추정할 수 있다.
알고리즘을 적용하기 위해 다음과 같이 행렬식으로 표현한다.
Figure 112010079216370-pat00006
... (10)
여기서,
Figure 112010079216370-pat00007
,
Figure 112010079216370-pat00008
,
Figure 112010079216370-pat00009
이다.
여기에서, wm은 모터의 속도,
Figure 112010079216370-pat00010
은 가속도 측정 데이터값을 나타내고, Va는 모터의 입력전압 데이터값을 나타낸다.
SVD(singular value decomposition)를 취하기 위해서 상기 행렬식 (10)은 첨가행렬의 형태로 다음의 식 (11)과 같이 표현되어진다.
Figure 112010079216370-pat00011
... (11)
상기 식 (11)에 대하여 SVD를 취한 이후에 그 결과값을 가지고, right singular vector의 완전최소자승법을 통해 최적의 추정된 파라미터 값을 얻을 수 있다.
Figure 112010079216370-pat00012
... (12)
추정된 파라미터 ξ1 및 ξ2를 이용하여 관성모멘트(J)와 마찰력상수(B)의 각 값을 얻을 수 있다.
즉,
Figure 112010079216370-pat00013
,
Figure 112010079216370-pat00014
.
이상과 같은 절차를 거쳐, DC 서보모터 시스템의 입력전압 값과 출력 각속도 및 각가속도 값의 데이터로부터 특이값 분해를 통해 완전최소자승법에 의한 추정된 파라미터 값인 관성모멘트와 마찰력상수를 얻을 수 있다.
한편, 상기에서 설명된 방법은 본 발명의 이해를 돕기 위한 일 실시예에 불과한 것이므로 본 발명의 권리범위 내지 기술적 범위가 상기 설명된 바에 한정되는 것으로 해석되어서는 곤란하다.
본 발명의 권리범위 내지 기술적 범위는 후술하는 특허청구범위 및 그 균등범위에 의해 정하여진다.

Claims (1)

  1. DC 서보모터 시스템의 관성모멘트(J)와 마찰계수(B)를 구하는 방법에 있어서,
    상기 DC 서보모터 시스템에 다양한 전압(Va)을 입력하고, 각각의 입력에 따른 상기 DC 서보모터 시스템의 출력으로서 각속도(wm)와 각가속도(
    Figure 112010079216370-pat00015
    )를 측정하는 단계와;
    상기 DC 서보모터 시스템에 적용되는 전기자 회로의 전압방정식, 유기기전력 관계식, 회전력 관계식 및 운동방정식으로부터 아래의 관계식 (Ⅰ)을 유도하는 단계와;
    Figure 112010079216370-pat00016
    ... (Ⅰ)
    (여기서, ξ1 = (Ra/ktφf)?J, ξ2 = (Ra/ktφf)?(B + ktkeφf 2/Ra); 이때, Ra는 전기자 권선의 저항, kt는 토크상수, φf는 계자 자속, ke는 역기전력 상수)
    상기 측정된 전압(Va)에 따른 각속도(wm)와 각가속도(
    Figure 112010079216370-pat00017
    )를 상기 관계식 (Ⅰ)에 대입하여 상기 ξ1과 ξ2를 구하되, 완전최소자승법을 이용하여 정확한 ξ1 값과 ξ2 값을 추정하는 단계; 및
    상기 추정된 ξ1 값과 ξ2 값으로부터 상기 관성모멘트(J)와 마찰계수(B)를 구하는 단계를 포함하는 것을 특징으로 하는 DC 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법.
KR1020100121366A 2010-12-01 2010-12-01 Dc 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법 KR101121474B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020100121366A KR101121474B1 (ko) 2010-12-01 2010-12-01 Dc 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020100121366A KR101121474B1 (ko) 2010-12-01 2010-12-01 Dc 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법

Publications (1)

Publication Number Publication Date
KR101121474B1 true KR101121474B1 (ko) 2012-03-13

Family

ID=46141622

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020100121366A KR101121474B1 (ko) 2010-12-01 2010-12-01 Dc 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법

Country Status (1)

Country Link
KR (1) KR101121474B1 (ko)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101378824B1 (ko) * 2012-04-20 2014-03-27 라인스트림 테크놀로지스 기계 시스템에서 자동적으로 관성을 추정하기 위한 방법 및 모션 프로파일을 생성하기 위한 방법
US8710777B2 (en) 2012-04-20 2014-04-29 Linestream Technologies Method for automatically estimating inertia in a mechanical system
WO2017023083A1 (ko) * 2015-08-05 2017-02-09 명지대학교 산학협력단 전동기 구동 시스템의 파라미터 추정 장치
US10061275B2 (en) 2014-07-29 2018-08-28 Linestream Technologies Optimized parameterization of active disturbance rejection control
US10126202B2 (en) 2015-09-11 2018-11-13 Linestream Technologies Method for automatically estimating inertia, coulomb friction, and viscous friction in a mechanical system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050054480A1 (en) 2003-09-10 2005-03-10 Ford Global Technologies, Llc Hybrid vehicle powertrain with a multiple-ratio power transmission mechanism
JP2009080092A (ja) 2007-01-19 2009-04-16 Toshihiro Ozasa 軸受摩擦測定方法及び軸受摩擦測定装置
KR20100023846A (ko) * 2007-05-24 2010-03-04 듀르 시스템스 게엠베하 탄성 로봇 구조를 위한 운동 제어기
KR20100076196A (ko) * 2008-12-26 2010-07-06 주식회사 포스코 워크롤의 마찰계수 측정방법

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050054480A1 (en) 2003-09-10 2005-03-10 Ford Global Technologies, Llc Hybrid vehicle powertrain with a multiple-ratio power transmission mechanism
JP2009080092A (ja) 2007-01-19 2009-04-16 Toshihiro Ozasa 軸受摩擦測定方法及び軸受摩擦測定装置
KR20100023846A (ko) * 2007-05-24 2010-03-04 듀르 시스템스 게엠베하 탄성 로봇 구조를 위한 운동 제어기
KR20100076196A (ko) * 2008-12-26 2010-07-06 주식회사 포스코 워크롤의 마찰계수 측정방법

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101378824B1 (ko) * 2012-04-20 2014-03-27 라인스트림 테크놀로지스 기계 시스템에서 자동적으로 관성을 추정하기 위한 방법 및 모션 프로파일을 생성하기 위한 방법
US8710777B2 (en) 2012-04-20 2014-04-29 Linestream Technologies Method for automatically estimating inertia in a mechanical system
US9268316B2 (en) 2012-04-20 2016-02-23 Linestream Technologies Method for automatically estimating a friction coefficient in a mechanical system
US10061275B2 (en) 2014-07-29 2018-08-28 Linestream Technologies Optimized parameterization of active disturbance rejection control
WO2017023083A1 (ko) * 2015-08-05 2017-02-09 명지대학교 산학협력단 전동기 구동 시스템의 파라미터 추정 장치
US10126202B2 (en) 2015-09-11 2018-11-13 Linestream Technologies Method for automatically estimating inertia, coulomb friction, and viscous friction in a mechanical system

Similar Documents

Publication Publication Date Title
JP5952332B2 (ja) 誘導電動機のセンサレスベクトル制御装置
KR101121474B1 (ko) Dc 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법
US20120007528A1 (en) Ac motor control device and ac motor driving system
JP5223109B2 (ja) 永久磁石形同期電動機の制御装置
JP6485810B2 (ja) 誘導モーターのローター角速度を制御するための方法及びシステム
CN103227600A (zh) 电机的无传感器控制设备及其控制方法
KR20160110685A (ko) 영구자석 동기 전동기의 회전자 초기 자극 위치 검출 방법
CN109194224B (zh) 基于扩张状态观测器的永磁同步电机无传感器控制方法
Kuehl et al. Bivariate polynomial approximation of cross-saturated flux curves in synchronous machine models
JP2008286779A (ja) Ipmモータのためのトルク推定器
JP2010035352A (ja) 同期電動機のロータ位置推定装置
KR20220016094A (ko) 동기식 전기 기계의 전자기 토크 추정 방법
CN110620534B (zh) 非线性柔化与变结构滑模控制永磁同步电机转速稳定的方法
Schubert et al. On the torque accuracy of stator flux observer based induction machine control
JP6108114B2 (ja) 永久磁石形同期電動機の制御装置
CN108649849A (zh) 一种简单的无传感器永磁同步电机速度估测方法
KR101590251B1 (ko) 정지형 회전자 시정수 추정 방법
JP2008219982A (ja) 同期モータの位置推定方法およびそのプログラム
KR100881360B1 (ko) 유도전동기의 센서리스 벡터 제어 장치 및 방법
CN106712620B (zh) 永磁同步电机转矩计算方法及装置
Sadoughi et al. Comparison and simulation of rotor flux estimators in stationary and rotating reference frame in the vector control of induction motors for low-speed applications
JP6573213B2 (ja) 永久磁石形同期電動機の制御装置
CN112688616B (zh) 开关磁阻电机相电流估测方法
Quang et al. A sensorless approach for tracking control problem of tubular linear synchronous motor
Zhou et al. High-order terminal sliding-mode observer for speed estimation of induction motors

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: 20150112

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20151217

Year of fee payment: 5

LAPS Lapse due to unpaid annual fee