JP4275426B2 - Evaluation method, evaluation program and evaluation system for demagnetization of permanent magnet - Google Patents

Evaluation method, evaluation program and evaluation system for demagnetization of permanent magnet Download PDF

Info

Publication number
JP4275426B2
JP4275426B2 JP2003049453A JP2003049453A JP4275426B2 JP 4275426 B2 JP4275426 B2 JP 4275426B2 JP 2003049453 A JP2003049453 A JP 2003049453A JP 2003049453 A JP2003049453 A JP 2003049453A JP 4275426 B2 JP4275426 B2 JP 4275426B2
Authority
JP
Japan
Prior art keywords
demagnetization
permanent magnet
factor
node
mesh
Prior art date
Legal status (The legal status 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 status listed.)
Expired - Fee Related
Application number
JP2003049453A
Other languages
Japanese (ja)
Other versions
JP2004257879A (en
Inventor
正嗣 中野
晴之 米谷
均 岡元
光弘 川村
浩二 宮田
信宏 斎藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shin Etsu Chemical Co Ltd
Mitsubishi Electric Corp
Toshiba Mitsubishi Electric Industrial Systems Corp
Original Assignee
Shin Etsu Chemical Co Ltd
Mitsubishi Electric Corp
Toshiba Mitsubishi Electric Industrial Systems Corp
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 Shin Etsu Chemical Co Ltd, Mitsubishi Electric Corp, Toshiba Mitsubishi Electric Industrial Systems Corp filed Critical Shin Etsu Chemical Co Ltd
Priority to JP2003049453A priority Critical patent/JP4275426B2/en
Publication of JP2004257879A publication Critical patent/JP2004257879A/en
Application granted granted Critical
Publication of JP4275426B2 publication Critical patent/JP4275426B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Description

【0001】
【発明の属する技術分野】
この発明は、永久磁石の減磁を評価する磁界解析の評価方法、評価プログラム及び評価システムに関するものである。特に、有限要素法あるいは積分要素法を用いた永久磁石の減磁の評価手法に関するものである。また、この評価手法を用いて永久磁石型モータや発電機など永久磁石を具備した電気機器の特性に減磁が及ぼす影響を評価するための評価方法に関するものである。
【0002】
【従来の技術】
従来、減磁曲線から永久磁石の減磁率を算出し、面積を用いた重み計算により永久磁石内部の平均的な減磁率を求める手法がある。この手法においては、平均的な減磁率を求めるには有効な手段である(例えば、非特許文献1参照)。
【0003】
また、B−Hカーブを再定義することにより減磁量を評価する手法がある(例えば、非特許文献2参照)。
【0004】
【非特許文献1】
石橋利之(安川電機)著「埋込磁石同期モータにおける減磁の条件とその影響」、MAG−01−181、2001年発行、第61頁−第65頁
【非特許文献2】
棗(なつめ)田充俊、高林博文(住友特殊金属株式会社)著「永久磁石の減磁量解析手法の検討」、JMAG Users Conference、2002年10月24日〜25日、第4−2頁−第4−7頁
【0005】
【発明が解決しようとする課題】
しかしながら、従来の手法では、発生した減磁によってモータをはじめとする電機機器の特性にどのような影響を与えるのかを知るには不充分である。特に、モータや発電機のコギングトルク、トルク脈動、誘起電圧波形、電磁加振力が減磁によってどのように変化するかは減磁率を反映した高精度な磁界解析を実施しなければ知ることはできないという問題点があった。
【0006】
また、従来のB−Hカーブを再定義する手法では、減磁率の経時変化を知ることはできないという問題点があった。
【0007】
この発明は、上述のような課題を解決するためになされたもので、第1の目的は、従来例より減磁の評価を高精度化することができ、減磁量のより正確な定量化を実現することができ、また、減磁による回転機の特性変化(コギングトルク、トルク脈動、誘起電圧波形、電磁加振力など)を予測することができる永久磁石の減磁の評価方法、評価プログラム及び評価システムを得るものである。
【0008】
第2の目的は、経時変化を考慮できるため回転機を始めとする電機機器の長期信頼性の向上を図ることができる永久磁石の減磁の評価方法、評価プログラム及び評価システムを得るものである。
【0009】
第3の目的は、このような評価手法を、永久磁石を具備した電気機器の設計に取り入れることにより設計の精度を向上することができる永久磁石の減磁の評価方法、評価プログラム及び評価システムを得るものである。
【0010】
【課題を解決するための手段】
この発明に係る永久磁石の減磁の評価方法は、永久磁石型モータの電子機巻線の電流、及び回転子と固定子の位置関係が変化して、前記永久磁石型モータの永久磁石が所定の逆磁界、所定の温度にさらされる各時点において、有限要素法または積分要素法を用いた磁界解析により、前記永久磁石をメッシュ分割し、メッシュの各要素または各節点の第1の減磁率あるいは第1の減磁量を求めるステップと、履歴である、前記各時点で求めたメッシュの各要素または各節点の第1の減磁率あるいは第1の減磁量のうち、メッシュの各要素または各節点の第1の減磁率あるいは第1の減磁量の最大値を、メッシュの各要素または各節点における第2の減磁率あるいは第2の減磁量として求めるステップと、メッシュの各要素または各節点における第2の減磁率あるいは第2の減磁量に応じて、前記永久磁石の残留磁束密度、保磁力、減磁曲線のうち少なくとも1つを求めるステップと、メッシュの各要素または各節点における第2の減磁率あるいは第2の減磁量に基づき、前記永久磁石型モータのコギングトルク、負荷時トルク、誘起電圧の評価例のうち少なくとも1つを求めるステップとを含むものである。
【0011】
【発明の実施の形態】
実施の形態1.
この発明の実施の形態1に係る永久磁石の減磁の評価システムについて図面を参照しながら説明する。図1は、この発明の実施の形態1に係る永久磁石の減磁の評価システムの構成を示す図である。なお、各図中、同一符号は同一又は相当部分を示す。
【0012】
図1において、評価システムは、記録装置1と、減磁評価手段2と、記憶装置3と、電機機器特性評価手段4とから構成されている。なお、破線で示すように、記録装置1は、必ずしも必要ではない。
【0013】
記録装置1は、例えばハードディスク装置などから構成され、磁束密度データなどの磁界解析結果や、減磁評価用パラメータΔと初期減磁率の特性データ、減磁評価用パラメータΔと経時変化率の特性データ、B−Hカーブデータなどをデータベースとして格納する。
【0014】
減磁評価手段2及び電機機器特性評価手段4は、例えばパソコン等の本体(CPU)で動作するプログラムである。また、記憶装置3は、例えばメモリなどから構成される。この減磁評価手段2は、永久磁石の減磁率あるいは減磁量を評価する。また、電機機器特性評価手段4は、コギングトルク、負荷時トルク、誘起電圧などの電機機器特性を評価する(詳細な説明は実施の形態5を参照)。
【0015】
減磁評価手段2は、磁界解析結果、及び減磁評価用パラメータΔと初期減磁率、経時変化率、B−Hカーブなどの減磁評価に必要なデータを記録装置1から呼び出す。そして、永久磁石各部の減磁率を算出する。その際、必要に応じて記憶装置3とデータをやり取りしてもよい。この減磁評価手段2により得られた結果は、記録装置1に記録される。
【0016】
また、減磁が電機機器にどのような影響を及ぼすかについては、減磁評価手段2は、減磁率計算結果を電機機器特性評価手段4に受け渡す。そして、この電機機器特性評価手段4は、コギングトルク、負荷時のトルク、誘起電圧などの特性を計算する。
【0017】
つぎに、この実施の形態1に係る永久磁石の減磁の評価システムの動作について図面を参照しながら説明する。
【0018】
図2は、この発明の実施の形態1に係る永久磁石の減磁の評価システムの動作を示すフローチャートである。また、図3は、この発明の実施の形態1に係る永久磁石の減磁の評価システムの減磁評価の具体例を示す図である。
【0019】
ここでは、減磁を反映した具体的な手法について説明する。
【0020】
まず、ステップ101において、減磁評価手段2は、磁界解析(既存のプログラム)により磁場条件における減磁を算出する。
【0021】
すなわち、磁界解析を行い、永久磁石内の減磁率あるいは減磁量(物理量)を算出する。この減磁率あるいは減磁量の算出の方法は、実施の形態3で説明するように減磁評価用パラメータΔを用いる手法でもよいし、実施の形態4で説明するように減磁曲線から幾何学的に算出する手法でもその他いかなる手法でもよい。また、永久磁石の各部における減磁率あるいは減磁量は、後述するような評価ステップ毎の履歴を反映させるものとする。ただし、評価ステップの数は1つでもよいし、2以上(複数)であってもよいことは言うまでもない。
【0022】
次に、ステップ102において、減磁が生じるかを判定する。
【0023】
次に、ステップ103において、減磁が生じる場合には、評価ステップ毎の履歴を追い、永久磁石各部の減磁率(減磁量)を評価条件に反映する。例として、残留磁束密度を変化させる。保磁力を変化させる。減磁曲線を変える。
【0024】
次に、ステップ104において、評価を実施し、減磁の影響を調べる。
【0025】
すなわち、減磁が生じる場合には、減磁率あるいは減磁量を次の磁界解析に反映する作業を行う。このとき、減磁率あるいは減磁量に応じて永久磁石の磁気特性を変化させる。例えば、残留磁束密度を変化させてもよいし、保磁力を変化させてもよい。さらに、減磁曲線を変えてもよい。例えば、減磁率が仮に1%であったとすれば、残留磁束密度を減磁する前の99%とすればよい。
【0026】
このように永久磁石内の減磁率を反映させて、磁界解析を行う。このようにすれば、減磁が生じたときの磁束密度の分布などを得ることができる。また、評価する対象が電気機器であった場合にはその特性の変化も知ることができる。
【0027】
もう少し詳しく説明する。図3は、直方体形状の永久磁石10の例を示す。なお、図3では軸方向から見ているため永久磁石10は長方形として示す。
【0028】
強い逆磁界が掛かったり、あるいは温度上昇が起こったりして永久磁石10の右肩部分を中心に減磁が生じたとする。また、その減磁率の分布は、図3(a)に示すように、右肩部分の減磁率の大きい部分(密の斜線)で1%、その周辺(粗の斜線)で0.5%、さらにその周辺では減磁は発生せずに減磁率は0%であったとする。ただし、実際、減磁率は連続的に変化して分布するが、ここでは簡単のため離散的に変化するとした。
【0029】
次に、この減磁率を次の評価に反映する。ここでは、図3(b)に示すように、永久磁石10内部の残留磁束密度Bを減磁率に応じて変化させる手法を取った。場合によっては、保磁力も変化させてもよいし、非線形性を考慮するため減磁曲線を減磁率に応じて変化させてもよい。
【0030】
このように永久磁石10内の減磁率を反映させて、磁界解析を再度行う。ここでは、図3(c)に示すように、2点X、Yを結ぶ直線上の磁束密度をグラフ(曲線)にした。減磁する前は左右対称なグラフAであったが、右肩部分を中心として減磁が発生したときには、従来例のグラフBに比べて、本実施の形態のグラフCは、右肩部分の磁束密度が低下しているのが分かる。
【0031】
従来の評価手法では局所的な減磁の影響を考慮できないため、正確な磁束密度分布を予測することはできない。一方、本実施の形態では減磁が生じた場合の磁束密度分布を高精度に知ることができる。
【0032】
以上により、本実施の形態に係る評価手法では、永久磁石において減磁が生じた場合の磁束密度分布を高精度に知ることができるという効果があり、さらに磁束密度分布を高精度に知ることで永久磁石を具備した電気機器の特性を高精度に知ることができるという効果もある。
【0033】
すなわち、この実施の形態1は、永久磁石の減磁を評価する手法において、永久磁石の各部の減磁率あるいは減磁量の導出を行い、導出された永久磁石の各部の減磁率あるいは減磁量を反映させる評価手法であって、永久磁石の各部の減磁率あるいは減磁量は評価ステップごとの履歴を追うことによって導出されることを特徴とするものである。
【0034】
実施の形態2.
この発明の実施の形態2に係る永久磁石の減磁の評価方法について図面を参照しながら説明する。図4は、この発明の実施の形態2に係る永久磁石型モータの断面構成を示す図である。
【0035】
図4において、永久磁石型モータ20は、固定子鉄心21と、電機子巻線22と、スロット23と、回転子シャフト24と、永久磁石25とを備える。
【0036】
上記の実施の形態1においては、非常に簡単なモデルを用いて説明した。しかし、例えば、図4に示すような永久磁石型モータ20の場合、運転中は電子機巻線22の電流が時々刻々と変化する上、さらに回転子と固定子の位置関係も時々刻々と変化する。すなわち、永久磁石25が受ける逆磁界が時々刻々と変化することになる。
【0037】
したがって、磁界解析を行って減磁率を算出する際、ある時刻あるいは回転子と固定子のある位置関係のみに注目して減磁率を算出しても、正確な減磁率の分布とは異なる場合が考えられる。このような場合には、永久磁石25の各部が経験する逆磁界や温度の最も悪い条件によって減磁率が決まると思われるので、磁界解析の各ステップの履歴を追って減磁率の算出を行う必要がある。上記の実施の形態1においても、評価ステップの履歴を追うと説明したが本実施の形態ではその方法を具体的に示すことにする。
【0038】
図5は、この発明の実施の形態2に係る永久磁石の減磁の評価方法の具体例を示す図である。
【0039】
図5(a)〜(e)は、評価ステップごとの履歴を追っていく様子を示す。ただし、実際、減磁率は連続的に変化して分布するが、ここでは簡単のため離散的に変化するとした。
【0040】
評価ステップ1においては、永久磁石25の左肩部分を中心に減磁が発生していたとする。それが評価ステップ2以降においては、電機子電流の値や回転子と固定子の位置関係が変化したりして、減磁の発生する位置が移動していき、そして、評価ステップNでは、永久磁石25の右肩部分に減磁の中心が移動したとする。
【0041】
以上の評価ステップ1から評価ステップNまでの履歴を追うと、図5(d)のような減磁率の分布を知ることができる。
【0042】
さらに、この減磁率の分布を次の評価に反映させる。ここでは、減磁率に応じて残留磁束密度Bを変化させた。上記実施の形態1において述べたように、場合によっては保磁力も変化させてもよいし、非線形性を考慮するため減磁曲線を減磁率に応じて変化させてもよい。
【0043】
以上のような履歴を追う手法を用いることによって、減磁率の分布を高精度に定量的に求めることができるという効果がある上、上記減磁率の分布によって永久磁石を具備する電気機器の特性がどのように変化するのか高精度に予測することができるという効果もある。
【0044】
図6は、この発明の実施の形態2に係る永久磁石の減磁の評価方法の別の具体例を示す図である。
【0045】
近年では、計算機の発達により有限要素法や積分要素法などを用いた磁界解析が電気機器の設計の場で活躍している。有限要素法や積分要素法などでは、磁気回路をモデル化する際、磁石部分はメッシュ分割される。
【0046】
図6(a)〜(d)は、メッシュ分割された永久磁石に関して、評価ステップの履歴を追っていく様子を示す。この場合は、永久磁石のメッシュの各要素において減磁率が定義される場合であり、節点に対して減磁率が定義される場合でもよいことは言うまでもない。
【0047】
評価ステップ1〜Nまでの履歴を追っていくことにより、この場合は、磁石表面に減磁率1%の領域が分布するという結果を得ることとなった。
【0048】
なお、この例では、メッシュの要素ごとに減磁率の履歴を追っていたが、複数の要素をまとめて扱うこともできる。このようにすれば、モデル作成等に掛かる手間が軽減されるという効果がある。
【0049】
履歴を追っていく際に、メッシュの各要素または各節点の減磁率あるいは減磁量の最大値をその要素または節点における減磁率あるいは減磁量とすれば、より高精度に永久磁石の減磁の分布を予測できるという効果がある。図6において、例えば、メッシュの左最上段の要素について着目すると、評価ステップ1では減磁率が1%、評価ステップ2では減磁率が0%、そして、評価ステップNでは減磁率が0%であるので、この要素の減磁率の最大値は1%となる。
【0050】
実施の形態3.
この発明の実施の形態3に係る永久磁石の減磁の評価方法について図面を参照しながら説明する。
【0051】
従来から、永久磁石の熱減磁について研究がなされている。例えば、文献『「長期安定性を考慮した永久磁石磁気回路の磁束条件」電気学会論文誌 1998』には、熱減磁を評価するためのパラメータを提案し、このパラメータの大小が熱減磁に大きく影響するということが述べられている。
【0052】
そのパラメータとは、永久磁石の残留磁束密度をB[T]とし、保磁力を[A/m]、リコイル透磁率をμ[H/m]とし、永久磁石の各部の動作点における磁界をH[A/m](ここでHは動作点における磁界の絶対値と定義する)としたとき、以下の式(1)で表される。
【0053】
Δ={μ−H)}/B (1)
【0054】
一般に、永久磁石の減磁は、永久磁石がある逆磁界、ある温度条件に暴露された時点で発生する初期減磁とその状態に暴露される時間によって変化する経時変化に分けられる。式(1)のパラメータΔと初期減磁率、経時変化率の特性は実測により求めることができる。
【0055】
図7は、この発明の実施の形態3に係るパラメータΔと初期減磁率の関係を示す図である。図7において、縦軸は初期減磁率[%]を表す。また、横軸はパラメータΔを表す。原点から離れる程、それらの値は大きくなる。
【0056】
また、図8は、この発明の実施の形態3に係るパラメータΔと経時変化率の関係を示す図である。図8において、縦軸は経時変化率[%]を表す。また、横軸はパラメータΔを表す。原点から離れる程、それらの値は大きくなる。
【0057】
図7、図8において、パラメータΔが小さいほど初期減磁率、経時変化率も大きい。上記実施の形態1及び2において、減磁率あるいは減磁量の算出法については触れなかったが、図7、図8に示すパラメータΔと初期減磁率、経時変化率の特性を用いる方法がある。例えば、有限要素法の評価において各要素あるいは各節点のパラメータΔを求めることができる。このパラメータΔと、図7、図8の特性から各要素あるいは各節点の減磁率を算出することができるのである。
【0058】
図7に示す初期減磁特性のみから減磁率あるいは減磁量を算出する手法は、モータや発電機などで突発短絡が生じて大電流が流れた場合の減磁の量や分布を予測することができるという効果がある。
【0059】
また、図7及び図8に示す初期減磁と経時変化から、ある一定時間暴露された後の減磁率あるいは減磁量を算出することができるため、モータや発電機などで長時間運転された後の減磁の量や分布を予測することができるという効果がある。
【0060】
また、式(1)のパラメータΔを用いた手法は、実測に基づいた手法であるので、減磁の量や分布を高精度に予測することができるという効果もある。
【0061】
さらに、B−Hカーブを再定義する従来例では、永久磁石の種類に応じてデータを用意しなければならないが、本実施の形態で用いるパラメータΔは例えば希土類磁石では永久磁石の種類が変わってもほぼ一定の特性を示すので、パラメータΔと初期減磁率の関係、パラメータΔと経時変化率の関係が1種類あるだけで様々な永久磁石についての減磁の評価が可能であるという効果もある。
【0062】
実施の形態4.
この発明の実施の形態4に係る永久磁石の減磁の評価方法について図面を参照しながら説明する。図9は、この発明の実施の形態4に係る作図によって減磁率の求める例を示す図である。
【0063】
上記の実施の形態3では減磁を評価するパラメータΔを用いる手法であったが、従来からよく知られているように幾何学的に求める手法もある。
【0064】
図9において、縦軸は磁束密度Bと磁化Jを表す。また、横軸は磁界Hを表す。動作点41にある状態からΔHの逆磁界を受けて動作点が変化し、再び磁界を取り除き元に戻ったときの減磁量ΔBは、動作点41→42→43→44→45のように求めることができる。
【0065】
この手法では初期減磁を減磁曲線から算出できるので、式(1)のパラメータΔを用いるような手法とは異なり、比較的に容易に減磁量を推定できるという効果がある。
【0066】
実施の形態5.
この発明の実施の形態5に係る永久磁石の減磁の評価方法について図面を参照しながら説明する。図10は、この発明の実施の形態5に係る永久磁石型モータの評価用モデルを示す図である。また、図11は、この発明の実施の形態5に係る減磁評価の結果の一例(減磁率とその分布)を示す図である。
【0067】
図10において、永久磁石型モータの評価用モデル50は、固定子鉄心51と、電機子巻線52と、回転子シャフト54と、永久磁石55の各部のモデルを備える。
【0068】
この永久磁石型モータの評価用モデル50を用いて磁界解析を行い、定格電流よりはるかに大きい電流が流れた場合を想定し、減磁率を求めたものが図11である。
【0069】
図11に示すように、永久磁石55に逆方向の起磁力が掛かる部分に減磁が確認された。なお、この結果は、回転子位置や電流値が変化する評価ステップ毎の履歴を追った結果である。
【0070】
さらに、この減磁率の分布を評価に反映させ、コギングトルク、負荷時のトルク波形、誘起電圧波形の評価を実施した。結果を図12〜14にそれぞれ示す。
【0071】
図12は、この発明の実施の形態5に係るコギングトルクの評価例を示す図である。図12において、縦軸はコギングトルク[Nm]を表す。また、横軸は回転子位置[度(電気角)]を表す。破線の波形は減磁発生前、実線の波形は減磁発生後をそれぞれ表す。
【0072】
また、図13は、この発明の実施の形態5に係る負荷時のトルクの評価例を示す図である。図13において、縦軸はトルク[Nm]を表す。また、横軸は回転子位置[度(電気角)]を表す。破線の波形は減磁発生前、実線の波形は減磁発生後をそれぞれ表す。
【0073】
さらに、図14は、この発明の実施の形態5に係る誘起電圧波形の評価例を示す図である。図14において、縦軸は無負荷誘起電圧[V]を表す。また、横軸は回転子位置[度(電気角)]を表す。破線の波形は減磁発生前、実線の波形は減磁発生後をそれぞれ表す。
【0074】
コギングトルクやトルク脈動を低減する設計にしても、何らかの原因で大電流が流れるなどして減磁が発生すればコギングトルクやトルク脈動が異なることが分かる。さらに、誘起電圧波形も若干変化することが分かる。
【0075】
また、音や振動に関係する電磁加振力の評価も可能なことは言うまでもない。
【0076】
このように、減磁が発生した場合、モータ特性にどの程度影響があるのか予測することが可能となるのである。
【0077】
このような評価手法を用いることにより、永久磁石型回転機の設計を高精度化できるという効果がある。さらに、永久磁石型回転機に大電流が流れたときや、高温になった場合における機器の特性を予測することができるので、これらを設計にフィードバックし、電気設計、磁気回路設計の信頼性の向上を図ることができるという効果がある。
【0078】
【発明の効果】
この発明に係る永久磁石の減磁の評価方法は、以上説明したとおり、永久磁石において減磁が生じた場合の磁束密度分布を高精度に知ることができ、さらに磁束密度分布を高精度に知ることで永久磁石を具備した電気機器の特性を高精度に知ることができるという効果を奏する。
【図面の簡単な説明】
【図1】 この発明の実施の形態1に係る永久磁石の減磁の評価システムの構成を示す図である。
【図2】 この発明の実施の形態1に係る永久磁石の減磁の評価システムの動作を示すフローチャートである。
【図3】 この発明の実施の形態1に係る永久磁石の減磁の評価システムの減磁評価の具体例を示す図である。
【図4】 この発明の実施の形態2に係る永久磁石型モータの断面構成を示す図である。
【図5】 この発明の実施の形態2に係る永久磁石の減磁の評価方法の具体例を示す図である。
【図6】 この発明の実施の形態2に係る永久磁石の減磁の評価方法の別の具体例を示す図である。
【図7】 この発明の実施の形態3に係るパラメータΔと初期減磁率の関係を示す図である。
【図8】 この発明の実施の形態3に係るパラメータΔと経時変化率の関係を示す図である。
【図9】 この発明の実施の形態4に係る作図によって減磁率の求める例を示す図である。
【図10】 この発明の実施の形態5に係る永久磁石型モータの評価用モデルを示す図である。
【図11】 この発明の実施の形態5に係る減磁評価の結果の一例(減磁率とその分布)を示す図である。
【図12】 この発明の実施の形態5に係るコギングトルクの評価例を示す図である。
【図13】 この発明の実施の形態5に係る負荷時のトルクの評価例を示す図である。
【図14】 この発明の実施の形態5に係る誘起電圧波形の評価例を示す図である。
【符号の説明】
1 記録装置、2 減磁評価手段、3 記憶装置、4 電機機器特性評価手段、10 永久磁石、20 永久磁石型モータ、21 固定子鉄心、22 電機子巻線、23 スロット、24 回転子シャフト、25 永久磁石、50 永久磁石型モータの評価用モデル。
[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an evaluation method, an evaluation program, and an evaluation system for magnetic field analysis for evaluating demagnetization of a permanent magnet. In particular, the present invention relates to a method for evaluating demagnetization of a permanent magnet using a finite element method or an integral element method. Further, the present invention relates to an evaluation method for evaluating the influence of demagnetization on the characteristics of an electric device equipped with a permanent magnet such as a permanent magnet type motor or a generator using this evaluation method.
[0002]
[Prior art]
Conventionally, there is a method of calculating a demagnetization factor of a permanent magnet from a demagnetization curve and obtaining an average demagnetization factor inside the permanent magnet by calculating a weight using an area. This method is an effective means for obtaining an average demagnetization factor (see Non-Patent Document 1, for example).
[0003]
There is also a method for evaluating the amount of demagnetization by redefining the BH curve (see, for example, Non-Patent Document 2).
[0004]
[Non-Patent Document 1]
Toshiyuki Ishibashi (Yaskawa Electric) “Demagnetizing Conditions and Effects on Embedded Magnet Synchronous Motors”, MAG-01-181, 2001, pp. 61-65 [Non-patent Document 2]
Natsume Tadatoshi, Takabayashi Hirofumi (Sumitomo Special Metals Co., Ltd.) “Examination of Demagnetization Analysis Method of Permanent Magnet”, JMAG Users Conference, October 24-25, 2002, 4-2- Page 4-7 [0005]
[Problems to be solved by the invention]
However, the conventional method is insufficient to know how the generated demagnetization affects the characteristics of electrical equipment such as a motor. In particular, it is not possible to know how cogging torque, torque pulsation, induced voltage waveform, and electromagnetic excitation force of motors and generators change due to demagnetization without performing high-precision magnetic field analysis that reflects the demagnetization factor. There was a problem that it was not possible.
[0006]
Further, the conventional method of redefining the BH curve has a problem that it is impossible to know the change in demagnetization with time.
[0007]
The present invention has been made in order to solve the above-described problems. The first object is to make it possible to improve the accuracy of demagnetization compared to the conventional example, and to more accurately quantify the amount of demagnetization. Evaluation method and evaluation of permanent magnet demagnetization capable of predicting changes in characteristics of rotating machine due to demagnetization (cogging torque, torque pulsation, induced voltage waveform, electromagnetic excitation force, etc.) Obtain programs and evaluation systems.
[0008]
The second object is to obtain an evaluation method, an evaluation program, and an evaluation system for demagnetization of a permanent magnet that can improve the long-term reliability of electrical equipment such as a rotating machine because changes with time can be taken into account. .
[0009]
The third object is to provide an evaluation method, an evaluation program, and an evaluation system for demagnetization of a permanent magnet that can improve the design accuracy by incorporating such an evaluation method into the design of an electrical device having a permanent magnet. To get.
[0010]
[Means for Solving the Problems]
According to the method for evaluating demagnetization of a permanent magnet according to the present invention, the current of the winding of the electric machine of the permanent magnet type motor and the positional relationship between the rotor and the stator are changed, and the permanent magnet of the permanent magnet type motor is predetermined. The permanent magnet is divided into meshes by a magnetic field analysis using a finite element method or an integral element method at each time point exposed to a predetermined magnetic field, and a first demagnetizing factor of each element or each node of the mesh. The step of obtaining the first demagnetization amount and the history, each element of the mesh or each element of the mesh among the first demagnetization factor or the first demagnetization amount of each node obtained at each time point. Obtaining a maximum value of the first demagnetization factor or first demagnetization amount of the node as the second demagnetization factor or second demagnetization amount of each element or each node of the mesh; and each element or each of the mesh At the node 2 to obtain at least one of a residual magnetic flux density, a coercive force, and a demagnetization curve of the permanent magnet according to a demagnetization factor of 2 or a second demagnetization amount, and a second element at each element or node of the mesh. Determining at least one of evaluation examples of cogging torque, on-load torque, and induced voltage of the permanent magnet motor based on the demagnetization factor or the second demagnetization amount .
[0011]
DETAILED DESCRIPTION OF THE INVENTION
Embodiment 1 FIG.
A permanent magnet demagnetization evaluation system according to Embodiment 1 of the present invention will be described with reference to the drawings. FIG. 1 is a diagram showing the configuration of a permanent magnet demagnetization evaluation system according to Embodiment 1 of the present invention. In addition, in each figure, the same code | symbol shows the same or equivalent part.
[0012]
In FIG. 1, the evaluation system includes a recording device 1, a demagnetization evaluation means 2, a storage device 3, and an electrical equipment characteristic evaluation means 4. Note that the recording apparatus 1 is not necessarily required as indicated by a broken line.
[0013]
The recording device 1 is composed of, for example, a hard disk device, etc., and includes magnetic field analysis results such as magnetic flux density data, demagnetization evaluation parameter Δ and initial demagnetization factor characteristic data, demagnetization evaluation parameter Δ and time-dependent change rate characteristic data. , BH curve data and the like are stored as a database.
[0014]
The demagnetization evaluation unit 2 and the electrical equipment characteristic evaluation unit 4 are programs that operate on a main body (CPU) such as a personal computer. Further, the storage device 3 is composed of a memory, for example. The demagnetization evaluation means 2 evaluates the demagnetization rate or demagnetization amount of the permanent magnet. Moreover, the electrical equipment characteristic evaluation means 4 evaluates electrical equipment characteristics such as cogging torque, load torque, and induced voltage (refer to Embodiment 5 for detailed description).
[0015]
The demagnetization evaluation means 2 calls the data necessary for demagnetization evaluation such as the magnetic field analysis result, the demagnetization evaluation parameter Δ, the initial demagnetization rate, the change rate with time, and the BH curve from the recording device 1. Then, the demagnetization factor of each part of the permanent magnet is calculated. At that time, data may be exchanged with the storage device 3 as necessary. The result obtained by the demagnetization evaluation means 2 is recorded in the recording device 1.
[0016]
Further, as to how demagnetization affects the electrical equipment, the demagnetization evaluation means 2 passes the demagnetization factor calculation result to the electrical equipment characteristic evaluation means 4. And this electrical equipment characteristic evaluation means 4 calculates characteristics, such as a cogging torque, the torque at the time of load, and an induced voltage.
[0017]
Next, the operation of the permanent magnet demagnetization evaluation system according to the first embodiment will be described with reference to the drawings.
[0018]
FIG. 2 is a flowchart showing the operation of the permanent magnet demagnetization evaluation system according to Embodiment 1 of the present invention. FIG. 3 is a diagram showing a specific example of demagnetization evaluation of the permanent magnet demagnetization evaluation system according to Embodiment 1 of the present invention.
[0019]
Here, a specific method reflecting demagnetization will be described.
[0020]
First, in step 101, the demagnetization evaluation means 2 calculates demagnetization under magnetic field conditions by magnetic field analysis (existing program).
[0021]
That is, magnetic field analysis is performed to calculate the demagnetization factor or demagnetization amount (physical quantity) in the permanent magnet. As a method for calculating the demagnetization factor or demagnetization amount, a method using the demagnetization evaluation parameter Δ as described in the third embodiment may be used, or a geometry from a demagnetization curve as described in the fourth embodiment. Any other method may be used. Further, the demagnetization rate or the amount of demagnetization in each part of the permanent magnet reflects a history for each evaluation step as will be described later. However, it goes without saying that the number of evaluation steps may be one or two or more (plural).
[0022]
Next, in step 102, it is determined whether demagnetization occurs.
[0023]
Next, when demagnetization occurs in step 103, the history for each evaluation step is followed and the demagnetization rate (demagnetization amount) of each part of the permanent magnet is reflected in the evaluation conditions. As an example, the residual magnetic flux density is changed. Change the coercivity. Change the demagnetization curve.
[0024]
Next, in step 104, evaluation is performed and the influence of demagnetization is examined.
[0025]
That is, when demagnetization occurs, an operation of reflecting the demagnetization rate or demagnetization amount in the next magnetic field analysis is performed. At this time, the magnetic characteristics of the permanent magnet are changed according to the demagnetization rate or the amount of demagnetization. For example, the residual magnetic flux density may be changed, or the coercive force may be changed. Further, the demagnetization curve may be changed. For example, if the demagnetization factor is 1%, the residual magnetic flux density may be 99% before demagnetization.
[0026]
Thus, the magnetic field analysis is performed by reflecting the demagnetization factor in the permanent magnet. In this way, the distribution of magnetic flux density when demagnetization occurs can be obtained. In addition, when the object to be evaluated is an electric device, the change in its characteristics can be known.
[0027]
I will explain a little more in detail. FIG. 3 shows an example of a rectangular parallelepiped permanent magnet 10. In addition, in FIG. 3, since it sees from an axial direction, the permanent magnet 10 is shown as a rectangle.
[0028]
It is assumed that demagnetization occurs around the right shoulder portion of the permanent magnet 10 due to a strong reverse magnetic field applied or a temperature rise. Further, as shown in FIG. 3 (a), the distribution of the demagnetization factor is 1% at the portion of the right shoulder where the demagnetization factor is large (the dense oblique line), 0.5% around it (the coarse oblique line), Further, it is assumed that no demagnetization occurs in the vicinity and the demagnetization factor is 0%. However, in practice, the demagnetization factor varies continuously and is distributed, but here it is assumed to change discretely for simplicity.
[0029]
Next, this demagnetization factor is reflected in the next evaluation. Here, as shown in FIG. 3 (b), we took the approach of changing depending on the demagnetizing factor of the residual magnetic flux density B r inside the permanent magnet 10. In some cases, the coercive force may be changed, or the demagnetization curve may be changed in accordance with the demagnetization factor in order to consider nonlinearity.
[0030]
In this way, the magnetic field analysis is performed again by reflecting the demagnetization factor in the permanent magnet 10. Here, as shown in FIG. 3C, the magnetic flux density on the straight line connecting the two points X and Y is shown as a graph (curve). Before the demagnetization, the graph A was symmetric. However, when demagnetization occurs around the right shoulder, the graph C of the present embodiment shows the right shoulder portion compared to the graph B of the conventional example. It can be seen that the magnetic flux density is decreasing.
[0031]
Since the conventional evaluation method cannot consider the influence of local demagnetization, an accurate magnetic flux density distribution cannot be predicted. On the other hand, in this embodiment, the magnetic flux density distribution when demagnetization occurs can be known with high accuracy.
[0032]
As described above, the evaluation method according to the present embodiment has an effect that the magnetic flux density distribution when demagnetization occurs in the permanent magnet can be known with high accuracy, and further, the magnetic flux density distribution can be known with high accuracy. There is also an effect that the characteristics of the electric device including the permanent magnet can be known with high accuracy.
[0033]
That is, in the first embodiment, in the method for evaluating the demagnetization of the permanent magnet, the demagnetization rate or demagnetization amount of each part of the permanent magnet is derived, and the demagnetization rate or demagnetization amount of each part of the derived permanent magnet is calculated. The demagnetization factor or demagnetization amount of each part of the permanent magnet is derived by following the history for each evaluation step.
[0034]
Embodiment 2. FIG.
A method for evaluating demagnetization of a permanent magnet according to Embodiment 2 of the present invention will be described with reference to the drawings. FIG. 4 is a diagram showing a cross-sectional configuration of a permanent magnet type motor according to Embodiment 2 of the present invention.
[0035]
In FIG. 4, the permanent magnet type motor 20 includes a stator core 21, an armature winding 22, a slot 23, a rotor shaft 24, and a permanent magnet 25.
[0036]
In the first embodiment, the description has been made using a very simple model. However, for example, in the case of a permanent magnet type motor 20 as shown in FIG. 4, the current of the electronic winding 22 changes every moment during operation, and the positional relationship between the rotor and the stator also changes every moment. To do. That is, the reverse magnetic field received by the permanent magnet 25 changes every moment.
[0037]
Therefore, when calculating the demagnetizing factor by performing a magnetic field analysis, even if the demagnetizing factor is calculated by paying attention only to a certain time or a certain positional relationship between the rotor and the stator, the distribution of the demagnetizing factor may be different. Conceivable. In such a case, the demagnetization factor is considered to be determined by the worst condition of the reverse magnetic field and temperature experienced by each part of the permanent magnet 25. Therefore, it is necessary to calculate the demagnetization factor following the history of each step of the magnetic field analysis. is there. In the first embodiment, it has been described that the history of the evaluation step is followed, but in this embodiment, the method will be specifically shown.
[0038]
FIG. 5 is a diagram showing a specific example of a method for evaluating demagnetization of a permanent magnet according to Embodiment 2 of the present invention.
[0039]
FIGS. 5A to 5E show how the history for each evaluation step is followed. However, in practice, the demagnetization factor varies continuously and is distributed, but here it is assumed to change discretely for simplicity.
[0040]
In the evaluation step 1, it is assumed that demagnetization has occurred around the left shoulder portion of the permanent magnet 25. In the evaluation step 2 and thereafter, the armature current value and the positional relationship between the rotor and the stator change, and the position where the demagnetization occurs is moved. It is assumed that the demagnetization center moves to the right shoulder of the magnet 25.
[0041]
By following the history from evaluation step 1 to evaluation step N, it is possible to know the demagnetization distribution as shown in FIG.
[0042]
Furthermore, this demagnetization factor distribution is reflected in the next evaluation. Here, the residual magnetic flux density Br was changed according to the demagnetization factor. As described in the first embodiment, the coercive force may be changed depending on the case, and the demagnetization curve may be changed according to the demagnetization factor in consideration of nonlinearity.
[0043]
By using the method of following the history as described above, there is an effect that the distribution of the demagnetization factor can be obtained quantitatively with high accuracy, and the characteristics of the electric device including the permanent magnet can be obtained by the distribution of the demagnetization factor. There is also an effect that the change can be predicted with high accuracy.
[0044]
FIG. 6 is a diagram showing another specific example of the method for evaluating demagnetization of the permanent magnet according to Embodiment 2 of the present invention.
[0045]
In recent years, with the development of computers, magnetic field analysis using the finite element method, the integral element method, etc. has been active in the design of electrical equipment. In the finite element method, the integral element method, and the like, when modeling a magnetic circuit, the magnet portion is divided into meshes.
[0046]
FIGS. 6A to 6D show a state in which the history of the evaluation step is followed for the mesh-divided permanent magnet. In this case, the demagnetization factor is defined in each element of the mesh of the permanent magnet, and it goes without saying that the demagnetization factor may be defined for the nodes.
[0047]
By following the history of evaluation steps 1 to N, in this case, a result that a region having a demagnetization rate of 1% is distributed on the magnet surface was obtained.
[0048]
In this example, the history of the demagnetization factor is tracked for each element of the mesh, but a plurality of elements can be handled together. In this way, there is an effect that the labor required for model creation and the like is reduced.
[0049]
When tracking the history, if the maximum value of the demagnetization factor or demagnetization amount of each element or node of the mesh is the demagnetization factor or demagnetization amount of that element or node, the demagnetization of the permanent magnet can be performed with higher accuracy. The effect is that the distribution can be predicted. In FIG. 6, for example, when attention is paid to the upper leftmost element of the mesh, the demagnetization factor is 1% in the evaluation step 1, the demagnetization factor is 0% in the evaluation step 2, and the demagnetization factor is 0% in the evaluation step N. Therefore, the maximum value of the demagnetization factor of this element is 1%.
[0050]
Embodiment 3 FIG.
A method for evaluating demagnetization of a permanent magnet according to Embodiment 3 of the present invention will be described with reference to the drawings.
[0051]
Conventionally, research has been made on thermal demagnetization of permanent magnets. For example, the document “Flux Conditions of Permanent Magnet Magnetic Circuit Considering Long-Term Stability”, IEEJ Transactions 1998 proposes a parameter for evaluating thermal demagnetization, and the magnitude of this parameter is related to thermal demagnetization. It is stated that it will have a big impact.
[0052]
The parameters are the residual magnetic flux density of the permanent magnet is B r [T], the coercive force is i H c [A / m], the recoil permeability is μ r [H / m], and the operation of each part of the permanent magnet When the magnetic field at the point is H [A / m] (where H is defined as the absolute value of the magnetic field at the operating point), it is expressed by the following equation (1).
[0053]
Δ = {μ r ( i H c −H)} / B r (1)
[0054]
In general, demagnetization of a permanent magnet is divided into an initial demagnetization that occurs when the permanent magnet is exposed to a certain reverse magnetic field and a certain temperature condition, and a change over time that varies depending on the time of exposure to the state. The characteristics of the parameter Δ, the initial demagnetization factor, and the change rate with time of the equation (1) can be obtained by actual measurement.
[0055]
FIG. 7 is a diagram showing the relationship between the parameter Δ and the initial demagnetization factor according to Embodiment 3 of the present invention. In FIG. 7, the vertical axis represents the initial demagnetization factor [%]. The horizontal axis represents the parameter Δ. The further away from the origin, the larger the values.
[0056]
FIG. 8 is a diagram showing the relationship between the parameter Δ and the rate of change with time according to Embodiment 3 of the present invention. In FIG. 8, the vertical axis represents the rate of change with time [%]. The horizontal axis represents the parameter Δ. The further away from the origin, the larger the values.
[0057]
7 and 8, the smaller the parameter Δ, the greater the initial demagnetization rate and the rate of change with time. In the first and second embodiments, the calculation method of the demagnetization factor or demagnetization amount has not been described, but there is a method using the characteristics of the parameter Δ, the initial demagnetization factor, and the change rate with time shown in FIGS. For example, the parameter Δ of each element or each node can be obtained in the finite element method evaluation. The demagnetizing factor of each element or each node can be calculated from this parameter Δ and the characteristics shown in FIGS.
[0058]
The method of calculating the demagnetization factor or demagnetization amount only from the initial demagnetization characteristics shown in FIG. 7 is to predict the amount and distribution of demagnetization when a large current flows due to a sudden short circuit in a motor or a generator. There is an effect that can be.
[0059]
Moreover, since the demagnetization rate or demagnetization amount after exposure for a certain period of time can be calculated from the initial demagnetization and change with time shown in FIGS. 7 and 8, the motor or generator was operated for a long time. There is an effect that the amount and distribution of subsequent demagnetization can be predicted.
[0060]
In addition, since the method using the parameter Δ in Expression (1) is a method based on actual measurement, there is an effect that the amount and distribution of demagnetization can be predicted with high accuracy.
[0061]
Furthermore, in the conventional example for redefining the BH curve, data must be prepared according to the type of permanent magnet. However, the parameter Δ used in the present embodiment is different for a rare earth magnet, for example. Since it has almost constant characteristics, there is an effect that it is possible to evaluate the demagnetization of various permanent magnets with only one kind of relationship between the parameter Δ and the initial demagnetization factor and the relationship between the parameter Δ and the change rate with time. .
[0062]
Embodiment 4 FIG.
A method for evaluating demagnetization of a permanent magnet according to Embodiment 4 of the present invention will be described with reference to the drawings. FIG. 9 is a diagram showing an example of obtaining the demagnetization factor by plotting according to the fourth embodiment of the present invention.
[0063]
In the above-described third embodiment, the method using the parameter Δ for evaluating demagnetization is used, but there is also a method for obtaining geometrically as is well known.
[0064]
In FIG. 9, the vertical axis represents the magnetic flux density B and the magnetization J. The horizontal axis represents the magnetic field H. The demagnetizing amount ΔB when the operating point is changed by receiving a reverse magnetic field of ΔH from the state at the operating point 41 and the magnetic field is removed again and returned to the original is as in the operating point 41 → 42 → 43 → 44 → 45. Can be sought.
[0065]
Since this method can calculate the initial demagnetization from the demagnetization curve, there is an effect that the demagnetization amount can be estimated relatively easily unlike the method using the parameter Δ in the equation (1).
[0066]
Embodiment 5 FIG.
A method for evaluating demagnetization of a permanent magnet according to Embodiment 5 of the present invention will be described with reference to the drawings. FIG. 10 is a diagram showing an evaluation model for a permanent magnet type motor according to Embodiment 5 of the present invention. Moreover, FIG. 11 is a figure which shows an example (demagnetization factor and its distribution) of the result of the demagnetization evaluation which concerns on Embodiment 5 of this invention.
[0067]
In FIG. 10, a permanent magnet motor evaluation model 50 includes models of respective parts of a stator core 51, an armature winding 52, a rotor shaft 54, and a permanent magnet 55.
[0068]
FIG. 11 shows the demagnetization factor obtained by conducting a magnetic field analysis using the evaluation model 50 for the permanent magnet motor and assuming that a current much larger than the rated current flows.
[0069]
As shown in FIG. 11, demagnetization was confirmed at the portion where the magnetomotive force in the reverse direction is applied to the permanent magnet 55. This result is a result of following a history for each evaluation step in which the rotor position and the current value change.
[0070]
Furthermore, the distribution of the demagnetization factor was reflected in the evaluation, and the cogging torque, the torque waveform at the time of loading, and the induced voltage waveform were evaluated. The results are shown in FIGS.
[0071]
FIG. 12 is a diagram showing an evaluation example of cogging torque according to the fifth embodiment of the present invention. In FIG. 12, the vertical axis represents cogging torque [Nm]. The horizontal axis represents the rotor position [degree (electrical angle)]. The broken line waveform represents before demagnetization, and the solid line waveform represents after demagnetization.
[0072]
Moreover, FIG. 13 is a figure which shows the example of evaluation of the torque at the time of the load which concerns on Embodiment 5 of this invention. In FIG. 13, the vertical axis represents torque [Nm]. The horizontal axis represents the rotor position [degree (electrical angle)]. The broken line waveform represents before demagnetization, and the solid line waveform represents after demagnetization.
[0073]
FIG. 14 is a diagram showing an evaluation example of the induced voltage waveform according to the fifth embodiment of the present invention. In FIG. 14, the vertical axis represents the no-load induced voltage [V]. The horizontal axis represents the rotor position [degree (electrical angle)]. The broken line waveform represents before demagnetization, and the solid line waveform represents after demagnetization.
[0074]
Even if the design reduces cogging torque and torque pulsation, it can be seen that cogging torque and torque pulsation differ if demagnetization occurs due to a large current flowing for some reason. Furthermore, it can be seen that the induced voltage waveform also changes slightly.
[0075]
Needless to say, it is possible to evaluate electromagnetic excitation force related to sound and vibration.
[0076]
Thus, when demagnetization occurs, it is possible to predict how much the motor characteristics are affected.
[0077]
By using such an evaluation method, there is an effect that the design of the permanent magnet type rotating machine can be made highly accurate. Furthermore, it is possible to predict the characteristics of the equipment when a large current flows through the permanent magnet type rotating machine or when the temperature becomes high, so these can be fed back to the design to improve the reliability of electrical design and magnetic circuit design. There is an effect that improvement can be achieved.
[0078]
【The invention's effect】
As described above, the evaluation method of demagnetization of the permanent magnet according to the present invention can know the magnetic flux density distribution when the demagnetization occurs in the permanent magnet with high accuracy, and further know the magnetic flux density distribution with high accuracy. Thus, there is an effect that the characteristics of the electric device including the permanent magnet can be known with high accuracy.
[Brief description of the drawings]
FIG. 1 is a diagram showing the configuration of a permanent magnet demagnetization evaluation system according to Embodiment 1 of the present invention;
FIG. 2 is a flowchart showing the operation of the permanent magnet demagnetization evaluation system according to Embodiment 1 of the present invention;
FIG. 3 is a diagram showing a specific example of demagnetization evaluation of the permanent magnet demagnetization evaluation system according to Embodiment 1 of the present invention;
FIG. 4 is a diagram showing a cross-sectional configuration of a permanent magnet type motor according to a second embodiment of the present invention.
FIG. 5 is a diagram showing a specific example of a method for evaluating demagnetization of a permanent magnet according to Embodiment 2 of the present invention.
FIG. 6 is a diagram showing another specific example of a method for evaluating demagnetization of a permanent magnet according to Embodiment 2 of the present invention.
FIG. 7 is a diagram showing a relationship between a parameter Δ and an initial demagnetization factor according to Embodiment 3 of the present invention.
FIG. 8 is a diagram showing a relationship between a parameter Δ and a rate of change with time according to Embodiment 3 of the present invention.
FIG. 9 is a diagram showing an example of obtaining a demagnetizing factor by plotting according to the fourth embodiment of the present invention.
FIG. 10 is a diagram showing an evaluation model for a permanent magnet type motor according to Embodiment 5 of the present invention;
FIG. 11 is a diagram showing an example of a demagnetization evaluation result (demagnetization factor and its distribution) according to Embodiment 5 of the present invention;
FIG. 12 is a diagram showing an evaluation example of cogging torque according to the fifth embodiment of the present invention.
FIG. 13 is a diagram showing an evaluation example of torque at the time of loading according to Embodiment 5 of the present invention.
FIG. 14 is a diagram showing an evaluation example of an induced voltage waveform according to the fifth embodiment of the present invention.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 1 Recording device, 2 Demagnetization evaluation means, 3 Memory | storage device, 4 Electric equipment characteristic evaluation means, 10 Permanent magnet, 20 Permanent magnet type motor, 21 Stator iron core, 22 Armature winding, 23 slots, 24 Rotor shaft, 25 Permanent magnet, 50 Evaluation model for permanent magnet motor.

Claims (5)

永久磁石型モータの電子機巻線の電流、及び回転子と固定子の位置関係が変化して、前記永久磁石型モータの永久磁石が所定の逆磁界、所定の温度にさらされる各時点において、有限要素法または積分要素法を用いた磁界解析により、前記永久磁石をメッシュ分割し、メッシュの各要素または各節点の第1の減磁率あるいは第1の減磁量を求めるステップと、
履歴である、前記各時点で求めたメッシュの各要素または各節点の第1の減磁率あるいは第1の減磁量のうち、メッシュの各要素または各節点の第1の減磁率あるいは第1の減磁量の最大値を、メッシュの各要素または各節点における第2の減磁率あるいは第2の減磁量として求めるステップと、
メッシュの各要素または各節点における第2の減磁率あるいは第2の減磁量に応じて、前記永久磁石の残留磁束密度、保磁力、減磁曲線のうち少なくとも1つを求めるステップと、
メッシュの各要素または各節点における第2の減磁率あるいは第2の減磁量に基づき、前記永久磁石型モータのコギングトルク、負荷時トルク、誘起電圧の評価例のうち少なくとも1つを求めるステップと
を含むことを特徴とする永久磁石の減磁の評価方法。
At each time when the current of the winding of the motor of the permanent magnet type motor and the positional relationship between the rotor and the stator are changed, the permanent magnet of the permanent magnet type motor is exposed to a predetermined reverse magnetic field and a predetermined temperature. Dividing the permanent magnet into meshes by magnetic field analysis using a finite element method or an integral element method, and obtaining a first demagnetizing factor or a first demagnetizing amount of each element or each node of the mesh;
Of the first demagnetization factor or the first demagnetization amount of each element or each node of the mesh obtained at each time point, which is a history, the first demagnetization factor or the first demagnetization factor of each element or each node of the mesh Obtaining a maximum value of demagnetization amount as a second demagnetization factor or a second demagnetization amount at each element or node of the mesh;
Determining at least one of residual magnetic flux density, coercivity, and demagnetization curve of the permanent magnet according to a second demagnetization factor or a second demagnetization amount at each element or each node of the mesh;
Obtaining at least one of evaluation examples of cogging torque, on-load torque, and induced voltage of the permanent magnet type motor based on a second demagnetization factor or a second demagnetization amount at each element or each node of the mesh; and A method for evaluating the demagnetization of a permanent magnet.
前記第1の減磁率あるいは前記第1の減磁量は、前記永久磁石の残留磁束密度をB 、保磁力を 、リコイル透磁率をμ 、メッシュの各要素または各節点における磁界の絶対値をHとしたとき、
Δ={μ −H)}/B
で表されるパラメータΔと、前記パラメータΔに対する初期減磁率とから求める
ことを特徴とする請求項1記載の永久磁石の減磁の評価方法。
The first demagnetization factor or the first demagnetization amount is defined as follows: the residual magnetic flux density of the permanent magnet is B r , the coercive force is i H c , the recoil permeability is μ r , each element of the mesh, or the magnetic field at each node. When the absolute value of is H,
Δ = {μ r ( i H c −H)} / B r
2. The method for evaluating demagnetization of a permanent magnet according to claim 1 , wherein the demagnetization is obtained from a parameter Δ represented by:
前記第1の減磁率あるいは前記第1の減磁量は、前記永久磁石の残留磁束密度をB 、保磁力を 、リコイル透磁率をμ 、メッシュの各要素または各節点における磁界の絶対値をHとしたとき、
Δ={μ −H)}/B
で表されるパラメータΔと、前記パラメータΔに対する初期減磁率と、前記パラメータΔに対する経時変化率とから求める
ことを特徴とする請求項記載の永久磁石の減磁の評価方法。
The first demagnetization factor or the first demagnetization amount is defined as follows: the residual magnetic flux density of the permanent magnet is B r , the coercive force is i H c , the recoil permeability is μ r , each element of the mesh, or the magnetic field at each node. When the absolute value of is H,
Δ = {μ r ( i H c −H)} / B r
A parameter delta represented in an initial flux loss for said parameter delta, the evaluation method of the demagnetization of the permanent magnet of claim 1, wherein the determining from the temporal change rate for the parameter delta.
コンピュータに
永久磁石型モータの電子機巻線の電流、及び回転子と固定子の位置関係が変化して、前記永久磁石型モータの永久磁石が所定の逆磁界、所定の温度にさらされる各時点において、有限要素法または積分要素法を用いた磁界解析により、前記永久磁石をメッシュ分割し、メッシュの各要素または各節点の第1の減磁率あるいは第1の減磁量を求める手順、
履歴である、前記各時点で求めたメッシュの各要素または各節点の第1の減磁率あるいは第1の減磁量のうち、メッシュの各要素または各節点の第1の減磁率あるいは第1の減磁量の最大値を、メッシュの各要素または各節点における第2の減磁率あるいは第2の減磁量として求める手順、
メッシュの各要素または各節点における第2の減磁率あるいは第2の減磁量に応じて、前記永久磁石の残留磁束密度、保磁力、減磁曲線のうち少なくとも1つを求める手順、
及び
メッシュの各要素または各節点における第2の減磁率あるいは第2の減磁量に基づき、前記永久磁石型モータのコギングトルク、負荷時トルク、誘起電圧の評価例のうち少なくとも1つを求める手順
を実行させるための永久磁石の減磁の評価プログラム
On the computer
At each time when the current of the winding of the motor of the permanent magnet type motor and the positional relationship between the rotor and the stator are changed, the permanent magnet of the permanent magnet type motor is exposed to a predetermined reverse magnetic field and a predetermined temperature. A step of dividing the permanent magnet into meshes by a magnetic field analysis using a finite element method or an integral element method, and obtaining a first demagnetizing factor or a first demagnetizing amount of each element or each node of the mesh;
Of the first demagnetization factor or the first demagnetization amount of each element or each node of the mesh obtained at each time point, which is a history, the first demagnetization factor or the first demagnetization factor of each element or each node of the mesh A procedure for obtaining the maximum value of the demagnetization amount as the second demagnetization factor or the second demagnetization amount at each element or each node of the mesh;
A procedure for obtaining at least one of a residual magnetic flux density, a coercive force, and a demagnetization curve of the permanent magnet according to a second demagnetization factor or a second demagnetization amount at each element or each node of the mesh;
as well as
Procedure for obtaining at least one of evaluation examples of cogging torque, load torque, and induced voltage of the permanent magnet type motor based on the second demagnetization factor or the second demagnetization amount at each element or each node of the mesh
An evaluation program for permanent magnet demagnetization .
永久磁石型モータの電子機巻線の電流、及び回転子と固定子の位置関係が変化して、前記永久磁石型モータの永久磁石が所定の逆磁界、所定の温度にさらされる各時点において、有限要素法または積分要素法を用いた磁界解析により、前記永久磁石をメッシュ分割し、メッシュの各要素または各節点の第1の減磁率あるいは第1の減磁量を求め、
履歴である、前記各時点で求めたメッシュの各要素または各節点の第1の減磁率あるいは第1の減磁量のうち、メッシュの各要素または各節点の第1の減磁率あるいは第1の減磁量の最大値を、メッシュの各要素または各節点における第2の減磁率あるいは第2の減磁量として求め、
メッシュの各要素または各節点における第2の減磁率あるいは第2の減磁量に応じて、前記永久磁石の残留磁束密度、保磁力、減磁曲線のうち少なくとも1つを求める減磁評価手段と、
前記減磁評価手段によって求めたメッシュの各要素または各節点における第2の減磁率あるいは第2の減磁量に基づき、前記永久磁石型モータのコギングトルク、負荷時トルク、誘起電圧の評価例のうち少なくとも1つを求める電機機器特性評価手段と
を備えたことを特徴とする永久磁石の減磁の評価システム
At each time when the current of the winding of the motor of the permanent magnet type motor and the positional relationship between the rotor and the stator are changed, the permanent magnet of the permanent magnet type motor is exposed to a predetermined reverse magnetic field and a predetermined temperature. The permanent magnet is divided into meshes by magnetic field analysis using a finite element method or an integral element method, and a first demagnetization factor or a first demagnetization amount of each element or each node of the mesh is obtained.
Of the first demagnetization factor or the first demagnetization amount of each element or each node of the mesh obtained at each time point, which is a history, the first demagnetization factor or the first demagnetization factor of each element or each node of the mesh The maximum value of the demagnetization amount is obtained as the second demagnetization factor or the second demagnetization amount at each element or each node of the mesh,
Demagnetization evaluation means for obtaining at least one of a residual magnetic flux density, a coercive force, and a demagnetization curve of the permanent magnet according to a second demagnetization factor or a second demagnetization amount at each element or each node of the mesh; ,
Based on the second demagnetization factor or the second demagnetization amount at each element or each node of the mesh obtained by the demagnetization evaluation means, an evaluation example of cogging torque, load torque, and induced voltage of the permanent magnet motor Electrical equipment characteristic evaluation means for obtaining at least one of them
A system for evaluating demagnetization of a permanent magnet .
JP2003049453A 2003-02-26 2003-02-26 Evaluation method, evaluation program and evaluation system for demagnetization of permanent magnet Expired - Fee Related JP4275426B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2003049453A JP4275426B2 (en) 2003-02-26 2003-02-26 Evaluation method, evaluation program and evaluation system for demagnetization of permanent magnet

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003049453A JP4275426B2 (en) 2003-02-26 2003-02-26 Evaluation method, evaluation program and evaluation system for demagnetization of permanent magnet

