Microsoft Word - 04 _ __262 박장현

Similar documents

歯522박병호.PDF

장완식.fm

슬라이드 1

정밀국문표지-5월

Microsoft Word - KSR2012A062.doc

Microsoft Word - 스프링백_KAIST_송정한.doc

<B4EBC7D0BCF6C7D02DBBEFB0A2C7D4BCF62E687770>

Microsoft PowerPoint - ch02-1.ppt

cat_data3.PDF

04 Çмú_±â¼ú±â»ç

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE. vol. 29, no. 6, Jun Rate). STAP(Space-Time Adaptive Processing)., -

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE. vol. 29, no. 10, Oct ,,. 0.5 %.., cm mm FR4 (ε r =4.4)

¼º¿øÁø Ãâ·Â-1

High Resolution Disparity Map Generation Using TOF Depth Camera In this paper, we propose a high-resolution disparity map generation method using a lo

미분기하학 II-16 복소평면의선형분수변환과쌍곡평면의등장사상 김영욱 (ÑñÁ) 강의양성덕 (zû ) 의강의록 Ø 'x! xxñ 2007 년 김영욱 (ÑñÁ) 강의양성덕 (zû ) 의강의록 (Ø 'x!) 미분기하 II 2007 년 1 / 26

Python과 함께 배우는 신호 해석 제 5 강. 복소수 연산 및 Python을 이용한 복소수 연산 (제 2 장. 복소수 기초)

À±½Â¿í Ãâ·Â

<4D F736F F D20BECBB1E220BDACBFEE20BAA4C5CD2C20C1C2C7A5B0E82C20C1C2C7A5BAAFC8AFC7E0B7C4>

소성해석

Microsoft Word - 15__ _ 임덕원

2 : 3 (Myeongah Cho et al.: Three-Dimensional Rotation Angle Preprocessing and Weighted Blending for Fast Panoramic Image Method) (Special Paper) 23 2

지능정보연구제 16 권제 1 호 2010 년 3 월 (pp.71~92),.,.,., Support Vector Machines,,., KOSPI200.,. * 지능정보연구제 16 권제 1 호 2010 년 3 월

05 목차(페이지 1,2).hwp

낙랑군

<30312DC1A4BAB8C5EBBDC5C7E0C1A4B9D7C1A4C3A52DC1A4BFB5C3B62E687770>

(Microsoft PowerPoint - Ch6_NumAnalysis.ppt [\310\243\310\257 \270\360\265\345])

Microsoft Word - 12__96-104_ 김진오

Microsoft Word - SDSw doc

PowerPoint 프레젠테이션

<3131BFF92D3828C6D0B3CEBFACB1B82DC0CCBBF3C8A D38302E687770>

DBPIA-NURIMEDIA

슬라이드 1

자기구성지도 기반 방법을 이용한 이상 탐지(Novelty Detection using SOM SOM-based Methods)

Chap 6: Graphs

<35335FBCDBC7D1C1A42DB8E2B8AEBDBAC5CDC0C720C0FCB1E2C0FB20C6AFBCBA20BAD0BCAE2E687770>

문제지 제시문 2 보이지 않는 영역에 대한 정보를 얻기 위하여 관측된 다른 정보를 분석하여 역으로 미 관측 영역 에 대한 정보를 얻을 수 있다. 가령 주어진 영역에 장애물이 있는 경우 한 끝 점에서 출발하여 다른 끝 점에 도달하는 최단 경로의 개수를 분석하여 장애물의

09권오설_ok.hwp

01

untitled

(JBE Vol. 21, No. 1, January 2016) (Regular Paper) 21 1, (JBE Vol. 21, No. 1, January 2016) ISSN 228

歯 PDF

Journal of Educational Innovation Research 2018, Vol. 28, No. 1, pp DOI: * A Analysis of

FGB-P 학번수학과권혁준 2008 년 5 월 19 일 Lemma 1 p 를 C([0, 1]) 에속하는음수가되지않는함수라하자. 이때 y C 2 (0, 1) C([0, 1]) 가미분방정식 y (t) + p(t)y(t) = 0, t (0, 1), y(0)

<353420B1C7B9CCB6F52DC1F5B0ADC7F6BDC7C0BB20C0CCBFEBC7D120BEC6B5BFB1B3C0B0C7C1B7CEB1D7B7A52E687770>


Chapter 1

<BAF9C7D8BFEEC7D7BCB1B9DA20C1F6C4A728B1B9B9AE292E687770>

Microsoft PowerPoint - LA_ch6_1 [호환 모드]

Microsoft Word - 4장_처짐각법.doc

PowerPoint 프레젠테이션

슬라이드 1

<31382D C3D6C0B1B1D42E666D>

벡터(0.6)-----.hwp

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Mar.; 28(3),

<313020BFA9BFB5B1B E666D>

고 학년도 9월고수학 1 전국연합학력평가영역문제지 1 1 제 2 교시 수학영역 5 지선다형 3. 두다항식, 에대하여 는? [ 점 ] 1. 의값은? ( 단, ) [ 점 ] 다항식 이 로인수분해될때, 의값은? ( 단,,

nonpara6.PDF

E-IC-D-1-065(수정).hwp

자기공명 임피던스 단층촬영 기술 연구센터 \(MREIT Research Center\)

08김현휘_ok.hwp

1 1 장. 함수와극한 1.1 함수를표현하는네가지방법 1.2 수학적모형 : 필수함수의목록 1.3 기존함수로부터새로운함수구하기 1.4 접선문제와속도문제 1.5 함수의극한 1.6 극한법칙을이용한극한계산 1.7 극한의엄밀한정의 1.8 연속

Journal of Educational Innovation Research 2016, Vol. 26, No. 3, pp DOI: * The Grounds and Cons

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Nov.; 26(11),

장연립방정식을풀기위한반복법 12.1 선형시스템 : Gauss-Seidel 12.2 비선형시스템 12.1 선형시스템 : Gauss-Seidel (1/10) 반복법은초기근을가정한후에더좋은근의값을추정하는체계적인절차를이용한다. G-S 방법은선형대수방정

2013unihangulchar {45380} 2unihangulchar {54617}unihangulchar {44592} unihangulchar {49328}unihangulchar {50629}unihangulchar {51312}unihangulchar {51

1 경영학을 위한 수학 Final Exam 2015/12/12(토) 13:00-15:00 풀이과정을 모두 명시하시오. 정리를 사용할 경우 명시하시오. 1. (각 6점) 다음 적분을 구하시오 Z 1 4 Z 1 (x + 1) dx (a) 1 (x 1)4 dx 1 Solut

<C7A5C1F620BEE7BDC4>

06_ÀÌÀçÈÆ¿Ü0926

<5B D B3E220C1A634B1C720C1A632C8A320B3EDB9AEC1F628C3D6C1BE292E687770>

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Jun.; 27(6),

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Dec.; 27(12),


REVIEW CHART 1

Microsoft PowerPoint - Ch2-1평형1.ppt [호환 모드]

<BFACBDC0B9AEC1A6C7AEC0CC5F F E687770>

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE. vol. 27, no. 8, Aug [3]. ±90,.,,,, 5,,., 0.01, 0.016, 99 %... 선형간섭

164

,,,.,,,, (, 2013).,.,, (,, 2011). (, 2007;, 2008), (, 2005;,, 2007).,, (,, 2010;, 2010), (2012),,,.. (, 2011:,, 2012). (2007) 26%., (,,, 2011;, 2006;

1_12-53(김동희)_.hwp

230 한국교육학연구 제20권 제3호 I. 서 론 청소년의 언어가 거칠어지고 있다. 개ㅅㄲ, ㅆㅂ놈(년), 미친ㅆㄲ, 닥쳐, 엠창, 뒤져 등과 같은 말은 주위에서 쉽게 들을 수 있다. 말과 글이 점차 된소리나 거센소리로 바뀌고, 외 국어 남용과 사이버 문화의 익명성 등

56(2)-15(PC-054).fm

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Jan.; 26(1),

조사연구 권 호 연구논문 한국노동패널조사자료의분석을위한패널가중치산출및사용방안사례연구 A Case Study on Construction and Use of Longitudinal Weights for Korea Labor Income Panel Survey 2)3) a

Microsoft Word - 03 _ _ 심종엽

인문사회과학기술융합학회

λx.x (λz.λx.x z) (λx.x)(λz.(λx.x)z) (λz.(λx.x) z) Call-by Name. Normal Order. (λz.z)


DBPIA-NURIMEDIA

PowerPoint Presentation

STATICS Page: 7-1 Tel: (02) Fax: (02) Instructor: Nam-Hoi, Park Date: / / Ch.7 트러스 (Truss) * 트러스의분류 트러스 ( 차원 ): 1. 평면트러스 (planar tru


< C7D0B3E2B5B520B4EBBCF6B4C920C7D8BCB3C1F628B1B9BEEE41C7FC20C8A6BCF6292E687770>

PDF

09È«¼®¿µ 5~152s

제 2 교시 2019 학년도 3 월고 1 전국연합학력평가문제지수학영역 1 5 지선다형 1. 의값은? [2점] 일차방정식 의해는? [2 점 ] 두수, 의최대공약수는? [2 점 ] 일차함수 의그래프에서

<4D F736F F F696E74202D2035BBF3C6F2C7FC5FBCF8BCF6B9B0C1FA2E BC8A3C8AF20B8F0B5E55D>

Microsoft Word - KSR2014A433

Chapter 5

<4D F736F F F696E74202D20C1A63134B0AD202D20BBE7BFF8BCF6BFCD20C8B8C0FC>

Transcription:

한국정밀공학회지제 3 권 3 호 pp. 53-61 J. Korean Soc. Precs. Eng., Vol. 3, No. 3, pp. 53-61 ISSN 15-9071(Prnt), ISSN 87-8769(Onlne) March 015 / 53 http://dx.do.org/10.7736/kspe.015.3.3.53 여유자유도를가진 3-SPS/S 병렬 메커니즘의 등각기하대수를 이용한기하학적특이점회피 Geometrc Sngularty Avodance of a 3-SPS/S Parallel Mechansm wth Redundancy usng Conformal Geometrc Algebra 김제석 1, 정진한 1, 박장현 1, Je Seok Km 1, Jn Han Jeong 1, and Jahng Hyon Park 1, 1 한양대학교미래자동차공학과 (Department of Automotve Engneerng, Hanyang Unversty) Correspondng author: jpark@hanyang.ac.kr, Tel: +8--0-0448 Manuscrpt receved: 014.8.7 / Revsed: 015..6 / Accepted: 015..16 A parallel mechansm wth redundancy can be regarded as a means for not only maxmzng the benefts of parallel mechansms but also overcomng ther drawbacks. We proposed a novel parallel mechansm by elmnatng an unnecessary degree of freedom of the confguraton space. Because of redundancy, however, the soluton for the nverse knematcs of the developed parallel mechansm s nfnte. Therefore, we defned a cost functon that can mnmze the movement tme to the target orentaton and found the soluton for the nverse knematcs by usng a numercal method. In addton, we proposed a method for determnng the boundary of the geometrc sngularty n order to avod sngulartes. Key Words: Inverse knematcs ( 역기구학 ), Geometrc sngularty ( 기하학적특이점 ), 3-SPS/S parallel mechansm (3- SPS/S 병렬메커니즘 ), Redundancy ( 여유자유도 ), Conformal geometrc algebra ( 등각기하대수 ) 기호설명 O = the poston of sphercal jont on the 1-S subchan B 0 = the center poston of movng platform O (X o, Y o, Z o ) = base coordnate system fxed on O B (X B, Y B, Z B ) = output coordnate system fxed on B 0 α = the offset angle between axs X o and axs X B β = the angle between axs X o and pont A γ = the angle between axs X B and pont B d A = the dstance between ponts O and A 0 d B = the dstance between ponts O and B 0 1. 서론 대부분의메커니즘은넓은작업공간 (workspace) 과쉬운조작성 (manpulablty) 이라는큰장점을바탕으로직렬형으로설계되었다. 하지만직렬메커니즘은낮은강성으로큰하중에취약하며정밀도가낮은단점을가지고있다. 병렬메커니즘은강성과정밀도, 그리고가반하중을크게향상시킬수있는장점이있다. 하지만제한된작업공간으로인해높은정밀도가요구되는작업에만제한적으로사용되어왔다. 1 Copyrght C The Korean Socety for Precson Engneerng Ths s an Open-Access artcle dstrbuted under the terms of the Creatve Commons Attrbuton Non-Commercal Lcense (http://creatvecommons.org/lcenses/by-nc/3.0) whch permts unrestrcted non-commercal use, dstrbuton, and reproducton n any medum, provded the orgnal work s properly cted.

한국정밀공학회지제 3 권 3 호 pp. 53-61 March 015 / 54 기존의병렬형메커니즘이가지는구조적단점에대해여유자유도를가지는병렬형메커니즘은대안이될수있다. 여유자유도는작업공간의자유도보다기구의자유도가더큰경우로정의되며, 추가된자유도로인하여말단장치를주어진경로로추종하면서부가적인작업을구현할수있게된다. Lee와 Km의연구 에따르면, 여유자유도를가지는병렬형메커니즘은아래와같이세가지로구분한다. 기존시스템의부속체인에추가적인구동기를추가한메커니즘 6 기존시스템의수동관절을능동관절로교체한메커니즘 3,4 기존의시스템에추가적인부속체인을추가한메커니즘 7,8 대표적으로 Wang 6 등은기존의부속체인에하나의구동기를추가하여해석학적특이점에대해회피할수있는방법을고안하였고, Kurtz 7 등과 Zhang 8 등은시스템이필요한자유도보다하나많은부속체인을추가하는방법을연구하였다. 하지만기존의연구들은부속체인에서추가적인자유도를부가하였기때문에, 자체운동 (self-moton) 은불가능하다. 자체운동은말단장치의운동에영향을주지않는관절운동으로정의되며, 특이점회피, 운동에너지최소화등다양한문제를해결할수있게된다. 이러한자체운동특성을이용하여, Park 5 등은 3개의회전자유도중에서하나의자유도를여유자유도로정의한새로운 3-SPS/S 병렬메커니즘을처음제안하였다. 이방법은미사일발사대와같이요구되는자유도가오직두개의회전자유도 (yaw and ptch) 만필요한경우, 목표자세로의이동시간을최소화할수있다. 본논문은기존의연구를발전시켜, 목표자세로의이동시간을최소화하는동시에기하학적특이점을회피할수있는역기구학해법을제안한다. 기하학적특이점은기계적구속조건에의해발생되며, 기계적구속조건은기하학적인방법을통해찾을수있다. 기존의연구 9 는주로 3차원유클리드공간 (Eucldean space) 상에서구속조건을찾고자하였으나, 복잡한해법으로인해작업공간을해석하기위한용도로제한적으로사용되었다. 이에반해본논문에서는 3차원이아닌 5차원의등각공간 (conformal space) 상에서기계적구속조건 Table 1 Lst of conformal geometrc enttes Entty IPNS OPNS Pont P = x+ 0.5x e + e 0 Sphere Plane S = P 0.5r e 의경계점을등각기하대수 (conformal geometrc algebra) 로손쉽게찾을수있는방법을제안한다. 이기계적구속조건의경계점을이용하면기하학적특이점이발생하지않은영역을찾을수있고, 이영역을좋은자세가능영역으로정의하여역기구학을위한해영역 (soluton regon) 으로제한한다. 그렇게함으로써기하학적특이점으로인한작업공간이제한되는것을보완한다. 장에서는문제를설정하고, 3 장에서는여유자유도역기구학문제의해법을기술하고, 4 장에서는제안하는특이점회피알고리즘을시뮬레이션을통하여효용성을보이도록하겠다.. 여유자유도를가진 3-SPS/S 병렬메커니즘.1 등각기하대수병렬메커니즘의복잡한형상으로인해유클리드공간상에서해석하는것은어렵다. 그러나등각기하대수를이용하면보다편하게해석할수있다. 등각기하대수는 Quaternon, Projectve Geometry, Le Algebra 등다양한수학적인이론을하나로통합하여이해하기쉽게만든수학적언어로서로봇공학, 10-1 컴퓨터그래픽스 13,14 등다양한분야에서새로운방법으로주목을받고있다. 등각기하대수는 3차원유클리드기하대수에서두개의기저벡터를추가하여 5차원으로확장시킨이론으로 5개의기저벡터들 ( e, e, e, e, ) 1 3 e 에의해 0 정의된다. 추가된기저벡터 e 는 3차원공간의원 0 점을표현하고, e 는무한한곳에서의점을표현한다. 등각기하대수의주요특징은 Table 1과같이직선, 평면, 구, 원등의다양한기하학적요소들을직관적으로표현할수있고, 이러한기하학적인요소들은내적, 외적, 그리고기하학적곱에의해하나의대수형식으로표현이가능한장점이 S P P P P = 1 3 4 π= n + de π = P P P e 1 3 Crcle Z = S S 1 Z = P P P 1 3 Lne = 1 L π π L P P = e 1 Pont-par PP = S S S 1 3 PP = P P 1

한국정밀공학회지제 3 권 3 호 pp. 53-61 March 015 / 55 있다. 등각 기하대수와 관련된 상세한 내용은 Hldenbrand의논문 15,16 을참조하기바란다.. 설계변수여유자유도를가진 3-SPS/S 병렬메커니즘은지면에고정된기저플랫폼 (base platform) 과직교공간상에서회전만이가능한이동플랫폼 (movng platform) 으로구성된다. 그리고이동플랫폼은 3- SPS 부속체인의직동관절 (prsmatc jont) 에의해구동되며, 1-S 부속체인에의해움직임이제한된다. Fg. 1 에서보듯이기저플랫폼에놓인세개의구형관절은반경이 r A 인원상에각 β 의간격으로놓여있으므로, 기저플랫폼상의구형관절의위치 A (=1,,3) 는다음과같다. a = r c e + r s e d e A 1 = a + a e e + 0 β 1 β 3 A A A 이동플랫폼에놓인세개의구형관절은반경이 r B 인원상에각 간격으로놓여있고, X o 축에대해각 α 만큼회전된위치에존재하므로, 이동플랫폼상의구형관절의위치 B ( =1,,3) 는다음과같다. b = rc e + rs e + d e ( α + γ ) ( α + ) γ 1 3 B B B = 1 + e + e 0 B b b l 는 번째 SPS 부속체인에서직동관절의관절변수로서점 A 에서부터점 B 까지의거리를나타내고, 등각기하대수에서두점사이의거리는다음과같이표현된다. (1) () l = ( A B ) (3).3 기구학적관계식기존의 3-SPS/S 병렬메커니즘에서요구되는출력공간의자유도가오직두개의회전자유도 (yaw and ptch) 만필요한경우, ZYZ 오일러각 (ϕ, θ, ψ) 중에서각 ψ 를여유자유도로정의할수있다. 그러므로여유자유도를가지는 3-SPS/S 병렬메커니즘의출력공간은일반적인출력공간의 6 개변수중 개만으로구성된다 : u =[ϕ, θ] T. base platform (a) 3D vew movng platform (b) Top vew Fg. 1 3-SPS/S parallel mechansm wth ts desgn parameters. S and P denote a sphercal jont, and a prsmatc jont, respectvely 기구학적여유자유도로인해 Z B 축을기준으로각 ψ의변화에따라자체운동이가능하므로점 B 는항상원 C B 상에존재한다. C = B B B (4) B 1 3 이동플랫폼의자세가주어질때, 출력좌표계의자세 R t 는각 ψ 에관한함수로등각기하대수의 Rotor 에의해표현된다. R ( ψ) = R R R ( ψ) (5) φ θ ψ t z y z,,, 여기서 R z, ϕ 는 Z B 축을기준으로각 ϕ 만큼의회전을, R y,θ 는 Z B 축을기준으로각 θ 만큼의회전을, R z,ψ 는 Z B 축을기준으로각 ψ 만큼의회전을표현하는 Rotor 이다. 그리고이동플랫폼의 번째구형관절까지의위치 B 는기하학적곱에의해표현된다. B ( ψ) = R ( ψ) B R ( ψ) (6) t t 식 (6) 을식 (3) 에대입하여관절변수 l 를각 ψ 에관한함수로정의하면다음과같다. l ( ψ) = d + d + r + r + d d A B A B A cos( θ) B + drsn( β B A φ θ) + drsn A B + + 0.5rr cos A B + + + rr cos A B + + / + rr cos + + + ( α γ ψ θ) ( α β γ ψ φ θ) ( α β γ ψ φ θ) ( α β γ ψ φ) A B (7)

한국정밀공학회지제 3 권 3 호 pp. 53-61 March 015 / 56 ( α β γ ψ φ) ( β φ θ) ( α γ ψ θ) ( α β γ ψ φ θ) ( α β γ ψ φ θ) rr cos + + + A B drsn B A drsn + + + A B rr cos + + + + / A B rr cos + + + + / A B 병렬메커니즘의각링크는관절좌표계로표현되기때문에메커니즘을제어하기위해서는로봇말단의작업공간을관절공간으로표현해야한다. m u R : Output varable (8) n l R : Jont varable ϕ, θ 는출력변수로서역기구학해석에앞서주어지고, 관절변수 l 는여유자유도 ψ 에의해결정된다. 그러므로각관절변수에서출력변수로사상시키는 f 를정의하면다음과같이기구학적관계식을정리할수있다. 3. 역기구학해석 n m f : R R (9) u = f l ( ) 3.1 비용함수의정의여유자유도를가지는로봇은다양한형태로운동이가능하다는장점으로인해작업자가원하는작업목적에부합하도록역기구학을해석할수있다. 본논문에서는무한의해가존재하는여유자유도문제에서최적해를찾기위해, 목표자세로의이동시간을최소화하도록목적함수를선정하였다. 이를통해높은강성등병렬형구조의장점을유지하면서직렬형구조에비해빠른이동성등의이점을가질수있기때문이다. 목표자세로의이동시간최소화는곧각체인의직동관절의변위량을최소화와마찬가지이므로비용함수 f cost 를다음과같이정의하였다. { ψ } f ( ψ ) = max Δ l ( ) cost (10) π ψ < π 여기서, Δ l ( ψ) = l ( ψ) l 이고, l,old,old 는현재의 번째링크길이이다. 만약임의의함수 f (x) 와 g(x) 가연속이면 Fg. Graph of the change n the cost functon f( x) + g( x) f( x) g( x) max ( f( x), g( x) ) = + (11) 이므로, 식 (10) 는다음과같은연속함수로표현된다. 17 ( ψ ) = max ( Δ, Δ, Δ ) f l l l cost π ψ < π 1 3 Δ l +Δ l +Δ l + Δl Δl = 1 3 3 4 Δl Δ l +Δ l + Δl Δl + 1 3 3 4 (1) 따라서링크의최대변위량으로표현되는비용함수 f cost 를최소화함으로써목표자세로의이동시간을감소시킬수있다. 3. 역기구학해 Fg. 에서보듯이비용함수 f cost 는비선형방정식이므로제안하는메커니즘의역기구학은비선형최적화문제이다. 일반적으로수치해법을이용하여근사적인최적해를찾을수있다. 본논문에서는황금분할탐색법과 차보간법을결합하여사용하였다. 황금분할탐색법은연산량은많지만신뢰도가높고, 차보간법은신뢰도는낮으나연산량이적은특징을가진다. 이두방법을결합함으로인해장점은부각시키고, 단점을상쇄시킬수있다. 18 4. 기하학적특이점회피 여유자유도로봇이가지는장점은일반적인역기구학과달리주어진작업에대해무수히많은관절공간의해를가지고있으므로기하학적특이

한국정밀공학회지제 3 권 3 호 pp. 53-61 March 015 / 57 점에대해서도회피가가능하다. 기하학적특이점은알고리즘적특이점과다르게다음과같은기계적구속조건에의해발생된다. 19,0 (1) 링크길이의제한 () 구형관절의이동범위제한 (3) 링크간의간섭기하학적특이점이발생하게되면, 기구에물리적인손상이가해지는등의피해가발생된다. 그렇기때문에이기구학적특이현상을회피하기위한방법이필요하다. 기하학적특이점을정확히알고있으면, 비용함수의해탐색범위를감소시켜기하학적특이점을회피하는방법을사용할수있다. 하지만이기하학적특이점에대한경계점을정확히알기는어렵다. 그러나기하학적인방법을통해근사적으로기하학적특이점의경계점을알수있다. 다만, 링크길이의변위량을최소화하기때문에, 기하학적특이점중에서 (3) 보다는 (1) 과 () 가발생할가능성이높다. 이번장에서는구속조건중 (1) 과 () 의경계점을기하학적으로찾을수있는방법을제안한다. 4.1 링크길이의제한직동구동기 (prsmatc actuator) 는기계적으로최대길이가제한될수밖에없다. 그러므로 3-SPS 부속체인의링크길이 l ( =1,,3) 는식과같은구속조건을만족해야한다. l l ( ψ ) l for = 1,,3 (13),mn,max 여기서 l,max 과 l,mn 는링크 의최대, 최소링크길이를의미한다. 이러한구속조건을회피하려면, 구속조건의경계점을정확히아는것이중요하다. 본구속조건의경계점은 Fg. 3(a) 과같이최대, 최소링크길이를의미하는구와이동플랫폼이만나는두점이다. 이는등각기하대수를이용하여기하학적인방법으로찾을수있다. 우선 S,max 는반지름이 l,max 이고, 중심이점 A 인구이고, S,mn 은반지름이 l,mn 이고, 중심이점 A 인구라고하자. e S = A 0.5 l,max max (14) Fg. 3 Geometrc sngularty usng geometrc approach 가교차할때생성되는쌍점 (pont-par) PP,max 혹은 PP,mn 은제한된링크길이의경계점을의미한다. PP = S C,max,max B (16) PP = S C,mn,mn B (17) 쌍점 PP 는다음과같이기하학적특이점의양경계점으로분해할수있다. 1 B j PP ± PP PP = e PP (18) 두기하학적요소가교차하는경우기하학적특이점이발생되는영역을위와같이찾을수있다. 하지만교차하지않는경우에는특이점이발생되지않는다. 따라서교차하는지에대해판별할필요성이있다. 두기하학적요소, 원 C B 와구 S 가교차하는지에대한판단은 B=(C B S ) 의평방놈 (squared norm) B = B 에의해아래와같이판별할수있다. 14 B >0 이면, C B 와 S 는두점에서교차한다. B =0 이면, C B 와 S 는한점에서접한다. B <0 이면, C B 와 S 는교차하지않는다. 만약 PP,max 0 이고 PP,mn 0 이라면, 여유자유도 ψ 가 [-π, π] 의범위내에서기하학적특이점이발생되지않음을의미한다. 반대의경우, Fg. 3(a) 와같이쌍점 PP 사이의영역은기하학적특이점이발생됨을의미하고, 쌍점내부에둘러싸인영역은기하학적특이점이발생되는영역이다. S = A 0.5 (15),mn mn l e 구 S,max 혹은 S,mn 과이동플랫폼상의원 C B 4. 구형관절의이동범위제한이동플랫폼상의점 B 와기저플랫폼상의점 A 는구형관절로이루어져있다. 이구형관절은

한국정밀공학회지제 3 권 3 호 pp. 53-61 March 015 / 58 기계적인구속조건에의해링크 와이동플랫폼사이에일정각도이상은움직일수없는영역이형성되므로다음과같은구속조건을만족해야한다. δ < δ for 1,,3 mn = (19) 여기서, δ mn 은구형관절의최소작동범위이고, δ 는링크 와이동플랫폼사이각이다. 마찬가지로구속조건을회피하려면, 구속조건의경계점을정확히아는것이중요하다. 본구속조건의경계점은 Fg. 3(b) 와같이이동플랫폼상의구형관절의최소작동범위를의미하는평면과이동플랫폼상의원이만나는두개의점이다. 이는등각기하대수를이용하여기하학적인방법으로찾을수있다. 이동플랫폼상의평면을 π B 라하고, 벡터 OB 0 와벡터 OA 가존재하는평면에수직한벡터를 n 라하자. 여유자유도 ψ 의모든구간에서기하학적특이점이발생되지않음을의미한다. 이외의경우, Fg. 3(b) 와같이쌍점 PP 사이의영역은식 (19) 의구속조건을만족하지않음을의미한다. 4.3 기하학적특이점회피우리는기하학적인방법을통해근사적으로기하학적특이점의경계점을찾았다. 이경계점들은비용함수의해탐색범위를감소시켜기하학적특이점을회피할수있다. 하지만해의탐색범위는점이아니라각 ψ 의값이므로, 점 B j 를각 ψ 로변환하여야한다. 각 ψ =0 일때직선을 L,0 라하고, L,j 를점 B 0 와기하학적특이점의양경계점 B,j 를지나는직선이라하자. L = B (0) B (0) (6) e,0 0 π = C ^ = B B B ^ (0) B B e e 0 1 3 n = OB OA (1) 이동플랫폼상의구형관절의최소작동범위를의미하는평면 π,mn 은이동플랫폼상의평면 π B 를벡터 n 를회전축으로각 δ mn 만큼회전시킨다음, 점 A 로평행이동시킨평면으로다음과같이정의된다. ( δ ) R = exp 0.5 n () B mn T = 1 0.5 A (3) B e π = R T π T R (4),mn B B B B B L = B (0) B (7) e, j 0, j 각 ψ j 는직선 L,0 와직선 L,j 가이루는각으로 OPNS 의내적에의해다음과같이계산된다. [9] ψ j = cos 1 L L,0,0 Lj L j (8) 우리는기하학적특이점이발생되지않는해의영역을알고있으므로비용함수 f cost 의범위를제한함으로써기하학적특이점을회피하여작업공간의확장이가능할뿐만아니라수치해법의단점인지역극소점 (local mnma) 에빠지는확률을감소시킬수있다. 그러므로평면 π,mn 과원 C B 가교차할때생성되는쌍점 PP,mn 은제한된구형관절의이동범위의경계점을나타낸다. PP = π C (5),mn,mn B 하지만원 C B 와평면 π,mn 는항상교차하지않는다. 여유자유도 ψ 가 [-π, π] 의범위내에서항상식 (19) 의구속조건을만족할수도, 만족하지않을수도있기때문이다. 그러므로원 C B 와구 S,j 에의해생성된쌍점 PP 의평방놈 PP 0 이라면 5. 시뮬레이션 병렬형기구는복잡한형상으로인해작업공간내에기하학적특이점이많이존재한다. 이러한이유로본논문에서는기하학적특이점을회피할수있는여유자유도병렬형구조의역기구학해석법을제안하였고, 이에대한시뮬레이션을통해제안하는메커니즘의효용성을검증하고자한다. 여유자유도를가지는 3-SPS/S 병렬메커니즘의시뮬레이션을위하여설계변수를 Table 와같이선정하였다.

한국정밀공학회지제 3 권 3 호 pp. 53-61 March 015 / 59 Table Desgn parameters r A r B d A d B α β,γ 0.35m 0.0m 0.33m 0.0m 30 10(-1) (a) (b) Fg. 4 Comparson of workspace: (a) non-sngularty avodance (b) geometrc sngularty avodance 업공간중 18.% 영역에서만특이점이발생한다. 이처럼제안하는특이점회피알고리즘을통해작업공간이확장되었지만본논문에서목표로삼은특이점회피와더불어목표자세로의이동시간최소화라는목적도동시에달성했는지확인이필요하다. 이를확인하기위해링크의최대변위량으로표현되는비용 f cost 을특이점을적용하지않은경우와적용한경우에대해서로비교하였다. Fg. 5 는 θ, ϕ 의목표자세가주어질경우해당자세에대한 f cost 값을표시한것으로, 제안하는특이점회피알고리즘을통해 Fg. 4 에서보듯이작업공간내의기하학적특이점을기존대비 41.1% 가량감소시켰음에도불구하고, Fg. 5 에서보듯이링크의최대변위량은크게증가하지않는것을보여주고있다. 6. 결론 (a) (b) Fg. 5 Comparson of optmal soluton: (a) nonsngularty avodance, (b) geometrc sngularty avodance 제안한기하학적특이점회피알고리즘의효용성을검증하기위해제안한메커니즘에기하학적특이점회피알고리즘을적용하지않을경우와적용한경우를비교하였다. 비교를위해출력변수의범위를 ϕ =[-π/, π/], θ =[-π/, π/] 로제한하고, 1 간격마다역기구학해석을수행하였다. 제안하는특이점회피알고리즘의효용성을알아보기위해 Fg. 4 와같이작업이가능한공간을붉은점으로표시하였다. Fg. 4(a) 는특이점회피알고리즘을적용하지않고, 직동관절의변위량을최소화만하였을경우작업공간을보여주고, Fg. 4(b) 는특이점회피와구동관절의변위량을모두최소화하였을경우의작업공간을보여준다. Fg. 4 에서보듯이, 기하학적특이점회피를적용하지않는경우전체작업공간중에 44.% 영역에서특이점이발생하지만, 본논문에서제안하는특이점회피알고리즘을사용하였을경우전체작 본논문에서는기존의여유자유도를가지는병렬형메커니즘과다르게구동기를추가하지않고출력공간의자유도를감소시켜기구학적으로자체운동이가능한병렬메커니즘을고안하였다. 그리고여유자유도를특성을최대한활용하여기존의직병렬형구조의단점을극복할수있는방향으로연구를수행하였다. 본논문에서는기구학적여유자유도로인해역기구학은무한한해를가지므로, 목표자세로의이동시간을최소화하도록비용함수를정의하고, 수치해법을통해역기구학의최적해를찾을수있는방법을제안하였다. 하지만병렬메커니즘은복잡한형상으로인해작업공간내에기하학적특이점많이존재하므로, 이를등각기하대수를이용하여기하학적특이점의경계점을찾고이를회피하기위한연구를수행하였다. 시뮬레이션을통해본논문에서제안하는여유자유도메커니즘의효용성을확인하였다. 후기 이논문은국방과학연구소생존성기술특화연구센터의사업으로지원받아연구되었음 ( 계약번호 UD090090GD). REFERENCES 1. Gosseln, C. and Angeles, J., The Optmum

한국정밀공학회지제 3 권 3 호 pp. 53-61 March 015 / 60 Knematc Desgn of a Planar Three-degree-of- Freedom Parallel Manpulator, Journal of Mechancal Desgn, Vol. 110, No. 1, pp. 35-41, 1988.. Lee, S. H. and Km, S. B., Knematc Feature Analyss of Parallel Manpulator Systems, Proc. of Internatonal Conference on Intellgent Robots and Systems, Vol., pp. 141-148, 1994. 3. Bonev, I. A. and Ryu, J., A New Method for Solvng the Drect Knematcs of General 6-6 Stewart Platforms usng Three Lnear Extra Sensors, Mechansm and Machne Theory, Vol. 35, No. 3, pp. 43-436, 000. 4. Lu, Y., Hu, B., and Lu, P.-L., Knematcs and Dynamcs Analyses of a Parallel Manpulator wth Three Actve Legs and One Passve Leg by a Vrtual Seral Mechansm, Multbody Syst. Dyn., Vol. 17, No. 4, pp. 9-41, 007. 5. Park, J. H., Km, J. S., and Jeong, J. H., Knematcs Optmzaton of a 3-SPS Parallel Redundant Moton Mechansm usng Conformal Geometrc Algebra, Int. Conf. on Mechancal Manufacturng and Control, 014. 6. Wang, J. and Gosseln, C., Knematc Analyss and Desgn of Knematcally Redundant Parallel Mechansms, Journal of Mechancal Desgn, Vol. 16, No. 1, pp. 109-118, 004. 7. Kurtz, R. and Hayward, V., Multple-goal Knematc Optmzaton of a Parallel Sphercal Mechansm wth Actuator Redundancy, IEEE Transactons on Robotcs and Automaton, Vol. 8, No. 5, pp. 644-651, 199. 8. Zhang, D. and Le, J., Knematc Analyss of a Novel 3-dof Actuaton Redundant Parallel Manpulator usng Artfcal Intellgence Approach, Robotcs and Computer-Integrated Manufacturng, Vol. 7, No. 1, pp. 157-163, 011. 9. Merlet, J. P., Determnaton of the Orentaton Workspace of Parallel Manpulators Journal of ntellgent and robotcs systems, Vol. 13, No., pp. 143-160, 1995. 10. Zamora, J. and Bayro-Corrochano, E., Inverse Knematcs, Fxaton and Graspng usng Conformal Geometrc Algebra, Proc. of IEEE/RSJ Int. Conf. on Inlelllgent Robots and Systems, Vol. 4, pp. 3841-3846, 004. 11. Arstdou, A., and Lasenby, J., Fabrk: A Fast, Iteratve Solver for the Inverse Knematcs Problem, Graphcal Models, Vol. 73, No. 5, pp. 43-60, 011. 1. Km, J. S., J, Y. K., and Park, J. H., Inverse Knematcs Analyss of 7-DOF Anthropomorphc Robot Arm usng Conformal Geometrc Algebra, J. Korean Soc. Precs. Eng., Vol. 9, No. 10, pp. 1119-117, 01. 13. Wareham, R., Cameron, J., and Lasenby, J., Applcatons of Conformal Geometrc Algebra n Computer Vson and Graphcs, n Computer Algebra and Geometrc Algebra wth Applcatons, L, H., Olver, P., and Sommer, G.(Eds.), Sprnger Berln Hedelberg, pp. 39-349, 005. 14. Roa, E., Theoktsto, V., Farén, M., and Navazo, I., GPU Collson Detecton n Conformal Geometrc Space, V Ibero-Amercan Symposum n Computer Graphcs SIACG, pp. 153-157, 011. 15. Hldenbrand, D., Fontjne, D., Perwass, C., and Dorst, L., Geometrc Algebra and ts Applcaton to Computer Graphcs, Tutoral notes of the EUROGRAPHICS conference, 004. 16. Hldenbrand, D., Foundatons of Geometrc Algebra Computng, Sprnger, Berln Hedelberg, pp. 17-4, 013. 17. Ross, K.A., Elementary Analyss, Sprnger-Verlag, pp. 19-130, 1980. 18. Press, W. H., Teukolsky, S. A., Vetterlng, W. T., and Flannery, B. P., Numercal Recpes n C++, Cambrdge UP, pp. 398-409, 00. 19. Merlet, J. P., Determnaton of 6D Workspaces of Gough-type Parallel Manpulator and Comparson between Dfferent Geometres, The Internatonal Journal of Robotcs Research, Vol. 18, No. 9, pp. 90-916, 1999. 0. Tsa, K. and Ln, J., Determnng the Compatble Orentaton Workspace of Stewart Gough Parallel Manpulators, Mechansm and machne theory, Vol. 41, No. 10, pp. 1168-1184, 006. 1. Ptt, J., Hldenbrand, D., Stelzer, M., and Koch, A., Inverse Knematcs of a Humanod Robot based on Conformal Geometrc Algebra usng Optmzed Code Generaton, Proc. of the 8th IEEE-RAS Internatonal Conference on Humanod Robot, pp. 681-686, 008.. Chaparro-Altamrano, D., Zavala-Yoe, R., and

한국정밀공학회지제 3 권 3 호 pp. 53-61 March 015 / 61 Ramrez-Mendoza, R., Knematc and Workspace Analyss of a Parallel Robot Used n Securty Applcatons, n Mechatroncs, Electroncs and Automotve Engneerng, pp. 3-8, 013.

한국정밀공학회지제 3 권 3 호 pp. 53-61 March 015 / 6