JP2019204269A - 電磁場解析システム、電磁場解析方法、およびプログラム - Google Patents
電磁場解析システム、電磁場解析方法、およびプログラム Download PDFInfo
- Publication number
- JP2019204269A JP2019204269A JP2018098756A JP2018098756A JP2019204269A JP 2019204269 A JP2019204269 A JP 2019204269A JP 2018098756 A JP2018098756 A JP 2018098756A JP 2018098756 A JP2018098756 A JP 2018098756A JP 2019204269 A JP2019204269 A JP 2019204269A
- Authority
- JP
- Japan
- Prior art keywords
- flux density
- magnetic flux
- iron loss
- derived
- relationship
- 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.)
- Granted
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 147
- 230000005672 electromagnetic field Effects 0.000 title claims abstract description 57
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims abstract description 344
- 229910052742 iron Inorganic materials 0.000 claims abstract description 172
- 230000004907 flux Effects 0.000 claims abstract description 146
- 230000005284 excitation Effects 0.000 claims abstract description 71
- 238000012360 testing method Methods 0.000 claims abstract description 64
- 239000000463 material Substances 0.000 claims abstract description 40
- 239000000696 magnetic material Substances 0.000 claims abstract description 29
- 238000005096 rolling process Methods 0.000 claims description 8
- 230000005415 magnetization Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 abstract description 21
- 238000000034 method Methods 0.000 description 34
- 229910000831 Steel Inorganic materials 0.000 description 17
- 239000010959 steel Substances 0.000 description 17
- 230000000052 comparative effect Effects 0.000 description 8
- 238000004141 dimensional analysis Methods 0.000 description 8
- 238000012986 modification Methods 0.000 description 8
- 230000004048 modification Effects 0.000 description 8
- 238000009795 derivation Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 229910000976 Electrical steel Inorganic materials 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000007429 general method Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000012916 structural analysis Methods 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 238000004870 electrical engineering Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010030 laminating Methods 0.000 description 1
- 238000003475 lamination Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/60—Other road transportation technologies with climate change mitigation effect
- Y02T10/64—Electric machine technologies in electromobility
Landscapes
- Measuring Magnetic Variables (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Abstract
Description
例えば、高効率のモータを実現するためには、モータの主要な損失である鉄損と銅損とを低減させることが、その具体的な方策となる。定常運転時のモータの損失には全損失に対する鉄損の占める割合が大きいことから、鉄損の低減が、モータ効率の改善に有効である。
特許文献1には、以下の技術が開示されている。まず、材料属性と励磁基本周波数とに対応する渦電流損係数特性(異常渦電流損係数と磁束密度の最大値との関係)を用意する。次に、コア(鉄心)の磁束密度ベクトルを数値解析(静磁場解析)により導出し、導出した磁束密度ベクトルから磁束密度の大きさの最大値を導出する。次に、このようにして導出した磁束密度の大きさに対応する異常渦電流損係数を用いて、コアを構成する電磁鋼板の導電率を補正した補正導電率を導出する。そして、このようにして導出した補正導電率を用いて、コアの磁束密度ベクトルおよび渦電流ベクトルを数値解析(動磁場解析(静磁場項に加えて渦電流項を考慮した解析))により導出し、導出した磁束密度ベクトルおよび渦電流ベクトルを用いてヒステリシス損と渦電流損とをそれぞれ導出する。このようにして導出したヒステリシス損と渦電流損の和がコアの鉄損になる。
インバータ電源により駆動するモータのコアにうず電流を発生させて鉄損を増加させているのは、モータに印加される励磁電圧の時間波形に含まれている高調波成分である。この高調波成分は、インバータ電源を使用することにより基本波成分に重畳されるものである。以下の説明では、このようにインバータ電源を用いることに起因して発生する高調波を必要に応じてインバータ高調波(図1のグラフ13、14の鋸歯状の部分を参照)と称する。また、純粋な(理想的な)正弦波を必要に応じて純正弦波(図1のグラフ11を参照)と称する。また、純正弦波の励磁電圧で磁性材料を励磁した場合、磁性材料における磁束密度の時間波形は歪む。このようにして歪んだ波形を必要に応じて歪み波(図1のグラフ12を参照)と称する。歪み波に含まれる高調波は、インバータ高調波に比べて低次の高調波である。
本発明者は、波高値が同じBmaxであれば、グラフ11に示すような、時間波形が純正弦波である磁束密度となるように励磁された磁性材料の鉄損と、グラフ12に示すような、時間波形が歪み波である磁束密度となるように励磁された磁性材料の鉄損と、グラフ13に示すような、時間波形が純正弦波にインバータ高調波が重畳された時間波形である磁束密度となるように励磁された磁性材料の鉄損と、グラフ14に示すような、時間波形が歪み波にインバータ高調波が重畳された時間波形である磁束密度となるように励磁された磁性材料の鉄損は、相互に相関関係があるという知見を得た。
以下、図面を参照しながら、本発明の一実施形態を説明する。本実施形態では、モータのコア(ステータコア)の鉄損を推定する場合(推定対象が、モータのコアの鉄損)を例に挙げて説明する。また、コアは、複数の電磁鋼板を積層することにより構成されるものとする。以下の説明では、このコアを必要に応じてモータコアと称する。
[磁気特性試験装置100]
磁気特性試験装置100は、モータコアを構成する電磁鋼板の試験材の鉄損を測定するための装置である。以下の説明では、このモータコアを構成する電磁鋼板の試験材を必要に応じて試験材と略称する。試験材の厚みおよび種類は、モータコアを構成する電磁鋼板の厚みおよび種類と同じである。
尚、図3では、変調波・搬送波の値を相対値で示す。
尚、図3では、励磁電圧360の値を相対値で示す。
以上のようなPWMインバータにおける変調動作(パルス幅変調の動作)を定めるパラメータには、変調率mとキャリア周波数でfcとが含まれる。PWMインバータでは、変調率mは、変調波(正相変調波310および負相変調波330)の振幅E0を、搬送波320の振幅Ecで割った値(m=E0÷Ec)で表される。キャリア周波数でfcは、搬送波320の周波数である。
例えば、非特許文献1には、鉄損を測定する際には磁束密度に対応する二次電圧の波形率が1.10[%]〜1.12[%]になるようにすることが記載されている。そこで、例えば、PWMインバータで試験材を励磁することにより磁気特性試験装置100で測定される二次電圧の基本波の波形率が1.10[%]〜1.12[%]になる状態にすることにより、試験材における磁束密度の時間波形が、純正弦波にインバータ高調波が重畳された時間波形となる状態にすることができる。
また、以下の説明では、磁束密度の波高値Bmax、磁化方向と磁化容易軸(圧延方向)とのなす角度θB、鉄損測定時に試験材に与えられる応力σを、それぞれ、磁束密度波高値Bmax、角度θB、応力σと略称する。
前述したようにして磁気特性試験装置100の測定を行うことにより、図3に示すように、磁束密度波高値Bmaxと鉄損wとの関係が、励磁基本周波数fs毎、応力σ毎、角度θB毎、キャリア周波数fc毎、および変調率m毎に得られる。以下の説明では、励磁基本周波数fs毎、応力σ毎、角度θB毎、キャリア周波数fc毎、および変調率m毎の、磁束密度波高値Bmaxと鉄損wとの関係を、必要に応じてBm−w特性と称する。
図2に、電磁場解析装置200の機能的な構成の一例を示す。電磁場解析装置200のハードウェアは、例えば、CPU、ROM、RAM、HDD、および各種のインターフェースを用いることにより実現される。
<データ入力部201>
データ入力部201は、Bm−w特性のデータを入力する。Bm−w特性のデータは、Bm−w特性を作成するために必要なデータであり、磁気特性試験装置100による試験材の鉄損の測定結果から得られる。データ入力部201は、例えば、外部装置から、Bm−w特性のデータを受信したり、電磁場解析装置200のユーザインターフェースの操作からBm−w特性のデータの内容を特定したり、可搬型の記憶媒体記憶されたBm−w特性のデータを読み出したりすることにより、Bm−w特性のデータを入力することができる。
鉄損データベース作成部202は、データ入力部201により入力されたBm−w特性のデータに基づいて、Bm−w特性を作成する。鉄損データベース作成部202は、Bm−w特性のデータに対し回帰分析等の手法を用いることにより、Bm−w特性を表す(少なくとも1つの)関数を導出し、当該関数をBm−w特性とすることができる。また、Bm−w特性励磁基本周波数fs毎、応力σ毎、角度θB毎、キャリア周波数fc毎、および変調率m毎に、磁束密度波高値Bmaxと、当該磁束密度波高値Bmaxに対応する鉄損wとを格納するテーブルを、Bm−w特性として作成してもよい。
鉄損データベース作成部202は、以上のようにして作成したBm−w特性を鉄損データベース記憶部203に記憶する。
解析条件入力部204は、モータコアの鉄損を導出する計算を行う際に必要なデータを入力する。以下の説明では、このデータを必要に応じて解析条件データと称する。解析条件データには、磁気解析用データと、鉄損用データと、応力解析用データとが含まれる。
応力解析部205は、非線形非定常2次元有限要素法を用いた構造解析を行うことにより、解析条件データ(応力解析用データ)に従って荷重が付加された場合のモータコアの応力(のx軸方向の成分およびy軸方向の成分)σx、σyを微小領域毎に導出する。ここでは、2次元解析(2次元の数値解析)を行うので、応力のz軸方向の成分σzについては導出しない。尚、x軸、y軸、z軸は、3次元直交座標の各軸であり、磁束密度の計算の際にも、応力の計算の際にも、共通の3次元直交座標を用いるものとする。また、z軸の方向は、例えば、電磁鋼板の積層方向である。
電磁場解析部206は、非線形非定常2次元有限要素法を用いた電磁場解析を行うことにより、微小領域(メッシュ)のそれぞれにおいて、解析条件入力部204により入力された解析条件データ(磁気解析用データ)に従って目標励磁電圧(目標駆動電圧)で励磁された場合のモータコアの磁束密度(のx軸方向の成分およびy軸方向の成分)Bx、Byを導出する。ここでは、2次元解析(2次元の数値解析)を行うので、磁束密度のz軸方向の成分Bzについては導出しない。尚、目標励磁電圧(目標駆動電圧)にはインバータ高調波は含まれない(2次元解析時にはインバータ高調波はないものとされる)。
図5は、電磁場解析部206により導出された、或る微小領域における2次元磁束密度ベクトルの時間波形の一例を示す図である。
微小領域内パラメータ導出部207は、微小領域の2次元磁束密度ベクトルの時間波形の波高値Bmaxを、微小領域のそれぞれについて導出する。図5に示す例では、時刻t1において微小領域の2次元磁束密度ベクトルの時間波形の波高値Bmaxが得られることを示す。以下の説明では、2次元磁束密度ベクトルの時間波形の波高値Bmaxを、必要に応じて磁束密度波高値Bmaxと称する。
また、微小領域内パラメータ導出部207は、微小領域の2次元応力ベクトルの大きさを応力σとして導出する。
微小領域内鉄損導出部208は、微小領域内パラメータ導出部207により導出された、各微小領域における角度θBおよび応力σと、解析条件入力部204により入力された解析条件データ(鉄損用データ)に含まれる、励磁基本周波数fs、キャリア周波数fc、および変調率mとに対応するBm−w特性を、鉄損データベース記憶部203から読み出す。そして、微小領域内鉄損導出部208は、読み出したBm−w特性から、微小領域内パラメータ導出部207により導出された磁束密度波高値Bmaxに対応する鉄損wを導出することを、全ての微小領域のそれぞれについて行う。
鉄損総和部209は、微小領域内鉄損導出部208により導出された各微小領域における鉄損wの総和をとる(積算する)ことにより、モータコア全体の鉄損Wを導出する。
<出力部210>
出力部210は、鉄損総和部209により導出されたモータコア全体の鉄損Wの情報を出力する。出力の形態としては、例えば、コンピュータディスプレイへの表示、電磁場解析装置200の内部または外部の記憶媒体への記憶、および外部装置への送信の少なくとも何れか1つを採用することができる。
図6のフローチャートを参照しながら、本実施形態の電磁場解析装置200による電磁場解析方法の一例を説明する。尚、ここでは、鉄損データベース記憶部203にBm−w特性が既に記憶されているものとする。
ステップS601において、解析条件入力部204は、解析条件データを入力する。
次に、ステップS602において、応力解析部205は、構造解析を行うことにより、解析条件データ(応力解析用データ)に従って荷重が付加された場合のモータコアの各微小領域における2次元応力ベクトルを導出する。尚、モータコアに荷重が付加されない場合には、ステップS602の処理は省略される。
次に、ステップS604において、微小領域内パラメータ導出部207は、2次元磁束密度ベクトルおよび2次元応力ベクトルを用いて、モータコアの各微小領域におけるパラメータ(磁束密度波高値Bmax、角度θB、および応力σ)を導出する。
次に、ステップS606において、鉄損総和部209は、微小領域における鉄損wの総和をとる(積算する)ことにより、モータコア全体の鉄損Wを導出する。
最後に、ステップS607において、出力部210は、鉄損総和部209により導出されたモータコア全体の鉄損Wの情報を出力する。
次に、実施例を説明する。本実施例では、同一のモータを同一の条件で動作させた場合のモータコアの鉄損を、本実施形態で説明した手法と、特許文献1に記載の手法とのそれぞれで導出した。本実施形態で説明した手法を発明例と表記し、特許文献1に記載の手法を比較例と表記する。本実施例では、モータに付加される応力σは0(ゼロ)とした。
図8に示すように、発明例では、比較例に比べ、微小領域の数(要素数)を大幅に削減することができ、波形一周期分の計算における時間分割の数も大幅に削減することができることが分かる。従って、メモリ容量を削減することができる。また、発明例では、比較例に比べ、計算時間を大幅に削減することができる(変調率mが0.45の場合には、2810分の1、変調率mが0.6の場合は1520分の1になった)。そして、このようにメモリ容量および計算時間を削減しても、厳密な3次元解析(3次元空間の全ての方向に電気伝導が生じる得るものとして行う3次元の数値解析)を行った比較例1、2と同程度の鉄損を得ることができた。具体的に、変調率mが0.45の場合には、比較例1に対する発明例の鉄損の誤差率は−8.4[%]、変調率mが0.6の場合には、比較例2に対する発明例の鉄損の誤差率は−2.3[%]になった。
以上のように本実施形態では、PWMインバータで励磁した場合の試験材の鉄損を磁気特性試験装置100で測定する。電磁場解析装置200は、その結果に基づいて、Bm−w特性(励磁基本周波数fs毎、応力σ毎、角度θB毎、キャリア周波数fc毎、および変調率m毎の、磁束密度波高値Bmaxと鉄損wとの関係)を作成して記憶する。Bm−w特性を記憶した後、電磁場解析装置200は、時間波形が純正弦波である励磁電圧(駆動電圧)をモータに印加した場合の、モータコアの各微小領域における磁束密度波高値Bmaxを2次元の数値解析により導出する。そして、電磁場解析装置200は、導出した磁束密度波高値Bmaxに対応する鉄損wを、Bm−w特性から導出し、各微小領域における鉄損wの総和をモータコア全体の鉄損Wとして導出する。
<変形例1>
本実施形態では、インバータ電源としてPWMインバータを用いる場合を例に挙げて説明した。しかしながら、インバータ電源は、PWMインバータに限定されない。インバータ電源の変調動作を定めるパラメータ(本実施形態では変調率mとキャリア周波数でfc)は、インバータ電源における変調方式に基づいて定められる。例えば、PAM(pulse Amplitude Modulation)インバータを用いる場合には、インバータ部に供給される直流電圧とモータへの出力電圧の比が変調動作を定めるパラメータに含まれる。
本実施形態では、非線形非定常2次元有限要素法を用いる場合を例に挙げて説明した。しかしながら、磁束密度の導出対象となる2次元平面を定める方向(3次元直交座標で3次元空間を表す場合には、例えば、x軸方向およびy軸方向)以外の方向に電気伝導が生じないものとして3次元解析を行い、3次元空間を定める方向(3次元直交座標で3次元空間を表す場合には、x軸方向、y軸方向、およびz軸方向)の成分の磁束密度Bx、By、Bzを導出してもよい。尚、磁束密度の導出対象となる2次元平面を定める方向以外の方向に電気伝導が生じないようにすることは、当該電気伝導が生じない方向における導電率σ´が0(ゼロ)になることに対応する。このようにしても、3次元空間の全ての方向において電気伝導が生じ得る(導電率が3次元空間の全ての方向で分布を持つ)ものとして3次元解析を行う場合に比べて計算負荷を軽減することができる。尚、このようにする場合であって、磁性体板を積み重ねた物体を数値解析する場合には、微小領域を磁性体板の厚みよりも大きくする。
本実施形態では、磁気特性試験装置100は、グラフ13に示すように、試験材における磁束密度の時間波形が、純正弦波にインバータ高調波が重畳された時間波形となる状態で試験材の鉄損を測定する場合を例に挙げて説明した。しかしながら、図1を参照しながら説明したように、必ずしもこのようにする必要はない。例えば、磁気特性試験装置100は、励磁電圧の時間波形が、純正弦波にインバータ高調波が重畳された時間波形(目標励磁電圧の時間波形が純正弦波)となる状態で試験材の鉄損を測定してもよい。
本実施形態では、磁気特性試験装置100による測定を行う場合を例に挙げて説明した。このようにすれば、モータ高調波が重畳された励磁電圧をモータコアに印加した場合の数値解析を行わないので、計算負荷をより軽減することができるので好ましい。しかしながら、必ずしもこのようにする必要はない。例えば、(3次元空間の全ての方向において電気伝導が生じ得るものとして行う)3次元解析の結果に基づいてBm−w特性を作成してもよい。試験材の形状は、モータコアの形状に比べて単純である。また、鉄損の導出対象の物体(本実施形態ではモータコア)の形状のみが異なる場合には、同じBm−w特性を流用することができる。従って、非特許文献1に記載の技術のように、複雑な形状の物体に対して、その都度3次元解析を行う場合に比べ、計算負荷を軽減することができる。
本実施形態では、磁束密度の時間波形を代表する値である磁束密度の代表値として、磁束密度最大値Bmaxを導出する。しかしながら、Bm−w特性における磁束密度(本実施形態では磁束密度最大値Bmax)と数値解析により導出された磁束密度(本実施形態では磁束密度最大値Bmax)とが相互に対応していればよいので、磁束密度代表値は、磁束密度最大値Bmaxでなくてもよい。例えば、磁束密度代表値として、磁束密度の一周期における大きさの平均値または実効値を用いてもよい。
<変形例6>
本実施形態では、目標励磁電圧の時間波形が純正弦波である場合を例に挙げて説明した。しかしながら、目標励磁電圧の時間波形は純正弦波でなくてもよい。
また、以上説明した本発明の実施形態は、何れも本発明を実施するにあたっての具体化の例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
以下に、請求項の記載と実施形態の記載との関係の一例を説明する。尚、請求項の記載が実施形態の記載に限定されないことは前述した通りである。
関係取得手段は、例えば、磁気特性試験装置100、データ入力部201、および鉄損データベース作成部202を用いることにより実現される。
鉄損の推定対象である物体に使用される磁性材料の試験材の磁束密度の大きさの代表値と鉄損との関係は、例えば、Bm−W特性を用いることにより実現される。
磁束密度導出手段は、例えば、電磁場解析部206および微小領域内パラメータ導出部207を用いることにより実現される。
マックスウェルの方程式に基づく数値解析は、例えば、(1)式〜(4)式に基づく有限要素法を用いた数値解析により実現される。
鉄損導出手段は、例えば、微小領域内鉄損導出部208および鉄損総和部209を用いることにより実現される。
インバータ電源は、例えば、PWMインバータを用いることにより実現される。
前記インバータ電源により前記物体を励磁する際の目標励磁電圧と同じ電圧であって、前記インバータ電源を用いることに起因して発生する高調波成分を含まない電圧は、例えば、電磁場解析部206が、モータコアの2次元磁束密度ベクトルを導出する際の励磁条件として、純正弦波の電圧(目標励磁電圧(目標駆動電圧))を用いることにより実現される。
前記磁束密度の導出対象となる2次元平面を定める方向以外の方向に電気伝導が生じない条件下で前記数値解析を行うことは、例えば、電磁場解析部206が、二次元磁束密度ベクトルを導出することと、<変形例2>により実現される。
前記インバータ電源における変調動作を定めるパラメータは、例えば、変調率mおよびキャリア周波数fcを用いることにより実現される。
励磁周波数の基本周波数は、例えば、励磁基本周波数fsを用いることにより実現される。
前記試験材の磁化方向と、前記試験材の圧延方向とのなす角度は、例えば、角度θBを用いることにより実現される。
前記試験材を励磁する際に前記試験材に与えられる応力は、例えば、応力σを用いることにより実現される。
磁気特性試験手段は、例えば、磁気特性試験装置100を用いることにより実現される。
Claims (12)
- 鉄損の推定対象である物体に使用される磁性材料の試験材の磁束密度の大きさの代表値と鉄損との関係を取得する関係取得手段と、
前記物体に使用される磁性材料の磁束密度の大きさの代表値を、マックスウェルの方程式に基づく数値解析により、複数の微小領域毎に導出する磁束密度導出手段と、
前記磁束密度導出手段により導出された磁束密度の大きさの代表値に対応する鉄損を、前記複数の微小領域毎の鉄損として、前記関係取得手段により取得された前記関係に基づいて導出し、導出した前記複数の微小領域毎の鉄損に基づいて、前記物体の鉄損を導出する鉄損導出手段と、を有することを特徴とする電磁場解析システム。 - 前記関係取得手段は、前記試験材をインバータ電源により励磁した場合の当該試験材における磁束密度の大きさの代表値と鉄損との関係を取得し、
前記磁束密度導出手段は、前記インバータ電源により前記物体を励磁する際の目標励磁電圧と同じ電圧であって、前記インバータ電源を用いることに起因して発生する高調波成分を含まない電圧で励磁した場合の前記物体における磁束密度の大きさの代表値を、マックスウェルの方程式に基づく数値解析により、複数の微小領域毎に導出し、
前記磁束密度導出手段は、前記磁束密度の導出対象となる2次元平面を定める方向以外の方向に電気伝導が生じない条件下で前記数値解析を行うことを特徴とする請求項1に記載の電磁場解析システム。 - 前記磁束密度導出手段は、マックスウェルの方程式に基づく数値解析により、前記磁束密度の導出対象となる2次元平面を定める方向の成分の値のみをもつ2次元磁束密度ベクトルを導出し、導出した2次元磁束密度ベクトルに基づいて、前記磁性材料の磁束密度の大きさの代表値を導出することを特徴とする請求項2に記載の電磁場解析システム。
- 前記関係取得手段は、前記インバータ電源における変調動作を定めるパラメータの少なくとも1つ毎に、前記関係を導出し、
前記鉄損導出手段は、前記物体を励磁する際の前記パラメータに対応する前記関係に基づいて、前記磁束密度導出手段により導出された磁束密度の大きさの代表値に対応する鉄損を、前記複数の微小領域毎の鉄損として導出することを特徴とする請求項2または3に記載の電磁場解析システム。 - 前記インバータは、パルス幅変調を行い、
前記インバータにおける変調動作を定めるパラメータは、変調率とキャリア周波数とを含むことを特徴とする請求項4に記載の電磁場解析システム。 - 前記関係取得手段は、前記試験材を励磁する際の励磁周波数の基本周波数ごとに、前記関係を導出し、
前記鉄損導出手段は、前記物体を励磁する際の励磁周波数の基本周波数に対応する前記関係に基づいて、前記磁束密度導出手段により導出された磁束密度の大きさの代表値に対応する鉄損を、前記複数の微小領域毎の鉄損として導出することを特徴とする請求項1〜5の何れか1項に記載の電磁場解析システム。 - 前記磁性材料は、磁性体板を含み、
前記関係取得手段は、前記試験材の磁化方向と、前記試験材の圧延方向とのなす角度毎に、前記関係を導出し、
前記磁束密度導出手段は、マックスウェルの方程式に基づく数値解析により導出した前記磁性体板の磁束密度ベクトルと、前記磁性体板の圧延方向とのなす角度を導出し、
前記鉄損導出手段は、前記磁束密度導出手段により導出された前記角度に対応する前記関係に基づいて、前記磁束密度導出手段により導出された磁束密度の大きさの代表値に対応する鉄損を、前記複数の微小領域毎の鉄損として導出することを特徴とする請求項1〜6の何れか1項に記載の電磁場解析システム。 - 前記関係取得手段は、前記試験材を励磁する際に前記試験材に与えられる応力ごとに、前記関係を導出し、
前記鉄損導出手段は、前記物体を励磁する際に前記物体に与えられる応力に対応する前記関係に基づいて、前記磁束密度導出手段により導出された磁束密度の大きさの代表値に対応する鉄損を、前記複数の微小領域毎の鉄損として導出することを特徴とする請求項1〜7の何れか1項に記載の電磁場解析システム。 - 前記関係取得手段は、インバータ電源により励磁した場合の前記試験材の磁束密度と鉄損を測定する磁気特性試験手段を更に有し、
前記磁気特性試験手段により測定された磁束密度と鉄損に基づいて前記関係を取得することを特徴とする請求項1〜8の何れか1項に記載の電磁場解析システム。 - 前記磁束密度の大きさの代表値は、前記磁束密度の波高値であることを特徴とする請求項1〜9の何れか1項に記載の電磁場解析システム。
- 鉄損の推定対象である物体に使用される磁性材料の試験材の磁束密度の大きさの代表値と鉄損との関係を取得する関係取得工程と、
前記物体に使用される磁性材料の磁束密度の大きさの代表値を、マックスウェルの方程式に基づく数値解析により、複数の微小領域毎に導出する磁束密度導出工程と、
前記磁束密度導出工程により導出された磁束密度の大きさの代表値に対応する鉄損を、前記複数の微小領域毎の鉄損として、前記関係取得工程により取得された前記関係に基づいて導出し、導出した前記複数の微小領域毎の鉄損に基づいて、前記物体の鉄損を導出する鉄損導出工程と、を有することを特徴とする電磁場解析方法。 - 鉄損の推定対象である物体に使用される磁性材料の試験材の磁束密度の大きさの代表値と鉄損との関係を取得する関係取得手段と、
前記物体に使用される磁性材料の磁束密度の大きさの代表値を、マックスウェルの方程式に基づく数値解析により、複数の微小領域毎に導出する磁束密度導出手段と、
前記磁束密度導出手段により導出された磁束密度の大きさの代表値に対応する鉄損を、前記複数の微小領域毎の鉄損として、前記関係取得手段により取得された前記関係に基づいて導出し、導出した前記複数の微小領域毎の鉄損に基づいて、前記物体の鉄損を導出する鉄損導出手段と、してコンピュータを機能させるためのプログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018098756A JP7124441B2 (ja) | 2018-05-23 | 2018-05-23 | 電磁場解析システム、電磁場解析方法、およびプログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018098756A JP7124441B2 (ja) | 2018-05-23 | 2018-05-23 | 電磁場解析システム、電磁場解析方法、およびプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2019204269A true JP2019204269A (ja) | 2019-11-28 |
JP7124441B2 JP7124441B2 (ja) | 2022-08-24 |
Family
ID=68726984
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2018098756A Active JP7124441B2 (ja) | 2018-05-23 | 2018-05-23 | 電磁場解析システム、電磁場解析方法、およびプログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7124441B2 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113704900A (zh) * | 2021-07-22 | 2021-11-26 | 无锡欧瑞京电机有限公司 | 基于磁路计算与电磁场校核的异步电机转子通风孔设计方法 |
CN114242425A (zh) * | 2021-10-22 | 2022-03-25 | 西南交通大学 | 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003098241A (ja) * | 2001-09-20 | 2003-04-03 | Nippon Steel Corp | 電磁場解析の鉄損評価システム |
JP2008170320A (ja) * | 2007-01-12 | 2008-07-24 | Nippon Steel Corp | 電磁場解析システム |
JP2017162187A (ja) * | 2016-03-09 | 2017-09-14 | 新日鐵住金株式会社 | 磁気特性解析装置、磁気特性解析方法、及びコンピュータプログラム |
JP6206608B1 (ja) * | 2017-01-06 | 2017-10-04 | 新日鐵住金株式会社 | 電磁場解析装置、電磁場解析方法、およびプログラム |
JP2018077682A (ja) * | 2016-11-09 | 2018-05-17 | 富士通株式会社 | 磁性材料シミュレーションプログラム、磁性材料シミュレーション方法および磁性材料シミュレーション装置 |
-
2018
- 2018-05-23 JP JP2018098756A patent/JP7124441B2/ja active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003098241A (ja) * | 2001-09-20 | 2003-04-03 | Nippon Steel Corp | 電磁場解析の鉄損評価システム |
JP2008170320A (ja) * | 2007-01-12 | 2008-07-24 | Nippon Steel Corp | 電磁場解析システム |
JP2017162187A (ja) * | 2016-03-09 | 2017-09-14 | 新日鐵住金株式会社 | 磁気特性解析装置、磁気特性解析方法、及びコンピュータプログラム |
JP2018077682A (ja) * | 2016-11-09 | 2018-05-17 | 富士通株式会社 | 磁性材料シミュレーションプログラム、磁性材料シミュレーション方法および磁性材料シミュレーション装置 |
JP6206608B1 (ja) * | 2017-01-06 | 2017-10-04 | 新日鐵住金株式会社 | 電磁場解析装置、電磁場解析方法、およびプログラム |
Non-Patent Citations (1)
Title |
---|
宮城大輔 外7名: "PWMインバータの変調率が無方向性電磁鋼板の鉄損に及ぼす影響", 平成22年 電気学会全国大会講演論文集[CD−ROM], JPN6022004797, 5 March 2010 (2010-03-05), pages 48 - 49, ISSN: 0004700829 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113704900A (zh) * | 2021-07-22 | 2021-11-26 | 无锡欧瑞京电机有限公司 | 基于磁路计算与电磁场校核的异步电机转子通风孔设计方法 |
CN113704900B (zh) * | 2021-07-22 | 2023-12-26 | 无锡欧瑞京电机有限公司 | 基于磁路计算与电磁场校核的异步电机转子通风孔设计方法 |
CN114242425A (zh) * | 2021-10-22 | 2022-03-25 | 西南交通大学 | 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 |
CN114242425B (zh) * | 2021-10-22 | 2023-06-02 | 西南交通大学 | 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 |
Also Published As
Publication number | Publication date |
---|---|
JP7124441B2 (ja) | 2022-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hargreaves et al. | Calculation of iron loss in electrical generators using finite-element analysis | |
Zarko et al. | Analytical calculation of magnetic field distribution in the slotted air gap of a surface permanent-magnet motor using complex relative air-gap permeance | |
Krings | Iron losses in electrical machines-influence of material properties, manufacturing processes, and inverter operation | |
Liu et al. | Nonlinear equivalent magnetic network of a linear permanent magnet vernier machine with end effect consideration | |
Gao et al. | Improved genetic optimization algorithm with subdomain model for multi-objective optimal design of SPMSM | |
Boughrara et al. | Analytical model of slotted air-gap surface mounted permanent-magnet synchronous motor with magnet bars magnetized in the shifting direction | |
JP6206608B1 (ja) | 電磁場解析装置、電磁場解析方法、およびプログラム | |
Bramerdorfer | Tolerance analysis for electric machine design optimization: Classification, modeling and evaluation, and example | |
JP6384208B2 (ja) | 電磁場解析装置、電磁場解析方法、およびプログラム | |
Krings et al. | Influence of the welding process on the performance of slotless PM motors with SiFe and NiFe stator laminations | |
Lin et al. | The effects of steel lamination core losses on 3D transient magnetic fields | |
Simao et al. | Evaluation of hysteresis losses in iron sheets under DC-biased inductions | |
JP7124441B2 (ja) | 電磁場解析システム、電磁場解析方法、およびプログラム | |
Gysen et al. | Analytical and numerical techniques for solving Laplace and Poisson equations in a tubular permanent-magnet actuator: Part I. Semi-analytical framework | |
Bramerdorfer et al. | Identification of a nonlinear PMSM model using symbolic regression and its application to current optimization scenarios | |
Pechlivanidou et al. | A comparative study on 2D and 3D magnetic field analysis of permanent magnet synchronous motor using FEM simulations | |
Zhu et al. | A general and accurate iron loss calculation method considering harmonics based on loss surface hysteresis model and finite-element method | |
Kim et al. | Simplified impedance modeling and analysis for inter-turn fault of IPM-type BLDC motor | |
JP5609376B2 (ja) | 鉄損の推定方法 | |
Herold et al. | System simulation of a PMSM servo drive using field-circuit coupling | |
CN106324505A (zh) | 一种电动机谐波分析方法 | |
JP6578845B2 (ja) | 電磁場解析装置、電磁場解析方法、およびプログラム | |
Pippuri et al. | Time-harmonic induction-machine model including hysteresis and eddy currents in steel laminations | |
Tekgun | Analysis, measurement and estimation of the core losses in electrical machines | |
Min | Inductance calculation of coreless-type linear PM machines based on analytical field projection and coil separation method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210112 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20220124 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220208 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220303 |
|
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: 20220712 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220725 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 7124441 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |