HOME >> 鉄道模型工学　>　モータ特性の測定とモデル化　モータ定数の推定方法

モータ特性の測定とモデル化　　モータ定数の推定方法

なぜ、モータの測定にこだわるのか？

　我が実験室では無謀にも鉄道模型の動力車のモデル化に挑戦しています。　そして、念願であったモータ単体の特性が測定出来るようになったので、その測定データを使って、トルク定数とか逆起電力定数などを推定し、モデル毎の特徴を見出そうと工夫しています。　その方法は、モータの外部から測定出来る特性を使用し、モータの機能を発揮している各定数を推定する方法です。　モータの研究をする訳では無いので、原理的な厳密性を追求するものではありませんが、モータに供給された電圧と電流値および回転数から、モータの出力トルクを計算出来れば良いと言う立場です。

　さらに、もう一つの目的は、動力車のブラックボックになっている定数の推定である。　「測定項目の追加　走行中のモータ回転数と端子電圧を測ろう」で説明したように、モータ単品状態での電圧と電流、および回転数とトルクを測定してモータを構成する定数をあらかじめ求めておけば、動力車に組み込まれた状態でも、走行中の動力車のモータ端子電圧と電流、およびモータ回転数よりモータトルクを算出することが可能となるのだ。　すると、減速機ブロックの入出力の４変数が判明することになり、伝達効率やウォームギヤの様子が解明できる、ブラックボックの中身を覗く有力な道具立てが完成するのである。

 

■　定数の推定の方法

　　.

　先回も同様なイラストを紹介したが、そのイラストを少し修正する。　変数の関係を示す線について、外から測定できる変数を赤線で示した。　そしてよく見ると、４端子ある変数は関係式をはさんで、二つのグループに別れることが分かる。　即ち、電圧・電流と回転数、および、電流・出力トルクと回転数のグループである。

　前者は定数 Ra、Ke、Eb が関係し、下の変数とは独立していることが分かる。　また、後者についても同様に、定数 Kt、λm、Rm とが関係して上の変数とは独立している。

　これは、前者が電圧に関する関係式であり、後者がトルクに関する関係式なのであるが、右に示したイラストの中の式を見ると、右辺が Nm に関する一次式であり、それぞれの勾配と Y 切片が求めたい定数を示している。

　また、左辺は二つの変数と一つの定数から成っているので、その定数を仮定すれば、左辺は測定データ毎に計算出来るのである。　その計算されたデータを、 Nm を横軸に取った散布図で示せば、その勾配と Y 切片について一次近似式より推定できるのである。　即ち、左辺の定数を仮定すれば、それに対応して、他の二つの求めたい定数を直接示唆しているのである。　その上に、二つのグループはお互いに独立しているので、それぞれのグループ単独で計算出来るのである。

　次に、仮定した一つの定数と求められた二つの定数のセットが、どれだけ実際と合致しているのかは、求めた定数を使って計算された計算値と実測データを比較して、その相関具合がチェックして行けば、確からしさが確認できるのである。　前者の電圧に関する関係では回転数を、後者のトルクに関する関係では出力トルクをその評価項目にして、計算値と実測値との比較を実施すれば、求める定数の確からしさが確認できるのである。

　この理屈を使って、EXCELで計算しながら、４個の変数より、６個の定数を逐次近似法もどきの手法で推定することが出来るのである

 

■　具体的な手順

　モータ単品で測定した特性データをもとに、モータの特性を決める６個の定数を容易に推定する方法を紹介しよう。

 

　　.

１）　トルクに関する関係式から定数 Kt 、Rm 、λm の最適値をそれぞれ推定する

　回転数 Nm 、電流 I 、出力トルクTm’の関係より、Kt 、Rm 、λm を推定をするものである。　それは、測定された電流 I と回転数Nm より、モータの関係式で説明した式（11）と式（12）を使用して出力トルクTm’を計算する。　この計算された値と Tm’の実測値との比較により、もっともよく合致する Kt 、Rm 、λm を求めるものである。　しかし、この３っの定数 Kt 、Rm 、λm を同時に推定しなければならないのであるが、その組み合わせを考えると膨大なケースとなってしまう。　しかし、 モータ軸損失トルクの計算式を活用すれば、その制約を回避できるのである。

 

　　

　定数推定の方法は、いわゆる 手動式逐次計算法 である。　ある定数を一つ仮定すれば、EXCELで自動的に計算してグラフ化してくれるので、評価項目により近い値はどちらの方向かを示してくれるのである。

 

２）　電圧に関する関係式から定数 Ke 、Ra 、Eb の最適値をそれぞれ推定する

　同じ考えを定数 Ke 、Ra 、Eb にも当てはめて、デルタ電圧δと勝手に名づけたが、この要素を計算することによって、これらの定数も同様に簡単に最適値が推定出来る。

　



　　　　　　

３）　推定した定数値の検証

　上記の方法で、６個の定数を推定したが、果たして本当に的を得た値であるのかどうを検証して置く必要がある。　そこで、実際の測定グラフの上に、推定した定数値から求めた線図を引いてみて、測定データとどれだけ合致しているのか見ておくと安心である。

◆　無負荷特性の比較

　まず、モータ単体での無負荷状態での測定データを使って、実測値と計算値の比較を実施する。

　　.

　

◆　モータのトルク特性の比較

　同様に負荷をかけた時のデータも比較する。

　

 

 

■　実際の計算方法

　実際の計算方法について説明する。　計算シートは Excel を使用して必要な欄に計算式を埋め込み、計算された値をグラフに表示する。

拡大図

　　.

　この計算シートは、測定データをペーストするゾーン、そのデータから必要項目を計算するゾーン、推定する定数を設定するゾーン、それらの値をメモして置く保管ゾーン、および決定した定数から実測値と比較するための計算ゾーン（図の下の方に隠れてしまったが・・・・・）、関係するグラフの表示ゾーンから成り立っている。

　次の、計算手順を説明する。

  1. 測定データを整理してデータ設定ゾーンに貼り付ける。
  2. 定数設定ゾーンに、Kt あるいは、Ra の値を適当に設定して関係するグラフの様子を見る。　これで大体の様子を把握したら、小さめの数値と大きめの数値を代入して定数保管ゾーンの両側に値だけをペーストして置く。　
  3. 設定する定数を中間の値にセットして評価項目の内容をチェックする。　そして定数保管ゾーンに値をペーストする。
  4. そして評価項目のグラフを見ながら、設定する定数の範囲を狭めていく。　手動の逐次近似法である！
  5. トルクと電圧の関係式は、それぞれ独立しているので別々に検討して行けば良いし、今回は相関係数をなどを自動的に計算するようにしたので、楽々である。
  6. ＯＫであれば、モータ特性や無負荷特性のグラフにて、実測値と計算値のマッチング具合をチェックする。
  7. もともと、モータ作動の不安定性や、やや怪しい測定方法の影響で、測定データ自身のバラツキや信頼性への疑問などを念頭において、ほどほどのマッチング具合で良しとすべきでしょう。

■　意外な収穫

　この方法で使用した「モータ軸損失トルク」と「モータのδ電圧」のグラフを見ていると、新たな発見がある事に気が付いた。

　例えば前者の「モータ軸損失トルク」グラフからは、モータ自身の損失トルクが意外と大きい事が分かる。　出力トルクが30～40 gf-mm に対して、損失トルクが 20 gf-mm も有るのである。　モータの回転子が発生したトルクの約1/3は自分の摩擦として消えている・・・・・・・？　その原因は整流子のブラシ部分の摩擦ではないかと推定するのであるが、不確実である。　また、２０，０００ rpm 以上では特性が変化しているので、摩擦現象が変化しているのではないかとも考えられる。

　一方の「モータのδ電圧」のグラフからは、高速回転部では非線形になっているのがわかる。　それまでは、綺麗な線形を示したいたので、理論どうりの逆起電力が発生していたものの、高速回転では何故非線形になるのだろうか？　もしかして、ブラシ部の電圧降下が影響しているのではないだろうか？　高速になるとだんだん電圧降下量が増えてくると考えれば合点がいくのであるが・・・・・・・・。

　また、モータ無負荷時の電流・電圧特性での非線形状態は、「モータ軸損失トルク」グラフの形態と類似しているので、原因は軸の摩擦トルクのバラツキが主要因であると考えることが出来ないだろうか？

　さらに、測定中に作動音の喧しいモータもあり、その音が途中で変化しているもので見受けられた。　これは明らかに軸の共振現象であり、軸のアンバランスが、ダイナミックバランスも含めた大きいのではないかと思われる。　そのようなモータは、手で持った状態で回転させると、ブルブルと手に感じるのですぐに分かるのである。

　などなど・・・・・・いろいろな現象について、あれこれ想定をめぐらすことが出来るが、やはりデータのN増しが必要であり、多くの現象を総合的に判断しないと間違った判断をすることになる。　

 

■　まとめ

　はじめは、ほどほどのマッチング具合で良しとすべきと思っていたが、モータ特性の測定データの向上と計算シートの改善によって、上記のグラフの例で示すように、意外とマッチする定数を推定することが出来るようになった。　まずは今回の取り組みは成功であると評価しよう。

　今後は、多くのモータのデータを集めて、モデル毎の特徴を見い出せたらと考えている。

 

************　　モータ特性のモデル化　改良版　（2016/10/26） を再編集　********