(57)【要約】【課題】実機ロボットを動作させることなくロボットの駆動源にかかる負荷状態を診断することができる診断方法を提供する。【解決手段】ロボットの動作プログラムにしたがってロボット動作をシミュレーションし（Ｓ２）、その間のロボット各自軸の角度変化量を取得して（Ｓ３）、角度変化量から角速度を求め（Ｓ４）、角速度と軸の慣性モーメントを掛けて負荷を求め（Ｓ５）、求めた負荷を加算して累積負荷を求めて（Ｓ６）、累積負荷をマハラノビスタグチ法における特徴量としてマハラノビス距離を求め（Ｓ７）、マハラノビス距離が基準値か以下か否かにより不具合の有無を判定する（Ｓ８〜１０）。【選択図】図２
2005293332.tif 000002
【特許請求の範囲】
【請求項１】
ロボットの動作を動作プログラムにしたがってシミュレーションし、当該シミュレーション中における前記ロボットの診断対象軸の動作に伴う角度変化量を任意の一定時間間隔ごとに取得する段階と、
前記角度変化量から、角速度、加速度、前記診断対象軸に作用する前記時間間隔ごとの負荷、および前記時間間隔ごとの前記負荷を全て加算した累積負荷の少なくともいずれか一つを求める段階と、
前記角速度、前記加速度、前記負荷、および前記累積負荷の少なくともいずれか一つの値から前記診断対象軸の負荷状態を診断する段階と、
を有することを特徴とするロボットの負荷状態診断方法。
【請求項２】
前記診断対象軸の負荷状態を診断する段階は、
前記角速度、前記加速度、前記負荷、および前記累積負荷のうち少なくともいずれか一つの値をマハラノビスタグチ法における特徴量としてマハラノビス距離を求める段階と、
前記マハラノビス距離と、あらかじめ決められた基準値を比較して前記診断対象軸の診断を行う段階と、
を有することを特徴とする請求項１記載のロボットの負荷状態診断方法。
【請求項３】
前記角速度は、前記角度変化量を微分することにより求めることを特徴とする請求項１または２記載のロボットの負荷状態診断方法。
【請求項４】
前記加速度は、前記角速度を微分することにより求めることを特徴とする請求項１または２記載のロボットの負荷状態診断方法。
【請求項５】
前記負荷は、前記角速度に前記診断対象軸の慣性モーメントを掛けることにより求めることを特徴とする請求項１または２記載のロボットの負荷状態診断方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ロボットの負荷状態診断方法に関する。
【背景技術】
【０００２】
ロボットの負荷状態を推定する方法として、従来から、ロボットや工作機械においては、その駆動源であるサーボモータへの電流量を測定して負荷状態を把握し、過負荷によるモータコイルの焼損やアンプなどの電気系統の破損を防止する方法がある。
【０００３】
たとえば、サーボモータを備えた制御系において、サーボモータの測定電流を用いて推測した発熱量により熱シミュレーションを行い、制御系の温度特性を推測し、過負荷を判定するアラーム特性と比較することによってサーボモータの過負荷を監視する方法において、過負荷と判定する条件を緩和してアラーム特性を設定し、測定電流に重み付けを行い、重み付けした電流を用いて温度特性を求め、求めた温度特性と緩和したアラーム特性とを比較してサーボモータの過負荷監視を行う方法がある（特許文献１）。
【特許文献１】特開平９−９３７９５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、従来の方法は、実際のロボットが動作中におけるサーボモータに流れる電流を測定して発熱量を推定して過負荷の監視を行っているため、実施にロボットが運用状態に供された後、または試験的に動作させなければ、負荷状態を診断することができない。また、電流の測定値から、負荷を求めているため、電流測定時の瞬間における負荷状態を推定しているに過ぎず、累積的にモータに加わる負荷量を知ることができないという問題がある。
【０００５】
そこで、本発明の目的は、ロボットを実際に動作させることなく、ロボットの駆動源における負荷状態を診断することのできるロボットの負荷状態診断方法を提供することである。
【課題を解決するための手段】
【０００６】
上記目的を達成するための本発明は、ロボットの動作を動作プログラムにしたがってシミュレーションし、当該シミュレーション中における前記ロボットの診断対象軸の動作に伴う角度変化量を任意の一定時間間隔ごとに取得する段階と、前記角度変化量から、角速度、加速度、前記診断対象軸に作用する前記時間間隔ごとの負荷、および前記時間間隔ごとの前記負荷を全て加算した累積負荷の少なくともいずれか一つを求める段階と、前記角速度、前記加速度、前記負荷、および前記累積負荷の少なくともいずれか一つの値から前記診断対象軸の負荷状態を診断する段階と、を有することを特徴とするロボットの負荷状態診断方法である。
【０００７】
また、本発明においては、前記診断対象軸の負荷状態を診断する段階は、前記角速度、前記加速度、前記負荷、および前記累積負荷のうち少なくともいずれか一つの値をマハラノビスタグチ法における特徴量としてマハラノビス距離を求める段階と、前記マハラノビス距離と、あらかじめ決められた基準値を比較して前記診断対象軸の診断を行う段階と、を有することを特徴とする。
【発明の効果】
【０００８】
本発明のロボットの負荷状態診断方法によれば、シミュレーションによって診断対象軸の角度変化量を取得し、その値から角速度、加速度、負荷、および累積負荷の少なくともいずれか一つを求めて負荷状態診断することとしたので、実際にロボットを動作させることなく、診断対象軸の負荷状態を診断することが可能となる。特に、マハラノビスタグチ法を用いて診断することで、確実に負荷状態を診断することができる。また、累積負荷を用いることで、累積された負荷による不具合の診断も可能となる。
【発明を実施するための最良の形態】
【０００９】
以下、図面を参照して本発明を実施するための最良の形態を説明する。
【００１０】
本実施の形態は、ロボットの動作シミュレーションを利用して評価対象となるロボットの負荷状態を診断する負荷診断方法である。
【００１１】
図１は評価対象となるロボットの一例を示す図である。
【００１２】
このロボットは６軸の多関節型ロボットであって、このロボット１は、ロボット全体を旋回させる第１軸１１、ロボット全体を傾斜させる第２軸１２、アーム２１を傾斜させる第３軸１３、アーム２１を回転させる第４軸１４、ハンド２２を傾斜させる第５軸１５、ハンド２２を回転させる第６軸１６よりなる６軸ロボットである。各軸には、サーボモータが内蔵されており、各軸を駆動する。ハンド２２には、エンドエフェクタとしてたとえば、スポット溶接機、把持装置などが接続される。
【００１３】
図２は負荷状態の診断方法を実施するための処理手順を示すフローチャートである。
【００１４】
まず、診断を行うロボットの動作プログラムをシミュレーション装置に取り込む（Ｓ１）。なお、ここで動作プログラムとは、ロボットを実際に動作させる際の動作軌跡とその動作軌跡上での動作速度を定めたものであり、具体的には、ロボットのティーチングプログラムまたは実際のロボットの動作軌跡と速度を示すものなどである。また、シミュレーションを行う装置は、周知のＣＡＤ装置などロボットの動作を、その動作プログラムにしたがってシミュレートできるものであれば特に限定されない。
【００１５】
次に、ロボットの動作シミュレーション装置によって取得した動作プログラムによりロボット動作をシミュレーションする（Ｓ２）。そして、シミュレーションにおけるロボットの動作に伴うロボット各軸の角度変化量を任意の一定時間間隔ごとに取得する（Ｓ３）。なお、ここで任意の一定時間間隔は、特に限定されず、ロボットの１サイクルの動作時間を分割できる値であればよい。たとえば、１サイクルと５等分、１０等分、または１００等分など、まったく任意の時間間隔でよいが、実際には使用するシミュレーション装置（コンピュータ）の処理能力に応じて決定することになる。
【００１６】
次に、取得した角度変化量から時間単位の角速度を算出する（Ｓ４）。これには、取得した角度変化量を時間微分すればよい。
【００１７】
次に、算出した角速度から各軸に作用する負荷を算出する（Ｓ５）。作用する負荷の算出には、あらかじめ入力された慣性モーメントを上記角速度に掛けることによって求める。
【００１８】
次に、各軸に作用する負荷の頻度から累積負荷を算出する（Ｓ６）。累積負荷の算出は、ロボットの１サイクルの動作におけるステップＳ５で求めた負荷の値を全て加算することによって算出する。
【００１９】
次に、ステップＳ６で求めた累積負荷の値と、その他の因子を、マハラノビスタグチ法（以下ＭＴ法と称する）で用いる特徴量として用いて、マハラノビス距離を求める（Ｓ７、詳細後述）。
【００２０】
次に、求めたマハラノビス距離をあらかじめ求められているマハラノビス距離による診断の基準値と比較する（Ｓ８）。このとき、ステップＳ７で求めたマハラノビス距離が基準値以下であれば、不具合の発生する可能性はないものと判断する（Ｓ９）。一方、ステップＳ７で求めたマハラノビス距離が基準値よりも大きければ不具合の発生する可能性が高いものと判断する（Ｓ１０）。不具合発生の可能性が高い場合には、その後、ロボットの動作（ティーチング）を見直したり、負荷（主にハンド２２に取り付けられているエンドエフェクタやハンド２２で把持するワークなど）を見直したりすることになる。
【００２１】
ＭＴ法は、周知のとおり、複数の特徴量から正常空間の距離を求める必要がある。そこで、本実施の形態では、特徴量として、１：ステップＳ６で求めた累積負荷値、２：ロボット型式、３：各軸における負荷重量、４：各軸の動作した角度最大値、５：各軸の負荷モーメント６：各軸の負荷イナーシャを用いている。
【００２２】
ここで、５：各軸の負荷モーメントは、Ｍｘ＝Σｍｉ・Ｌｘｉで求める。ただし、Ｍｘは各軸の負荷モーメント、ｍｉはエンドエフェクタやワークなど各要素の質量、Ｌｘｉは各軸中心からエンドエフェクタおよびワークなどの各要素の重心位置までの距離であり、添え字ｘは各軸をあらわし（６軸ロボットの場合１〜６となる）、ｉは各要素をあらわす（要素の数だけ変化する数となる）。また、６：各軸の負荷イナーシャは、Ｊｘ＝Σｍｉ・Ｌｘｉ２＋ΣＪｘｉで求める。ただし、Ｊｘは各軸の負荷イナーシャ、ｍｉはエンドエフェクタやワークなど各要素の質量、Ｌｘｉは各軸中心からエンドエフェクタおよびワークなどの各要素の重心位置までの距離、Ｊｘｉは各軸方向のエンドエフェクタやワークなどの要素の自分自身の慣性モーメントであり、添え字ｘは各軸をあらわし、ｉは要素をあらわす。
【００２３】
ＭＴ法では、これらの特徴量を用いて正常空間の距離を求める。したがって、実際に稼動しており、不具合の発生していないロボットについて、その動作軌跡と動作速度、または動作プログラムを用いたシミュレーション結果から上記各特徴量を取得して上記同様にしてマハラノビス距離を算出し、これによって算出されたマハラノビス距離の最大値が、上記ステップＳ４で用いる基準値となる。なお、ＭＴ法の実行には、市販のソフトウェアを利用すればよい。ＭＴ法を実行するための市販ソフトとしては、たとえば、株式会社オーケン製ＭＴＳ  ｆｏｒ  ｗｉｎｄｏｗｓ  Ｅｘｃｅｌ版などがある。
【００２４】
なお、上記特徴量のうち、特徴量２：ロボット形式は、評価対象となるロボット形式が同じものであれば、特徴量としての入力を省いてもよい。また、特徴量３〜６についても、ロボットが同じで負荷重量も同じ（すなわち動作経路および速度のみが異なる）場合には、省略してもよい。一方、本実施形態のように、ロボット形式をも特徴量とすることで、形式の異なるさまざまなロボットを用いて正常空間を規定すれば、ロボットの形式が異なるロボットを対象に診断を行う場合でも、負荷状態の診断を的確に行うことが可能となる。
【００２５】
ここで、実際の診断結果の一例を説明する。
【００２６】
図３はロボットの第６軸における動作角度を０．１秒ごとに示すグラフであり、（ａ）は過負荷を与えたロボット、（ｂ）は過負荷のないロボットである。図４は図３に示した動作角度を時間微分した結果を示すグラフである（（ａ）は過負荷を与えたロボット、（ｂ）は過負荷のないロボットである。図５も同様）。図５は、図４に示した角速度の結果に慣性モーメントを掛けて負荷を求めた結果を示すグラフである。図６は図５に示した負荷の結果を時間間隔ごとに過負荷を与えたロボットと過負荷のないロボットについて示したグラフである。図７は図６に示した各負荷の値を全て加算した累積負荷を示すグラフである。各図はそれぞれ上記ステップＳ３〜６の処理結果に対応するものである。
【００２７】
図３〜５の結果から、ここでは、過負荷を与えたロボットでは、過負荷を与えていないロボットに対していずれの結果も激しく振動していることがわかり、図６において、負荷が一時的に高くなっていることがわかる。そして、図７に示したように、累積負荷も過負荷を与えたロボットのほうが高くなっている。
【００２８】
さらに、本実施形態では、ロボットの累積負荷の値と、ロボットの諸特性を示す上述したその他の因子とからマハラノビス距離を求めてＭＴ法により判定している。
【００２９】
その結果、上記過負荷をかけたロボットについては、マハラノビス法における正常空間（基準値）から外に出た位置になり、一方、過負荷を与えていないロボットは正常空間内となった。
【００３０】
次に、さらに多くの実機ロボットについて、用いる特徴量のみを変えて、上述した手順にしたがってシミュレーションによりマハラノビス距離を求め、実機ロボットの実際のガタつきを比較した。
【００３１】
用いた特徴量は、上述した特徴量のうち１：累積負荷の値を後述する各図の説明のとおり変えたものであり、その他の特徴量２〜８については、同じである。なお、以下の各図においては、いずれもロボットの第６軸についての結果である。
【００３２】
図８は、特徴量として、角度変化量を用いた場合のマハラノビス距離とガタつきの関係を示すグラフである。実際の入力特徴量は、上記特徴量２〜８の他に、角速度変化量（一つ前の時点の角度からある時点の角度に変化した量、以下同様）の累計、角速度の累計、角度変化の３σ（３σは統計処理におけるばらつきの値、以下同様）、角度変化量の３σである。
【００３３】
なお、角度変化は回転方向が変化するまでの角度であり、角度変化量は１つ前の角度からある時点の角度まで変化した量である。
【００３４】
図から、マハラノビス距離とガタ量に比例関係があることがわかる。
【００３５】
図９は、特徴量として、角度変化量と角速度を用いた場合のマハラノビス距離とガタつきの関係を示すグラフである。実際の入力特徴量は、上記特徴量２〜８の他に、角速度変化量の累計、角速度の累計、角度変化の３σ、角度変化量の３σ、角速度の最大値、角速度の平均値、角速度の３σである。
【００３６】
図から、マハラノビス距離とガタ量に比例関係があり、比例線に実測値が近くなっていることがわかる。
【００３７】
図１０は、特徴量として、角度変化量と角加速度を用いた場合のマハラノビス距離とガタつきの関係を示すグラフである。実際の入力特徴量は、上記特徴量２〜８の他に、角速度変化量の累計、角加速度の累計、角度変化の３σ、角度変化量の３σ、角加速度の最大値、角加速度の平均値、角加速度の３σである。
【００３８】
図から、マハラノビス距離とガタ量に比例関係があり、比例線に実測値が良く乗っていることがわかる。
【００３９】
図１１は、特徴量として、角度変化量と累積負荷を用いた場合のマハラノビス距離とガタつきの関係を示すグラフである。実際の入力特徴量は、上記特徴量２〜８の他に、角速度変化量の累計、累積負荷の累計、角度変化の３σ、角度変化量の３σ、累積負荷の最大値、累積負荷の平均値、累積負荷の３σである。なお、ここで累積負荷とは、
図から、マハラノビス距離とガタ量に比例関係があり、最も良く比例線の周囲に実測値が集まっていることがわかる。
【００４０】
以上の結果から、シミュレーションで算出したマハラノビス距離と、実機におけるロボットの過負荷が原因となって現れるガタつきとが比例しており、マハラノビスタグチ法により不良判定が可能であることがわかる。特に図１１に示した累積負荷を特徴量として用いる場合が最も比例関係が良く、正確に不良判定を行えることがわかる。
【００４１】
なお、上記の各グラフからわかるように、角度変化量、角速度、各加速度、負荷、および累積負荷の値を、いずれか一つ、またはさまざまに組み合わせて特徴量とすることでロボット駆動源の負荷状態をＭＴ法を用いて診断することが可能である。
【００４２】
また、本実施形態によれば、シミュレーションとマハラノビスタグチ法を組み合わせて診断することで、これから動作させるロボットの諸特性（負荷モーメントや負荷イナーシャ、ロボットの型式などの特徴量とする値）と、そのロボットの動作プログラムさえあれば、実際のロボット（実機）を動かすことなく、負荷の状態を予測して診断することが可能となる。このため、実地運用前の段階で、過負荷状態や累積負荷の状態などから、動作プログラムや負荷重量の見直しができる。
【００４３】
以上説明した実施形態においては、処理手順として、シミュレーション中の各軸の角度変化量を取得することとしたが、必ずしもロボットの全ての軸についての角度変化量を取得する必要はなく、過負荷が起こる可能性のある軸（診断対象軸）についての角度変化量のみを取得し、その軸についてのみ負荷の状態を診断するようにしてもよい。
【００４４】
これは、たとえば、先の具体例のように第６軸についての結果のみを求めたい場合には、第６軸についてのみ角度変化量を取得すればよいことである。なお、先の具体例として第６軸を診断対象軸として示したのは、ハンドに持たせるエンドフェクタやワークなどの荷重が変化し、また、ハンドの移動量が変化した場合に、この第６軸にかかる負荷が最も大きく変化するためである。なお、第６軸以外でも、当然に、本発明を適用して過負荷状態の診断に利用できることは言うまでもない。たとえば、過負荷の状態は、第６軸の他に、第４軸、第５軸なども第１〜第３軸に比較して負荷が過剰にかかりやすくなる傾向があるので、その場合には、第４〜６軸の３軸を診断対象軸としてそれらの角度変化量を取得するようにするとよい。
【００４５】
このように診断対象軸を決めて診断することで、どの軸に過負荷や累積負荷などによる不具合が生じるかを個別に診断することができる。
【００４６】
また、基準値については、ＭＴ法における正常空間の距離の値を使用すればよいが、これに限らず正常空間の距離の値に、さらに危険率（正常範囲内でも不具合の発生する可能性のある率）を考慮して、正常空間の距離よりも小さな値を基準値としたり、逆に正常空間よりも広い大きな値としてもよい。
【００４７】
以上説明した実施形態においては、角度変化量を微分することにより角速度を、同様に、角速度を微分することにより角加速度を求めているが、簡易的には、たとえば、図１２に示すように、角度変化量（または角速度）の時間変化の値を所定の値で区切った線（イ〜チ）を引き、時間変化の値がイ〜チの各線を横切った値を区間ごとに求めてその横切った回数を角速度（または角加速度）としてもよい。
【００４８】
このような簡易的な微分値を用いても累積的な過負荷の状態を推定して、ＭＴ法を用いて診断することができる。
【産業上の利用可能性】
【００４９】
本発明は、新たに設置するロボットの動作に伴う過負荷の状態を検証できるほか、既存のロボットにおいても新たに動作軌跡をティーチングしなおす場合、また、既存のロボットにおいて、現状の見た目では、不具合が発生していないものでも、今後不具合の発生する可能性があるか否かを検証することもできる。また、本発明は、サーボモータを工作機械などにおける駆動源にかかる過負荷の状態を診断することもできる。
【図面の簡単な説明】
【００５０】
【図１】診断対象となるロボットの一例を示す図面である。
【図２】負荷状態診断方法の処理手順を示すフローチャートである。
【図３】ロボットの第６軸における動作角度を０．１秒ごとに示すグラフである。
【図４】図３に示した動作角度を時間微分した結果を示すグラフである。
【図５】図４に示した角速度の結果に慣性モーメントをかけて負荷を求めた結果を示すグラフである。
【図６】図５に示した負荷の結果を時間間隔ごとに過負荷を与えたロボットと過負荷のないロボットについて示したグラフである。
【図７】図６に示した各負荷の値を全て加算した累積負荷を示すグラフである。
【図８】特徴量として、角度変化量を用いた場合のマハラノビス距離とガタつきの関係を示すグラフである。
【図９】特徴量として、角度変化量と角速度を用いた場合のマハラノビス距離とガタつきの関係を示すグラフである。
【図１０】特徴量として、角度変化量と角加速度を用いた場合のマハラノビス距離とガタつきの関係を示すグラフである。
【図１１】特徴量として、角度変化量と累積負荷を用いた場合のマハラノビス距離とガタつきの関係を示すグラフである。
【図１２】簡易的な微分、積分を説明するための特徴量の時間変化を示すグラフである。
【符号の説明】
【００５１】
１…ロボット、
１１…第１軸、
１２…第２軸、
１３…第３軸、
１４…第４軸、
１５…第５軸、
１６…第６軸、
２１…アーム、
２２…ハンド。
【図１】
2005293332.tif 000003
【図２】
2005293332.tif 000004
【図６】
2005293332.tif 000005
【図７】
2005293332.tif 000006
【図８】
2005293332.tif 000007
【図９】
2005293332.tif 000008
【図１０】
2005293332.tif 000009
【図１１】
2005293332.tif 000010
【図１２】
2005293332.tif 000011
【図３】
2005293332.tif 000012
【図４】
2005293332.tif 000013
【図５】
2005293332.tif 000014