JP7087875B2 - Processing system, processing method, and program - Google Patents
Processing system, processing method, and program Download PDFInfo
- Publication number
- JP7087875B2 JP7087875B2 JP2018177461A JP2018177461A JP7087875B2 JP 7087875 B2 JP7087875 B2 JP 7087875B2 JP 2018177461 A JP2018177461 A JP 2018177461A JP 2018177461 A JP2018177461 A JP 2018177461A JP 7087875 B2 JP7087875 B2 JP 7087875B2
- Authority
- JP
- Japan
- Prior art keywords
- target signal
- signal waveform
- iron loss
- magnetic flux
- eddy current
- 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.)
- Active
Links
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
- Control Of Electric Motors In General (AREA)
Description
本発明は、処理システム、処理方法、およびプログラムに関し、特に、鉄心を励磁することによる動作する電気機器に用いて好適なものである。 The present invention relates to processing systems, processing methods, and programs, and is particularly suitable for use in electrical equipment that operates by exciting an iron core.
鉄心を励磁することにより動作する電気機器においては、その損失を低減することが求められる。例えば、モータは、一般に、インバータ電源を用いて、ステータコイルに励磁電圧を印加することで駆動する。このため、インバータ電源から出力される電圧波形を制御することでモータの損失（鉄損）を低減することが行われる。この種の技術として、特許文献１では、励磁電圧の波形をパルス波形とすると共に、電気一周期のパルス数を５として、モータにおける電力損失が最小になるパルス幅を導出する。このとき、特許文献１では、モータを励磁した際のモータの損失（鉄損）を、電磁場解析を行うことにより導出する。また、特許文献１では、パルス波形が擬似正弦波になるようにする。 In electrical equipment that operates by exciting an iron core, it is required to reduce the loss. For example, a motor is generally driven by applying an exciting voltage to a stator coil using an inverter power supply. Therefore, the loss (iron loss) of the motor can be reduced by controlling the voltage waveform output from the inverter power supply. As this kind of technique, in Patent Document 1, the waveform of the exciting voltage is used as a pulse waveform, and the number of pulses in one electric cycle is set to 5, and the pulse width that minimizes the power loss in the motor is derived. At this time, in Patent Document 1, the loss (iron loss) of the motor when the motor is excited is derived by performing electromagnetic field analysis. Further, in Patent Document 1, the pulse waveform is set to be a pseudo sine wave.
しかしながら、特許文献１に記載の技術では、電気一周期のパルス数を５としている。実際にモータを駆動する場合には、電気一周期のパルス数が５である場合は少ない。このため、特許文献１に記載の技術を、実際のモータの駆動に適用する場合には、電気一周期のパルス数を多くしなければならず、決定すべき変数が多くなる。更に、特許文献１に記載の技術では、モータの鉄心（ステータコア）そのものを電磁場解析する。従って、計算負荷が非常に大きくなる。また、特許文献１に記載の技術では、電流の高調波成分に着目しているが、どのような電磁場解析を行うのかについて、具体的に示していない。また、特許文献１に記載の技術では、パルス幅のみを変数としているため、パルスの信号レベルを考慮することができない。従って、モータの鉄損を十分に低減することが容易ではない。 However, in the technique described in Patent Document 1, the number of pulses in one electric cycle is set to 5. When actually driving the motor, the number of pulses in one electric cycle is small. Therefore, when the technique described in Patent Document 1 is applied to drive an actual motor, the number of pulses in one electric cycle must be increased, and the number of variables to be determined increases. Further, in the technique described in Patent Document 1, the iron core (stator core) of the motor itself is subjected to electromagnetic field analysis. Therefore, the calculation load becomes very large. Further, in the technique described in Patent Document 1, the harmonic component of the current is focused on, but what kind of electromagnetic field analysis is performed is not specifically shown. Further, in the technique described in Patent Document 1, since only the pulse width is used as a variable, the signal level of the pulse cannot be considered. Therefore, it is not easy to sufficiently reduce the iron loss of the motor.
本発明は、以上のような問題点に鑑みてなされたものであり、鉄心を励磁することによる動作する電気機器における鉄損を低減することができる励磁信号の時間波形を短時間に決定することができるようにすることを目的とする。 The present invention has been made in view of the above problems, and it is intended to determine a time waveform of an excitation signal in a short time, which can reduce iron loss in an operating electric device by exciting an iron core. The purpose is to be able to do.
本発明の処理システムは、鉄心を有する電気機器を動作させるための処理を行う処理システムであって、前記電気機器に印加される励磁信号の時間波形の目標値である目標信号波形を決定する目標信号波形決定手段を有し、前記目標信号波形決定手段は、前記鉄心の鉄損解析用試料が励磁された際の磁束密度および渦電流密度の分布を、マクスウェルの方程式に基づいて導出し、当該磁束密度および渦電流密度の分布に基づいて、当該鉄損解析用試料の鉄損を導出し、当該鉄損解析用試料の鉄損に基づいて、前記目標信号波形を決定し、前記磁束密度および渦電流密度の分布は、前記鉄損解析用試料の内部において磁束および渦電流が浸透する深さ方向における分布を含むことを特徴とする。 The processing system of the present invention is a processing system that performs processing for operating an electric device having an iron core, and is a target for determining a target signal waveform which is a target value of a time waveform of an excitation signal applied to the electric device. The target signal waveform determining means has a signal waveform determining means, and the target signal waveform determining means derives the distribution of the magnetic flux density and the eddy current density when the iron loss analysis sample of the iron core is excited based on Maxwell's equation. Based on the distribution of the magnetic flux density and the eddy current density, the iron loss of the iron loss analysis sample is derived, and the target signal waveform is determined based on the iron loss of the iron loss analysis sample, and the magnetic flux density and the magnetic flux density and the iron loss are determined. The distribution of the eddy current density is characterized by including the distribution in the depth direction in which the magnetic flux and the eddy current permeate inside the sample for iron loss analysis.
本発明の処理方法は、鉄心を有する電気機器を動作させるための処理を行う処理方法であって、前記電気機器に印加される励磁信号の時間波形の目標値である目標信号波形を決定する目標信号波形決定工程を有し、前記目標信号波形決定工程は、前記鉄心の鉄損解析用試料が励磁された際の磁束密度および渦電流密度の分布を、マクスウェルの方程式に基づいて導出し、当該磁束密度および渦電流密度の分布に基づいて、当該鉄損解析用試料の鉄損を導出し、当該鉄損解析用試料の鉄損に基づいて、前記目標信号波形を決定し、前記磁束密度および渦電流密度の分布は、前記鉄損解析用試料の内部において磁束および渦電流が浸透する深さ方向における分布を含むことを特徴とする。 The processing method of the present invention is a processing method for operating an electric device having an iron core, and is a target for determining a target signal waveform which is a target value of a time waveform of an excitation signal applied to the electric device. It has a signal waveform determination step, and the target signal waveform determination step derives the distribution of the magnetic flux density and the eddy current density when the sample for iron loss analysis of the iron core is excited based on Maxwell's equation. Based on the distribution of the magnetic flux density and the eddy current density, the iron loss of the iron loss analysis sample is derived, and the target signal waveform is determined based on the iron loss of the iron loss analysis sample, and the magnetic flux density and the magnetic flux density and the iron loss are determined. The distribution of the eddy current density is characterized by including the distribution in the depth direction in which the magnetic flux and the eddy current permeate inside the sample for iron loss analysis.
本発明のプログラムは、前記処理システムの各手段としてコンピュータを機能させるためのものである。 The program of the present invention is for making a computer function as each means of the processing system.
本発明によれば、鉄心を励磁することによる動作する電気機器における鉄損を低減することができる励磁信号の時間波形を短時間に決定することができる。 According to the present invention, it is possible to determine the time waveform of the excitation signal in a short time, which can reduce the iron loss in the electric device operating by exciting the iron core.
以下、図面を参照しながら、本発明の一実施形態を説明する。
本実施形態の処理システムは、目標信号波形決定装置１００と、駆動装置２００と、モータＭとを有する。このように本実施形態では、鉄心を励磁することにより動作する電気機器がモータＭである場合を例に挙げて説明する。また、本実施形態では、ＰＷＭ（Pulse Width Modulation）制御により得られるパルス信号（パルス電圧）をモータＭに印加する場合を例に挙げて説明する。
Hereinafter, an embodiment of the present invention will be described with reference to the drawings.
The processing system of the present embodiment includes a target signal
［目標信号波形決定装置１００］
目標信号波形決定装置１００は、目標信号波形を生成する。目標信号波形は、駆動装置２００から電気機器に印加される励磁信号の時間波形の目標値である。本実施形態では、モータＭのステータコイルに印加する励磁電圧の時間波形の目標値として、ＰＷＭ制御により得られるパルス信号の時間波形の目標値が、目標信号波形になる。目標信号波形は、電気一周期の各時刻における値が特定されればよい。例えば、目標信号波形を、電気一周期の各時刻における値とし、当該値から、各周期の各時刻における値を導出することができる。また、例えば、電気半周期の各時刻における値から、残りの半周期の各時刻における値を導出する場合、目標信号波形は、電気半周期の各時刻における値とすればよい。以下の説明では、電気一周期の各時刻における電圧値を、目標信号波形とする場合を例に挙げて説明する。
[Target signal waveform determining device 100]
The target signal
図１は、目標信号波形決定装置１００の機能的な構成の一例を示す図である。目標信号波形決定装置１００のハードウェアは、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ、ＨＤＤ、および各種のインターフェースを備える情報処理装置、または、専用のハードウェアを用いることにより実現される。本実施形態では、目標信号波形決定装置１００は、遺伝的アルゴリズム（ＧＡ：Genetic Algorithm）を用いて、目標信号波形の最適値を導出する。その一例を以下に説明する。
FIG. 1 is a diagram showing an example of a functional configuration of the target signal
＜候補解設定部１０１＞
候補解設定部１０１は、モータＭの速度指令値およびトルク指令値として想定される値を入力する。以下の説明では、モータＭの速度指令値として想定される値を、必要に応じて、モータＭの速度指令の想定値と称する。また、モータＭのトルク指令値として想定される値を、必要に応じて、モータＭのトルク指令の想定値と称する。モータＭの速度指令の想定値およびトルク指令の想定値の入力の形態としては、例えば、目標信号波形決定装置１００のユーザインタフェースのオペレータによる入力操作、外部装置からの受信、または、可搬型記憶媒体からの読み出しが挙げられる。
<Candidate solution setting unit 101>
The candidate solution setting unit 101 inputs values assumed as the speed command value and the torque command value of the motor M. In the following description, the value assumed as the speed command value of the motor M is referred to as the assumed value of the speed command of the motor M, if necessary. Further, the value assumed as the torque command value of the motor M is referred to as the assumed value of the torque command of the motor M, if necessary. The input form of the assumed value of the speed command and the assumed value of the torque command of the motor M is, for example, an input operation by the operator of the user interface of the target signal
候補解設定部１０１は、遺伝的アルゴリズムに従って、それぞれが、目標信号波形を特定する情報からなる個体を候補解として複数設定する。このとき、候補解設定部１０１は、電気周期（電気周波数）が、モータＭの速度指令の想定値に応じた電気周期（電気周波数）になるように、複数の候補解（目標信号波形）を設定する。目標信号波形を特定する情報は、例えば、電気一周期の各時刻における目標信号波形の値である。 The candidate solution setting unit 101 sets, as candidate solutions, a plurality of individuals each consisting of information for specifying a target signal waveform according to a genetic algorithm. At this time, the candidate solution setting unit 101 sets a plurality of candidate solutions (target signal waveforms) so that the electric cycle (electric frequency) becomes the electric cycle (electric frequency) corresponding to the assumed value of the speed command of the motor M. Set. The information that identifies the target signal waveform is, for example, the value of the target signal waveform at each time of the electric cycle.
候補解設定部１０１は、初期の候補解を設定する際には、例えば、電気周期（電気周波数）および振幅以外のパラメータについてはランダムに設定して初期の候補解を設定する。
また、候補解設定部１０１は、２回目以降に候補解を設定する際には、例えば、複数の候補解のうち、後述する電磁場解析部１０２により導出される鉄心試料の鉄損が小さくなるものから順に所定数の候補解を選択する。このとき、候補解設定部１０１は、後述する電磁場解析部１０２により導出される鉄心試料の磁束密度の実効値が、トルク指令の想定値に対応する値以上になる候補解のみを選択する。即ち、候補解設定部１０１は、後述する電磁場解析部１０２により導出される鉄心試料の磁束密度の実効値が、トルク指令の想定値に対応する値以上になる候補解を、鉄心が小さくなるものから順番に並べ、並べた候補解を鉄損が小さいものから順に所定数だけ選択する。
When setting the initial candidate solution, the candidate solution setting unit 101 sets the initial candidate solution by randomly setting, for example, parameters other than the electric cycle (electric frequency) and the amplitude.
Further, when the candidate solution setting unit 101 sets the candidate solution from the second time onward, for example, among a plurality of candidate solutions, the iron loss of the iron core sample derived by the electromagnetic field analysis unit 102, which will be described later, becomes small. Select a predetermined number of candidate solutions in order from. At this time, the candidate solution setting unit 101 selects only candidate solutions in which the effective value of the magnetic flux density of the iron core sample derived by the electromagnetic field analysis unit 102, which will be described later, is equal to or greater than the value corresponding to the assumed value of the torque command. That is, the candidate solution setting unit 101 reduces the size of the candidate solution in which the effective value of the magnetic flux density of the iron core sample derived by the electromagnetic field analysis unit 102, which will be described later, is equal to or greater than the value corresponding to the assumed value of the torque command. Arrange in order from the beginning, and select a predetermined number of the arranged candidate solutions in order from the one with the smallest iron loss.
そして、候補解設定部１０１は、交叉や突然変異を行い、新たな候補解を設定する。新たな候補解（目標信号波形）も、電気周期（電気周波数）が、モータＭの速度指令の想定値に応じた電気周期（電気周波数）になるように設定される。尚、トルク指令の想定値に対応する磁束密度は、トルクと磁束密度との関係式から得られる。例えば、トルクをマクスウェル応力とし、マクスウェル応力と磁束密度との公知の関係式から、トルク（トルク指令の想定値）に対応する磁束密度（の実効値）を得ることができる。後述するように、本実施形態では、鉄心試料の鉄損が最小となる解を探索する。従って、トルク指令の想定値に対応する磁束密度の実効値の下限値を定めておけば、トルク指令の想定値に近い（好ましくは一致する）磁束密度の実効値が得られる。 Then, the candidate solution setting unit 101 crosses or mutates and sets a new candidate solution. The new candidate solution (target signal waveform) is also set so that the electric cycle (electric frequency) becomes the electric cycle (electric frequency) corresponding to the assumed value of the speed command of the motor M. The magnetic flux density corresponding to the assumed value of the torque command is obtained from the relational expression between the torque and the magnetic flux density. For example, the torque is defined as Maxwell stress, and the magnetic flux density (effective value) corresponding to the torque (assumed value of the torque command) can be obtained from the known relational expression between the Maxwell stress and the magnetic flux density. As will be described later, in this embodiment, a solution that minimizes the iron loss of the iron core sample is searched for. Therefore, if the lower limit of the effective value of the magnetic flux density corresponding to the assumed value of the torque command is set, the effective value of the magnetic flux density close to (preferably matching) the assumed value of the torque command can be obtained.
＜電磁場解析部１０２＞
電磁場解析部１０２は、候補解設定部１０１により設定された複数の候補解（目標信号波形）のそれぞれで鉄心試料を励磁した際の鉄損を、マクスウェル（Maxwell）の方程式に基づく電磁場解析を行うことにより導出する。
<Electromagnetic field analysis unit 102>
The electromagnetic field analysis unit 102 performs electromagnetic field analysis based on Maxwell's equations for the iron loss when the iron core sample is excited by each of the plurality of candidate solutions (target signal waveforms) set by the candidate solution setting unit 101. Derived by
<<着想>>
一般に、電気自動車やハイブリッド車等の駆動用モータのように、可変速、可変トルクで運転するモータは、インバータによって駆動し、モータの損失はインバータの駆動条件によって異なる。そのため、インバータの電圧パルスパターンを適切に制御することで、モータ効率を向上できる可能性がある。モータでの損失は、銅損、機械損、鉄損に大別される。これらのうち、銅損、機械損は、インバータ条件によって大きく変動しないのに対し、鉄損は、インバータ条件によって大きく変動する。鉄損は、渦電流損およびヒステリシス損に分類され、渦電流損は電磁鋼板などの磁性体内部の渦電流密度によって決定され、ヒステリシス損は磁束密度によって決定される。
<< Idea >>
Generally, a motor that operates at a variable speed and a variable torque, such as a drive motor for an electric vehicle or a hybrid vehicle, is driven by an inverter, and the loss of the motor differs depending on the drive conditions of the inverter. Therefore, there is a possibility that the motor efficiency can be improved by appropriately controlling the voltage pulse pattern of the inverter. Losses in motors are roughly classified into copper loss, mechanical loss, and iron loss. Of these, copper loss and mechanical loss do not fluctuate significantly depending on the inverter conditions, whereas iron loss fluctuates greatly depending on the inverter conditions. Iron loss is classified into eddy current loss and hysteresis loss, eddy current loss is determined by the eddy current density inside a magnetic material such as an electromagnetic steel plate, and hysteresis loss is determined by the magnetic flux density.
ここで、本発明者は、渦電流密度および磁束密度は、表皮効果により磁性体内部で減衰するため、渦電流損およびヒステリシス損を精度良く見積もるためには、この減衰の影響を含めて電磁場解析をする必要があることを着想した。図２は、渦電流密度が磁性体内部で減衰する様子を概念的に示す図である。図２において、ｚ軸方向は、鉄心試料において、磁束および渦電流が浸透する深さ方向（表皮深さの方向）である。本実施形態では、鉄心試料は、磁性体板（電磁鋼板）により構成されるので、ｚ軸方向は、磁性体板の板面に対して垂直な方向（板厚方向）である。また、図２において、ｚ軸に垂直な方向を向いている実線の矢印線は、渦電流を示す。当該矢印線の長さは、渦電流密度の大きさに対応する。図２に示すように、表皮効果により、磁性体板の内部の位置であるほど、渦電流密度は小さくなる。このことは、磁束密度についても同じである。 Here, since the eddy current density and the magnetic flux density are attenuated inside the magnetic material due to the skin effect, the present inventor analyzes the electromagnetic field including the influence of the attenuation in order to accurately estimate the eddy current loss and the hysteresis loss. I came up with the idea that I need to do. FIG. 2 is a diagram conceptually showing how the eddy current density is attenuated inside the magnetic material. In FIG. 2, the z-axis direction is the depth direction (direction of the skin depth) in which the magnetic flux and the eddy current permeate in the iron core sample. In the present embodiment, since the iron core sample is composed of a magnetic material plate (electromagnetic steel plate), the z-axis direction is a direction perpendicular to the plate surface of the magnetic material plate (plate thickness direction). Further, in FIG. 2, a solid arrow line pointing in a direction perpendicular to the z-axis indicates an eddy current. The length of the arrow line corresponds to the magnitude of the eddy current density. As shown in FIG. 2, due to the skin effect, the eddy current density becomes smaller as the position is inside the magnetic plate. This also applies to the magnetic flux density.
本発明者は、渦電流の表皮深さは、磁性体の透磁率によって変化するため、高透磁率領域では渦電流の減衰量が大きいために鉄損が低減されるという現象から、鉄損を低減できるような励磁信号の時間波形は、擬似正弦波に限られないと考えた。また、前述したように、このような渦電流密度および磁束密度の減衰は、表皮深さの方向（板厚方向）で生じる。本発明者は、このような減衰を考慮することができれば、モータＭのステータコアそのものに対して電磁場解析を行わなくても、実用上要求される精度で鉄損を導出することができると考えた。
このような着想の下、本実施形態では、モータＭのステータコアを簡略化した形状の鉄心試料における磁束密度および渦電流密度の板厚方向の分布を導出し、当該磁束密度および渦電流密度に基づいて鉄損を導出し、当該鉄損をモータＭの鉄損として推定する。このようにすることにより、大きく精度を低下させることなく短時間（実用上実現できる時間内）に、モータＭの鉄損を推定することができる。
According to the present inventor, since the skin depth of the eddy current changes depending on the magnetic permeability of the magnetic material, the iron loss is reduced due to the large attenuation of the eddy current in the high magnetic permeability region. It is considered that the time waveform of the excitation signal that can be reduced is not limited to the pseudo sine wave. Further, as described above, such attenuation of the eddy current density and the magnetic flux density occurs in the direction of the skin depth (the plate thickness direction). The present inventor thinks that if such attenuation can be taken into consideration, iron loss can be derived with the accuracy required for practical use without performing electromagnetic field analysis on the stator core itself of the motor M. ..
Based on such an idea, in the present embodiment, the distribution of the magnetic flux density and the eddy current density in the plate thickness direction in the iron core sample having a simplified shape of the stator core of the motor M is derived, and the distribution is based on the magnetic flux density and the eddy current density. The iron loss is derived, and the iron loss is estimated as the iron loss of the motor M. By doing so, it is possible to estimate the iron loss of the motor M in a short time (within a practically feasible time) without significantly reducing the accuracy.
<<鉄心試料>>
鉄心試料の一例について説明する。
本実施形態では、鉄心試料は、モータＭのステータコアの鉄損解析用試料である。鉄心試料は、モータＭのステータコアよりも簡略化した形状を有し、モータＭのステータコアに使用される材料と同じ材質の軟質磁性材料（電磁鋼板等）で構成されるものとするのが好ましい。図３は、鉄心試料の一例を示す図である。図３（ａ）は、単板（１枚の矩形状の板）である鉄心試料ＣＳ１を示し、図３（ｂ）は、リング状（１枚のリング状の板）である鉄心試料ＣＳ２を示す。図３（ａ）および図３（ｂ）において、鉄心試料ＣＳ１、ＣＳ２に対して巻き回されているコイルＣ１、Ｃ２は、鉄心試料ＣＳ１、ＣＳ２を励磁する際に用いる励磁コイルを概念的に示す図である。目標信号波形に対応する励磁電流がコイルＣ１、Ｃ２に流れることにより、図３（ａ）および図３（ｂ）において鉄心試料ＣＳ１、ＣＳ２の内側に両矢印線で示すように磁束が発生する。鉄心試料ＣＳ１、ＣＳ２の板厚は、モータＭのステータコアに使用される軟磁性体板の板厚と同じであるのが好ましいが、異なっていてもよい。ただし、後述する電磁場解析において、板厚方向に複数の離散化領域（複数のメッシュ）を設定することができる厚みを有する必要がある。
<< Iron core sample >>
An example of an iron core sample will be described.
In the present embodiment, the iron core sample is a sample for iron loss analysis of the stator core of the motor M. The iron core sample has a shape simplified more than that of the stator core of the motor M, and is preferably made of a soft magnetic material (electromagnetic steel sheet or the like) made of the same material as the material used for the stator core of the motor M. FIG. 3 is a diagram showing an example of an iron core sample. FIG. 3 (a) shows the iron core sample CS1 which is a single plate (one rectangular plate), and FIG. 3 (b) shows the iron core sample CS2 which is a ring-shaped (one ring-shaped plate). show. In FIGS. 3 (a) and 3 (b), the coils C1 and C2 wound around the iron core samples CS1 and CS2 conceptually show the excitation coils used when exciting the iron core samples CS1 and CS2. It is a figure. When the exciting current corresponding to the target signal waveform flows through the coils C1 and C2, magnetic flux is generated inside the iron core samples CS1 and CS2 in FIGS. 3A and 3B as shown by double-headed arrows. The plate thicknesses of the iron core samples CS1 and CS2 are preferably the same as the plate thickness of the soft magnetic material plate used for the stator core of the motor M, but may be different. However, in the electromagnetic field analysis described later, it is necessary to have a thickness that allows a plurality of discretized regions (multiple meshes) to be set in the plate thickness direction.
図３（ａ）および図３（ｂ）では、鉄心試料ＣＳ１、ＣＳ２が、１枚の板である場合を例に挙げて説明した。しかしながら、鉄心試料は、必ずしも１枚の板である必要はない。例えば、複数枚の板を積み重ねたものを鉄心試料としてもよい。複数枚の板は、同じ形状および大きさを有していてもいなくてもよい。また、同じ形状および大きさを有する複数枚の板を積み重ねる際に、縁の全ての部分を一致させなくてもよい。同じ形状および大きさを有する複数枚の板を積み重ねる際に、縁の全ての部分を一致させないようにする例として、公知のエプスタイン試験器における試験片のように、端部が一枚ずつ交互に重なり合うようにして正方形に組み、長さと断面積の等しい４個の辺が形成されるようにする構成が挙げられる。 In FIGS. 3A and 3B, the case where the iron core samples CS1 and CS2 are one plate has been described as an example. However, the iron core sample does not necessarily have to be a single plate. For example, a stack of a plurality of plates may be used as an iron core sample. The plurality of plates may or may not have the same shape and size. Also, when stacking a plurality of boards having the same shape and size, it is not necessary to match all the edges. As an example of preventing all edges from matching when stacking multiple plates of the same shape and size, the edges alternate one by one, as in a test piece in a known Epstein tester. There is a configuration in which the squares are assembled so as to overlap each other so that four sides having the same length and cross-sectional area are formed.
<<電磁場解析>>
次に、電磁場解析の手法について説明する。
本実施形態では、電磁場解析部１０２は、非線形非定常有限要素法を用いた電磁場解析を行うことにより、鉄心試料に対して設定した要素（メッシュ）のそれぞれにおいて、目標信号波形に従って励磁された場合の鉄心試料の磁束密度Ｂと渦電流密度Ｊｅを導出する場合を例に挙げて説明する。
<< Electromagnetic field analysis >>
Next, the method of electromagnetic field analysis will be described.
In the present embodiment, when the electromagnetic field analysis unit 102 is excited according to the target signal waveform in each of the elements (mesh) set for the iron core sample by performing the electromagnetic field analysis using the non-stationary unsteady finite element method. The case of deriving the magnetic flux density B and the eddy current density Je of the iron core sample will be described as an example.
有限要素法を用いた電磁場解析の手法としてＡ－φ法を用いる手法がある。この場合、電磁場解析を行うための基礎方程式は、マクスウェルの方程式に基づき以下の（１）式～（４）式で与えられる。尚、各式において、→は、ベクトルであることを表す。 There is a method using the A-φ method as a method of electromagnetic field analysis using the finite element method. In this case, the basic equations for performing electromagnetic field analysis are given by the following equations (1) to (4) based on Maxwell's equations. In each equation, → indicates that it is a vector.
（１）式～（４）式において、μは、透磁率であり、Ａは、ベクトルポテンシャルであり、σは、導電率であり、Ｊ0は、励磁電流密度であり、Ｊｅは、渦電流密度であり、Ｂは、磁束密度である。（１）式および（２）式を連立して解いて、ベクトルポテンシャルＡとスカラーポテンシャルφを求めた後、（３）式および（４）式から磁束密度Ｂと、渦電流密度Ｊｅを要素のそれぞれに対して求める。尚、（１）式では、表記を簡素化するため、透磁率のｘ成分μx、ｙ成分μy、ｚ成分μzが等しい場合（μx＝μy＝μzの場合）の式を示す。 In equations (1) to (4), μ is the magnetic permeability, A is the vector potential, σ is the conductivity, J 0 is the exciting current density, and Je is the eddy current. It is a density, and B is a magnetic flux density. After solving equations (1) and (2) simultaneously to obtain the vector potential A and the scalar potential φ, the magnetic flux density B and the eddy current density Je are the elements from the equations (3) and (4). Ask for each. In equation (1), in order to simplify the notation, the equation when the x component μ x , y component μ y , and z component μ z of magnetic permeability are equal (when μ x = μ y = μ z ) is used. show.
鉄心試料の磁束密度Ｂ、渦電流密度Ｊｅの板厚方向における分布を導出していれば、必ずしも３次元解析を行う必要はない（ｘ成分、ｙ成分、ｚ成分の全てを導出する必要はない）。例えば、図２に示すように、鉄心試料の板面に垂直な断面（鉄心試料を板厚方向（ｚ軸方向）に沿って切った断面）において、板面方向の位置（ｘ－ｙ平面の位置）が同じ位置（ｘ軸およびｙ軸の値が同じ位置）での、磁束密度Ｂおよび渦電流密度Ｊｅの板厚方向の分布を導出してもよい。また、板面方向の複数の位置（ｘ－ｙ平面の複数の位置）での、磁束密度Ｂおよび渦電流密度Ｊｅの板厚方向の分布を導出してもよい。以上のようにする場合、例えば、（１）式において、透磁率μをｘ成分、ｙ成分、ｚ成分毎に表記した上で、ベクトルポテンシャルのｚ成分Ａz、∂／∂ｘ、∂／∂ｙをそれぞれ０（ゼロ）とする（Ａz＝０、∂／∂ｘ＝０、∂／∂ｙ＝０）。尚、電磁場解析を行う手法は、非特許文献１等に記載されているように一般的な手法であるので、その詳細な説明を省略する。 If the distributions of the magnetic flux density B and the eddy current density Je of the iron core sample in the plate thickness direction are derived, it is not always necessary to perform three-dimensional analysis (it is not necessary to derive all of the x component, y component, and z component). ). For example, as shown in FIG. 2, in a cross section perpendicular to the plate surface of the iron core sample (a cross section obtained by cutting the iron core sample along the plate thickness direction (z-axis direction)), the position in the plate surface direction (of the xy plane). The distribution of the magnetic flux density B and the eddy current density Je in the plate thickness direction may be derived at the same position (position) at the same position (the value of the x-axis and the y-axis is the same). Further, the distributions of the magnetic flux density B and the eddy current density Je in the plate thickness direction at a plurality of positions in the plate surface direction (plural positions in the xy plane) may be derived. In the above case, for example, in equation (1), the magnetic permeability μ is expressed for each of the x component, y component, and z component, and then the z component A z , ∂ / ∂x, and ∂ / ∂ of the vector potential. Let y be 0 (zero), respectively (A z = 0, ∂ / ∂x = 0, ∂ / ∂y = 0). Since the method for performing electromagnetic field analysis is a general method as described in Non-Patent Document 1 and the like, detailed description thereof will be omitted.
電磁場解析部１０２は、以上のようにして、目標信号波形の電気一周期の各時間ステップｔにおいて、鉄心試料の各要素における磁束密度Ｂおよび渦電流密度Ｊｅを導出する。例えば、電磁場解析部１０２は、目標信号波形の電気一周期の各時間ステップｔにおいて、鉄心試料の各要素における、鉄心試料の磁束密度（のｘ成分およびｙ成分）Ｂx、Ｂyおよび渦電流密度（のｘ成分、ｙ成分）Ｊｅx、Ｊｅyを導出する。
そして、電磁場解析部１０２は、鉄心試料の各要素における磁束密度Ｂおよび渦電流密度Ｊｅを用いて、鉄心試料のヒステリシス損および渦電流損（古典的渦電流損）を導出する。
As described above, the electromagnetic field analysis unit 102 derives the magnetic flux density B and the eddy current density Je in each element of the iron core sample at each time step t of the electric cycle of the target signal waveform. For example, the electromagnetic field analysis unit 102 may use the magnetic flux densities (x and y components) of the iron core sample (x and y components) B x , By and eddy currents in each element of the iron core sample at each time step t of the electric cycle of the target signal waveform. The densities ( x component, y component) Jex and Jey are derived.
Then, the electromagnetic field analysis unit 102 derives the hysteresis loss and the eddy current loss (classical eddy current loss) of the iron core sample by using the magnetic flux density B and the eddy current density Je in each element of the iron core sample.
例えば、電磁場解析部１０２は、各要素の磁束密度から磁束密度ベクトルを導出し、当該磁束密度ベクトルの大きさを導出することを、目標信号波形における電気一周期の各時刻ステップｔのそれぞれにおいて行う。電磁場解析部１０２は、その結果から、磁束密度ベクトルの大きさと時間との関係を示す波形の最大値と最小値の差Ｂmを導出し、以下の（５）式により、各要素におけるヒステリシス損Ｗｈを導出する。 For example, the electromagnetic field analysis unit 102 derives the magnetic flux density vector from the magnetic flux density of each element and derives the magnitude of the magnetic flux density vector at each time step t of one electric cycle in the target signal waveform. .. From the result, the electromagnetic field analysis unit 102 derives the difference B m between the maximum value and the minimum value of the waveform showing the relationship between the magnitude and time of the magnetic flux density vector, and the hysteresis loss in each element is derived by the following equation (5). Wh is derived.
（５）式において、ｆは、励磁周波数（モータＭの速度指令の想定値に応じた電気周波数）であり、Ｋｈは、ヒステリシス損係数であり、βは定数（例えば、１．６または２）である。ヒステリシス損係数Ｋｈは、例えば、二周波法により予め求められるものである。
電磁場解析部１０２は、以上のヒステリシス損Ｗｈの導出を全ての要素に対して行い、全ての要素におけるヒステリシス損Ｗｈの総和を鉄心試料のヒステリシス損として導出する。尚、ヒステリシス損は、公知の方法で導出することができ、各要素の磁束密度を用いて導出する方法であれば、どのような方法で導出してもよい。
In equation (5), f is the excitation frequency (electrical frequency according to the assumed value of the speed command of the motor M), Kh is the hysteresis loss coefficient, and β is a constant (for example, 1.6 or 2). Is. The hysteresis loss coefficient Kh is obtained in advance by, for example, the dual frequency method.
The electromagnetic field analysis unit 102 derives the above-mentioned hysteresis loss Wh for all the elements, and derives the sum of the hysteresis loss Wh in all the elements as the hysteresis loss of the iron core sample. The hysteresis loss can be derived by a known method, and may be derived by any method as long as it is derived by using the magnetic flux density of each element.
また、電磁場解析部１０２は、各要素の渦電流密度Ｊｅから渦電流密度ベクトルを導出し、当該渦電流密度ベクトルの大きさを導出することを、目標信号波形の電気一周期の各時刻ステップｔのそれぞれにおいて行う。そして、電磁場解析部１０２は、以下の（６）式により、各要素における渦電流損Ｗｅ（古典的渦電流損）を導出する。 Further, the electromagnetic field analysis unit 102 derives the eddy current density vector from the eddy current density Je of each element and derives the magnitude of the eddy current density vector at each time step t of the electric cycle of the target signal waveform. Do in each of. Then, the electromagnetic field analysis unit 102 derives the eddy current loss We (classical eddy current loss) in each element by the following equation (6).
（６）式において、ｖは、要素の大きさであり、Ｔは、目標信号波形における電気一周期に相当する時間である。要素の大きさは、１次元解析を行う場合には長さ、２次元解析を行う場合には面積、３次元解析を行う場合には体積の大きさになる。
電磁場解析部１０２は、以上の渦電流損Ｗｅの導出を全ての要素に対して行い、全ての要素における渦電流損Ｗｅの総和を電磁鋼板の渦電流損として導出する。尚、渦電流損は、公知の方法で導出することができ、各要素の渦電流密度を用いて導出する方法であれば、どのような方法で導出してもよい。
そして、電磁場解析部１０２は、鉄心試料のヒステリシス損と渦電流損の和を鉄心試料の鉄損として導出する。
In equation (6), v is the size of the element, and T is the time corresponding to one electrical cycle in the target signal waveform. The size of the element is the size when performing one-dimensional analysis, the area when performing two-dimensional analysis, and the size of volume when performing three-dimensional analysis.
The electromagnetic field analysis unit 102 derives the above eddy current loss We for all the elements, and derives the sum of the eddy current loss We in all the elements as the eddy current loss of the electromagnetic steel sheet. The eddy current loss can be derived by a known method, and any method may be used as long as it is derived by using the eddy current density of each element.
Then, the electromagnetic field analysis unit 102 derives the sum of the hysteresis loss and the eddy current loss of the iron core sample as the iron loss of the iron core sample.
以上のようにして、１つの目標信号波形に対して、モータＭの速度指令の想定値およびトルク指令の想定値に対応する鉄心試料の鉄損が１つ導出される。電磁場解析部１０２は、候補解設定部１０１で設定された複数の候補解（目標信号波形）のそれぞれに対して、モータＭの速度指令の想定値およびトルク指令の想定値に対応する鉄心試料の鉄損を導出する。 As described above, for one target signal waveform, one iron loss of the iron core sample corresponding to the assumed value of the speed command of the motor M and the assumed value of the torque command is derived. The electromagnetic field analysis unit 102 sets the iron core sample corresponding to the assumed value of the speed command and the assumed value of the torque command of the motor M for each of the plurality of candidate solutions (target signal waveforms) set by the candidate solution setting unit 101. Derive iron loss.
そして、前述したように候補解設定部１０１は、電磁場解析部１０２により導出された鉄心試料の磁束密度の実効値が、トルク指令の想定値に対応する値以上になる候補解（目標信号波形）を、鉄心が小さくなるものから順番に並べ、並べた候補解を、以上のようにして導出した鉄損が小さいものから順に所定数だけ選択する。そして、候補解設定部１０１は、交叉や突然変異を行い、新たな候補解を設定する。 Then, as described above, the candidate solution setting unit 101 has a candidate solution (target signal waveform) in which the effective value of the magnetic flux density of the iron core sample derived by the electromagnetic field analysis unit 102 is equal to or higher than the value corresponding to the assumed value of the torque command. Are arranged in order from the one with the smallest iron core, and a predetermined number of the arranged candidate solutions are selected in order from the one with the smallest iron loss derived as described above. Then, the candidate solution setting unit 101 crosses or mutates and sets a new candidate solution.
＜終了判定部１０３＞
終了判定部１０３は、電磁場解析部１０２における鉄損の導出の終了条件を満足したか否かを判定する。終了条件としては、例えば、所定数の世代交代が行われたことや、前回と今回とでの候補解の差が所定の条件になったこと（例えば、前回と今回とでの候補解が変わらなくなったこと）等、遺伝的アルゴリズムの手法で一般的に採用されている条件を用いることができる。
尚、遺伝的アルゴリズム自体は公知の技術で実現することができるので、ここでは、その詳細な説明を省略する。
前述した候補解設定部１０１における候補解の設定と、電磁場解析部１０２における鉄心試料の鉄損の導出は、終了判定部１０３により終了条件を満足すると判定されるまで繰り返し実行される。
<End determination unit 103>
The end determination unit 103 determines whether or not the end condition for deriving the iron loss in the electromagnetic field analysis unit 102 is satisfied. As the end conditions, for example, a predetermined number of generations have been changed, and the difference between the candidate solutions of the previous time and this time has become a predetermined condition (for example, the candidate solutions have changed between the previous time and this time). Conditions generally adopted in the method of genetic algorithms such as (disappeared) can be used.
Since the genetic algorithm itself can be realized by a known technique, detailed description thereof will be omitted here.
The setting of the candidate solution in the candidate solution setting unit 101 and the derivation of the iron loss of the iron core sample in the electromagnetic field analysis unit 102 are repeatedly executed until the end determination unit 103 determines that the end condition is satisfied.
＜出力部１０４＞
出力部１０４は、終了判定部１０３により終了条件を満足したと判定された直前に電磁場解析部１０２により導出された鉄損のうち、最小の鉄損に対応する候補解（目標信号波形）を最適解（最適な目標信号波形）として出力する。出力の形態としては、例えば、外部装置への送信、コンピュータディスプレイへの表示、または目標信号波形決定装置１００の外部または内部の記憶媒体への記憶が挙げられる。
以上のようにして、モータＭの速度指令の想定値およびトルク指令の想定値に対応する目標信号波形の最適値が１つ導出される。モータＭの速度指令の想定値およびトルク指令の想定値の組として想定される組のそれぞれに対して候補解設定部１０１、電磁場解析部１０２、終了判定部１０３、および出力部１０４の処理を行う。これにより、モータＭの速度指令値およびトルク指令値に対応する目標信号波形の最適値が得られる。
<
The
As described above, one optimum value of the target signal waveform corresponding to the assumed value of the speed command of the motor M and the assumed value of the torque command is derived. The candidate solution setting unit 101, the electromagnetic field analysis unit 102, the end determination unit 103, and the
［駆動装置２００］
駆動装置２００は、速度指令値およびトルク指令値に対応する目標信号波形を取得して、モータＭに印加される励磁電圧が目標信号波形になるようにインバータの動作を制御する。
図４は、駆動装置２００の機能的な構成の一例を示す図である。駆動装置２００のハードウェアは、例えば、ＣＰＵ、ＲＯＭ、ＲＡＭ、ＨＤＤ、各種のインターフェース、およびインバータ回路、または、専用のハードウェアを用いることにより実現される。
＜目標信号波形記憶部２０１＞
目標信号波形記憶部２０１は、モータＭの速度指令値およびトルク指令値と目標信号波形の情報とを相互に関連付けて記憶する。モータＭの速度指令値およびトルク指令値に関連付けられて記憶される目標信号波形の情報は、出力部１０４から出力される、モータＭの当該速度指令値および当該トルク指令値に対応する目標信号波形の最適値を特定する情報である。目標信号波形記憶部２０１は、モータＭの速度指令の想定値およびトルク指令の想定値の組として想定される組のそれぞれに対する目標信号波形の情報を記憶する。このようにして、モータＭの速度指令の想定値およびトルク指令の想定値の組として想定される組のそれぞれに対する目標信号波形の情報が記憶された後に、目標信号波形取得部２０２およびインバータ部２０３の動作が開始する。
[Drive device 200]
The
FIG. 4 is a diagram showing an example of a functional configuration of the
<Target signal
The target signal
＜目標信号波形取得部２０２＞
目標信号波形取得部２０２は、モータＭの速度指令値およびトルク指令値の現在値に対応する目標信号波形の情報を、目標信号波形記憶部２０１から取得する。モータＭの速度指令値およびトルク指令値の現在値は、モータＭの運転中に駆動装置２００の外部で発せられるものである。モータＭが三相モータである場合、目標信号波形取得部２０２は、例えば、モータＭの速度指令値およびトルク指令値の現在値に対応する目標信号波形の位相をずらすことにより、各相の目標信号波形を導出することができる。
<Target signal
The target signal
＜インバータ部２０３＞
インバータ部２０３は、インバータ回路と、当該インバータ回路を制御する制御回路とを有する。本実施形態では、目標信号波形は、パルス信号の時間波形である。インバータ回路は、スイッチング素子を有する。モータＭが三相モータである場合、スイッチング素子を、例えば、上アームおよび下アームを有する三相フルブリッジ型回路とすることができる。制御装置は、目標信号波形取得部２０２により取得された目標信号波形がインバータ回路から出力されるように各スイッチング素子の開閉動作を制御する。これにより、目標信号波形（または目標信号波形に近い時間波形）の励磁電圧がモータＭ（のステータコイル）に印加される。
<
The
［動作フローチャート］
次に、図５のフローチャートを参照しながら、目標信号波形決定装置１００による目標信号波形の決定方法の一例を説明する。図５のフローチャートは、モータＭの（１つの）速度指令の想定値および（１つの）トルク指令の想定値の組ごとに繰り返し実行される。ここでは、モータＭの（１つの）速度指令の想定値および（１つの）トルク指令の想定値の組を運転条件と称する。
[Operation flow chart]
Next, an example of a method of determining the target signal waveform by the target signal
ステップＳ５０１において、候補解設定部１０１は、運転条件を１つ入力する。
次に、ステップＳ５０２において、候補解設定部１０１は、初期の候補解を設定する。このとき、候補解設定部１０１は、ステップＳ５０１で入力した運転条件に含まれるモータＭの速度指令の想定値に応じた電気周期（電気周波数）になるように、複数の候補解（目標信号波形）を初期の候補解として設定する。
In step S501, the candidate solution setting unit 101 inputs one operating condition.
Next, in step S502, the candidate solution setting unit 101 sets an initial candidate solution. At this time, the candidate solution setting unit 101 has a plurality of candidate solutions (target signal waveforms) so that the electric cycle (electrical frequency) corresponds to the assumed value of the speed command of the motor M included in the operating conditions input in step S501. ) Is set as the initial candidate solution.
次に、ステップＳ５０３において、電磁場解析部１０２は、ステップＳ５０２で設定された複数の候補解のそれぞれについて、鉄心試料の鉄損を導出する。このとき、鉄心試料の磁束密度の実効値も導出される。
次に、ステップＳ５０４において、候補解設定部１０１は、ステップＳ５０３で導出された鉄心試料の鉄損および磁束密度の実効値に基づいて、ステップＳ５０２で設定された複数の候補解から、所定数の候補解を選択すると共に交叉や突然変異を行い、新たな候補解を設定する（複数の候補解を更新する）。
Next, in step S503, the electromagnetic field analysis unit 102 derives the iron loss of the iron core sample for each of the plurality of candidate solutions set in step S502. At this time, the effective value of the magnetic flux density of the iron core sample is also derived.
Next, in step S504, the candidate solution setting unit 101 has a predetermined number of candidate solutions from the plurality of candidate solutions set in step S502 based on the effective values of the iron loss and the magnetic flux density of the iron core sample derived in step S503. Select a candidate solution, cross or mutate it, and set a new candidate solution (update multiple candidate solutions).
次に、ステップＳ５０５において、終了判定部１０３は、終了条件を満足するか否かを判定する。この判定の結果、終了条件を満足しない場合、処理はステップＳ５０３に戻る。この場合、ステップＳ５０３において、電磁場解析部１０２は、ステップＳ５０４で更新された複数の候補解のそれぞれについて、鉄心試料の鉄損を導出する。そして、ステップＳ５０５において終了条件を満足すると判定されるまで、ステップＳ５０３～Ｓ５０５の処理が繰り返し実行される。 Next, in step S505, the end determination unit 103 determines whether or not the end condition is satisfied. As a result of this determination, if the end condition is not satisfied, the process returns to step S503. In this case, in step S503, the electromagnetic field analysis unit 102 derives the iron loss of the iron core sample for each of the plurality of candidate solutions updated in step S504. Then, the processes of steps S503 to S505 are repeatedly executed until it is determined in step S505 that the end condition is satisfied.
ステップＳ５０５において、終了条件を満足すると判定されると、処理はステップＳ５０６に進む。処理がステップＳ５０６に進むと、出力部１０４は、ステップＳ５０５で終了条件を満足すると判定された直前にステップＳ５０４で更新された複数の候補解のうち、最小の鉄損に対応する候補解を、ステップＳ５０１で入力した運転条件に対応する最適解（最適な目標信号波形）として出力する。そして、図５のフローチャートによる処理が終了する。尚、出力部１０４は、想定される全ての運転条件に対応する最適解が導出された時点で、各運転条件に対応する最適解（最適な目標信号波形）を出力してもよい。
If it is determined in step S505 that the end condition is satisfied, the process proceeds to step S506. When the process proceeds to step S506, the
次に、図６のフローチャートを参照しながら、駆動装置２００によるモータＭの駆動方法の一例を説明する。尚、図６のフローチャートは、想定される全ての運転条件に対応する最適解（最適な目標信号波形）の情報が、目標信号波形記憶部２０１に記憶された後に開始される。
Next, an example of a method of driving the motor M by the
ステップＳ６０１において、目標信号波形取得部２０２は、モータＭの速度指令値およびトルク指令値を、駆動装置２００の外部（外部装置）から取得するまで待機する。モータＭの速度指令値およびトルク指令値を、外部装置から取得すると、処理はステップＳ６０２に進む。
処理がステップＳ６０２に進むと、目標信号波形取得部２０２は、ステップＳ６０１で取得したモータＭの速度指令値およびトルク指令値に対応する目標信号波形の情報を、目標信号波形記憶部２０１から取得する。
In step S601, the target signal
When the process proceeds to step S602, the target signal
次に、ステップＳ６０３において、インバータ部２０３は、モータＭのステータコイルに印加される励磁電圧の時間波形が、ステップＳ６０２で取得された目標信号波形になるようにインバータ回路を動作させて目標信号波形を生成し、モータＭのステータコイルに印加（出力）する。そして、図６のフローチャートによる処理が終了する。尚、図６のフローチャートにおいて、ステップＳ６０２、Ｓ６０３は、ステップＳ６０１において、モータＭの速度指令値およびトルク指令値が、駆動装置２００の外部から取得されるたびに繰り返し実行される。
Next, in step S603, the
［計算例］
次に、計算例を説明する。本計算例では、解析対象のモータを同一のモータ（集中巻のＩＰＭＳＭ（Interior Permanent Magnet Synchronous Motor）とし、本実施形態の手法で求めた目標信号波形（電圧パルスパターン）をＩＰＭＳＭに印加した場合（発明例）の鉄損と、擬似正弦波の目標信号波形をＩＰＭＳＭに印加した場合（比較例）の鉄損とを電磁場解析により求めた。このとき、何れの手法においても、トルク指令値を３［Ｎｍ］、速度指令値を１００００［ｒｐｍ］とした。
[Calculation example]
Next, a calculation example will be described. In this calculation example, the motor to be analyzed is the same motor (IPMSM (Interior Permanent Magnet Synchronous Motor) with centralized winding, and the target signal waveform (voltage pulse pattern) obtained by the method of this embodiment is applied to the IPMSM ( The iron loss of the invention example) and the iron loss when the target signal waveform of the pseudo sine wave was applied to the IPMSM (comparative example) were obtained by electromagnetic field analysis. At this time, the torque command value was set to 3 in any of the methods. [Nm] and the speed command value were set to 10000 [rpm].
図７は、目標信号波形の一例を示す図である。図７（ａ）は、擬似正弦波の目標信号波形（比較例（正弦波ＰＷＭ））を示し、図７（ｂ）は、本実施形態の手法で求めた目標信号波形（発明例（最適化電圧波形））を示す。尚、図７（ａ）および図７（ｂ）では、電気一周期分の目標信号波形を示す。また、図７（ａ）および図７（ｂ）における目標信号波形の振幅および電気一周期は同じ値である。図８は、発明例と比較例における鉄損、渦電流損、およびヒステリシス損を示す図である。図８に示すように、本実施形態の手法を用いることで、目標信号波形を擬似正弦波とする場合に比べ、鉄損を８［％］程度低減できることが分かる。また、図７（ａ）および図７（ｂ）に示すように、モータの鉄損を低減する目標信号波形は、擬似正弦波とは大きく異なることが分かる。 FIG. 7 is a diagram showing an example of a target signal waveform. FIG. 7A shows a target signal waveform of a pseudo sine wave (comparative example (sine wave PWM)), and FIG. 7B shows a target signal waveform obtained by the method of the present embodiment (invention example (optimization)). Voltage waveform)) is shown. Note that FIGS. 7 (a) and 7 (b) show target signal waveforms for one electrical cycle. Further, the amplitude and the electric cycle of the target signal waveforms in FIGS. 7 (a) and 7 (b) have the same values. FIG. 8 is a diagram showing iron loss, eddy current loss, and hysteresis loss in the invention example and the comparative example. As shown in FIG. 8, it can be seen that by using the method of the present embodiment, the iron loss can be reduced by about 8 [%] as compared with the case where the target signal waveform is a pseudo sine wave. Further, as shown in FIGS. 7 (a) and 7 (b), it can be seen that the target signal waveform for reducing the iron loss of the motor is significantly different from the pseudo sine wave.
［まとめ］
以上のように本実施形態では、目標信号波形決定装置１００は、目標信号波形で鉄心試料を励磁した場合の、当該鉄心試料の板厚方向における磁束密度および渦電流密度の分布を導出することにより、当該鉄心試料の鉄損が最小になる目標信号波形を目標信号波形の最適値として導出する。駆動装置２００は、この目標信号波形の最適値の励磁電圧がモータＭのステータコイルに印加されるようにインバータ回路を動作させる。したがって、モータＭの鉄損を低減することができる励磁電圧の時間波形を短時間（実用的な計算時間内）に決定することができる。
[summary]
As described above, in the present embodiment, the target signal
また、本実施形態では、目標信号波形（パルスパターン）自体を最適化の設計変数とし、遺伝的アルゴリズム（最適化手法）を用いて設計変数の最適値を導出するため、解空間を限定しない探索が可能となる。したがって、特許文献１のような擬似正弦波に限定されず、目標信号波形（パルスパターン）を敢えて擬似正弦波（正弦波ＰＷＭ）とは異なる形にすることで、渦電流密度の減衰などを考慮して、モータＭの鉄損を低減することができる目標信号波形を設定することができる。
また、本実施形態では、マルチレベルインバータにおける多段階の電圧レベルを有する候補解（目標信号波形）を設定することにより、シングルレベルインバータに限らず、マルチレベルインバータに対しても本実施形態の手法を適用することができる。
Further, in the present embodiment, the target signal waveform (pulse pattern) itself is used as the optimization design variable, and the optimum value of the design variable is derived using the genetic algorithm (optimization method), so that the search does not limit the solution space. Is possible. Therefore, it is not limited to the pseudo sine wave as in Patent Document 1, and the target signal waveform (pulse pattern) is intentionally made different from the pseudo sine wave (sine wave PWM) to consider the attenuation of the eddy current density. Then, a target signal waveform that can reduce the iron loss of the motor M can be set.
Further, in the present embodiment, by setting a candidate solution (target signal waveform) having a multi-step voltage level in the multi-level inverter, the method of the present embodiment is applied not only to the single-level inverter but also to the multi-level inverter. Can be applied.
［変形例］
＜変形例１＞
本実施形態では、遺伝的アルゴリズムを用いて目標信号波形の最適値を導出する場合を例に挙げて説明した。しかしながら、必ずしも遺伝的アルゴリズムを用いて目標信号波形の最適値を導出する必要はない。例えば、遺伝的アルゴリズム以外のメタヒューリスティック手法による最適化計算を実行することにより、目標信号波形の最適値を導出してもよい。
[Modification example]
<Modification 1>
In this embodiment, a case where the optimum value of the target signal waveform is derived by using a genetic algorithm has been described as an example. However, it is not always necessary to derive the optimum value of the target signal waveform using a genetic algorithm. For example, the optimum value of the target signal waveform may be derived by executing the optimization calculation by a metaheuristic method other than the genetic algorithm.
＜変形例２＞
本実施形態では、有限要素法を用いて電磁場解析を行う場合を例に挙げて説明した。しかしながら、必ずしも有限要素法を用いて電磁場解析を行う必要はない。例えば、有限要素法以外の数値解析の手法（離散化手法）を用いて電磁場解析を行ってもよい。
<Modification 2>
In this embodiment, a case where electromagnetic field analysis is performed using the finite element method has been described as an example. However, it is not always necessary to perform electromagnetic field analysis using the finite element method. For example, electromagnetic field analysis may be performed using a numerical analysis method (discretization method) other than the finite element method.
＜変形例３＞
本実施形態では、ＰＷＭ制御を実行する場合を例に挙げて説明した。しかしながら、必ずしもこのようにする必要はない。例えば、ＰＡＭ（Pulse Amplitude Modulation）制御を実行するようにしてもよい。
＜変形例４＞
本実施形態では、目標信号波形決定装置１００と駆動装置２００とが別の装置である場合を例に挙げて示した。しかしながら、必ずしもこのようにする必要はない。例えば、目標信号波形決定装置１００の機能を駆動装置２００に含めてもよい。
<Modification 3>
In this embodiment, a case where PWM control is executed has been described as an example. However, it is not always necessary to do this. For example, PAM (Pulse Amplitude Modulation) control may be executed.
<Modification example 4>
In the present embodiment, the case where the target signal
＜その他の変形例＞
尚、以上説明した本発明の実施形態のうち、目標信号波形決定装置１００および駆動装置２００のインバータ回路を除く部分は、コンピュータがプログラムを実行することによって実現することができる。また、前記プログラムを記録したコンピュータ読み取り可能な記録媒体および前記プログラム等のコンピュータプログラムプロダクトも本発明の実施形態として適用することができる。記録媒体としては、例えば、フレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ－ＲＯＭ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
また、以上説明した本発明の実施形態は、何れも本発明を実施するにあたっての具体化の例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
<Other variants>
In the embodiment of the present invention described above, the portion of the target signal
In addition, the embodiments of the present invention described above are merely examples of embodiment of the present invention, and the technical scope of the present invention should not be construed in a limited manner by these. It is a thing. That is, the present invention can be implemented in various forms without departing from the technical idea or its main features.
１００：目標信号波形決定装置、１０１：候補解設定部、１０２：電磁場解析部、１０３：終了判定部、１０４：出力部、２００：駆動装置、２０１：目標信号波形記憶部、２０２：目標信号波形取得部、２０３：インバータ部、Ｍ：モータ 100: Target signal waveform determination device, 101: Candidate solution setting unit, 102: Electromagnetic field analysis unit, 103: End determination unit, 104: Output unit, 200: Drive device, 201: Target signal waveform storage unit, 202: Target signal waveform Acquisition part, 203: Inverter part, M: Motor
Claims (10)
前記電気機器に印加される励磁信号の時間波形の目標値である目標信号波形を決定する目標信号波形決定手段を有し、
前記目標信号波形決定手段は、前記鉄心の鉄損解析用試料が励磁された際の磁束密度および渦電流密度の分布を、マクスウェルの方程式に基づいて導出し、当該磁束密度および渦電流密度の分布に基づいて、当該鉄損解析用試料の鉄損を導出し、当該鉄損解析用試料の鉄損に基づいて、前記目標信号波形を決定し、
前記磁束密度および渦電流密度の分布は、前記鉄損解析用試料の内部において磁束および渦電流が浸透する深さ方向における分布を含むことを特徴とする処理システム。 It is a processing system that performs processing to operate electrical equipment with an iron core.
It has a target signal waveform determining means for determining a target signal waveform which is a target value of a time waveform of an excitation signal applied to the electric device.
The target signal waveform determining means derives the distribution of the magnetic flux density and the eddy current density when the iron loss analysis sample of the iron core is excited based on Maxwell's equations, and the distribution of the magnetic flux density and the eddy current density. Based on, the iron loss of the iron loss analysis sample is derived, and the target signal waveform is determined based on the iron loss of the iron loss analysis sample.
The processing system characterized in that the distribution of the magnetic flux density and the eddy current density includes a distribution in the depth direction in which the magnetic flux and the eddy current permeate inside the sample for iron loss analysis.
前記目標信号波形の候補は、前記モータの速度指令値に対応する周期を有し、
前記目標信号波形決定手段は、前記モータのトルク指令値に対応する磁束密度を満足する磁束密度が得られる範囲内で前記鉄損解析用試料の鉄損が最小になる前記目標信号波形の候補を、前記目標信号波形として決定することを特徴とする請求項３に記載の処理システム。 The electric device is a motor and
The candidate of the target signal waveform has a period corresponding to the speed command value of the motor, and has a period corresponding to the speed command value of the motor.
The target signal waveform determining means selects a candidate for the target signal waveform that minimizes the iron loss of the iron loss analysis sample within a range in which a magnetic flux density satisfying the magnetic flux density corresponding to the torque command value of the motor can be obtained. The processing system according to claim 3, wherein the target signal waveform is determined.
前記駆動手段は、前記モータの速度指令値およびトルク指令値と、前記目標信号波形を特定するための情報とを相互に関連付けて記憶する目標信号波形記憶手段と、
前記モータの速度指令値およびトルク指令値に対応する前記目標信号波形を、前記目標信号波形記憶手段により記憶された前記情報に基づいて取得する目標信号波形取得手段と、を有し、
前記目標信号波形取得手段により取得された前記目標信号波形を目標値として励磁信号の時間波形を生成して前記モータに印加することを特徴とする請求項４に記載の処理システム。 Further having a driving means for generating a time waveform of an excitation signal with the target signal waveform determined by the target signal waveform determining means as a target value and applying the time waveform to the electric device.
The drive means includes a target signal waveform storage means for storing the speed command value and the torque command value of the motor and information for specifying the target signal waveform in association with each other.
It has a target signal waveform acquisition means for acquiring the target signal waveform corresponding to the speed command value and the torque command value of the motor based on the information stored by the target signal waveform storage means.
The processing system according to claim 4, wherein a time waveform of an excitation signal is generated using the target signal waveform acquired by the target signal waveform acquisition means as a target value and applied to the motor.
前記磁束密度および渦電流密度の分布は、前記板の板厚方向の分布であることを特徴とする請求項１～５の何れか１項に記載の処理システム。 The iron loss analysis sample has one or more plates and has one or more plates.
The processing system according to any one of claims 1 to 5, wherein the distribution of the magnetic flux density and the eddy current density is a distribution in the plate thickness direction of the plate.
前記励磁信号は、パルス信号であることを特徴とする請求項１～６の何れか１項に記載の処理システム。 The target signal waveform is a target signal waveform for the inverter.
The processing system according to any one of claims 1 to 6, wherein the excitation signal is a pulse signal.
前記電気機器に印加される励磁信号の時間波形の目標値である目標信号波形を決定する目標信号波形決定工程を有し、
前記目標信号波形決定工程は、前記鉄心の鉄損解析用試料が励磁された際の磁束密度および渦電流密度の分布を、マクスウェルの方程式に基づいて導出し、当該磁束密度および渦電流密度の分布に基づいて、当該鉄損解析用試料の鉄損を導出し、当該鉄損解析用試料の鉄損に基づいて、前記目標信号波形を決定し、
前記磁束密度および渦電流密度の分布は、前記鉄損解析用試料の内部において磁束および渦電流が浸透する深さ方向における分布を含むことを特徴とする処理方法。 It is a processing method for operating an electric device having an iron core.
It has a target signal waveform determination step of determining a target signal waveform which is a target value of a time waveform of an excitation signal applied to the electric device.
In the target signal waveform determination step, the distribution of the magnetic flux density and the eddy current density when the iron loss analysis sample of the iron core is excited is derived based on Maxwell's equations, and the distribution of the magnetic flux density and the eddy current density is derived. Based on, the iron loss of the iron loss analysis sample is derived, and the target signal waveform is determined based on the iron loss of the iron loss analysis sample.
The processing method characterized in that the distribution of the magnetic flux density and the eddy current density includes a distribution in the depth direction in which the magnetic flux and the eddy current permeate inside the sample for iron loss analysis.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018177461A JP7087875B2 (en) | 2018-09-21 | 2018-09-21 | Processing system, processing method, and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018177461A JP7087875B2 (en) | 2018-09-21 | 2018-09-21 | Processing system, processing method, and program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2020048384A JP2020048384A (en) | 2020-03-26 |
| JP7087875B2 true JP7087875B2 (en) | 2022-06-21 |
Family
ID=69901842
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2018177461A Active JP7087875B2 (en) | 2018-09-21 | 2018-09-21 | Processing system, processing method, and program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7087875B2 (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118216082A (en) * | 2021-11-15 | 2024-06-18 | 日本制铁株式会社 | Excitation waveform determining device, motor driving device, excitation waveform determining method, motor driving method, and program |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008220035A (en) | 2007-03-02 | 2008-09-18 | Ricoh Co Ltd | Brushless motor starter |
| JP2013131071A (en) | 2011-12-21 | 2013-07-04 | Fujitsu Ltd | Program, apparatus and method for magnetic body property analysis |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4995518B2 (en) * | 2006-09-14 | 2012-08-08 | 三菱電機株式会社 | AC motor control device and method for calculating superimposed current for suppressing iron loss |
| JP4932640B2 (en) * | 2007-08-23 | 2012-05-16 | 新日本製鐵株式会社 | Iron loss optimization system |
| JP6384208B2 (en) * | 2014-09-01 | 2018-09-05 | 新日鐵住金株式会社 | Electromagnetic field analysis apparatus, electromagnetic field analysis method, and program |
-
2018
- 2018-09-21 JP JP2018177461A patent/JP7087875B2/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008220035A (en) | 2007-03-02 | 2008-09-18 | Ricoh Co Ltd | Brushless motor starter |
| JP2013131071A (en) | 2011-12-21 | 2013-07-04 | Fujitsu Ltd | Program, apparatus and method for magnetic body property analysis |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2020048384A (en) | 2020-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Min et al. | 3-D performance analysis and multiobjective optimization of coreless-type PM linear synchronous motors | |
| JP6206608B1 (en) | Electromagnetic field analysis apparatus, electromagnetic field analysis method, and program | |
| JP6384208B2 (en) | Electromagnetic field analysis apparatus, electromagnetic field analysis method, and program | |
| JP5948766B2 (en) | Simulation method for permanent magnet rotating electrical machine, loss calculation method using the same, simulation program, loss calculation program, simulation system, and loss calculation system | |
| JP6834670B2 (en) | Magnetic property measurement system and magnetic property measurement method | |
| JP7087875B2 (en) | Processing system, processing method, and program | |
| Pinhal et al. | Performance map calculation of a salient-pole synchronous motor with hairpin winding | |
| Li et al. | Fast and accurate analytical calculation of the unsaturated phase inductance profile of 6/4 switched reluctance machines | |
| JP6798385B2 (en) | Magnetic property analysis system and magnetic property analysis method | |
| Pries et al. | Magnetic and thermal scaling of electric machines | |
| CN114580235A (en) | Design method of vibration scaling model of converter transformer | |
| JP6380072B2 (en) | Magnetic characteristic measuring method, magnetic characteristic measuring system, and program | |
| CN112514242B (en) | Excitation system, excitation method, program, and modulation operation setting device for iron core in electric device, and modulation operation setting device for inverter power supply | |
| JP7436778B2 (en) | Processing systems, processing methods, and programs | |
| JP6578845B2 (en) | Electromagnetic field analysis apparatus, electromagnetic field analysis method, and program | |
| JP2006258480A (en) | System and method for analyzing magnetic property | |
| Geldhof et al. | Analysis of hysteresis in resonance-based position estimation of switched reluctance drives | |
| JP6984426B2 (en) | Electromagnetic field analyzer, electromagnetic field analysis method, and program | |
| JP6984427B2 (en) | Electromagnetic field analyzer, electromagnetic field analysis method, and program | |
| Pinhal et al. | Discussion on the control of wound-rotor synchronous motors for variable speed drives | |
| JP2019204269A (en) | Electromagnetic field analysis system, electromagnetic field analysis method, and program | |
| Kushwaha et al. | Response surface-assisted trust-region method for optimization of electromagnetic devices | |
| Gao et al. | Losses analysis of soft magnetic ring core under sinusoidal pulse width modulation (SPWM) and space vector pulse width modulation (SVPWM) excitations | |
| Yaich et al. | Control strategy for switched reluctance motor based on embedded system | |
| Çakal et al. | Design of double sided linear motor with easy to manufacture hairpin plate winding |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210514 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20220425 |
|
| 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: 20220510 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220523 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 7087875 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |