ポテンシャル法によるロボット製品の障害物回避技術の開発,三菱重工技報 Vol.51 No.1(2014)

Similar documents
Microsoft PowerPoint - 6.PID制御.pptx

相対性理論入門 1 Lorentz 変換 光がどのような座標系に対しても同一の速さ c で進むことから導かれる座標の一次変換である. (x, y, z, t ) の座標系が (x, y, z, t) の座標系に対して x 軸方向に w の速度で進んでいる場合, 座標系が一次変換で関係づけられるとする

数値計算で学ぶ物理学 4 放物運動と惑星運動 地上のように下向きに重力がはたらいているような場においては 物体を投げると放物運動をする 一方 中心星のまわりの重力場中では 惑星は 円 だ円 放物線または双曲線を描きながら運動する ここでは 放物運動と惑星運動を 運動方程式を導出したうえで 数値シミュ

パナソニック技報

Microsoft PowerPoint - chap8.ppt

Microsoft PowerPoint - 【最終提出版】 MATLAB_EXPO2014講演資料_ルネサス菅原.pptx

untitled

次世代過給システムを支える高速モデル予測制御技術,三菱重工技報 Vol.54 No.1(2017)

平成20年2月10日号

航空機の運動方程式

目次 1: 安全性とソフトウェア 2: 宇宙機ソフトウェアにおける 安全 とは 3:CBCS 安全要求とは 4: 宇宙機ソフトウェアの実装例 5: 安全設計から得た新たな知見 6: 今後 2

パソコンシミュレータの現状

Microsoft PowerPoint - 三次元座標測定 ppt

Microsoft Word - NumericalComputation.docx

77

Microsoft PowerPoint - 知財報告会H20kobayakawa.ppt [互換モード]

2009 年 11 月 16 日版 ( 久家 ) 遠地 P 波の変位波形の作成 遠地 P 波の変位波形 ( 変位の時間関数 ) は 波線理論をもとに P U () t = S()* t E()* t P() t で近似的に計算できる * は畳み込み積分 (convolution) を表す ( 付録

クラス図とシーケンス図の整合性確保 マニュアル

線積分.indd

0 21 カラー反射率 slope aspect 図 2.9: 復元結果例 2.4 画像生成技術としての計算フォトグラフィ 3 次元情報を復元することにより, 画像生成 ( レンダリング ) に応用することが可能である. 近年, コンピュータにより, カメラで直接得られない画像を生成する技術分野が生

航空機の縦系モデルに対する、非線形制御の適用例

OCW-iダランベールの原理

コンピュータグラフィックス第6回

例 e 指数関数的に減衰する信号を h( a < + a a すると, それらのラプラス変換は, H ( ) { e } e インパルス応答が h( a < ( ただし a >, U( ) { } となるシステムにステップ信号 ( y( のラプラス変換 Y () は, Y ( ) H ( ) X (

人間の視野と同等の広視野画像を取得・提示する簡易な装置

構造力学Ⅰ第12回

Microsoft PowerPoint - 電装研_2波長赤外線センサを用いた2波長融合処理について

<4D F736F F D208D5C91A297CD8A7793FC96E591E631308FCD2E646F63>

efit+ for TCU 測量計算

3.3 導入ポテンシャルの簡易シミュレーション シミュレーションの流れ導入ポテンシャルの簡易シミュレーションは 河川上の任意の水路 100m セグメント地点で取水し 導水管を任意の箇所に設定して 任意の放水地点に発電施設を設置して発電する場合の導入ポテンシャル値 ( 設備容量 ) 及び概

画像類似度測定の初歩的な手法の検証

Chap2.key

.( 斜面上の放物運動 ) 目的 : 放物運動の方向の分け方は, 鉛直と水平だけではない 図のように, 水平面から角 だけ傾いた固定した滑らかな斜面 と, 質量 の小球を用意する 原点 から斜面に垂直な向きに, 速さ V で小球を投げ上げた 重力の加速度を g として, 次の問い に答えよ () 小

News Release 国立研究開発法人新エネルギー 産業技術総合開発機構 福島県 南相馬市 株式会社 SUBARU 日本無線株式会社 日本アビオニクス株式会社 三菱電機株式会社 株式会社自律制御システム研究所 世界初 無人航空機に搭載した衝突回避システムの探知性能試験を実施

0 部分的最小二乗回帰 Partial Least Squares Regression PLS 明治大学理 学部応用化学科 データ化学 学研究室 弘昌

ベクトル公式.rtf

コンテンツセントリックネットワーク技術を用いた ストリームデータ配信システムの設計と実装

SIP自動走行システム推進委員会

ビッグデータ分析を高速化する 分散処理技術を開発 日本電気株式会社

風力発電インデックスの算出方法について 1. 風力発電インデックスについて風力発電インデックスは 気象庁 GPV(RSM) 1 局地気象モデル 2 (ANEMOS:LAWEPS-1 次領域モデル ) マスコンモデル 3 により 1km メッシュの地上高 70m における 24 時間の毎時風速を予測し

リードタイムが変動する在庫管理モデルの安定性解析 西平 w(k)= u(k-l(k)) ( 2 ) となる このモデルに対して, メモリーレスフィードバック u(k)= Kx(k) (3) を施すことを考える また, 本稿では内部安定性を考えるため, 外生信号 d(k)= 0とすると, システム (

Kumamoto University Center for Multimedia and Information Technologies Lab. 熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI 宮崎県美郷

スライド 1

2/17 目次 I. はじめに... 3 II. 操作手順 (Controlの場合) 断面の作成 寸法測定 異なる断面間の寸法測定 繰り返し処理...11 III. 操作手順 (Verifyの場合) 断面の作成... 1

Microsoft Word - thesis.doc

2018年度 岡山大・理系数学

Functional Programming

Presentation Title Arial 28pt Bold Agilent Blue

図 5 一次微分 図 6 コントラスト変化に伴う微分プロファイルの変化 価し, 合否判定を行う. 3. エッジ検出の原理ここでは, 一般的なエッジ検出の処理内容と, それぞれの処理におけるパラメータについて述べる. 3.1 濃度投影検出線と直交する方向に各画素をスキャンし, その濃度平均値を検出線上

8. 自由曲線と曲面の概要 陽関数 陰関数 f x f x x y y y f f x y z g x y z パラメータ表現された 次元曲線 パラメータ表現は xyx 毎のパラメータによる陽関数表現 形状普遍性 座標独立性 曲線上の点を直接に計算可能 多価の曲線も表現可能 gx 低次の多項式は 計

Mobile robotics KUKA Navigation Solution JP


1/17 平成 29 年 3 月 25 日 ( 土 ) 午前 11 時 1 分量子力学とクライン ゴルドン方程式 ( 学部 3 年次秋学期向 ) 量子力学とクライン ゴルドン方程式 素粒子の満たす場 y ( x,t) の運動方程式 : クライン ゴルドン方程式 : æ 3 ö ç å è m= 0

組込みシステムにおける UMLモデルカタログの実践研究

Microsoft Word - NDIS1204意見受付用_表紙.docx

Microsoft PowerPoint - H22制御工学I-2回.ppt

<4D F736F F D20342E899E D2091E52D81848FAC82D682CC88F8897A2E646F6378>

概論 : 人工の爆発と自然地震の違い ~ 波形の違いを調べる前に ~ 人為起源の爆発が起こり得ない場所がある 震源決定の結果から 人為起源の爆発ではない事象が ある程度ふるい分けられる 1 深い場所 ( 深さ約 2km 以上での爆発は困難 ) 2 海底下 ( 海底下での爆発は技術的に困難 ) 海中や

第20回AIチャレンジ

<4D F736F F D2091E E8FDB C588ECE926E816A2E646F63>

高分解能衛星データによる地形図作成手法に関する調査研究 ( 第 2 年次 ) 実施期間平成 18 年度 ~ 測図部測図技術開発室水田良幸小井土今朝巳田中宏明 佐藤壮紀大野裕幸 1. はじめに国土地理院では, 平成 18 年 1 月に打ち上げられた陸域観測技術衛星 ALOS に関して, 宇宙航空研究開

vecrot

周期時系列の統計解析 (3) 移動平均とフーリエ変換 nino 2017 年 12 月 18 日 移動平均は, 周期時系列における特定の周期成分の消去や不規則変動 ( ノイズ ) の低減に汎用されている統計手法である. ここでは, 周期時系列をコサイン関数で近似し, その移動平均により周期成分の振幅

スライド 1

<4D F736F F F696E74202D208CB48E7197CD8A7789EF F4882CC91E589EF8AE989E A2E B8CDD8AB B83685D>

モデリングとは

アプリケーション インスペクションの特別なアクション(インスペクション ポリシー マップ)

CANapeを用いたラピッドコントロールプロトタイピングのバイパス手法による制御モデル開発

スライド 1

Probit , Mixed logit

技術資料 JARI Research Journal OpenFOAM を用いた沿道大気質モデルの開発 Development of a Roadside Air Quality Model with OpenFOAM 木村真 *1 Shin KIMURA 伊藤晃佳 *2 Akiy

<4D F736F F F696E74202D B F08A84816A955D89BF B838B91808DEC8EE88F878F F312E B8CDD8AB B83685D>

DVIOUT-SS_Ma

SICE東北支部研究集会資料(2011年)

卒業研究発表会  メタバースアバタの属性が パーソナルスペースの形状に及ぼす効果分析

<4D F736F F D208CF68BA48C6F8DCF8A C30342C CFA90B68C6F8DCF8A7782CC8AEE967B92E8979D32288F4390B394C529332E646F63>

OS バージョンアップ実行後のご注意 OS バージョンアップ後 更新完了通知が自動的にNECカシオモバイルコミュニケーションズ株式会社の運用するサーバへ送信されます なお NECカシオモバイルコミュニケーションズ株式会社は送信された情報を OS バージョンアップ以外の目的には利用いたしません また

バリデーション基準 1. 医薬品 医薬部外品 GMP 省令に規定するバリデーションについては 品質リスクを考慮し 以下の バリデーション基準 に基づいて実施すること 2. バリデーション基準 (1) バリデーションの目的バリデーションは 製造所の構造設備並びに手順 工程その他の製造管理及び品質管理の

(fnirs: Functional Near-Infrared Spectroscopy) [3] fnirs (oxyhb) Bulling [4] Kunze [5] [6] 2. 2 [7] [8] fnirs 3. 1 fnirs fnirs fnirs 1

多変量解析 ~ 重回帰分析 ~ 2006 年 4 月 21 日 ( 金 ) 南慶典

ISO ISO ISO ISO ISO ISO ISO ISO/TR 機械類の安全性 機械類への常設接近手段 第 2 部 : 作業用プラットフォーム及び通路機械類の安全性 機械類への常

延命セキュリティ製品 製品名お客様の想定対象 OS McAfee Embedded Control 特定の業務で利用する物理 PC 仮想 PC や Server 2003 Server 2003 ホワイトリスト型 Trend Micro Safe Lock 特定の業務で利用するスタンドアロン PC

P01

Robot Platform Project(RPP) "Spur" "YP-Spur" rev. 4 [ ] Robot Platform Project(RPP) WATANABE Atsushi 1.,,., Fig. 1.,,,,,.,,,..,,..,,..,,,,. "

1. はじめに 本書は スプリット演算器 MFS2 用コンフィギュレータソフトウェア の取扱方法 操作手順 注意事項などを説明したものです Windows の操作や用語を理解している方を前提にしています Windows の操作や用語については それぞれのマニュアルを参照してください 1.1. MFS

2014 年電子情報通信学会総合大会ネットワークシステム B DNS ラウンドロビンと OpenFlow スイッチを用いた省電力法 Electric Power Reduc8on by DNS round- robin with OpenFlow switches 池田賢斗, 後藤滋樹

1.民営化

2-2 需要予測モデルの全体構造交通需要予測の方法としては,1950 年代より四段階推定法が開発され, 広く実務的に適用されてきた 四段階推定法とは, 以下の4つの手順によって交通需要を予測する方法である 四段階推定法将来人口を出発点に, 1 発生集中交通量 ( 交通が, どこで発生し, どこへ集中

Microsoft Word - 01.docx

Microsoft PowerPoint - 口頭発表_折り畳み自転車

名古屋大学.indd

OS バージョンアップ実行中のご注意 OS バージョンアップ中は 故障の原因になりますので 絶対に N-03E 本体の電源を切ったり 電池パックを外したりしないでください OS バージョンアップ中は 電話の発着信を含めすべての機能がご利用になれません OS バージョンアップ中は 他のアプリケーション

機能紹介:コンテキスト分析エンジン

航空機の運動方程式

インターリーブADCでのタイミングスキュー影響のデジタル補正技術

, COMPUTATION OF SHALLOW WATER EQUATION WITH HIERARCHICAL QUADTREE GRID SYSTEM 1 2 Hiroyasu YASUDA and Tsuyoshi HOSHINO

Using VectorCAST/C++ with Test Driven Development

RMS(Root Mean Square value 実効値 ) 実効値は AC の電圧と電流両方の値を規定する 最も一般的で便利な値です AC 波形の実効値はその波形から得られる パワーのレベルを示すものであり AC 信号の最も重要な属性となります 実効値の計算は AC の電流波形と それによって

4Kチューナー(録画)_TU-BUHD100

Microsoft Word - 92.doc

Transcription:

三菱重工技報 Vol.51 No.1 (2014) 新製品 新技術特集 技 術 論 文 40 ポテンシャル法によるロボット製品の障害物回避技術の開発 Development of Obstacle Avoidance Control for Robot Products Using Potential Method 彌 城 祐 亮 *1 江 口 和 樹 *2 Yusuke Yashiro Kazuki Eguchi 岩 崎 聡 *3 山 内 由 章 *4 Satoshi Iwasaki Yoshiaki Yamauchi 中 田 昌 宏 *5 Masahiro Nakata 深海 宇宙空間 放射線発生領域等 人間が進入することが困難な場所で活動が可能なロボッ ト製品開発において 人間の遠隔操作無しで動作する自律製品が求められている より一般的な 問題への適用を考えると 自律動作には三次元空間での 未知の障害物を回避する機能が重要 である そこでポテンシャル法を用いて制御ロジックを構築し,その有効性をシミュレーションにより 検証した 具体的には 水平方向 鉛直方向それぞれの障害物回避経路を 障害物座標のみか ら適切に導出し 三次元空間での回避動作を効率化できることを確認した 1. はじめに 当社は従来 深海 宇宙空間 原子力発電所内等 人間の進入が難しい過酷環境下での運 用を前提とするロボット製品を複数開発している 近年ではコンピュータ センサ性能の発達に伴 い 自律動作を行う製品の開発が可能となった ロボットの自律化は 運用距離延長や人間の長 時間操作による疲労の低減が見込めることから お客様ニーズの一つとして市場からも要求され ている 一方 ロボット製品投入が想定される領域では 災害等の想定外事象により 事前に用意され たマップの情報に含まれない障害物が存在する可能性がある このとき 障害物回避機能を有す るロボットは障害物回避用の動作に切り替わるが 例えば事前に設定したシーケンス動作による 回避を行う場合 障害物形状によっては 同じ動作の繰り返しや過剰な回避動作となり 非効率 的な場合がある したがって ロボット製品の回避動作効率化には 障害物形状に応じた回避経路を導出する手 法が必要である また ロボットを投入する上でより一般的な問題への適用を考えると 手法は三 次元空間での移動に適用可能であることが望ましい 経路導出手法の一つに ポテンシャル法(1)と呼ばれる手法がある これは 目標位置と障害物 位置にポテンシャル関数を定義し その関数の勾配を進行方向とすることで 障害物を回避しつ つ 目標位置へと移動する経路を導出する手法である また ポテンシャル関数の計算は三次元 *1 技術統括本部 長崎研究所 *2 技術統括本部 長崎研究所 工博 *3 技術統括本部 長崎研究所 主席研究員 *4 交通 輸送ドメイン 船舶 海洋事業部 技術統括室 工博 *5 防衛 宇宙ドメイン 艦艇事業部 潜水艦部

三菱重工技報 Vol.51 No.1 (2014) 41 に拡張可能であり その時点までの観測情報を基に リアルタイムに経路を変化させられるため 未知の障害物にも対応可能となる 本稿では このポテンシャル法を用いてロボット製品の障害物回避経路を導出する制御ロジッ クを構成し その有効性をシミュレーションにより検証した結果を紹介する 制御ロジックでは セ ンサによって検知された障害物の座標を一定数保存 更新し その障害物座標と目標位置座標 からポテンシャル関数と勾配を計算し 目標方向を導出する ポテンシャル関数には重み係数を 設定し これを状況に応じて増減させる機能を追加することで 回避経路の効率化や 目標方向 がうまく導出されず ある場所に留まってしまう課題 停留 を解決している 制御ロジックの検証で は パラメータを変更することなく 異なる形状の障害物によるシミュレーションを行うことで セン サによる検知結果のみで適切な経路を導出し 効率的な障害物回避が可能であることを確認し ている 2. ポテンシャル法 ポテンシャル法は 障害物と目標位置の座標にポテンシャル関数と呼ばれる関数を定義し そ の関数の勾配 座標成分ごとの偏微分 から進むべき方向を導出する手法である 図1に2次元 平面上に設計したポテンシャル場 ポテンシャル関数の重ね合わせ のイメージを示す 障害物 に対しては正無限大 目標位置に対しては負無限大に発散するような滑らかな関数を定義する 障害物に定義した関数の勾配は 障害物から離れる方向であり 目標位置に定義した関数の勾 配は 逆に接近する方向となる この関数群の重ね合わせがポテンシャル場となり ポテンシャル 場の勾配 微分 は 障害物から離れ かつ目標位置へ接近する方向となる 障害物を回避しな がら目標位置へ向かうには この勾配をその時点でのロボットの目標方向とすれば良い ポテン シャル法の特徴として 次の3つが挙げられる その時点での情報に従って目標方向を導出するため 地形マップ等の事前情報が不要 である 障害物を発見した時点でポテンシャル場が変化するため リアルタイムに目標経路を変更 できる x y 等の座標成分別に勾配を導出するため 三次元空間への適用が容易に行える 図1 ポテンシャル法による経路導出例 二次元平面 二次元平面状に定義したポテンシャル関数の勾配から 目標方向 経路が導出 される一連の流れを示している

式 (1),(2) にポテンシャル関数を,(3) にポテンシャル場の計算を,(4) に勾配の計算をそれぞれ 示す,, 1 1 2 2 2 42 1,, 2 2 2 2,, 3 grad,,,, 4,, : ロボット位置,,, : 障害物位置,,, : 目標位置,, : 重み 3. 制御ロジックの構築 ポテンシャル法による制御ロジックの概要を図 2に示す このロジックにより導出した目標方向に従ってアクチュエータの制御を行うことで, ロボットが移動する 図 2 ポテンシャル法による障害物回避制御ロジックポテンシャル法を用いたロジックのフローと, その実機とアクチュエータに対する位置付けを示す ロジックには, より効率的に障害物回避と目標位置への接近を行うために, 以下の機能を追加している 障害物座標を保存するリストを更新する機能 ポテンシャル場の重みを調整する機能障害物座標リストの更新により, 回避済みの障害物を削除し, 考慮すべき障害物のみ対象としたポテンシャル場を作成できる また重みを調整することで, より迅速な障害物の回避や, 回避優先度の低い障害物の経路への影響を小さくすることも可能となる 障害物座標の保存と, 重みを調整する機能により, 停留を回避する手順を図 3に示す すべてのポテンシャル関数からポテンシャル勾配を計算したとき, 目標位置と障害物によるポテンシャルが釣り合う場合がある このとき, 目標方向は小さい異符号の値が繰り返し導出され, その結果ロボットがその場に停止するような経路となる これは停留と呼ばれ, ポテンシャル法を用いる際の課題となる この問題に対して, 一度接近した障害物を保存したまま, 目標位置の重みを減らすことで一時的に障害物から離れる動作を行う この後に再度目標位置の重みを元に戻すと, 記憶している障害物を回避しつつ目標位置へ向かうため, 停留を回避する経路となる

43 図 3 停留回避ロジック重みを変動させることで, 障害物と目標位置のポテンシャル勾配が釣り合う停留状態を抜け出していくロジックの流れを示す 4. シミュレーション 制御ロジックの有効性をシミュレーションによって確認した まず二次元空間での, 制御ロジックの機能確認の結果を図 4に示す 障害物と目標位置は, 停留が生じる位置関係とした また, ここでは制御ロジック単体検証のため, 障害物は全方向検知可能とし, ロボットは動特性を持たず, 目標方向に遅れなく移動できるものとした その結果, 障害物と停留を回避しながら目標位置へ進んでいることから, 制御ロジックの妥当性が確認できた 図 4 二次元平面でのシミュレーション結果開発した制御ロジックにより, 障害物と停留を回避しながら目標位置へ向かうことが可能であることを示している 次に, 鉛直方向を加えた三次元空間でのシミュレーションの結果を示す このシミュレーションでは, ロボットの動特性を考慮しており, 実際のロボットと同等の特性となっている 障害物は前方と, 左右斜め下 2 方の, 計 3 方向のみ検知可能とした 水平方向の障害物回避の結果を図 5に, 鉛直方向の障害物回避の結果を図 6にそれぞれ示す 回避動作効率化の比較のため, 方向転換 上昇 直進 下降からなる, シーケンス動作での回避結果を併せて示す それぞれのシミュレーションにおいて, 制御ロジックの設計パラメータは共通である また, 障害物座標は事前に与えず, 設定したセンサ範囲に存在する障害物の座標を取得できるとした

44 図 5 三次元空間での鉛直方向障害物回避シミュレーション結果 図 6 三次元空間での水平方向障害物回避シミュレーション結果 まず鉛直方向の障害物回避では, シーケンス動作による回避の場合, 一度上昇を行った後障害物を上方に避け, 次に下降して目標位置へ到達する ポテンシャル法による制御ロジックの場合, 衝突を避ける程度に旋回しながら上昇し, 上方回避後に下降して目標位置へ到達する 次に水平方向回避においては, シーケンス動作による回避の場合, 不要な旋回動作を含んだ経路を通り, 目標位置へ移動している これはそれぞれの動作の動作時間を, 鉛直方向の障害物回避に合わせてチューニングしているためであり, シーケンス動作では異なる障害物形状に対して, 非効率な回避動作を含んでしまうことを示している 一方でポテンシャル法による回避経路には, 不要な旋回は含まれない このことからポテンシャル法では, 障害物形状に合わせたポテンシャル関数が定義されるため, 同一のパラメータでも適切な回避経路が導出されることを示している なお, ポテンシャル法による経路が波打つような形状となる理由は, ロボットの動特性を考慮すると, 瞬時に方向転換ができないためである 回避動作効率を数値的に比較したものを図 7,8に示す 鉛直方向, 水平方向いずれもポテンシャル法による制御ロジックの所要時間が短く, アクチュエータ動作量から求めた消費エネルギー指標も減少している 以上より, ポテンシャル法による制御ロジックによって, 障害物回避動作が効率化されているといえる 図 7 障害物回避所要時間比較 図 8 障害物回避消費エネルギー比較

45 5. まとめ ポテンシャル法を用いて三次元空間に適用可能な障害物回避ロジックを構築し, シミュレーションによって, 障害物形状に応じた経路が導出され, 回避動作が効率化されることを確認した 本手法を用いることで, 未知の障害物を効率的に回避することを可能にし, 複雑地形領域でのロボットの自律移動につながると考えられる 今後はセンサ範囲やロボットの動特性に合わせた重みの設計手法等, 実機適用に必要な要素の研究を進めていき, 運用につなげたい 参考文献 (1) Jen-Hui Chuang et al,an analytically tractable potential field model of free space and its application in obstacle avoidance,ieee TRANSACTIONS ON SYSTEMS,MAN,AND CYBERNETICS-PART B: CYBERNETICS,VOL. 28,NO. 5 (1998) p.729-736