Publications (2)

Publication Number Publication Date
JP2004257879A JP2004257879A (en) 2004-09-16
JP4275426B2 true JP4275426B2 (en) 2009-06-10

Family

ID=33115168

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003049453A Expired - Fee Related JP4275426B2 (en) 2003-02-26 2003-02-26 Evaluation method, evaluation program and evaluation system for demagnetization of permanent magnet

Country Status (1)

Country Link
JP (1) JP4275426B2 (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101019123B1 (en) * 2008-10-24 2011-03-07 고려대학교 산학협력단 Apparatus and method for monitoring permanent magnet demagnetization of permanent magnet synchronous motors, and a medium having computer readable program for executing the method
EP2535907B1 (en) 2010-02-10 2017-04-05 Hitachi Metals, Ltd. Magnetic force characteristics computing method, magnetic force characteristics computing device, and computer program
JP5175312B2 (en) * 2010-03-19 2013-04-03 高麗大學校産學協力團 Demagnetization state diagnosis device for permanent magnet of permanent magnet synchronous motor, diagnosis method, recording medium recording program for executing the diagnosis method, and drive device for permanent magnet synchronous motor
EP2518456A1 (en) * 2011-04-29 2012-10-31 ABB Technology AG Method for monitoring demagnetization
US9547051B2 (en) 2011-05-17 2017-01-17 Hitachi Metals, Ltd. Calculating method of magnetic force characteristic, and magnetic force characteristic computing device
JP5994269B2 (en) * 2012-02-09 2016-09-21 マツダ株式会社 Motor torque setting method and setting device
KR101912495B1 (en) * 2013-01-02 2018-10-26 트레인 인터내셔날 인코포레이티드 Permanent magnet motor degradation diagnostics system
CN103986397B (en) * 2014-05-19 2016-08-24 安徽中家智锐科技有限公司 Brshless DC motor permanent magnet fault detection method
JP7283967B2 (en) * 2019-04-25 2023-05-30 株式会社日立製作所 Permanent magnet demagnetization tolerance analysis method, calculation device, calculation system

Also Published As

Publication number Publication date
JP2004257879A (en) 2004-09-16

Similar Documents

Publication Publication Date Title
Ruoho et al. Comparison of demagnetization models for finite-element analysis of permanent-magnet synchronous machines
Tariq et al. Iron and magnet losses and torque calculation of interior permanent magnet synchronous machines using magnetic equivalent circuit
Ajily et al. Magnet defect and rotor eccentricity modeling in axial-flux permanent-magnet machines via 3-D field reconstruction method
Williamson et al. Finite-element models for cage induction motor analysis
Roshen Iron loss model for permanent-magnet synchronous motors
Nasiri-Zarandi et al. Analysis, optimization, and prototyping of a brushless DC limited-angle torque-motor with segmented rotor pole tip structure
Gao et al. Robust design for torque optimization using response surface methodology
JP4275426B2 (en) Evaluation method, evaluation program and evaluation system for demagnetization of permanent magnet
Lee et al. Permanent magnet demagnetization characteristic analysis of a variable flux memory motor using coupled Preisach modeling and FEM
JP6213407B2 (en) Method and apparatus for estimating motor torque
Zhu et al. Cogging torque minimisation in FSPM machines by right‐angle‐based tooth chamfering technique
Jung et al. A novel approach for 2-D electromagnetic field analysis of surface mounted permanent magnet synchronous motor taking into account axial end leakage flux
Taran et al. A comparative study of methods for calculating AC winding losses in permanent magnet machines
Golzarzadeh et al. Analytical modelling of the linear switched reluctance motor with segmental translator
Kim et al. Core loss analysis of permanent magnet synchronous generator with slotless stator
Egorov et al. Model-based hysteresis loss assessment in PMSMs with ferrite magnets
Gotovac et al. Analytical model of permeance variation losses in permanent magnets of the multipole synchronous machine
Winterborne et al. Estimation of AC copper loss in electrical machine windings with consideration of end effects
Duan et al. Method for multi-objective optimized designs of surface mount permanent magnet motors with concentrated or distributed stator windings
Wang et al. Influence of the permanent magnet overhang on the performance of the brushless dc motor
Li et al. Braking force characteristics of permanent magnet eddy current brake with high acceleration
Li et al. Parametric study for the design of the end region of large synchronous generators based on three-dimensional transient finite element analysis
Gholizad et al. Dynamic analysis of highly saturated switched reluctance motors using coupled magnetic equivalent circuit and the analytical solution
Gómez et al. Methodology to study demagnetization risk in permanent magnet machines by finite element method
Cao et al. Magnetic field analytical solution and electromagnetic force calculation of coreless stator axial-flux permanent-magnet machines

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040617

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20051227

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20081202

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090129

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20090303

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20090304

R150 Certificate of patent or registration of utility model

Ref document number: 4275426

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120313

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120313

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130313

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130313

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140313

Year of fee payment: 5

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees